JP2013015406A

JP2013015406A - 超音波計測装置

Info

Publication number: JP2013015406A
Application number: JP2011148342A
Authority: JP
Inventors: Kazutoshi Okamoto; 和年岡本
Original assignee: Yokogawa Electric Corp
Current assignee: Yokogawa Electric Corp
Priority date: 2011-07-04
Filing date: 2011-07-04
Publication date: 2013-01-24

Abstract

【課題】相関係数の多峰性による誤流速検出を減少させ、特に、大口径や高流流速条件下での測定においても安定した正しい流速測定を可能にする。
【解決手段】超音波を流体に照射して得られる複数の反射波に相関演算処理を施して流体に関する所定の物理量を計測する超音波計測装置１Ａであって、相関演算処理によって得られる複数の相関係数を所定回数積算演算処理するかまたはフィルタリング演算処理することによって、所定の物理量を演算する物理量演算手段（相関係数演算積算処理部１５２、積算メモリ１５３、流量演算部１５４）、を有する。
【選択図】図１

Description

本発明は、超音波を流体に照射して得られる反射波に相関演算処理を施し、流量や流速等の所定の物理量を計測する、超音波計測装置に関する。

流量計の一つに反射相関法を利用した反射相関型超音波流量計がある。この反射相関型超音波流量計は、流体に照射した超音波バルスが流体に含まれる気泡や粒子によって反射して得られる反射波に相互相関処理を施すことにより流速を計測し、この流速と流路断面積とに基づき流量を演算するものである。

２つの反射波の相関係数を求めると、例えば、図６に示すようになる。図６によれば、反射波（信号）の類似性が高ければ、正しい流速に相当する時間差（横軸）において相関係数値（縦軸）が最大になるが、実際には、種々の要因によって波形の類似性が崩れ、あるいは正弦波状の波形が連続していたりするため、誤った流速に相当する時間差においても複数のピークが存在し、場合によっては正しい流速のピークよりも大きな相関係数値になってしまうことがあり、誤流速検出の原因になっている。これは相関係数の多峰性の問題と呼ばれている。

上記した相関係数の多峰性の問題に対し、従来、インタバルの異なる２組の反射信号の相関の組み合わせから流速を求める方法（マルチインタバル法）や、発振周波数を変えて異なる周波数の２組の反射信号の相関の組み合わせから流速を求める方法が提案されている（例えば、特許文献１参照）。

上記したマルチインタバル法によれば、図６に両方向矢印で示した、流速に対する相関係数値のピークの周期Ｗは、Ｗ＝Ｃ／（２＊ｓｉｎθ＊ｔ＊ｆ）で表現される。ここで、ｆは信号の周波数、ｔは発振インタバル、θは超音波入射角、Ｃは流体の音速である。これらの方法は、ピークの周期Ｗが、発振インタバルｔや信号波形の周波数ｆによって異なるため、正しい流速に相当するピークの位置は、発振インタバルや周波数に依存することなく一致するが、隣接するピークではその位置がずれることを利用した方法である。

特開２０１０−１５１４５２号公報

上記した特許文献１に開示された技術によれば、誤流速検出を減らすことはできても、特に、大口径や高流速条件下での測定においては十分に誤流速検出を減らすことはできない。その原因の一つとして、流速に相当する相関係数のピーク位置の算出結果は、波形類似性の低下や実際の反射体の運動の乱れ等により大きなばらつきを持つことが挙げられ、特に大口径や高流速条件下でばらつきが大きくなる。すなわち、１回の相関演算結果（ピーク位置）はそれぞればらつきを持っており、マルチインタバル法を使用しても各インタバルの相関演算結果がばらつきを持つため、図７（ａ）に示すように、図中、点線で示す正しい流速付近に相関係数のピークが一致して現れず、図７（ｂ）に示すように、ばらつきによって正しい流速から外れたところに相関係数のピークが現れ、誤った流速付近にある相関係数のピークを検出してしまう。

