JP6682500B2 - 信号伝搬時間差式流量計 - Google Patents

信号伝搬時間差式流量計 Download PDF

Info

Publication number
JP6682500B2
JP6682500B2 JP2017503164A JP2017503164A JP6682500B2 JP 6682500 B2 JP6682500 B2 JP 6682500B2 JP 2017503164 A JP2017503164 A JP 2017503164A JP 2017503164 A JP2017503164 A JP 2017503164A JP 6682500 B2 JP6682500 B2 JP 6682500B2
Authority
JP
Japan
Prior art keywords
signal
measurement
connector
response signal
response
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.)
Active
Application number
JP2017503164A
Other languages
English (en)
Other versions
JP2017525948A (ja
Inventor
トーマス ワーナー ヒース,
トーマス ワーナー ヒース,
トラン ダン ルオン,
トラン ダン ルオン,
クラウス−ディーター オール,
クラウス−ディーター オール,
ユルゲン ハインツ−フリードリヒ スクリパレ,
ユルゲン ハインツ−フリードリヒ スクリパレ,
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.)
GWF MessSysteme AG
Original Assignee
GWF MessSysteme AG
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
Application filed by GWF MessSysteme AG filed Critical GWF MessSysteme AG
Publication of JP2017525948A publication Critical patent/JP2017525948A/ja
Application granted granted Critical
Publication of JP6682500B2 publication Critical patent/JP6682500B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/245Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by measuring transit time of acoustical waves
    • 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/66Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • 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/66Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
    • G01F1/667Arrangements of transducers for ultrasonic flowmeters; Circuits for operating ultrasonic flowmeters
    • 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
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/241Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • G01P5/24Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
    • G01P5/241Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect
    • G01P5/244Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave by using reflection of acoustical waves, i.e. Doppler-effect involving pulsed waves
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Description