本発明は上記した課題を解決するためになされたものであり、相関係数の多峰性による誤流速検出を減少させ、特に、大口径や高流流速条件下での測定においても安定した正しい流速測定を可能にした超音波計測装置を提供することを目的とする。また、マルチインタバル法を改良することによって誤流速検出を減少させる超音波計測装置を提供することも目的とする。

上記した課題を解決するために本発明は、超音波を流体に照射して得られる複数の反射波に相関演算処理を施して前記流体に関する所定の物理量を計測する超音波計測装置であって、前記相関演算処理によって得られる複数の相関係数を所定回数積算演算処理するかまたはフィルタリング演算処理することによって、前記所定の物理量を演算する物理量演算手段、を有することを特徴とする。

本発明によれば、物理量演算手段は、相関係数演算結果を積算演算処理し、あるいはフィルタリング演算処理することで、相関係数がピークとなる流速がばらつくのを抑制する。このため、相関係数の多峰性による誤流速検出を減少させ、特に、大口径や高流流速条件下での測定においても安定した正しい流速測定を可能にした超音波計測装置を提供することができる。

本発明において、異なるタイムインタバルで間欠する少なくとも２組の前記超音波を前記流体に照射する超音波照射手段、を有し、前記物理量演算手段は、前記流体内で反射して得られる少なくとも２組の反射信号を受信すると、前記相関演算処理によって得られる少なくとも２組の相関係数を所定回数積算演算処理するかまたはフィルタリング演算処理することによって前記所定の物理量を演算することを特徴とする。

本発明によれば、物理量演算手段がインタバルの異なる２組の反射信号の相関に積算演算処理を施すことにより、積算回数に相当する時間内の、いわば平均的な流速に相当するところに相関係数のピークが現れるため、各々が持つばらつきは小さくなり、ばらつきの小さくなったピーク位置同士の組み合わせでマルチインタバル処理を行うことにより、マルチインタバル法を改良することによって誤流速検出を減少させる超音波計測装置を提供することができる。また、積算演算処理によらず、ローパスフィルタ処理を施すことでばらつきを減らすことができ、例えば、ＩＩＲ（Infinite Impulse Response）フィルタを用いる場合、上記したばらつきを減らす他に、相関係数の演算結果を蓄積するメモリの容量が少なくて済むという効果も期待できる。

本発明によれば、相関係数の多峰性による誤流速検出を減少させ、特に、大口径や高流流速条件下での測定においても安定した正しい流速測定を可能にした超音波計測装置を提供することができる。また、マルチインタバル法を改良することによって誤流速検出を減少させる超音波計測装置を提供することもできる。

本発明の実施の形態１に係る超音波計測装置の構成を示すブロック図である。本発明の実施の形態１に係る超音波計測装置の動作を示すフローチャートである。本発明の実施の形態２に係る超音波計測装置の構成を示すブロック図である。本発明の実施の形態２に係る超音波計測装置の動作を示すフローチャートである。本発明の実施の形態２に係る超音波計測装置の相関係数演算結果を時間軸上に示した図である。相関係数の多峰性を説明するために示した図である。従来の超音波計測装置による相関係数演算結果を時間軸上に示した図である。

以下、添付図面を参照して本発明を実施するための実施の形態（以下、単に本実施形態という）について詳細に説明する。

（実施形態１の構成）
図１は、本実施形態１に係る超音波計測装置の構成を示すブロック図である。図１によれば、本実施形態に係る超音波計測装置１Ａは、バースト信号送信回路部１１と、増幅器１２ａ、１２ｂと、トランスデューサ１３と、ＡＤＣ（Analog-Digital Converter）１４と、演算処理部１５と、制御部１６と、により構成される。

バースト信号送信回路部１１は、超音波の周波数帯域である、例えば５００ｋＨｚ〜２ＭＨｚの周波数帯域の何れかの周波数ｆを有する複数のバースト信号を増幅器１２ａに出力する電子回路である。バースト信号送信回路部１１は、正弦波発振回路と、この正弦波発振回路が発振した周波数ｆの正弦波信号を連続して繰り返すバースト信号に強度変調する強度変調回路とから構成される。増幅器１２ａは、バースト信号送信回路部１１から入力された送信バースト信号を所定の増幅度で増幅してトランスデューサ１３に出力する。なお、バースト信号の波形は、正弦波の他に、矩形波を用いてもよい。