本出願は流量計に関し、特に伝搬時間差式超音波流量計に関する。
現在、各種の流量計が、パイプを流れる液体や気体など流体の体積流量を測定するために使われている。超音波流量計としては、音響ドップラー効果を利用したドップラー式流量計や、超音波源と媒体の相対移動によって生じる伝搬時間差を利用した、伝送式流量計と呼ばれることもある、信号伝搬時間差式流量計が存在する。伝搬時間(travel time)は、タイムオブフライト(Time of flight)あるいはトランジットタイム(transit time)とも称される。
伝搬時間差式超音波流量計は、流れ方向及びその反対方向に伝搬する超音波パルスの伝搬時間差を評価する。各種の超音波流量計が、侵入(intrusive)式又は湿式流量計として知られるインライン流量計、または、非侵入式流量計として知られるクランプオン流量計として提供されている。その他の流量計として、ベンチュリ管方式、オーバフローシル(overflow sill)方式、レーダ式流量計、コリオリ式流量計、差圧式流量計、電磁誘導式流量計などの種類がある。
不規則な流量分布が存在する場合や開水路(open channel)では、平均流速を求めるために、1より多い伝搬路が必要になることがある。特にマルチパス方式の手順は、IEC41やEN ISO 6416などの流速測定規格に記載されている。更に別の応用として、超音波流量計は、例えば超音波ドップラー式流速計(ADCP)によって流量分布を測定するのにも使われている。ADCPは、河川や開水域(open water)において水の速度や流出を測定するのにも適している。
本明細書の目的は、一般には流体、特に水、油などの液体又は気体の平均流速あるいは流量分布を測定する改良型の伝搬時間差式流量計及びそれに対応した方法を提供することにある。
本明細書による流量測定装置では、テスト信号及び測定信号を発信及び受信するために、圧電トランスデューサとしても知られる、圧電素子などの音響トランスデューサが用いられる。
また、音波送信器として、金属膜を励振して振動させるレーザや、単純なスピーカを備えてもよい。また、別の方法で圧力波を発生させることもできる。受信器側は、超音波を検出できれば、圧電トランスデューサと異なる他の手段を用いてもよい。
本明細書では「圧電トランスデューサ」という用語を多く用いるが、この用語は超音波を発生又は検出するその他の音波トランスデューサも意味する。
本明細書における測定信号は、整合フィルタによってモデル形成が可能である。プローブ又はテスト信号としてピークがシャープなインパルスが使われる場合、トランスデューサにおける受信信号は、流体の導管や流路でのインパルス応答となる。本出願によれば、インパルス応答の時間的な反転が、測定信号として、同じ流路を通して、逆方向又は同一方向に送信される。その結果、ピークを有する信号が、本来の発生源が存在する起点あるいは元の受信器において得られる。
時間的な反転は、いくつかの方法で実現できる。応答信号を記録するのにアナログ手段が使われる場合には、記録された応答信号を逆モードで再生すればよい。応答信号の各サンプルを記録するのにデジタル手段が使われる場合には、記録された各サンプルの順序を逆にすれば、反転信号が得られる。後者は、記録された各サンプルの時間スタンプの値を、それぞれの時間値に (−1) を掛けて反転することによって行える。得られた時間スタンプの各値を昇順で再生すれば、記録された各サンプルが逆の順序で再生される。すなわち、反転応答信号は、記録された応答信号が逆向きに再生されたものである。
本明細書による超音波流量計は、上記の反転信号又は同様の形状の信号を用いることで、超音波流量計に集束性を与え、空間及び時間両方において集束された応答信号を形成する。その結果、受信側の圧電素子でより大きい振幅が得られ、より優れたS/N比が得られる。
本明細書による超音波流量計では、上記の集束性は極めて一般的な条件下において得られる。例えば、集束性は、1つの超音波発振器だけを励起して、反転信号を振幅範囲で粗くデジタル化された信号にした場合でも、反転信号の時間分解能が充分であれば得られる。更に、本明細書による超音波流量計は、パイプに取り付け容易で、パイプの改造を必要としないクランプオン式トランスデューサでも使用可能である。
本明細書による流量測定方法では、測定信号の振幅に関するビット分解能を調整できる。特に、大きい振幅の応答信号が得られるように、ビット分解能を調整可能である。
一実施形態によれば、測定信号に対する応答信号の振幅を大きくするために、ビット分解能が高められる。一実施形態において、ビット分解能は予め決められた工程により高められ、最大の振幅を有する応答信号をもたらすビット分解能が選択され、対応する測定信号の表記がコンピュータのメモリに記憶される。
別の実施形態によれば、測定信号に対する応答信号の振幅を大きくするために、ビット分解能が減少される。一実施形態において、ビット分解能は予め決められた工程により減少され、最大の振幅を有する応答信号をもたらすビット分解能が選択され、対応する測定信号の表記がコンピュータのメモリに記憶される。
具体的には、ビット分解能は、1桁又は2桁などで記憶される分解能、特に、1個又は2個の2進数などで記憶される低ビット分解能でもよい。他の実施形態によれば、低ビット分解能は、最小で1ビット分解能、最大で64ビット分解能とすることができる。
別の実施形態によれば、導管の壁厚の変化を求めるため、あるいは縦波及び横波の音波特性を測定して導管壁の素材特性を求め又は導くために、第1の応答信号が処理される。例えば、縦波及び横波の音波特性は、音波の異なる到着時間に対応する受信信号つまり応答信号の対応部分から導出し得る。
この実施形態では、流速を求めるのと、上記の素材特性を検出するのに、同じ応答信号を用いることができる。従って、汚れや素材欠陥などの影響を検出するために、別個の信号又は別個の装置を用いる必要がない。但し、別個の信号又は別個の装置を用いることもできる。更に、導き出された伝搬路特性は、流速をより正確に推定するために使うこともできる。
本明細書による超音波流量計においては、クランプオン式トランスデューサの連結性や方向性を良好にするための、あるいは、散乱を減少させるための技術的特徴が必要とされない。逆に、散乱が集束性を向上させることもある。散乱を増大させるために、液体やトランスデューサの屈折率に適合した連結材料を選択したり、より多くのせん断波を与えるトランスデューサ連結具を用いてもよい。
本明細書による流量計で使われる音波の周波数は、0.5 マイクロ秒(μs)の振動周期に対応する 20kHz超と2MHzとの間が好ましいが、800MHz程度の高い周波数でもよい。
多くの場合、超音波流量計は、可聴しきい値よりはるかに高い数百kHz以上の周波数で作動する。伝搬時間差式超音波流量計の周波数は通常、kHz又はMHzの単位範囲である。
1の観点からすれば、本明細書は、流体導管又は流路中の流体の流速を、伝搬時間差式超音波流量計を用いて求めるための、コンピュータで実行される方法を開示する。特に本方法は、パイプ(管)やチューブだけでなく、排水路や灌漑路への応用など、開水路に関する用途にも使える。好ましい実施形態において、「コンピュータで実行される」とは、一般的にはワークステーションやメインフレームコンピュータより小さいサイズの、流体パイプに沿った必要箇所に設置可能な、可搬式もしくはコンパクトな固定デジタル信号処理装置で使用できるマイクロプロセッサ、ASICs、FPGAsなどの小規模電子構成要素において実行されることを意味する。
以下の説明において、用語「流路」、「導管」、「通路」、「パイプ」等は同義語として使われる。本出願の主題は、あらゆる種類の流体用導管に対して、それらの導管の形状や、それらの導管が開放されているか閉じられているか、もしくは全体に満たされているか部分的に満たされているかにかかわりなく、適用可能である。また本出願の主題は、あらゆる種類の流体又は気体に対して、それらが気体であるか液体であるか、もしくは両方の混合物であるかにかかわりなく適用可能である。
測定信号の生成フェーズでは、流体導管に対して所定の速度を有する流体、特に流体導管に対して実質上静止している流体が、流体導管に供給される。測定信号は、最初に与えられるインパルス信号に対して伝搬路が返信する応答信号から得られる。
インパルス信号は、例えば圧電トランスデューサである第1の超音波トランスデューサに与えられる。ここで、インパルス信号とは、特に短い時間に集中された信号エネルギーを有する信号を意味する。ある特定の実施形態において、インパルス信号は、10〜20以下の振動周期など、搬送波のわずかな振動周期の間だけ持続する。具体的には、インパルス信号の包絡線は矩形状でもよいが、その他の形状も可能である。例えば、インパルス信号は、1つの時間ピークつまり単一のインパルス、短い矩形バースト、又は鋸波状、矩形波状、チャープ(chirp) 、サイン波などその他任意の信号形状、もしくは1/fノイズとしても知られるホワイトノイズやピンクノイズなど所定のノイズバーストに対応するものであってもよい。本方法は、インパルス信号がほぼどんな信号形状でも機能する。
信号生成フェーズは、各測定毎に繰り返す必要はない。例えば、信号生成フェーズは、最初の測定の前と、それより後であって、沈殿、腐食及び熱応力などにより、流体導管の条件が変化した時点に実施することができる。
「較正フェーズ」という用語が、測定信号生成フェーズを指して使われることがある。しかし、これは全く正しくない。流量測定の場合、流量計は、流量の測定値と目標値が比較される較正装置に配置されるのが一般的である。それら2つの値間の関連因子が較正因子と呼ばれ、較正因子は、流量測定において特定できないハードウェア及びソフトウェア両方の誤差を含んでいる。本出願の主題においては、測定信号生成フェーズと較正フェーズを区別する方がより適切である。測定信号生成フェーズは、測定信号が使われた際、測定信号に対してシャープなピークを有する応答信号を生じさせる、測定信号を提供する一方、較正フェーズは、正確な流量測定をもたらす流量計を提供する。
本明細書による方法の下記の各ステップ:
−インパルス信号を、流体導管の第1の箇所に位置された第1の超音波トランスデューサに与えること、
−前記インパルス信号に対する応答信号を、前記流体導管の第2の箇所に位置された第2の超音波トランスデューサで得ること、
−前記応答信号から測定信号を導出すること、且つ該測定信号の導出は応答信号の一部分又は応答信号から導出され信号の一部分を選択すること、及び前記選択された信号部分を時間反転することを含むこと、
は、実際の信号を現実に実際の導管に与え、測定することで得られる。但し、インパルス信号がデジタル又はアナログ信号により与えられれば、インパルス信号に対する応答信号を第2の超音波トランスデューサで得るステップ、及び測定信号を導出するステップは、数値シミュレーション又はアナログシミュレーションによって得られることが判明している。それには、有限要素法ソフトウェアが使用可能である。
複数の圧電トランスデューサは、流体導管に配置される。すなわち、複数の圧電トランスデューサは、それぞれ流体導管に取り付けて配置可能である。第1の圧電トランスデューサは、流体導管周囲の第1の位置に取り付けて配置される。特定の一実施形態において、第1の圧電トランスデューサは流体導管の周囲にクランプ止めされる。インパルス信号に対する応答信号は、第2の圧電トランスデューサで受信される。
例えば圧電トランスデューサである第2の超音波トランスデューサは、流体導管周囲の第2の位置に取り付けて配置され、第2の位置は、第1の位置に対し導管の中心軸を通る断面に沿って、流体導管の長さ方向にずれている。ここで、流体導管の長さ方向は、流路を流れる平均流の方向に対応する。流体表面における反射及びその他の影響が望ましくないときは、流体導管が流体で完全に満たされる。
測定信号は、最初のインパルス信号に対する伝搬路の応答である応答信号から、アナログ的手段で又はデジタル的手段により導出される。測定信号の導出は、応答信号の一部分又は応答信号から導出された信号の一部分を選択すること、及び前記選択された信号部分を時間反転することを含む。またその導出は、測定信号をその後における使用のため、例えばデジタル化した形で、コンピュータで読み取り可能なメモリに記憶するステップを含んでもよい。その際、方法の各ステップは、異なる順序で実施可能である。例えば、信号の時間反転は、信号を記憶した後に行ってもよい。
流体が圧力、重力、パイプの傾斜などの外的条件に従い流体導管に対して移動する測定フェーズの間、測定信号が、圧電トランスデューサなどの第1及び第2の超音波トランスデューサの一方に与えられる。つまり、記憶した測定信号から導き出された電気信号が、上記トランスデューサに与えられる。
測定信号に対する第1の応答信号が、圧電トランスデューサなどの他方の超音波トランスデューサで測定され、流体の流速が、少なくとも第1の応答信号から導出される。具体的には、この測定は、上流又は下流両方向の伝搬時間を測定することを含む。速度の推定は、測定された伝搬時間と、例えば流体の温度を測定して得た現在の条件下における音速を考慮して較正した伝搬時間とを比較することによって得られる。別のステップにおいて、体積流量や質量流量を、流速や流速の分布から導出してもよい。
より正確な推定を得るため、測定は両方向、つまり各々が例えば圧電トランスデューサである第1の超音波トランスデューサから第2の超音波トランスデューサへの方向とその逆方向に行ってもよい。そのようにすれば特に、伝搬時間を音速として推定することが可能となり、また現時点における音速の信頼度の高い推定が可能になる。
本明細書による流量測定は、2つのトランスデューサを備えた構成だけでなく、図43及び図44の構成又は図4及び図5の構成などのように、2つよりも多いマルチトランスデューサの構成でも行える。具体的には、相互に対向して配置された一対のトランスデューサ又はマルチトランスデューサの構成によって流量測定を実施可能である。対のトランスデューサは、図43に示すように導管の中心軸を通る平面に配置されてもよいし、図44に示すように導管の中心軸に対してずれた平面に配置されてもよい。図44の構成は、中心軸から所定の距離に位置する流体層における流速を求めるのに使える。
このように、測定信号を与えるステップと応答信号を測定するステップは、逆方向にも繰り返される。言い換えれば、先の受信器が送信器として使われ、先の送信器が受信器として使われ、各々が例えば圧電トランスデューサである2つの超音波トランスデューサのうち最初の他方の超音波トランスデューサから、最初の一方の超音波トランスデューサへと信号が送られ、第2の応答信号が得られる。流体の流速は、第1の応答信号と第2の応答信号から導出される。具体的には、この導出は下流及び上流両方向の伝搬時間を導くことを含む。
例えば圧電トランスデューサである1つの超音波トランスデューサから、例えば圧電トランスデューサである別の超音波トランスデューサへ測定信号を送ることができるが、流速又は流量測定を行う場合に、上記の送信を順方向と逆方向両方で行うのが有利である。
すなわち、その手順は次の通りである。
順方向:
−インパルス信号を、第1の超音波トランスデューサから第2の超音波トランスデューサへ送る、
−インパルス信号に対する応答信号を、第2の超音波トランスデューサで受信する、
−第2の超音波トランスデューサで受信した応答信号を時間反転し、測定信号を得る、
−測定信号を、第1の超音波トランスデューサから第2の超音波トランスデューサへ送る、
−測定信号に対する応答信号を、第2の超音波トランスデューサで受信する。
逆方向:
−インパルス信号を、例えば各々圧電トランスデューサである第2の超音波トランスデューサから第1の超音波トランスデューサへ送る、
−インパルス信号に対する応答信号を、第1の超音波トランスデューサで受信する、
−第1の超音波トランスデューサで受信した、インパルス信号に対する応答信号を時間反転し測定信号を得る、
−測定信号を、第2の超音波トランスデューサから第1の超音波トランスデューサへ送る、
−測定信号に対する応答信号を、第1の超音波トランスデューサで受信する、
−第2の超音波トランスデューサと第1の超音波トランスデューサで受信した応答信号間の時間差を測定する。この時間差は、例えば各々圧電トランスデューサである2つの超音波トランスデューサ間での流速に比例する。
なお、順方向での測定信号と逆方向での測定信号は異なるものであってもよい。簡単な構成とする場合には、同一の測定信号を使用できるが、通常、測定信号は各伝搬方向において特有の形状を有する。
本出願全体を通して、「コンピュータ」という用語がよく用いられている。コンピュータはラップトップやデスクトップコンピュータなどの装置を備えるものであるが、マイクロコントローラ、ACIDs、FPGAsなどによって、信号の送信および受信も行える。
所定の層における流速を得るため、トランスデューサ間における想定接続線は流体導管の中心から幾何学的にずれていてもよく、また一対より多いトランスデューサを設けてもよい。さらに、測定信号を1より多いトランスデューサから与えること、及び測定信号に対する応答信号を1より多いトランスデューサで測定することの両方又は一方を行うようにしてもよい。
簡単な実施形態では、複数の受信トランスデューサからの応答信号を線形重畳して平均の応答信号が生成され、この平均の応答信号に上述した信号処理ステップを施して、測定信号が得られる。
さらに別の実施形態では、例えばN個の等しい数の送信および受信トランスデューサが設けられ、複数の送信トランスデューサの相対的な配置が、複数の受信トランスデューサの相対的な配置と同等にされる。N個の受信された応答信号がそれぞれ上述の信号処理ステップに従って処理され、N個の測定信号が各々得られる。
上記N個のトランスデューサは通例、例えばクランプオン式トランスデューサあるいは挿入式又は内部装着式トランスデューサとして配置される。一例として、図43は8個のクランプオン式トランスデューサの配置を示し、図44は8個の挿入式トランスデューサの配置を示す。図43における8個のトランスデューサは、それぞれ導管の軸中心を通る4つの平面に配置されている。
図44における8個の挿入式トランスデューサは、4つの平行な平面に配置されている。
トランスデューサ間における接続線は、トランスデューサの動作モードを示している。図43の動作モードでは、各信号が第1のトランスデューサから、水導管の中心軸上の中心点を挟んで第1のトランスデューサと対向する第2のトランスデューサに送られる。
図44の動作モードでは、各信号が第1のトランスデューサから、4つの平行な平面の1つ内で各方形配置の中心に位置する中心点を挟んで対向する第2のトランスデューサに送られる。
一実施形態によれば、応答信号のうち測定信号を導出するのに使われる信号部分は、応答信号の最大振幅の周囲の第1部分と、後続する信号部分とを含み、後続信号部分は最大振幅の到着時間より遅い時間側に延びている。後続信号部分は、ダイレクト信号の近傍における信号とは別の、他の反射から生じた信号も含み、さらに良好な集束をもたらすことができる。
より改善された測定信号を生成するために、インパルス信号を与えるステップ及び対応する応答信号を受信するステップは、一回行うだけでなく、複数回、少なくとも2回繰り返してもよい。繰り返すことによって、複数の応答信号が得られる。これら複数の受信された応答信号の平均から、測定信号が導出される。
一実施形態においては、測定が複数回繰り返されるが、超音波信号は一方向にだけ伝搬する。別の実施形態においては、超音波信号を両方向に伝搬させて、測定が複数回繰り返される。さらに別の実施形態においては、測定が両方向で複数回繰り返され、両方向についてそれぞれ別の平均が導出される。
更なる実施形態では、1つより多い受信信号からの測定信号の導出が、応答信号の包絡線又は応答信号から導出した信号の包絡線を求めることを含む。この場合、包絡線に基づいて振幅変調された信号である振幅変調振動信号が得られる。信号サンプルの代わりに、もしくはそれに加えて包絡線を用いることは、記憶スペース及び演算速度の点で有利である。
具体的には、変調振幅が、測定信号又はその一部分について求められた包絡線の形状を有してもよい。搬送波振動の振動周期は、少なくとも20kHzである。別の実施形態では、周波数が少なくとも100kHz、少なくとも500kHz、あるいは少なくとも1MHzである。周波数の選択は散乱過程に影響を及ぼし、周波数が高いほど、導管壁のサンプリングがよりきめ細かくなり、ひいては超音波信号のより正確な形成を可能にすることがある。
更なる実施形態では、応答信号又は応答信号から導出した信号が、特に1ビットから8ビットの間の分解能で、振幅に関してデジタル化される。本明細書は、振幅を粗くデジタル化した場合でも、超音波信号の充分な集束が得られることを示している。低い分解能の使用が、演算時間及びメモリスペースを節約させる一方、より高い分解能は、導管を流れる流体の流量についてより正確な測定結果を必ずしも与えるものではない。応答信号又は測定信号の分解能の増減が、S/N比及び時間測定の精度を向上させるのに役立つことも判明している。分解能の減少は、測定信号に対する応答信号において、よりシャープな、つまり、より特徴的なピークをもたらす。これは、S/N比を高くするのに、測定信号の送信出力を大きくする代わりに、測定信号又は測定信号に対する応答信号の分解能を減少させてもよいことを意味している。
本明細書の別の観点によれば、流体導管つまりパイプ内の流体の流速を求める一部の方法において、伝搬時間差式超音波流量計の振幅変調測定信号又は振幅変調応答信号を用いてもよい。この方法では、信号生成フェーズのステップを必ずしも含まなくてもよいが、測定信号を得るために1回の信号生成フェーズのステップを用いてもよい。例えば、上記方法では、工場現場において予め生成した測定信号に依拠してもよい。その場合、例えば圧電トランスデューサである他方の超音波トランスデューサから送信された一連の発振を受信した、例えば圧電トランスデューサである一方の超音波トランスデューサからの時間反転された受信信号を、測定信号として生成することができる。
第1のステップでは、圧力、重力、パイプの傾斜などの外的条件下で流体導管に対して移動する流体が、流体導管に供給される。
第1の超音波トランスデューサは、流体導管の第1の位置に設けられる。例えば圧電トランスデューサである第2の超音波トランスデューサは、流体導管の第2の位置に設けられる。第2の位置は第1の位置に対して流体導管の長さ方向に沿ってずれており、流体導管の長さ方向は流体路の流体流方向に対応している。
測定信号は,例えば各々圧電トランスデューサである第1又は第2の超音波トランスデューサに与えられる。具体的には、振幅変調信号から導出された電気信号が、当該トランスデューサに送信可能である。
測定信号に対する第1の応答信号が、例えば圧電トランスデューサである他方の超音波トランスデューサで測定され、流体の流速が第1の応答信号から導出される。具体的には、この導出は下流又は上流方向の伝搬時間を導くことを含む。
前述の方法と同じく、測定を逆方向にも繰り返し下流及び上流両方向の伝搬時間を得ることで、より高い精度が達せられる。図43及び図44に示したように、例えば、平均フローをより正確に推定したり、流体導管の中心軸から所定の距離にある平面におけるフローを推定するために、N対のトランスデューサを用いることもできる。
具体的には、測定信号を与えるステップと応答信号を測定するステップとが逆方向に繰り返されて第2の応答信号が得られ、流体の流速が第1の応答信号および第2の応答信号から導出される。ここで、その導出は、下流及び上流両方向の伝搬時間を導くことを含む。
これらのステップは、各測定の前に装置の調整を行わずに複数の測定がなされる点を除けば、前述した方法のステップと非常に類似している。
下記の特徴は、各測定前における信号生成フェーズを含む方法および信号生成フェーズを含まない方法の、両方の方法に当てはまる。
更なる実施形態によれば、測定信号の振幅又は応答信号の振幅が、所定数の振動、例えば搬送波の5以上の振動周期にわたって、最大の振幅に高められる。振幅が所定の時間にわたって高められると、例えば圧電トランスデューサである超音波トランスデューサの反応時間における慣性が測定に及ぼす影響を減らすことができる。
1つの特定の実施形態では、測定信号又は応答信号が、搬送波信号の少なくとも5個の振動周期にわたり、最大振幅まで指数関数的に増大する。
更なる実施形態によれば、測定信号は先行部分を含み、この先行部分は測定信号の信号最大値の半値幅の複数倍の時間にわたって延びており、かつ、信号最大値の半値幅の少なくとも1つ分の領域だけ時間が先行している。
さらに別の実施形態によれば、測定信号が先行部分を含み、この先行部分は、受信信号のうち受信信号の信号最大値から時間的に後に続く後続部分から導出されたものである。
上記先行部分は、受信信号の信号最大値周囲において信号最大値の半値幅の少なくとも3倍にわたって延びている。
更なる実施形態によれば、上記先行部分は、測定信号の信号エネルギーの少なくとも10%もしくは少なくとも50%を含む。
信号S(t)のある時間間隔における信号エネルギーEは、下記式1
又は、その離散的表現である下記式2
によって定義できる。ここで、時間間隔は、[T1,T2]又は[−m×Δt,n×Δt]によってそれぞれ与えられる。
測定信号の先行部分は、空間及び時間においてピークとなる信号の生成に顕著に寄与することがある。
一部の特定の実施形態において、測定信号又は応答信号は、例えば1ビットから8ビットの間の分解能で振幅に関してデジタル化された振幅変調振動信号とすることができる。これは演算速度及びメモリスペースの点で有利であり、ひいては信号ピークを高めることにもなる。
更なる実施形態では、トランスデューサに与えられる測定信号が、所定の振幅かあるいは無振幅つまりゼロ振幅を与える0−1変調に従って変調された振動信号であってもよい。
具体的には、振幅変調された測定信号は、流体導管に対して所定の速度を有する流体又は流体導管に対して実質上静止している流体が流体導管に供給される信号生成フェーズに基づいて測定された応答信号から導出してもよい。
インパルス信号が、例えば圧電トランスデューサである第1の超音波トランスデューサに与えられ、そのインパルス信号に対する応答信号が、例えば圧電トランスデューサである第2の超音波トランスデューサで受信される。
測定信号は、応答信号から導出される。測定信号の導出は、応答信号の一部分又は応答信号から導出された信号の一部分を選択すること、及びその選択された信号部分を時間反転することを含み、デジタル化された測定信号は、後で使用するためコンピュータで読み取り可能なメモリに記憶可能である。
1つの特定の実施形態においては、測定信号又は応答信号の包絡線の振幅を、測定信号の先行信号部分から最大の振幅まで少なくとも1桁倍大きくすることができる。先行信号部分は、信号の最大値よりも時間的に先行する部分である。すなわち、先行信号部分は先に送信されたものとなる。更なる実施形態では、振幅が少なくとも2桁倍さらには少なくとも3桁倍増大される。
別の視点によれば、伝搬時間差式超音波流量計において流速を測定する装置が開示される。この装置は、第1の圧電素子と接続する第1コネクタと、第2の圧電素子と接続する第2コネクタと、第1コネクタと接続される任意のデジタル−アナログ変換器(DAC)と、第2コネクタと接続される任意のアナログ−デジタル変換器(ADC)とを備える。
また上記装置は、コンピュータで読み取り可能なメモリと、電子タイマーまたは発振器と、インパルス信号を第1コネクタに送る送信器と、インパルス信号に対する応答信号を第2コネクタから受信する受信器とを備える。
さらに上記装置は、受信した応答信号から測定信号を生成する手段、例えば、受信した応答信号の一部分又は受信した応答信号から導き出した信号の一部分を選択する選択部や、応答信号の上記選択部分を時間反転して反転信号を得る反転部を備える。任意の構成要素として、望ましくない信号成分を取り除く帯域フィルタを設けてもよい。少なくとも上記反転信号から測定信号を導出し、その測定信号をコンピュータで読み取り可能なメモリに記憶する処理部も設けられる。
さらに上記装置は、流速を測定する手段を備える。送信側には、第1コネクタ又は第2コネクタに接続可能な測定信号生成器と、測定信号を第1コネクタに送るDACや両コネクタなどの送信手段とが設けられる。受信側には、測定信号に対応する応答信号を第2コネクタから受信する受信器と、受信した応答信号から流速を導出する速度処理部とが設けられる。尚、本明細書において、「流れの速度(velocity of flow)」、「フロー速度(flow velocity)」、及び「流速(flow speed)」という用語は同義で使われる。
上記装置は、A/D変換器、D/A変換器、コンピュータ読み取り可能メモリを含まないアナログ装置としてもよいが、デジタルコンピュータシステムで上記装置あるいはその各部を構成することも可能である。
すなわち、速度処理部、選択部、及び反転部など各種の信号処理部は、全体的に又は部分的に、特定用途向けの電子部品もしくはコンピュータで読み取り可能な命令セットを有するプログラムメモリとしてもよい。同様に、送信器側の測定信号生成器及びインパルス信号生成器も、全体的に又は部分的に、特定用途向け電子部品でもよいし、もしくはコンピュータで読み取り可能な命令セットを有するプログラムメモリに含めてもよい。
別の実施形態によれば、上記装置は、上記ADCを含むダイレクトデジタル信号シンセサイザ(DDS)を備える。DDSは、周波数制御レジスタ、標準発振器、数値制御発振器、及び再構成ローパスフィルタを備える。さらに、上記ADCは再構成ローパスフィルタを経て第1及び第2コネクタに接続可能である。
特にデジタル信号シンセサイザは、コンピュータで読み取り可能なメモリを含むメモリ部に記憶された所定のアルゴリズム又は所定の値を用いて、測定信号などの信号を合成するように構成可能である。信号は例えば、直接的な信号生成又はDDS(ダイレクトデジタル合成)によって生成可能である。
また本明細書は、第1コネクタに接続された第1の圧電トランスデューサと、第2コネクタに接続された、例えば圧電トランスデューサである第2の超音波トランスデューサとを備えた流量測定装置を開示する。具体的には、例えば圧電トランスデューサである両超音波トランスデューサは、クランプ止め機構など、それらをパイプに取り付けるための取り付け領域を備えてもよい。
さらに本明細書は、パイプ部を含む流量測定装置を開示する。例えば圧電トランスデューサである第1の超音波トランスデューサがパイプ部の第1の位置に取り付けられ、例えば圧電トランスデューサである第2の超音波トランスデューサがパイプ部の第2の位置に取り付けられる。具体的には、両トランスデューサがパイプ部にクランプ止めされてもよい。パイプ部を含む装置を提供することは、装置をパイプ部に関して予め較正しておくことができるため有利となる。
本装置は、コンパクトかつ可搬式に構成できる。クランプオン式トランスデューサなど表面装着可能なトランスデューサを備えた本発明による可搬装置は、アクセス可能などんな箇所にあるパイプをチェックするのにも使える。一般的に、本装置は固定式又は可搬式のどちらでもよい。本装置は、必要な箇所に設置できるよう充分にコンパクトであるとともに、湿度、熱及び腐食性物質などの環境条件に対して充分に保護されることが好ましい。
さらに本明細書は、本明細書による流量測定方法を実行するためのコンピュータで読み取り可能なコード、このコンピュータで読み取り可能なコードを記憶するコンピュータ読み取り可能なメモリ、及び本明細書による方法の各ステップを実行するように動作可能な特定用途向け電子部品を開示する。
具体的には、特定用途向け電子部品は、例えばEPROM、EEPROM、フラッシュメモリなどの上記コンピュータで読み取り可能なメモリを含む電子部品でもよい。他の実施形態によれば、特定用途向け電子部品が、例えば特定用途向け集積回路(ASIC)やフィールドプログラマブルゲートアレイ(FPGA)など、配線接続された部品又は再構成可能な回路を備えた部品を備えてもよい。
別の実施形態において、本明細書による特定用途向け電子部品は、例えば、マルチダイ(チップ)構成で適切にプログラムされたEPROMに接続されているFPGAなど、複数の相互接続された電子部品を備える。特定用途向け電子部品の他の例は、プログラマブル論理アレイ(PLA)やコンプレックスプログラマブル論理装置(GPLD)などのプログラマブル集積回路である。
既製品のテスト装置が、本出願に従って流体導管内における流体の流速を測定しているかどうかを判定することが有用である。そのためには、流体導管に対して所定の速度を有する流体を、流体導管に供給する。テストインパルス信号を、例えばテスト装置の圧電トランスデューサであって、流体導管の第1の位置に取り付けられた第1の超音波トランスデューサに与えた後、そのテストインパルス信号に対するテスト応答信号が、例えば流体導管の第2の位置に取り付けられた圧電トランスデューサである、テスト装置の第2の超音波トランスデューサによって受信される。
次いで、テスト測定信号が上記応答信号から導出される。このテスト測定信号の導出は、信号を時間反転することを含む。その後、テスト測定信号が、例えば各々圧電トランスデューサである第1及び第2の超音波トランスデューサの他方から発信される測定信号と比較される。この測定信号は、メーカーからの供給時に、テスト装置をしばしば管に取り付けた状態で製造した後、工場で1回生成させた測定信号に基づく、テスト装置が有する信号である。
テスト装置が本出願に従って流体導管内における流体の流速を求める方法を使用している場合、上記テスト測定信号と上記測定信号は同じである。言い換えれば、テスト測定信号と測定信号が同じになるまで、テスト信号を選択し、本出願の信号生成フェーズを繰り返すことで、本出願の主題が逆行分析(リバースエンジニアリング)される。ここで「同じ」という用語は、テスト測定信号と測定信号との間に有意な相関性が存在することを意味する。
本方法は、テスト応答信号の一部分又はテスト応答信号から導出した信号の一部分を選択すること、及び後で使用するためテスト測定信号を記憶することを含んでもよい。
従って、伝搬時間差式超音波流量計において流速を測定する装置は、機能的特徴によって表現すれば、第1の圧電素子用の第1コネクタと、第2の圧電素子用の第2コネクタと、インパルス信号を前記第1コネクタに送る送信器と、前記インパルス信号に対する応答信号を、前記第2コネクタから受信する受信器と、前記応答信号を時間反転して反転信号を得る反転部と、前記反転信号から測定信号を導出する処理部とを備える。流体導管における流体の流速を求める上記装置を使用する際は、流体導管に対して速度を有する流体を、流体導管に供給する。その後、例えば各々圧電トランスデューサである第1及び第2の超音波トランスデューサの一方に測定信号を与え、例えば各々圧電トランスデューサである第1及び第2の超音波トランスデューサの他方によって、上記測定信号に対する第1の応答信号を受信する。その後、第1の応答信号から流体の流速を導出することができる。本装置の逆行分析は、テストインパルス信号を、例えばテスト装置の圧電トランスデューサである第1の超音波トランスデューサに与え、そのテストインパルス信号に対するテスト応答信号を、例えば流体導管の第2の位置に取り付けられた圧電トランスデューサである、テスト装置の第2の超音波トランスデューサで受信し、更に信号を時間反転することを含め上記応答信号からテスト測定信号を導出した際、そのテスト測定信号と例えば各々圧電トランスデューサである第1及び第2の超音波トランスデューサの他方から発せられる測定信号とが同じであることを示す。かかる機能的記述は、本出願の装置を、発信される信号の構造及び形状を説明せずに特徴付けるのに有効である。
本装置が、前記第1コネクタに接続されたD/A変換器と、前記第2コネクタに接続されたA/D変換器と、コンピュータ読み取り可能メモリとを備えてもよいことは明らかである。本装置はさらに、前記受信応答信号の一部分又は前記受信応答信号から導出された信号の一部分を選択する選択部を備えてもよい。その場合、前述の評価は、前記受信応答信号の選択された一部分又は前記受信応答信号から導出された信号の選択された一部分について行われる。
次に、本明細書の主題を、以下の図面を参照してさらに詳しく説明する。
図1は、2つの圧電素子を備えた第1の流量計装置を示す。 図2は、図1の流量計装置における1つのダイレクト信号を示す。 図3は、流れの方向から見た図1の流量計装置を示す。 図4は、4つの圧電素子と4つのダイレクト信号を備えた第2の流量計装置を示す。 図5は、流れの方向から見た図4の流量計装置を示す。 図6は、テスト信号の概略図を示す。 図7は、テスト信号の応答の概略図を示す。 図8は、反転信号の概略図を示す。 図9は、反転信号の応答の概略図を示す。 図10は、高分解能の第1反転信号を示す。 図11は、図10の反転信号の応答を示す。 図12は、高分解能の別の反転信号を示す。 図13は、図12の反転信号の応答を示す。 図14は、高分解能の別の反転信号を示す。 図15は、図14の反転信号の応答を示す。 図16は、高分解能の別の反転信号を示す。 図17は、図16の反転信号の応答を示す。 図18は、高分解能の別の反転信号を示す。 図19は、図18の反転信号の応答を示す。 図20は、高分解能の別の反転信号を示す。 図21は、図20の反転信号の応答を示す。 図22は、高分解能の別の反転信号を示す。 図23は、図22の反転信号の応答を示す。 図24は、高分解能の別の反転信号を示す。 図25は、図24の反転信号の応答を示す。 図26は、高分解能の別の反転信号を示す。 図27は、図26の反転信号の応答を示す。 図28は、12ビット分解能の別の反転信号を示す。 図29は、図28の信号の応答を示す。 図30は、3ビット分解能の別の反転信号を示す。 図31は、図30の信号の応答を示す。 図32は、2ビット分解能の別の反転信号を示す。 図33は、図32の信号の応答を示す。 図34は、1ビット分解能の別の反転信号を示す。 図35は、図34の信号の応答を示す。 図36は、図1の流量計の圧電素子における短いインパルスを示す。 図37は、図36の信号の反転応答から導出した、図1の流量計における圧電素子の信号を示す。 図38は、図37の信号の応答を示す。 図39は、上流側及び下流側の相互相関関数を示す。 図40は、図39の拡大断面を示す。 図41は、本明細書に基づき流速を測定する装置を示す。 図42は、図41の装置で使われるダイレクトデジタルシンセサイザを示す。 図43は、第1のマルチトランスデューサ配置を示す。 図44は、第2のマルチトランスデューサ配置を示す。 図45は、クランプオン式トランスデューサのZ形構成を示す。 図46は、クランプオン式トランスデューサのV形構成を示す。 図47は、クランプオン式トランスデューサのW形構成を示す。 図48は、1サイクルの送信信号を示す。 図49は、10サイクルの送信信号を示す。 図50は、TRA送信信号を示す。 図51は、図48の1サイクルの送信信号に対する応答信号を示す。 図52は、図49の10サイクルの送信信号に対する応答信号を示す。 図53は、図50のTRA送信信号に対する応答信号を示す。 図54は、TRA送信信号の圧力曲線とTRA送信信号に対する応答信号を示す。 図55は、TRA送信信号の圧力曲線とTRA送信信号に対する応答信号を示す。 図56は、図55の信号入力を生成するのに用いられるインパルス信号を示す。 図57は、チャネル特性を示す第1の応答信号を示す。 図58は、チャネル特性を表す第2の応答信号を示す。 図59は、別の応答信号を示す。 図60は、別の応答信号を示す。
以下の説明で、本明細書の実施形態を詳述する。但し、以下の詳細に拠らずとも各種の実施形態を実施可能なことは、当業者にとって自明である。
図1は、第1の流量計装置10を示す。本流量計装置においては、第1の圧電素子11が、管12とも称されるパイプ12の外壁に配置されている。第2の圧電素子13がパイプ12の反対側に配置され、第1の圧電素子11と下流側の第2の圧電素子13との間の直線が、パイプ12の対称軸方向でもある平均流の方向14に対して角度βを成している。角度βは図1の例では約45度に選定されているが、その傾斜角度は、例えば60度などもっと急でも、あるいは例えば30度などもっと緩やかでもよい。
図1に示した圧電素子11、13等の圧電素子は一般的に、音響送信器及び音響センサとして作動可能である。音響送信器及び音響センサは、同じ圧電素子で構成してもよいし、同じ圧電素子の異なる領域によって構成してもよい。その場合、圧電素子又はトランスデューサは、送信器又は音響源として作動されるとき圧電送信器とも称され、また音響センサとして作動されるとき音響センサ又は受信器とも称される。
流れ方向が図1に示す向きの場合、第1の圧電素子11は「上流側」圧電素子とも称され、第2の圧電素子13は「下流側」圧電素子とも称される。本明細書による流量計は、両方向の流れに対して実質上同様に作動し、図1の流れ方向は例示としてのみ示してある。
図1は、上流側圧電素子11が圧電トランスデューサとして作動され、下流側圧電素子13が音響センサとして作動される構成の場合の、図1における電気信号の流れも示している。分かり易くするため、本出願は上流側と下流側で同様に作動する。すなわち、両圧電素子の位置は交換可能である。
第1の演算部15が上流側圧電素子11に接続され、第2の演算部16が下流側圧電素子13に接続されている。第1の演算部15は、第1のデジタル信号プロセッサ、第1のデジタル-アナログ変換器(DAC)、及び第1のアナログ-デジタル変換器(ADC)を備える。同じく、第2の演算部16は、第2のデジタル信号プロセッサ、第2のデジタルアナログ変換器(DAC)、及び第2のアナログデジタル変換器(ADC)を備える。第1の演算部15は、第2の演算部16に接続されている。
図1に示した2つの演算部15、16を備えた構成は、例示としてのみ示したものである。他の実施形態として、図示と異なる数及び構成の複数の演算部を備えてもよい。例えば、1つの中央演算部だけでもよく、2つの AD/DC変換器と1つの中央演算部でもよく、もしくは各トランスデューサに設けた2つの小型演算部と1つの比較的大型の演算部でもよい。
1つ又は複数の演算部は、例えばマイクロコントローラ、特定用途向け集積回路(ASIC)、ACID、フィールドプログラマブルゲートアレイ(FPGA)によって得られる。特に、記憶されたデジタル信号からの電気信号の合成は、デジタル-アナログ変換器(DA,DAC)からなるダイレクトデジタルシンセサイザ(DDS)で行ってもよい。
本明細書に基づき測定信号を生成する方法は、以下のステップを含む。
第1の演算部15のデジタル信号プロセッサで音響信号を合成することにより、所定のデジタルテスト信号が生成される。このデジタル信号は、信号路17を介して、第1の演算部15から圧電素子11へと送られる。圧電素子11は、対応する超音波テスト信号を生成する。両演算部15と16を、1つの単一ユニットとして設けることもできる。
テスト信号は、例えば単一の1MHz振動または10個の同様の振動により、短いパルスとして与えられる。具体的には、テスト信号は一定の振幅の少数の振動によって与えることで、矩形信号に近似させてもよい。単一又は複数の振動は、サイン波状、三角状、矩形状、又はその他の形状とし得る。
超音波テスト信号は、パイプ12内の液体を通り圧電素子13まで伝搬する。図1において、超音波信号のダイレクト信号路が、矢印18で示してある。同じく、逆方向における超音波信号のダイレクト信号路が、矢印19で示してある。応答信号は、圧電素子13によって検出され、信号路20を介して第2の演算部16に送られ、第2の演算部16によってデジタル化される。
別のステップでは、デジタル測定信号がデジタル化された応答信号から導出される。この測定信号の導出は、デジタル化応答信号の時間反転により行われる。別の実施形態によれば、上記測定信号の導出は、振幅範囲における低分解能への変換や、低周波ノイズ及び高周波ノイズなどのノイズを除去する、信号の帯域幅フィルタリングなどの更なるステップを含む。特に、帯域幅フィルタリングのステップは、信号を時間反転するステップの前に行ってもよい。
信号の反転は、例えば、メモリ領域を逆方向に読み出したり、あるいはフーリエ表記のサイン成分の符号を逆にするなど、各種の方法で実施し得る。
1の実施形態では、デジタル化応答信号のうち、ダイレクト信号からの応答を含む適切な一部分が選択される。次いで、その応答信号の一部分が、時間転回もしくは時間反転される。言い換えれば、応答信号のうち後で受信された信号部分が、反転測定信号においては先に送出される。例えば、信号が時系列の振幅サンプルによって表わされている場合、上記の信号反転では、振幅信号の順序が反転もしくは逆転される。
時間の方向もしくは符号が反転されて得られた信号は、「反転信号」とも称される。ここでの「反転」という表現は、時間の方向に関する反転を意味し、振幅値など、値に関する反転を意味するものではない。
図10−19は、本明細書によるデジタル信号を例示している。
本明細書の一実施形態による流量計においては、上流及び下流方向の両方向18、19で同じ測定信号が使われ、単純且つ効率的な構成が得られる。他の実施形態によれば、異なる信号が両方向で使われる。具体的には、測定信号がテスト信号の本来の受信器に与えられてもよい。このような構成は、非対称な条件及びパイプ形状の場合に有利である。
上記の反転信号を測定信号として用いて、パイプを流れる液体の流速を測定する方法は、以下のステップを含む。
上記の測定信号が、第1の演算部15から圧電トランスデューサ11へ信号路17を介して送られる。圧電トランスデューサ11は、対応する超音波測定信号を生成する。このような測定信号の各例が、図10、12、14,16,18,20,22,24,26,28,30,32,34,37及び38に示してある。
超音波測定信号は、パイプ12内の液体を通って圧電センサ13まで伝搬する。応答信号が圧電センサ13によって検出され、信号路20に沿って第2の演算部16に送られ、第2の演算部16によりデジタル化される。
第2の演算部16は、デジタル化応答信号を第1の演算部15に送る。第1の演算部15は、例えば以下に詳しく説明する方法の1つを用い、受信信号の伝搬時間を求める。
同様の処理が、逆方向に伝搬する信号に対しても行われる。つまり、上記の測定信号が下流側圧電トランスデューサ13に与えられ、応答信号が上流側圧電トランスデューサ11により測定されて、逆方向における上流伝搬時間TOFUPが得られる。第1の演算部15は、例えば次式3に基づいて流れの速度を求める。
ここで、Lは両圧電素子11と13間のダイレクト路の長さ、βは両圧電素子11と13間のダイレクト路と平均流の方向との間の傾斜角、cは所定の圧力及び温度条件下における液体中の音速である。
音速の二乗cは、次式4よって2次近似可能である。
この式から、次式5が得られる。
これにより、温度や圧力を求めて、更に流体密度及び音速を求める必要がなくなり、もしくは音速や流体密度を直接測定する必要もなくなる。一方、1次の誤差は、1つの測定方向だけでは相殺されない。
因子 (2・L・cosβ)を用いる代わりに、既知の流速での較正測定から比例定数を得ることもできる。較正の比例定数では、流量分布、及び散乱し直線に沿って伝搬しない音波による寄与など別の影響も考慮する。
別の実施形態によれば、インパルス信号を生成し、応答信号を記録し、応答信号から反転した測定信号を導出する処理が、コンピュータでシミュレートされる。パイプ12の管径や各センサの位置など、関連のパラメータがシミュレーションのための入力パラメータとして与えられる。
更に別の実施形態によれば、送信側圧電素子に供給される測定信号が、図37及び図38に示す信号形状のような、インパルス信号に対する典型的な応答信号の形状を用いて合成される。例えば、1MHzのサイン波振動を、半値幅10msecのガウス確率関数に基づき包絡線で振幅変調して測定信号を得てもよい。半値幅は、管径及びセンサ位置など実際の構成に応じた入力パラメータとして選択し得る。
本明細書による流量計は、特に流量計がパイプ部分と一緒に提供される場合、工場でのテスト実施時に測定信号を生成させた事前設定した流量計としてもよい。
本明細書の簡単な実施形態では、上流及び下流方向の伝搬時間が、受信信号のピーク振幅の時間を測定信号の送信時間に対して比較評価することで求められる。より高い精度を得るため、受信信号の包絡線を用いて最大値を求めてもよい。別の実施形態では、測定を複数回繰り返し、伝搬の平均時間を用いる。
本明細書の別の実施形態によれば、信号の伝搬時間が、相互相関技法を用いて評価される。すなわち、各時間シフトは次式6に基づき、受信された上流又は下流信号を流速0における受信信号と相互相関させることによって評価できる。
ここで、SigFlowは、流体がパイプを流れる場合の測定条件下における上流又は下流信号を表し、SigNoFlowは、流速が0の場合の較正条件下における信号を表す。無限和の限界範囲は、充分に長い時間窓 [−T1,+T2]を表す。より一般的に言えば、−T1と+T2は同じである必要はなく、この点は流量計にとって実用上有利である。
時間シフト(TOFUP-TOFDown) は、上流側相関関数の最大値が得られた時刻と、下流側相関関数の最大値が得られた時刻とを比較することによって得られる。最大値の位置をより正確に求めるため、相関関数の包絡線を用いてもよい。
別の実施形態においては、第1の演算部15と第2の演算部16との間に別個の評価装置が設けられ、各信号到着時間及び流速の計算を行う。
一般に、音響センサの測定信号は、散乱した信号とダイレクト信号が重畳して得られる。散乱信号は、パイプ壁で1回又は複数回散乱される。この点が、図2及び図3に例示されている。
図1のトランスデューサ構成は、直線もしくはZ形構成である。V形及びW形構成など、パイプの反対側での反射を利用したその他の構成も可能である。V形及びW形構成は、パイプ壁での反射に基づいて動作し、Z形よりも多く散乱を生じる。本出願の主題は、経路が適切に考慮されている限り、上記のような構成でも効果が得られる。
V形構成の場合、2つのトランスデューサが、パイプの同じ側に取り付けられる。45度の反射を記録するため、2つのトランスデューサは、流れの方向に沿ってほぼ管径の距離だけ離れて配置される。W形構成では、3回の反射を利用する。V形構成の場合と同様、2つのトランスデューサがパイプの同じ側に取り付けられる。2回の45度の反射後の信号を記録するため、2つのトランスデューサは、流れの方向に沿って管径の2倍の距離だけ離れて配置される。
図2は、圧電素子11から圧電素子13へダイレクトに伝搬する第1の音響信号を例示している。
単純化のため、図2-5では散乱事象を反射として示しているが、実際の散乱過程はもっと複雑である。具体的には、最も関連のある散乱は通常、パイプ壁や各圧電トランスデューサの前方に配置された素材において生じる。受信される散乱も、センサの配置に依存する。図45,図46及び図47は、それぞれZ、V及びW形のセンサ配置を例示している。図3は、図1を流れ方向に沿って矢視方向A−Aから見た図である。
図4及び図5は、別の圧電素子22が圧電素子11に対して45度の角度に配置され、さらに別の圧電素子23が圧電素子13に対して45度の角度に配置された第2のセンサ配置を示す。
また図4及び図5は、圧電素子11、22が圧電トランスデューサとして作動され、圧電素子13,23が音響センサとして作動される場合におけるダイレクト(つまり直線の)音響信号路を示している。図4中パイプ12の後ろ側にある圧電素子23は、図4において点線で示されている。
図6-9は、テスト信号の応答から測定信号を生成する方法を単純化して示す。図6-9において、散乱による損失が、信号のうち斜線を施した部分と矢印で示してある。
図6-9においては、音響信号が、直線路、遅延時間Δt後に第1の散乱路、遅延時間2Δt後に第2の散乱路を沿って伝搬するものと仮定する。経路における信号の減衰は、考慮しない。
矩形状スパイクの形のテスト信号が、圧電素子11に与えられる。散乱により、信号振幅の第1部分が第1の散乱路で失われて時間Δt後に現れ、信号振幅の第2部分が第2の散乱路で失われて時間2Δt後に現れる。その結果、図7に白い柱で示す信号が生じ、この信号が圧電素子13で記録される。
信号プロセッサが上記の記録信号を時間反転し、反転信号を圧電素子11に与える。ここで、上記と同じ散乱過程が、3つの信号成分全てに加えられる。その結果、図9に示す信号が圧電素子13で記録され、この信号はほぼ対称形である。
実際には、受信信号は時間軸上で分布しており、パイプの素材中を伝搬し、ダイレクト信号よりも早く到着する「表面波」もしばしば存在する。この表面波は、反転測定信号を生成するための時間窓を、適切に選択することで排除される。同じように、多重反射によって生じ遅く到着する信号も、時間窓を限定する及び/又は信号の特定部分を選択することによって排除できる。
以下の表は、ダイレクト配置の場合、つまりDN250パイプの長さ方向に直角な平面においてDN250パイプに配置されたクランプオン式圧電素子間を直線でつないだ場合の測定信号の遅延を示す。流量は、DN250パイプを通る水の流量を表す。
表中、「TOF 1サイクル」は、圧電素子で発生された図36に示すようなインパルスが、1μmの周期を有する1発振の電気信号により励起されたものを示す。「TOF 10サイクル」は、圧電素子で発生された信号が、1μmの周期を有する定振幅の10個のサイン波発振の電気信号により励起されたものを示す。
図10-27は、高分解能の反転信号とそれぞれの応答信号を示す。マイクロ秒単位の時間に対して、電圧が任意単位で示してある。
各対のうち上方の図(図10-27のうちの偶数の図)の時間軸は、反転信号の送信時間を示す。送信時間は、受信信号を記録するのに用いられる時間窓に限定される。図10-27の例において、時間窓はダイレクト信号から得られる最大値の始まりよりもわずか前にスタートし、その後100マイクロ秒で終了する。
下方の図(図10-27のうちの奇数の図)の時間軸は、応答信号の最大値付近を中心として、反転信号用の時間窓のサイズである100マイクロ秒の範囲で、応答信号の最大値の前後に延びている。
図28-35は、高分解能の反転デジタル信号で、振幅範囲が12,3,2及び1ビットの分解能の信号と、それぞれの応答信号を示す。マイクロ秒単位の時間に対し、電圧はボルトで示してある。図28〜35の各信号は、水を満たしたDN250パイプについて得られたものである。
反転信号用の時間窓の長さは、450マイクロ秒である。従って、図28-35の時間窓は、前述した図8〜27の時間窓の4倍以上である。
図28〜35から分かるように、1ビット分解能のデジタル化でも、シャープなスパイクが得られる。また、分解能が低いほど、スパイクがよりシャープになることも分かる。この効果は恐らく、図28〜35の例においては、応答信号を時間的に集中させたまま、振幅範囲のデジタル化を次第に粗くしたことにより、入力信号の合計エネルギーが増加したためと思われる。
図36は、3.57MHzの周波数と等価な約0.56秒間持続した電気信号を受信した後、圧電素子によって発生された信号を示す。圧電素子の慣性のため、負の電圧における最大振幅は正の電圧における最大振幅よりも小さく、圧電素子が静止状態になる前に、複数の残響が生じている。
図37は、図1の上流側圧電素子11などの圧電素子に与えられる電気信号を示す。図37の信号は、図36に示したような信号に対する10個のデジタル化応答信号の平均を形成し、その信号を時間反転することによって得られ、各応答信号は図1の下流側圧電素子13などの圧電素子によって受信される。
図37の例では、応答信号のうち、応答信号の包絡線の開始より約10マイクロ秒前に始まり、応答信号の包絡線より約55マイクロ秒後に終了する時間部分を切り出すことによって、デジタル化信号が得られる。図37の応答信号の包絡線形状は、ガウス確率分布の形状、すなわち exp(−x)を適切にシフト及び倍率調整した形状に類似する。
図38は、図37に示した信号に対する応答信号の一部を示し、図37の信号が上流側圧電素子11などの第1の圧電素子に与えられ、図1の下流側圧電素子13などの第2の圧電素子で受信された場合に相当する。
図39は、図1の配置における上流側信号及び下流側信号を、流量0で得られた信号とそれぞれ相互相関させて得られた上流側相互相関関数と下流側相互相関関数を示す。
図40は、図39の拡大断面を示す。2つの位置マーカが、上流側及び下流側相互相関関数それぞれの最大値の位置を示している。2つの最大値間の時間差が、上流側及び下流側信号間の時間差を示す尺度である。
図48,図49及び図50は3つの異なる送信信号を示す。図48は通常のパルス(1サイクル)を示し、図49は、上記のように生成される図50の信号等の測定信号との比較のための、10サイクルのパルスを示す。各トランスデューサは、DN250パイプにクランプ止めされている。
図51,図52及び図53は、それぞれ図48,図49及び図50に示した各信号を送信後に受信された対応する信号を示す。これらの比較から、測定信号はエネルギーを集中させ、図48及び図49の通常の信号(例えば1又は10サイクル)に応答した受信信号と比べ、振幅が2倍以上大きい受信信号を発生させることが容易にわかる。
図41は、図1に示した構成あるいは本明細書によるその他の構成において、流れの測定を行う流量測定装置60を例示する。図1の構成において、フロー測定装置60は、第1及び第2の演算部15、16を備える。
流量測定装置60は、第1の圧電素子を接続する第1コネクタ61と、第2の圧電素子を接続する第2コネクタ62とを有する。第1コネクタ61は、マルチプレクサ63を経てデジタル−アナログ変換器(DAC)64に接続される。第2コネクタ62は、デマルチプレクサ66を経てアナログ−デジタル変換器(ADC)65に接続される。
ADC65は信号選択部67に接続され、信号選択部67は信号反転部68に接続され、信号反転部68は帯域フィルタ69に接続され、帯域フィルタ69はコンピュータ読み取り可能メモリ70に接続されている。さらに、ADC65は速度演算部71にも接続されている。
DAC64は、インパルス信号生成器72と測定信号生成器73に接続される。測定信号生成器73は、指令ライン74を経てインパルス信号生成器72に接続される。速度演算部71は、別の指令ライン75を経て測定信号生成器73に接続される。
一般に、インパルス信号生成器72と測定信号生成器73は、発振器などのハードウェア要素と、インパルス信号生成モジュール及び測定信号生成モジュールなど複数のソフトウェア要素を備える。この場合、指令ライン74,75は各モジュール間のソフトウェアインタフェースにより提供されてもよい。
信号生成フェーズにおいて、インパルス信号生成器が信号をDAC64に送り、選択部67が対応する入力信号をADC65を経て受信し、入力信号の一部分を選択する。反転部68が選択された信号部分を時間反転し、オプションの帯域フィルタ69が範囲外の低周波及び高周波をフィルターして除去し、得られた測定信号がコンピュータメモリ70に記憶される。信号操作のステップに関連して「信号」という用語が使われる場合、その信号は特にコンピュータメモリ中における信号を示す。
具体的には、信号の表記は、デジタル化された振幅と、それに関連した離散的な時間との、一対の値によって定義が可能である。その他の表記としては、特にフーリエ係数、ウェーブレット係数、及び信号を振幅変調する際の包絡線がある。
図42は、図1に示した構成あるいは本明細書によるその他の構成において、流れの測定を行う流量測定装置60’の第2の実施形態を示す。流量測定装置60’は、ダイレクトデジタルシンセサイザ(DDS)76を備える。単純化のため、DDS76の構成要素だけが示してある。DDS76は、任意波形生成器(AWG)とも称される。
DDS76は標準発振器77を備え、標準発振器77は周波数制御レジスタ78、数値制御式発振器(NCO)79、及びDAC64に接続される。Nチャネル用NOC79の入力が、周波数制御レジスタ78の出力に接続される。Mチャネル用DAC64の入力がNOC79に接続され、再構成ローパスフィルタ80の入力が DAC64に接続される。一例として、100MHzのクロック周波数を有するダイレクト数値制御式発振器79を用いて、振幅変調された1MHzの信号を発生してもよい。
再構成ローパスフィルタ80の出力が、図1の圧電素子11,13に接続される。
例えばダイレクトデジタルシンセサイザを用いて所定の振幅変調信号を得るためには、発振器結晶の慣性のため、搬送波の周波数より高い周波数の発振器を使用することが有利な場合が多い。
図45,図46及び図47は、それぞれ前述したZ、V及びW形の流量測定構成を示す。図45,図46及び図47の例では、クランプオン式トランスデューサがそれぞれの連結具を介して、導管に取り付けられている。
図54及び図55は、時間反転手順を用いなかった場合と時間反転手順を用いた場合に発生された各送信信号に対する、それぞれの受信つまり応答信号を比較して示す。
図54の例では、ガウス型の包絡線を有する変調サイン波が、送信信号として使われている。送信信号の信号エネルギーは1.3×10−7(Pa/m)sに比例し、信号の振幅は0.1Paである。この値は、単位長さ当たりの圧力の二乗を時間積分することによって得られる。応答信号は、極大値から極大値までの振幅が約0.09Paとなっている。
図55の例では、図56のインパルス信号に対する応答信号から得られた時間反転信号が送信信号として使われている。この送信信号は、図54の送信信号と同じ1.3×10−7(Pa/m)sの信号エネルギーを有するように調整されている。これにより、極大値から極大値までの振幅が約0.375Paである受信信号が得られる。
図55の受信信号の振幅は、図54の受信信号の振幅の4倍以上である。受信側における振幅の増大は、より容易で安定な信号認識を可能とする。振幅の増大は、時間反転信号の振幅のビット分解能を調整することで、特により大きい振幅を得るようにビット分解能を増減することで、調整可能である。
図56及び図57は、伝送路に関する情報、特に導管の壁厚や壁への堆積物に関する情報を得るのに、受信信号をどのように使用できるかを示している。本明細書によれば、時間反転信号である測定信号に対する応答を分析して、クラックや付着物などのような、パイプ素材の性状変化を求めることができる。本明細書の一実施形態による流量測定においては、伝搬時間の測定に使われたのと同じ受信信号を分析することで、上記のような性状変化が求められる。
図57は、第1の伝送路に関する情報を含む第1の応答信号を示す。
図58は、第2の伝送路に関する情報を含む第2の応答信号を示す。中央の主ローブ(main lobe)に付した水平の矢印の長さは、主ローブに対してそれぞれ左右にある左側ローブと右側ローブとの間に延びている。図46に従って信号が発生された場合、上記矢印の長さがパイプ壁の厚さを表す。測定される壁厚は、図46におけるパイプの下方部で音波が反射された箇所で求められたものである。パイプ壁に堆積物があると、測定壁厚が大きくなる。
図59は別の応答信号を示す。図59の信号を得るための試験設定は、複数のクランプオン式アングルトランスデューサ、アクリル製のトランスデューサ連結ヘッド、音速 c=2370m/s、連結角度40°、ステンレス鋼製の壁、横波速度c=3230m/s、61.17°、流体は水、流体中の音速c=1480m/s、横波の軸角度23.67°、及びフロー角66.33°であり、図59から得られた。
図60は別の応答信号を示す。図60の信号を得るための試験設定は、アクリル製のトランスデューサ連結ヘッド、音速c=2370m/s、連結角度20°、ステンレス鋼製の壁、縦波速度c=5790m/s、56.68°、横波速度c=3230m/s、流体は水、流体中の音速c=1480m/s、縦波の軸角度12.33°、横波の軸角度12.33°、及びフロー角66.33°であり、図60から得られた。
図59及び図60の代替の試験構成は、図45,図46及び図47に示してある。
本明細書の一実施形態によれば、図57−60に示した信号などの受信信号を分析することで、チャネル特性が導出される。
図59及び図60の例は、パイプ素材中における縦波及び横波の存在に応じた受信信号の差を示している。それら両波の存在は、選択される素材及び幾何学形状にとって一般的であり、素材の分析に使用可能である。こうした超音波テスト波に基づく素材分析は、非破壊的テスト(NDT)の応用分野で使われる。受信信号が伝送路と素材環境を含む測定系のインパルス応答を含んでいるため、本明細書は流量及び配管素材など両方の同時分析を可能とすることができる。
受信信号の分析は、例えば壁厚を求めるために、受信信号を先に受信したインパルス応答と比較したり、インパルス応答を直接評価したりなど、各種の方法で行える。
上記の説明は個別の限定を多く含んでいるが、それらは実施形態の範囲を制限するものでなく、予見可能な実施形態を例示するに過ぎないものと解されるべきである。方法の各ステップは説明した実施形態における順序と異なる順序で行われてもよく、また各測定装置の複数の処理ユニットへの再分割や、それらの間の相互接続は、説明した実施形態と異なってもよい。
すなわち、上記方法のステップ中、信号のデジタル表記を記憶するステップ、信号部分の選択などの操作を行うステップ、信号を時間反転するステップ、及び信号をフィルターするステップは、相互に入れ替えてもよい。例えば、信号を時間反転した形で記憶したり、あるいは信号を逆の順序に読み出して時間反転信号を得てもよい。
本発明を円形の DN250パイプに関連して説明したが、本発明はその他のパイプサイズやその他のパイプ形状にも適用可能である。また、実施形態をクランプオン式トランスデューサに関連して説明したが、パイプ内に突出させたり、開放路内に設置される湿式トランスデューサにも用い得る。
特に、各実施形態の上述した利点は、実施形態の範囲を制限するものでなく、実施形態を実施した場合に考えられる利点を説明したにすぎないものであることが解されるべきである。つまり、実施形態の範囲は、上述した各例によってではなく、請求の範囲の記載及びそれらの均等物によって定まる。
本明細書の実施形態は、各要素を項目別にまとめた下記のリストとして記述することもできる。項目リストに開示される特徴の各組み合わせは、それぞれ独立した主題と見なされ、本出願のその他の特徴と組み合わせることもできる。
1.流体導管内の流体の流速を求める方法であって、
−前記流体導管に対して所定の速度を有する流体を、前記流体導管に供給すること、
−インパルス信号を、前記流体導管の第1の箇所に取り付けられた第1の超音波トランスデューサに与えること、
−前記インパルス信号に対する応答信号を、前記流体導管の第2の箇所に位置している第2の超音波トランスデューサで受信すること、
−前記応答信号から測定信号を導出すること、且つ該測定信号の導出は応答信号の一部分又は応答信号から得た信号の一部分を選択すること、及び前記選択された信号部分を時間反転することを含むこと、
−前記測定信号を後で使用するために記憶すること、
−前記流体導管に対して移動する前記流体を、前記流体導管に供給すること、
−前記測定信号を前記第1及び第2の超音波トランスデューサの一方に与えること、
−前記測定信号に対する第1の応答信号を、前記第1及び第2の超音波トランスデューサの他方で測定すること、及び
−前記第1の応答信号から前記流体の流速を導出すること、
を含む方法。なお、上記方法のうち、次のステップは、随意に選択され、前記測定信号が予め定められている場合には省略可能である。
−インパルス信号を、前記流体導管の第1の箇所に取り付けられた第1の超音波トランスデューサに与えること、
−前記インパルス信号に対する応答信号を、前記流体導管の第2の箇所に位置している第2の超音波トランスデューサで受信すること、
−前記応答信号から測定信号を導出すること、且つ該測定信号の導出は応答信号の一部分又は応答信号から得た信号の一部分を選択すること、及び前記選択された信号部分を時間反転することを含むこと、及び
−前記測定信号を後で使用するために記憶すること。
2.項目1に記載の方法において、
−前記測定信号を与えるステップと、前記応答信号を測定するステップとを逆方向に繰り返し、第2の応答信号を得ること、及び
−前記第1の応答信号及び前記第2の応答信号から前記流体の流速を導出することを含む方法。
3.項目1又は2に記載の方法において、前記測定信号を得るのに使われる前記信号部分は、前記応答信号の最大振幅周囲の第1部分と、それに後続する信号部分とを含み、該後続信号部分は、前記最大振幅の到着時間より遅い時間に属している方法。
4.いずれかの前記項目に記載の方法において、
−インパルス信号を与えるステップと、対応する応答信号を受信するステップとを複数回繰り返し、複数の応答信号を得ること、及び
−前記受信された複数の応答信号の平均から測定信号を導出すること、
を含む方法。
5.いずれかの前記項目に記載の方法において、前記測定信号の導出は、前記応答信号又は前記応答信号から導出された信号の振幅をデジタル化することを含む方法。
6.項目5に記載の方法において、前記デジタル化信号のビット分解能を高めて、前記測定信号に対する応答信号の振幅を増大させることを含む方法。
7.項目5に記載の方法において、前記デジタル化信号のビット分解能を減少させて、前記測定信号に対する応答信号の振幅を増大させることを含む方法。
8.前記項目5−7の1つに記載の方法において、前記振幅に関するデジタル化信号の分解能が低ビット分解能である方法。
9.前記各項目の1つに記載の方法において、複数の応答信号のうち少なくとも1つを処理することで、導管の壁厚の変化を求め、あるいは、縦波及び横波の音波特性を測定して前記導管壁の素材特性を求めることを含む方法。
10.伝搬時間差式超音波流量計で流速を測定する装置であって、
−第1の超音波素子用の第1コネクタ、
−第2の超音波素子用の第2コネクタ、
−インパルス信号を前記第1コネクタに送る送信器、
−前記インパルス信号に対する応答信号を、前記第2コネクタから受信する受信器、
−前記応答信号を時間反転し、反転信号を得る反転部、及び
−前記反転信号から測定信号を導出し、該測定信号を記憶する処理部、
を備えており、これらの構成要素のうち、下記の構成は、随意に選択され、前記測定信号が予め定められすぐに利用できる場合には省略可能である。
−インパルス信号を前記第1コネクタに送る送信器、
−前記インパルス信号に対する応答信号を、前記第2コネクタから受信する受信器、
−前記応答信号を時間に関して反転し、反転信号を得る反転部、及び
−前記反転信号から測定信号を導出し、該測定信号を記憶する処理部。
11.項目10に記載の装置において、
−前記第1コネクタに接続されたD/A変換器、
−前記第2コネクタに接続されたA/D変換器、及び
−前記測定信号を記憶するコンピュータで読み取り可能なメモリ、
を更に備えた装置。
12.項目10又は11に記載の装置において、前記受信した応答信号の一部分又は前記受信した応答信号から導出された信号の一部分を選択する選択部を更に備え、前記反転部は、前記応答信号の選択部分を時間反転し、前記反転信号を得る装置。
13.項目10−12のいずれかに記載の装置において、
−前記第1コネクタ又は前記第2コネクタに接続可能な測定信号生成器、
−前記測定信号を前記第1コネクタに送信する送信手段、
−前記測定信号に対する応答信号を前記第2コネクタから受信する受信器、及び
−前記受信された応答信号から流速を導出する第2の処理部、
を更に備えた装置。
14.項目10−13のいずれかに記載の装置において、
−前記A/D変換器を備えたダイレクトデジタル信号シンセサイザ、及び
−周波数制御レジスタ、標準発振器、数値制御式発振器、及び再構成ローパスフィルタを備え、前記A/D変換器は再構成ローパスフィルタを経て前記第1及び第2コネクタに接続可能である装置。
15.項目10−14のいずれかに記載の装置において、
−前記第1コネクタに接続された第1の超音波トランスデューサ、及び
−前記第2コネクタに接続された第2の超音波トランスデューサ、
を更に備えた装置。
16.項目10−15のいずれかに記載の装置において、パイプの一部分を備え、前記第1の超音波トランスデューサが前記パイプの一部分の第1の箇所に取り付けられており、
前記第2の超音波トランスデューサが前記パイプの一部分の第2の箇所に取り付けられている装置。
17.項目1−9のいずれかに記載の方法を実行する複数のコンピュータで読み取り可能な命令を含むコンピュータ読み取り可能プログラムコード。
18.項目17に記載のコンピュータ読み取り可能プログラムコードを含むコンピュータ読み取り可能メモリ。
19.項目1−9のいずれかに記載の方法を実行するように動作可能な特定用途向け電子部品。
20.テスト装置が項目1−5のいずれかに基づいて流体導管内の流体の流速を測定しているかどうかを判定する方法であって、
−前記流体導管に対して所定の速度を有する流体を、前記流体導管に供給すること、
−テストインパルス信号を、前記流体導管の第1の箇所に取り付けられた、前記テスト装置の第1の超音波トランスデューサに与えること、
−前記テストインパルス信号に対するテスト応答信号を、前記流体導管の第2の箇所に取り付けられた、前記テスト装置の第2の超音波トランスデューサで受信すること、
−前記応答信号からテスト測定信号を導出すること、且つ該テスト測定信号の導出は応答信号を時間反転することを含むこと、及び
−前記テスト測定信号を、前記第1及び第2の超音波トランスデューサのいずれか他方で発せられる測定信号と比較することを含み、
前記テスト測定信号と前記測定信号が同じである場合、前記テスト装置は項目1−5のいずれかに基づいて流体導管内の流体の流速を求める方法を用いていると判定する方法。
21.伝搬時間差式超音波流量計で流速を測定する装置であって、
−第1の超音波素子用の第1コネクタ、
−第2の超音波素子用の第2コネクタ、
−インパルス信号を前記第1コネクタに送る送信器、
−前記インパルス信号に対する応答信号を、前記第2コネクタから受信する受信器、
−前記応答信号のうち選択された一部分を時間反転し、反転信号を得る反転部、及び
−前記反転信号から測定信号を導出し、該測定信号をコンピュータ読み取り可能メモリに記憶する処理部を備えており、
流体導管内の流体の流速を求める前記装置が、
−前記流体導管に対して速度を有する流体を、前記流体導管に供給すること、
−測定信号を前記第1及び第2の超音波素子の一方に与えること、
−前記測定信号に対する第1の応答信号を、前記第1及び第2の超音波素子の他方で測定すること、及び
−前記第1の応答信号から前記流体の流速を導出すること
によって用いられ、
テストインパルス信号をテスト装置の第1の超音波素子に与える場合には、
−前記流体導管の第2の箇所に取り付けられた、前記テスト装置の第2の超音波素子で、前記テストインパルス信号に対するテスト応答信号を受信すること、及び
−前記応答信号からテスト測定信号を導出すること、且つ該テスト測定信号の導出は応答信号を時間反転することを含み、
−前記テスト測定信号と前記第1又は第2の超音波素子で発せられる測定信号が同じである装置。
10 流量計装置
11 上流側圧電素子
12 パイプ
13 下流側圧電素子
14 平均流の方向
15 第1の演算部
16 第2の演算部
17 信号路
20 信号路
22 圧電素子
23 圧電素子
31−52 圧電素子
60,60’ フロー測定装置
61 第1コネクタ
62 第2コネクタ
63 マルチプレクサ
64 DAC
65 ADC
66 デマルチプレクサ
67 信号選択部
68 信号反転部
69 帯域フィルタ
70 メモリ
71 速度演算部
72 インパルス信号生成器
73 測定信号生成器
74 指令ライン
75 指令ライン
76 DDS
77 標準発振器
78 周波数制御レジスタ
79 数値制御式発振器
80 ローパスフィルタ