トランスデューサ１３は、図示するように、管路１０の外表面に設けられており、増幅器１２ａから入力された送信バースト信号を計測用超音波バースト信号に変換し、流体Ｘの流れ方向に直交する方向に対して角度θで管路１０内に出射すると共に、計測用超音波バースト信号が管路１０内の気泡等の粒子で反射することにより管路１０から入射した評価用超音波バースト信号を電気信号（受信信号）に変換する電気／音響変換器である。

増幅器１２ｂは、トランスデューサ１３から入力された受信信号を所定の増幅度で増幅してＡＤＣ１４に出力する。ＡＤＣ１４は、増幅器１２ｂから入力された受信信号を、サンプリング定理を満たす所定のサンプリングレート、つまり、上記した周波数ｆの２倍以上の繰り返し周期でサンプリングすることによりデジタル信号に変換して演算処理部１５（ウォールフィルタ）へ出力する。

演算処理部１５は、制御部１６による制御の下で、超音波を流体Ｘに照射して得られる複数の反射波から相関演算処理によって得られる複数の相関関数を所定回数積算演算処理することにより、流体Ｘの流量や流速等の物理量を演算する物理量演算手段として機能する。このため、演算処理部１５は、ウォールフィルタ１５１と、相関係数積算演算処理部１５２と、積算メモリ１５３と、流量演算部１５４と、に機能的に分解して表現することができる。

ウォールフィルタ１５１は、反射信号からなる評価用超音波バースト信号から残響振動成分を除去した受信データを相関係数積算演算処理部１５２に出力する。相関係数積算演算処理部１５２は、ウォールフィルタ１５１から入力される受信データ、つまり反射信号からなる評価用超音波バースト信号から上記した残響振動成分が除去されたものに相当するバースト信号について複数の時間領域毎に相関演算処理を施す。相関係数積算演算処理部１５２は、更に、相関演算処理によって得られる相関係数に所定回数積算演算処理を施し、過去複数回の相関係数を累積加算した結果を積算メモリ１５３に格納する。

ここで、相関係数積算演算処理とは、例えば、相関係数の積算回数をｍ回とすれば、積算相関係数Ｆ（ｖ）＝Σ_{ｋ（ｋ＝１〜ｍ）}ｆ_ｋ（ｖ）を求めることにより、積算回数ｍに相当する時間内の、いわば平均的な流速に相当するところに相関係数のピークを出現させることである。そして、気泡や粒子の移動速度の分布、つまり、管路１０内における流体Ｘの流速の分布（流速分布）を求め、積算メモリ１５３に格納するとともに流量演算部１５４に出力する。流量演算部１５４は、相関係数積算演算処理部１５２で求められた流速分布と、流速を算出した位置近傍での管路１０の断面積Ｓとに基づき流体Ｘの流量を演算する。

なお、制御部１６は、バースト信号送信回路部１１、ＡＤＣ１４、および上記したウォールフィルタ１５１、相関係数積算演算処理部１５２、積算メモリ１５３、流量演算部１５４からなる演算処理部１５、とを同期制御することにより、バースト信号送信回路部１１における送信バースト信号の発生タイミングに同期して受信信号をＡＤＣ１４で受信データに変換させると共に、この受信データを演算処理部１５で演算処理させる。

（実施形態１の動作）
以下、図２のフローチャートを参照しながら、図１に示す本実施形態１に係る超音波計測装置１Ａの動作について詳細に説明する。

まず、制御部１６による制御の下、バースト信号送信回路部１１は、例えば、５００ｋＨｚ〜２ＭＨｚの超音波周波数帯域の何れかの周波数ｆを有する複数のバースト信号を、増幅器１２ａ経由でトランスデューサ１３に出力する。これを受けてトランスデューサ１３は、増幅器１２ａから入力された送信バースト信号を計測用超音波バースト信号に変換し、流体Ｘの流れ方向に直交する方向に対して角度θで管路１０内に出射する（ステップＳ１０１）。また、トランスデューサ１３は、計測用超音波バースト信号が管路１０内の気泡等の粒子で反射することにより管路１０から入射した評価用超音波バースト信号を受信すると（ステップＳ１０２“ＹＥＳ”）、受信信号に変換して増幅器１２ｂに出力する（ステップＳ１０３）。