Claims (5)

  1. 伝搬時間差式超音波流量計で流速を測定する装置であって、
    第1の超音波素子用の第1コネクタ、
    第2の超音波素子用の第2コネクタ、
    インパルス信号を前記第1コネクタに送る送信器、
    前記第1コネクタに接続されたD/A変換器、
    前記第2コネクタに接続されたA/D変換器、
    測定信号を記憶するコンピュータで読み取り可能なメモリ、
    前記インパルス信号に対する応答信号を、前記第2コネクタから受信する受信器、
    前記受信した応答信号の一部分又は前記受信した応答信号から導出された信号の一部分を選択する選択部、
    前記応答信号を時間反転し、反転信号を得る反転部、且つ、該反転部は、受信された応答信号を記録したサンプルの順序を逆にすることで、前記応答信号の選択部分を時間反転し、前記反転信号を得ることを含み、
    前記応答信号又は前記応答信号から導出された信号の振幅をデジタル化することで、前記反転信号から測定信号を導出し、該測定信号を前記メモリに記憶する処理部、を備えた装置。
  2. 請求項1に記載の装置において、
    前記第1コネクタ又は前記第2コネクタに接続可能な、前記測定信号を導出する測定信号生成器、及び、前記受信された応答信号から流速を導出する第2の処理部、を更に備え
    前記受信器は、前記第2コネクタからの前記測定信号に対する応答信号を受信するように構成されており、且つ、前記送信器は、さらに前記測定信号を前記第1コネクタに送信するように構成されている、装置。
  3. 請求項1に記載の装置において、
    D/A変換器を備えたダイレクトデジタル信号シンセサイザ、及び
    周波数制御レジスタ、標準発振器、数値制御式発振器、及び再構成ローパスフィルタを備え、前記D/A変換器は再構成ローパスフィルタを経て前記第1及び第2コネクタに接続可能である装置。
  4. 請求項1に記載の装置において、
    前記第1コネクタに接続された第1の超音波トランスデューサ、及び
    前記第2コネクタに接続された第2の超音波トランスデューサ、
    を更に備えた装置。
  5. 請求項1に記載の装置において、パイプの一部分を備え、第1の超音波トランスデューサが前記パイプの一部分の第1の箇所に取り付けられており、
    第2の超音波トランスデューサが前記パイプの一部分の第2の箇所に取り付けられている装置。