増幅器１２ｂは、トランスデューサ１３から入力された受信信号を所定の増幅度で増幅してＡＤＣ１４に出力し、ＡＤＣ１４は、増幅器１２ｂから入力された受信信号を、サンプリングすることによりデジタル信号に変換して演算処理部１５のウォールフィルタ１５１へ出力する。ウォールフィルタ１５１は、反射信号からなる評価用超音波バースト信号から残響振動成分を除去した受信データを相関係数積算演算処理部１５２に出力する。これを受けて相関係数積算演算処理部１５２は、ウォールフィルタ１５１から入力される受信データ、つまり反射信号からなる評価用超音波バースト信号から上記した残響振動成分が除去されたものに相当するバースト信号について複数の時間領域毎に相関演算処理を施して相関係数ｆ（ｖ）を算出する（ステップＳ１０４）。

相関係数積算演算処理部１５２は、更に、この相関演算処理によって得られる相関係数ｆ（ｖ）に積算演算処理をｍ回繰り返し（ステップＳ１０５“ＹＥＳ”）、積算相関係数Ｆ（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ）を求めて積算メモリ１５３に格納する（ステップＳ１０６）。ここで、相関演算処理を１回行うと、例えば、図６や図７に示したような周期的なピークを持つ相関係数ｆ（ｖ）が得られ、この相関演算を繰り返し行うことにより、ｆｋ（ｖ）が求まるが、相関係数がピークになる流速が都度ばらつき、従来のマルチインタバル処理を用いても誤流速を検出してしまう原因になっている。しかしながら、上記した積算演算処理によって積算相関係数Ｆ（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ）を求め、相関係数がピークになる流速のばらつきを小さくすることができる。すなわち、測定対象の平均流速から大きく外れた流速において相関係数がピークとなることが少なくなる。このようにして、流量演算部１５４は、過去複数回の流速結果をもとにヒストグラムを作成して、ヒストグラム分布の山から明らかに外れた誤流速を排除する（ステップＳ１０７：ヒストグラム処理）。

最後に、流量演算部１５４は、積算メモリ１５３に格納された相関係数積算演算処理部１５２で求めた流速分布を読み出し、流速を算出した位置近傍での管路１０の断面積Ｓに基づき流体Ｘの流量を演算する（ステップＳ１０８：流速プロファイル算出）。

（実施形態１の効果）
以上説明のように、本実施形態１に係る超音波計測装置１Ａによれば、相関係数積算演算処理部１５２、積算メモリ１５３、流量演算部１５４からなる物理量演算手段が、相関係数演算結果を積算演算処理することで、相関係数がピークとなる流速がばらつくのを抑制する。このため、相関係数の多峰性による誤流速検出を減少させ、特に、大口径や高流流速条件下での測定においても安定した正しい流速測定が可能になる。すなわち、相関演算に積算演算処理を付加するで、積算回数に相当する時間内の、いわば平均的な流速に相当するところに相関係数のピークが現れるため、各々が持つ流速のばらつきは小さくなる。また、本実施形態１に係る超音波計測装置１Ａによれば、ハードウェアの変更や追加の必要はなく、ソフトウェアによる演算処理の追加のみで誤流速の検出を減らすことができるため、廉価構成で新規な超音波計測装置１Ａを設計することができる。

なお、本実施形態１に係る超音波計測装置１Ａによれば、相関係数の積算にあたり、ｍ回単純加算するものとして説明したが、積算演算処理に代わって、ＦＩＲ（有限インパルス応答）フィルタやＩＩＲ（無限インパルス応答）フィルタ等を使用した重み付き加算で代替しても良い。要するに、従来のマルチインタバル法によれば、１回の相関演算毎のピークとなる流速結果がばらつくため、ＬＰＦ（Low Pass Filter）演算処理を施すことで、上記したばらつきを減らすことができる。特に、ＩＩＲフィルタを用いる場合では、相関係数演算の結果を蓄積するメモリの容量を減らすことができるといった効果も合わせ持つ。

（実施形態２の構成）
図３は、本実施形態２に係る超音波計測装置１Ｂの構成を示すブロック図である。図３において、図１に示す実施形態１との構成上の差異は、バースト信号送信回路部１１に代わってマルチインタバル送信回路部１７が設けられ、更に、演算処理部１５に、マルチインタバル処理部１５５が付加された点にある。他の構成は、図１に示した実施形態１と同様である。

マルチインタバル送信回路部１７は、超音波の周波数帯域、例えば５００ｋＨｚ〜２ＭＨｚの周波数帯域の何れかの周波数ｆを有すると共に、異なる２つのタイムインタバル（第１インタバルＴ１、第２インタバルＴ２）でパルス状に間欠する複数の送信バースト信号を増幅器１２ａに出力する電子回路である。マルチインタバル送信回路部１７は、正弦波発振回路と、この正弦波発振回路が発振した周波数ｆの正弦波信号（連続信号）を上記したタイムインタバルで繰り返すバースト信号に強度変調する強度変調回路から構成される。増幅器１２ａは、このマルチインタバル送信回路部１７から入力された送信バースト信号を所定の増幅度で増幅してトランスデューサ１３に出力する。

マルチインタバル処理部１５５は、相関係数のピークが一致する流速を算出するものであり、相関係数積算演算処理部１５２から、第１インタバルで得られる積算相関係数Ｆ１（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ）と、第２インタバルで得られる積算相関係数Ｆ２（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ）とから算出する。

（実施形態２の動作）
以下、図４のフローチャートを参照しながら図３に示す本実施形態２に係る超音波計測装置１Ｂの動作について、図２に示す実施の形態１の動作との差異に着目して説明する。

まず、制御部１６による制御の下、マルチインタバル送信回路部１７は、超音波の周波数帯域、例えば５００ｋＨｚ〜２ＭＨｚの周波数帯域の何れかの周波数ｆを有すると共に、第１インタバルＴ１（１ｓｔ−２ｎｄ），第２インタバルＴ２（２ｎｄ−３ｒｄ）でパルス状に間欠する複数の送信バースト信号を増幅器１２ａに出力する（ステップＳ２０１）。

これを受けて、演算処理部１５では、相関係数積算演算処理部１５２が、第１インタバルでは（ステップＳ２０２“ＹＥＳ”）、トランスデューサ１３で受信し（ステップＳ２０４）、ＡＤＣ１４、ウォールフィルタ１５１を介して取り込まれる受信データから残響振動成分が除去されたものに相当するバースト信号について、複数の時間領域毎に相関演算処理を施し（ステップＳ２０５）、更に、この相関演算処理によって得られる相関係数にｍ回の積算演算処理（積算相関係数関数Ｆ（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ））を施すことにより、気泡や粒子の移動速度の分布、つまり管路１０内における流体Ｘの流速の分布（流速分布）を求め、積算メモリ１５３を介してマルチインタバル処理部１５５に出力する（ステップＳ２０６）。

また、第２インタバルでも同様（ステップＳ２０３“ＹＥＳ”）、トランスデューサ１３で受信し（ステップＳ２０４）、ＡＤＣ１４、ウォールフィルタ１５１を介して取り込まれる受信データから残響振動成分が除去されたものに相当するバースト信号について、複数の時間領域毎に相関演算処理を施し（ステップＳ２０５）、更に、この相関演算処理によって得られる相関係数にｍ回の積算演算処理（積算相関係数関数Ｆ（ｖ）＝Σ（ｋ＝１〜ｍ）ｆｋ（ｖ））を施すことにより、気泡や粒子の移動速度の分布、つまり管路１０内における流体Ｘの流速の分布（流速分布）を求め、積算メモリ１５３を介してマルチインタバル処理部１５５に出力する（ステップＳ２０６）。