JP2017503164A 2014-07-29 2015-07-29 信号伝搬時間差式流量計 Active JP6682500B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
IB2014063502 2014-07-29
IBPCT/IB2014/063502 2014-07-29
PCT/IB2015/055724 WO2016016818A1 (en) 2014-07-29 2015-07-29 Improved signal travel time flow meter

Publications (2)

Publication Number Publication Date
JP2017525948A JP2017525948A (ja) 2017-09-07
JP6682500B2 true JP6682500B2 (ja) 2020-04-15

Family

ID=55216831

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017503164A Active JP6682500B2 (ja) 2014-07-29 2015-07-29 信号伝搬時間差式流量計

Country Status (19)

Country Link
US (3) US10551405B2 (ja)
EP (2) EP3273206B1 (ja)
JP (1) JP6682500B2 (ja)
CN (2) CN112098676A (ja)
AU (2) AU2015294885B2 (ja)
BR (1) BR112017001145B1 (ja)
CA (1) CA2955165C (ja)
CL (1) CL2017000181A1 (ja)
DE (1) DE202015009558U1 (ja)
DK (2) DK3273206T3 (ja)
ES (2) ES2792857T3 (ja)
IL (1) IL248650B (ja)
MX (1) MX370992B (ja)
PL (2) PL3273206T3 (ja)
RU (1) RU2657343C2 (ja)
SA (1) SA517380780B1 (ja)
SG (1) SG11201609000YA (ja)
UA (1) UA122209C2 (ja)
WO (1) WO2016016818A1 (ja)