続いて、マルチインタバル処理部１５５は、第１インタバルで求めた相関係数演算結果（相関係数値のピークの周期Ｗ１＝Ｃ／（２＊ｓｉｎθ＊ｔ１＊ｆ）と、第２インタバルで求めた相関係数演算結果（相関係数値のピークの周期Ｗ２＝Ｃ／（２＊ｓｉｎθ＊ｔ２＊ｆ））とから、相関係数のピークが一致する流速を算出する（ステップＳ２０７）。ここでは、図５（ａ）（ｂ）に双方向矢印で示すように、第１インタバル、第２インタバルでの積算演算処理により、（ａ）に示す第１インタバルで得られた相関係数演算結果と、（ｂ）に示す第２インタバルで得られた相関係数演算結果のいずれにおいてもばらつきの範囲が縮小されるため、測定対象の平均流速から大きく外れた流速において相関係数がピークとなることが少なくなる。このため、マルチインタバル処理部１５５は、ばらつきが縮小したピーク位置（流速に相当）同士の組み合わせでマルチインタバル処理を行うことで誤流速検出の機会が減少する。このことは、従来例として示した図７（ｂ）と比較して明らかである。

続いて、マルチインタバル処理部１５５は、過去複数回の流速の結果に基づき誤流速を排除し（ヒストグラム処理：ステップＳ２０８）、最後に、流量演算部１５４は、積算メモリ１５３に格納された相関係数積算演算処理部１５２で求めた流速分布を読み出し、流速を算出した位置近傍での管路１０の断面積Ｓに基づき流体Ｘの流量（流速プロファイル）を演算する処理は（ステップＳ２０９）、実施形態１と同様である。

（実施形態２の効果）
以上説明のように、本実施形態２に係る超音波計測装置１Ｂによれば、相関係数積算演算処理部１５２と、積算メモリ１５３と、流量演算部１５４と、マルチインタバル処理部１５５とからなる物理量演算手段が、インタバルの異なる２組の反射信号の相関に積算演算処理を施すことにより、積算回数に相当する時間内の、いわば平均的な流速に相当するところに相関係数のピークが現れるため、各々が持つばらつきは小さくなり、ばらつきの小さくなったピーク位置同士の組み合わせでマルチインタバル処理を行うことにより、マルチインタバル法を改良することによって誤流速検出を減少させる超音波計測装置１Ｂを提供することができる。

また、実施形態１同様、積算演算処理によらず、ローパスフィルタ処理を施すことでばらつきを減らすことができ、例えば、ＩＩＲフィルタを用いる場合、上記したばらつきを減らす他に、相関係数の演算結果を蓄積するメモリの容量が少なくて済むという効果も期待できる。

以上、本発明の好ましい実施形態について詳述したが、本発明の技術的範囲は上記実施形態に記載の範囲には限定されないことは言うまでもない。上記実施形態に、多様な変更または改良を加えることが可能であることが当業者に明らかである。またその様な変更または改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲の記載から明らかである。

１Ａ、１Ｂ・・・超音波計測装置、１１・・・バースト信号送信回路部、
１２ａ、１２ｂ・・・増幅器、１３・・・トランスデューサ、１４・・・ＡＤＣ、１５・・・演算処理部、１６・・・制御部、１７・・・マルチインタバル送信回路部、１５１・・・ウォールフィルタ、１５２・・・相関係数積算演算処理部、１５３・・・積算メモリ、１５４・・・流量演算部、１５５・・・マルチインタバル処理部

Claims

超音波を流体に照射して得られる複数の反射波に相関演算処理を施して前記流体に関する所定の物理量を計測する超音波計測装置であって、
前記相関演算処理によって得られる複数の相関係数を所定回数積算演算処理するか、またはフィルタリング演算処理することによって前記所定の物理量を演算する物理量演算手段、
を有することを特徴とする超音波計測装置。
異なるタイムインタバルで間欠する少なくとも２組の前記超音波を前記流体に照射する超音波照射手段、を有し、
前記物理量演算手段は、
前記流体内で反射して得られる少なくとも２組の反射信号を受信すると、前記相関演算処理によって得られる少なくとも２組の相関係数を所定回数積算演算処理するか、またはフィルタリング演算処理することによって前記所定の物理量を演算することを特徴とする請求項１記載の超音波計測装置。