Families Citing this family (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10801868B2 (en) * 2014-06-10 2020-10-13 Texas Instruments Incorporated Extended range ADC flow meter
SG11201609000YA (en) 2014-07-29 2016-11-29 Hydrovision Asia Pte Ltd Improved signal travel time flow meter
DE102015107750A1 (de) * 2015-05-18 2016-11-24 Endress + Hauser Flowtec Ag Meßsystem zum Messen wenigstens eines Parameters eines Fluids
CA2895361C (en) * 2015-06-19 2023-08-01 Accutron Instruments Inc. Method and system for ultrasonic airflow measurements
ES2763077T3 (es) * 2015-08-28 2020-05-27 Reliance Worldwide Corp Configuración y calibración de medidor de flujo
UA122259C2 (uk) 2016-01-18 2020-10-12 Гевееф Месзюстеме Аг Удосконалений часопролітний витратомір із формуванням променя звукового сигналу
ES2901159T3 (es) * 2016-04-25 2022-03-21 Gwf Messsysteme Ag Transductor acústico compacto de ángulo amplio
WO2018006912A1 (en) * 2016-07-08 2018-01-11 Apator Miitors Aps Ultrasonic flow meter with improved adc arrangement
BR112019014821A2 (pt) * 2017-01-20 2020-02-27 Gilbarco Inc. Analisador, dispensador e ambiente de combustível
FR3063815B1 (fr) * 2017-03-10 2019-03-22 Sagemcom Energy & Telecom Sas Procede de mesure d’une vitesse d’un fluide
FR3063814B1 (fr) * 2017-03-10 2019-03-22 Sagemcom Energy & Telecom Sas Procede de mesure d’une vitesse d’un fluide
US10830735B2 (en) 2017-03-20 2020-11-10 Triad National Security, Llc Simultaneous real-time measurement of composition, flow, attenuation, density, and pipe-wallthickness in multiphase fluids
FR3065106B1 (fr) * 2017-04-10 2019-08-02 Efs Sa Procede de commande d'un appareil de mesure a ultrasons du debit d'un fluide dans un canal de mesure
FR3065070B1 (fr) 2017-04-10 2019-06-14 Efs Sa Structure de gestion des ondes parasites d’un debitmetre a ultrason
US10473499B2 (en) * 2017-11-16 2019-11-12 Cameron International Corporation System and method of metering with array of transducers
FR3074913B1 (fr) * 2017-12-08 2019-11-22 Sagemcom Energy & Telecom Sas Procede de mesure d'une vitesse d'un fluide
EP3521773B1 (de) 2018-02-06 2021-09-29 SICK Engineering GmbH Ultraschall-durchflussmessvorrichtung und verfahren zum bestimmen einer strömungsgeschwindigkeit
CN108387279A (zh) * 2018-04-09 2018-08-10 河北科鼎机电科技有限公司 一种基于声波矩阵的液体流量测量装置及测量方法
RU2695282C1 (ru) * 2018-11-08 2019-07-22 Открытое акционерное общество Омское производственное объединение "Радиозавод им. А.С. Попова" (РЕЛЕРО) Ультразвуковой расходомер
WO2020186473A1 (zh) * 2019-03-20 2020-09-24 深圳市汇顶科技股份有限公司 飞行时间产生电路以及相关芯片、流量计及方法
DE102019108189A1 (de) * 2019-03-29 2020-10-01 Krohne Ag Ultraschalldurchflussmessgerät, Verfahren zum Betreiben eines Ultraschall-Durchflussmessgeräts, Messverbund und Verfahren zum Betreiben eines Messverbunds
US11231311B2 (en) 2019-05-31 2022-01-25 Perceptive Sensor Technologies Llc Non-linear ultrasound method and apparatus for quantitative detection of materials
BR112021025099A2 (pt) * 2019-05-31 2022-04-26 Perceptive Sensor Tech Llc Método de ultrassom não-linear e aparato para detecção quantivativa de materiais (líquidos, gás, plasma)
CN211452465U (zh) * 2019-09-30 2020-09-08 霍尼韦尔(天津)有限公司 超声波流量计和流体管路
CN110887549A (zh) * 2019-10-17 2020-03-17 浙江威星智能仪表股份有限公司 一种超声波燃气表超声波飞行时间校准方法
EP4043838A4 (en) * 2019-11-15 2022-11-02 Shenzhen Goodix Technology Co., Ltd. FLOW VELOCITY MEASUREMENT CIRCUIT, ASSOCIATED CHIP, AND FLOW METER
US11359950B2 (en) 2019-12-10 2022-06-14 Johnson Controls Tyco IP Holdings LLP Reduced length valve assembly with ultrasonic flow sensor
US11732927B2 (en) * 2020-04-09 2023-08-22 Rheem Manufacturing Company Systems and methods for preventing and removing chemical deposits in a fluid heating device
GB2595224A (en) * 2020-05-18 2021-11-24 Univ Warwick Ultrasonic flow measurement
CN111649791B (zh) * 2020-06-18 2022-03-11 浙江省水利河口研究院 一种水面上测量多流态水体流量的装置及方法
CN112415223B (zh) * 2020-11-19 2021-09-24 中国科学院大学 一种液态金属内部的速度测量方法、装置及存储介质
WO2022120074A1 (en) 2020-12-02 2022-06-09 Perceptive Sensor Technologies Llc Variable angle transducer interface block
US11536696B2 (en) 2020-12-04 2022-12-27 Perceptive Sensor Technologies, Inc. In-wall multi-bounce material property detection and acoustic signal amplification
US11585690B2 (en) 2020-12-04 2023-02-21 Perceptive Sensor Technologies, Inc. Multi-path acoustic signal improvement for material detection
CN116829915A (zh) 2020-12-04 2023-09-29 感知传感器技术股份有限公司 在分层环境中的声学温度测量
US11788904B2 (en) 2020-12-04 2023-10-17 Perceptive Sensor Technologies, Inc. Acoustic temperature measurement in layered environments
CA3201086A1 (en) 2020-12-04 2022-06-09 Lazar Bivolarsky Apparatus, system, and method for the detection of objects and activity within a container
US11549839B2 (en) 2020-12-04 2023-01-10 Perceptive Sensor Technologies, Inc. Systems and methods for determining floating roof level tilt and characterizing runoff
US11604294B2 (en) 2020-12-04 2023-03-14 Perceptive Sensor Technologies, Inc. Determining layer characteristics in multi-layered environments
WO2022147234A1 (en) 2020-12-30 2022-07-07 Perceptive Sensor Technologies, Inc. Evaluation of fluid quality with signals
DK4067833T3 (da) * 2021-03-31 2024-03-18 Gwf Messsysteme Ag Fulddigital løbetidsflowmåler med tidsrevers akustik
WO2023154514A1 (en) 2022-02-11 2023-08-17 Perceptive Sensor Technologies, Inc. Acoustic signal detection of material composition in static and dynamic conditions
US11940420B2 (en) 2022-07-19 2024-03-26 Perceptive Sensor Technologies, Inc. Acoustic signal material identification with nanotube couplant

Family Cites Families (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6018005B2 (ja) * 1979-12-16 1985-05-08 株式会社荏原製作所 透過形測定モ−ドと反射形測定モ−ドとを自動切換可能な超音波流速流量計
US4462261A (en) * 1982-04-27 1984-07-31 The Babcock & Wilcox Company Mass and velocity flowmeter
US4882934A (en) * 1986-03-12 1989-11-28 Charles B. Leffert Ultrasonic instrument to measure the gas velocity and/or the solids loading in a flowing gas stream
US4787252A (en) * 1987-09-30 1988-11-29 Panametrics, Inc. Differential correlation analyzer
US5117698A (en) * 1988-12-07 1992-06-02 Joseph Baumoel Pulse train detection in transit time flowmeter
US5650571A (en) * 1995-03-13 1997-07-22 Freud; Paul J. Low power signal processing and measurement apparatus
US6386018B1 (en) * 1996-03-11 2002-05-14 Daniel Industries, Inc. Ultrasonic 2-phase flow apparatus and stratified level detector
FR2781047B1 (fr) * 1998-07-10 2000-09-01 Faure Herman Debitmetre a ultrasons multicorde
AU2003902318A0 (en) * 2003-05-14 2003-05-29 Vision Fire And Security Pty Ltd Improved Sensing Apparatus And Method
EP1719980B1 (en) * 2004-02-26 2017-05-03 Fuji Electric Co., Ltd. Ultrasonic flowmeter and ultrasonic flow rate measurement method
US7845242B2 (en) * 2006-07-28 2010-12-07 Micro Motion, Inc. Three pickoff sensor flow meter
DE102009046871A1 (de) * 2008-11-21 2010-06-17 Flexim Flexible Industriemesstechnik Gmbh Verfahren und Vorrichtung zur Kalibrierung von Messumformern von Ultraschall-Durchflussmessgeräten
WO2011069131A1 (en) * 2009-12-04 2011-06-09 Google Inc. Presenting real time search results
US8141434B2 (en) * 2009-12-21 2012-03-27 Tecom As Flow measuring apparatus
CN102116777A (zh) * 2009-12-30 2011-07-06 上海申瑞电力科技股份有限公司 有压管道液体流动速度及流量的超声波测量方法
US8820147B2 (en) * 2010-09-03 2014-09-02 Los Alamos National Security, Llc Multiphase fluid characterization system
US20120271569A1 (en) * 2011-04-20 2012-10-25 Neptune Technology Group, Inc. Ultrasonic flow meter
US8700344B2 (en) * 2011-04-20 2014-04-15 Neptune Technology Group Inc. Ultrasonic flow meter
EP2739942B1 (en) * 2011-08-04 2016-04-27 SP Technical Research Institute Of Sweden Fluid visualisation and characterisation system and method
WO2013165314A1 (en) * 2012-05-04 2013-11-07 Dhi Water & Environment (S) Pte Ltd A flow meter system
US9689726B2 (en) * 2012-07-09 2017-06-27 Texas Instruments Incorporated Flow meter
DE102012212901A1 (de) * 2012-07-24 2014-01-30 Continental Automotive Gmbh Verfahren und Vorrichtung zum Betreiben einer akustischen Messvorrichtung
US9160325B2 (en) * 2013-01-22 2015-10-13 General Electric Company Systems and methods for fast kilovolt switching in an X-ray system
CN103336145B (zh) * 2013-06-18 2015-11-18 清华大学 一种针对管道内流体轴向流场的超声成像方法及设备
SG11201609000YA (en) 2014-07-29 2016-11-29 Hydrovision Asia Pte Ltd Improved signal travel time flow meter

Also Published As

Publication number Publication date
EP3273206A1 (en) 2018-01-24
DK3273206T3 (da) 2020-08-31
US9903882B2 (en) 2018-02-27
AU2015294885B2 (en) 2020-06-18
SG11201609000YA (en) 2016-11-29
US20200150146A1 (en) 2020-05-14
US20170082650A1 (en) 2017-03-23
AU2020233686B2 (en) 2022-04-21
BR112017001145B1 (pt) 2021-05-04
MX370992B (es) 2020-01-13
CA2955165A1 (en) 2016-02-04
DE202015009558U1 (de) 2018-04-03
WO2016016818A1 (en) 2016-02-04
EP3175205A1 (en) 2017-06-07
DK202000051U1 (en) 2020-06-10
DK202000051Y3 (da) 2020-09-10
MX2017001126A (es) 2017-07-04
RU2016145626A3 (ja) 2018-05-23
UA122209C2 (uk) 2020-10-12
CL2017000181A1 (es) 2017-09-15
CN112098676A (zh) 2020-12-18
EP3175205B1 (en) 2020-01-08
RU2657343C2 (ru) 2018-06-13
IL248650A0 (en) 2017-01-31
AU2015294885A1 (en) 2017-02-02
PL3175205T3 (pl) 2020-09-07
PL3273206T3 (pl) 2020-11-16
US10928414B2 (en) 2021-02-23
AU2020233686A1 (en) 2020-10-08
JP2017525948A (ja) 2017-09-07
ES2799527T3 (es) 2020-12-18
CN107003332B (zh) 2020-10-23
ES2792857T3 (es) 2020-11-12
CA2955165C (en) 2023-07-04
CN107003332A (zh) 2017-08-01
IL248650B (en) 2021-02-28
BR112017001145A2 (pt) 2018-01-30
RU2016145626A (ru) 2018-05-23
EP3175205A4 (en) 2017-07-26
US20170227568A1 (en) 2017-08-10
EP3273206B1 (en) 2020-06-10
SA517380780B1 (ar) 2020-08-25
US10551405B2 (en) 2020-02-04

Similar Documents

Publication Publication Date Title
JP6682500B2 (ja) 信号伝搬時間差式流量計
JP6727308B2 (ja) 改良型ビーム整形音響信号伝搬時間差式流量計
JP4983787B2 (ja) 超音波計測器
Rahiman et al. Design and development of ultrasonic process tomography
EA043606B1 (ru) Усовершенствованный времяпролетный расходомер с формированием луча звукового сигнала

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170317

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180523

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180524

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190308

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190604

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190829

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191203

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200225

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200304

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200325

R150 Certificate of patent or registration of utility model

Ref document number: 6682500

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250