JPH10221141A - 超音波ドップラー流速・流量計 - Google Patents
超音波ドップラー流速・流量計Info
- Publication number
- JPH10221141A JPH10221141A JP9024662A JP2466297A JPH10221141A JP H10221141 A JPH10221141 A JP H10221141A JP 9024662 A JP9024662 A JP 9024662A JP 2466297 A JP2466297 A JP 2466297A JP H10221141 A JPH10221141 A JP H10221141A
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- flow velocity
- peak
- signal
- fft
- 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.)
- Granted
Links
Landscapes
- Measuring Volume Flow (AREA)
Abstract
(57)【要約】
【課題】 水路の流速分布の乱れや泡による突発的なデ
ータによる誤動作をなくす。表面波による誤動作をなく
す。センサの前面に堆積するごみや泥による計測誤差を
なくす。 【解決手段】 送信信号と受信信号を混合部１１で混合
し、検波回路５で検波する。フィルタ部６で両信号の差
信号としての変移周波数信号を取り出し、ＦＦＴ部７で
スペクトル変換する。変換結果を演算部８で移動平均す
る（（ｂ）図）。ピーク山の検出のための比較を周波数
の最大値Ｎｍａｘより始め、第１のピーク山Ｎｐ1 を検
出した後、パワー値Ｍｅａｎ・ｆｆｔ（Ｎ）がピーク山
Ｎｐ1 のパワー値Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）に対する
比率で１／２以下になると比較を中止する。第１のピー
ク山Ｎｐ1 の周波数を平均流速に換算し流量を演算す
る。
ータによる誤動作をなくす。表面波による誤動作をなく
す。センサの前面に堆積するごみや泥による計測誤差を
なくす。 【解決手段】 送信信号と受信信号を混合部１１で混合
し、検波回路５で検波する。フィルタ部６で両信号の差
信号としての変移周波数信号を取り出し、ＦＦＴ部７で
スペクトル変換する。変換結果を演算部８で移動平均す
る（（ｂ）図）。ピーク山の検出のための比較を周波数
の最大値Ｎｍａｘより始め、第１のピーク山Ｎｐ1 を検
出した後、パワー値Ｍｅａｎ・ｆｆｔ（Ｎ）がピーク山
Ｎｐ1 のパワー値Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）に対する
比率で１／２以下になると比較を中止する。第１のピー
ク山Ｎｐ1 の周波数を平均流速に換算し流量を演算す
る。
Description
【０００１】
【発明の属する技術分野】本発明は超音波ドップラー効
果を利用して流体の流速や流量を計測する流速・流量計
に係わり、特に開水路に好適な超音波流速・流量計に関
する。
果を利用して流体の流速や流量を計測する流速・流量計
に係わり、特に開水路に好適な超音波流速・流量計に関
する。
【０００２】
【従来の技術】超音波ドップラー方式の流体速度測定方
法およびその装置として、特公平７−３３５０号の従来
技術が公知である。以下これを第１の従来技術と言う。
法およびその装置として、特公平７−３３５０号の従来
技術が公知である。以下これを第１の従来技術と言う。
【０００３】このものは、信号発生器において発生され
た特定周波数の電気信号を発信用トランスデューサ（発
信子）を介して流路内を流れる流体に発信し、流体中の
物体から発射される反射波を受信用トランスデューサ
（受信子）で受信し、発信用トランスデューサから発信
信号の周波数と受信用トランスデューサで受信された受
信信号の周波数との差信号（シフト信号）に基づいて前
記流体の速度を測定するドップラーシフトを利用した流
体速度測定法において、発信用トランスデューサおよび
受信用トランスデューサを互いに近接した状態で流路の
底部に設けると共に、前記差信号を、一定の周波数幅を
もった所定の周波数帯でスキャンすることにより各周波
数帯の強度を求めた後、加重平均法を適用することによ
り前記流体の平均流速を求めると共に、前記特定周波数
の電気信号を９０°移相して得られる９０°移相信号の
周波数と前記受信信号の周波数との差信号を前記差信号
と比較することにより前記流体の流れる方向を検出する
ようにしている。
た特定周波数の電気信号を発信用トランスデューサ（発
信子）を介して流路内を流れる流体に発信し、流体中の
物体から発射される反射波を受信用トランスデューサ
（受信子）で受信し、発信用トランスデューサから発信
信号の周波数と受信用トランスデューサで受信された受
信信号の周波数との差信号（シフト信号）に基づいて前
記流体の速度を測定するドップラーシフトを利用した流
体速度測定法において、発信用トランスデューサおよび
受信用トランスデューサを互いに近接した状態で流路の
底部に設けると共に、前記差信号を、一定の周波数幅を
もった所定の周波数帯でスキャンすることにより各周波
数帯の強度を求めた後、加重平均法を適用することによ
り前記流体の平均流速を求めると共に、前記特定周波数
の電気信号を９０°移相して得られる９０°移相信号の
周波数と前記受信信号の周波数との差信号を前記差信号
と比較することにより前記流体の流れる方向を検出する
ようにしている。
【０００４】また、この第１の従来技術では、上記流体
速度測定装置を水位センサと組み合わせることにより、
流体の流量を測定できることを示唆している。この第１
の従来技術では、流量を求めるには水位センサを必要と
し、構造が複雑になる。そこで本願出願人は、水位セン
サを要しない超音波ドップラー流量計を特願平８−３９
２４３号で提案した。以下これを第２の従来技術と言
う。
速度測定装置を水位センサと組み合わせることにより、
流体の流量を測定できることを示唆している。この第１
の従来技術では、流量を求めるには水位センサを必要と
し、構造が複雑になる。そこで本願出願人は、水位セン
サを要しない超音波ドップラー流量計を特願平８−３９
２４３号で提案した。以下これを第２の従来技術と言
う。
【０００５】この第２の従来技術の超音波ドップラー流
量計は、開水路の底面中央部に流体の流れ方向に対して
一定の仰角θで超音波信号を送信する送信素子と、該送
信素子に隣接してほぼ同じ方向に向けて配設され、流体
中の固体粒子又は気泡等からの超音波の反射信号を受信
して電気信号に変換する受信素子と、送信素子に連続し
て高周波信号を供給する送信回路部と、前記送信素子の
送信信号の周波数と受信素子の受信信号の周波数の差の
周波数をとるヘテロダイン検波部と、該ヘテロダイン検
波部で得た差信号をデジタル信号に変換するＡＤコンバ
ーターと、その信号を周波数スペクトルに変換する高速
フーリエ変換部（ＦＦＴ部）と、該高速フーリエ変換部
で得た周波数スペクトルのピーク山の周波数を求めるピ
ーク山検出部と、該ピーク山検出部で求めたピーク山の
周波数に基づいて流量を演算する流量演算部とを有する
受信演算部とを備えたものである。
量計は、開水路の底面中央部に流体の流れ方向に対して
一定の仰角θで超音波信号を送信する送信素子と、該送
信素子に隣接してほぼ同じ方向に向けて配設され、流体
中の固体粒子又は気泡等からの超音波の反射信号を受信
して電気信号に変換する受信素子と、送信素子に連続し
て高周波信号を供給する送信回路部と、前記送信素子の
送信信号の周波数と受信素子の受信信号の周波数の差の
周波数をとるヘテロダイン検波部と、該ヘテロダイン検
波部で得た差信号をデジタル信号に変換するＡＤコンバ
ーターと、その信号を周波数スペクトルに変換する高速
フーリエ変換部（ＦＦＴ部）と、該高速フーリエ変換部
で得た周波数スペクトルのピーク山の周波数を求めるピ
ーク山検出部と、該ピーク山検出部で求めたピーク山の
周波数に基づいて流量を演算する流量演算部とを有する
受信演算部とを備えたものである。
【０００６】そして、この第２の従来技術では、受信信
号をヘテロダイン検波する際、混合する信号を送信信号
より一定だけ高い周波数又は低い周波数でヘテロダイン
検波を行い、高速フーリエ変換部で得た周波数スペクト
ルのピーク山を検出することによって、受信信号の周波
数が送信信号の周波数より大きいか小さいかにより流体
の流れの方向を判別するようにしている。
号をヘテロダイン検波する際、混合する信号を送信信号
より一定だけ高い周波数又は低い周波数でヘテロダイン
検波を行い、高速フーリエ変換部で得た周波数スペクト
ルのピーク山を検出することによって、受信信号の周波
数が送信信号の周波数より大きいか小さいかにより流体
の流れの方向を判別するようにしている。
【０００７】送信信号の周波数をｆ0 、仰角をθ、流体
の流速をＶ、流体中の音速をＣとすると、ドップラーシ
フトΔｆは、 ｜Δｆ｜＝（Ｖ／Ｃ）・２ｆ0 ｃｏｓθ ・・・（１） であらわされ、流速Ｖに比例する。
の流速をＶ、流体中の音速をＣとすると、ドップラーシ
フトΔｆは、 ｜Δｆ｜＝（Ｖ／Ｃ）・２ｆ0 ｃｏｓθ ・・・（１） であらわされ、流速Ｖに比例する。
【０００８】周波数スペクトルのピーク山の周波数をピ
ーク山検出部で求め、このピーク山の周波数を数式
（１）を変形した下記数式（２）のΔｆに代入して流速
Ｖを求める。
ーク山検出部で求め、このピーク山の周波数を数式
（１）を変形した下記数式（２）のΔｆに代入して流速
Ｖを求める。
【０００９】 Ｖ＝Ｃ・Δｆ／（２ｆ0 ｃｏｓθ） ・・・（２） こうして求めた流速Ｖが、いわゆるマニングの式で求め
た流路の平均流速Ｖに良く対応していることから、開水
路の断面形状寸法、水面勾配、壁面の粗度係数などが既
知であれば、平均流速に関するマニングの式から水位が
一つに定まり、流量を求められる。こうして流量を求め
る演算は、流量演算部で行っている。
た流路の平均流速Ｖに良く対応していることから、開水
路の断面形状寸法、水面勾配、壁面の粗度係数などが既
知であれば、平均流速に関するマニングの式から水位が
一つに定まり、流量を求められる。こうして流量を求め
る演算は、流量演算部で行っている。
【００１０】マニング（Ｍａｎｎｉｎｇ）の式は、開水
路に対する平均流速を求める式として、中小河川や水路
に対し、かなり良く一致するものとしてよく使用されて
いる次の数式（３）として周知である（日刊工業新聞
社、昭和５４年発行、流量計測ハンドブック、４０７
頁）。
路に対する平均流速を求める式として、中小河川や水路
に対し、かなり良く一致するものとしてよく使用されて
いる次の数式（３）として周知である（日刊工業新聞
社、昭和５４年発行、流量計測ハンドブック、４０７
頁）。
【００１１】 Ｖ＝（１／ｎ）Ｒ2/3 ／Ｉ 1/2 ・・・（３） Ｖ：平均流速（ｍ／ｓ）。 Ｒ：径深（ｍ）、水流断面積Ａ／潤辺長Ｐで定義され
る。ただし潤辺長は、水流に接している水路壁の長さを
言う。
る。ただし潤辺長は、水流に接している水路壁の長さを
言う。
【００１２】Ｉ：水面勾配。 ｎ：マニングの粗度係数。 なお、第２の従来技術では、周波数スペクトルのピーク
山の検出は、周波数スペクトルの低周波域から最大周波
数までの各点数におけるパワー値を１つずつ比較して、
単に最大パワー値をもつ点の周波数をピーク山の周波数
として求めていた。
山の検出は、周波数スペクトルの低周波域から最大周波
数までの各点数におけるパワー値を１つずつ比較して、
単に最大パワー値をもつ点の周波数をピーク山の周波数
として求めていた。
【００１３】
【発明が解決しようとする課題】出願人は、第２の従来
技術を提案した後、更に研究を進めた結果、次のような
問題点を見出した。
技術を提案した後、更に研究を進めた結果、次のような
問題点を見出した。
【００１４】流れによるドップラー変移周波数Δｆは、
流路による流れの分布の乱れとか、泡やごみからの反射
が離散的であることなどから、ＦＦＴで得られた結果を
単純に加算平均するだけでは、突発的なパワー値を持つ
周波数が発生して、図４に示すように周波数スペクトル
の分布がなめらかな曲線ではない。そのためにドップラ
ー変移周波数のピーク山を判定しにくく、流速や流量を
測定するのに誤動作を生じるという問題点があった。
流路による流れの分布の乱れとか、泡やごみからの反射
が離散的であることなどから、ＦＦＴで得られた結果を
単純に加算平均するだけでは、突発的なパワー値を持つ
周波数が発生して、図４に示すように周波数スペクトル
の分布がなめらかな曲線ではない。そのためにドップラ
ー変移周波数のピーク山を判定しにくく、流速や流量を
測定するのに誤動作を生じるという問題点があった。
【００１５】また、下水排水路のように多量のごみや泥
が含まれている場合は、送信子、受信子などの超音波振
動子の前面にごみや泥が堆積したり、高濃度のものが超
音波振動子の近くをゆっくり移動するため、近くからの
反射波は強く、遠くからの受信波は強度が相当に低下し
て受信される。従って、図５に点線で示すように、本来
の全体の流体からの反射波によるピーク山ｆD が、同図
に実線で示すように、流速の低い信号Ｐ（ｆ0 ）より小
さな値ｆD ′になり、ピーク山検出が低い方になってし
まい、本来のピーク山ｆD を正しく検出できなくて、流
速や流量測定に誤動作を生じることがあるという問題点
があった。
が含まれている場合は、送信子、受信子などの超音波振
動子の前面にごみや泥が堆積したり、高濃度のものが超
音波振動子の近くをゆっくり移動するため、近くからの
反射波は強く、遠くからの受信波は強度が相当に低下し
て受信される。従って、図５に点線で示すように、本来
の全体の流体からの反射波によるピーク山ｆD が、同図
に実線で示すように、流速の低い信号Ｐ（ｆ0 ）より小
さな値ｆD ′になり、ピーク山検出が低い方になってし
まい、本来のピーク山ｆD を正しく検出できなくて、流
速や流量測定に誤動作を生じることがあるという問題点
があった。
【００１６】このことは、超音波振動子の前面に付着堆
積したごみを除去すると直ちに誤動作がなくなって正常
な動作に回復するという現象が確認されていることから
も明らかである。
積したごみを除去すると直ちに誤動作がなくなって正常
な動作に回復するという現象が確認されていることから
も明らかである。
【００１７】また、前記第２の従来技術において、開水
路又は非満管路では、ドップラー変移（ドップラーシフ
ト）のパワースペクトルのピーク周波数（ピーク山の周
波数）は、実測データでは最大流速の周波数より２０〜
２５％低い値、即ちほぼ平均流速に対応する値になって
いるが、水位が非常に低くなった場合には、流れによっ
て生じる表面波によるドップラーシフトが発生し、これ
は、本来の平均流速を示すピーク山ｆD よりも低速側に
発生し、その周波数ｆｗは零から平均流速のピーク山ｆ
D の間で不定であり、パワーのピーク値も平均流速のピ
ーク山ｆD のピーク値（ピークパワー値）を越えること
がある。
路又は非満管路では、ドップラー変移（ドップラーシフ
ト）のパワースペクトルのピーク周波数（ピーク山の周
波数）は、実測データでは最大流速の周波数より２０〜
２５％低い値、即ちほぼ平均流速に対応する値になって
いるが、水位が非常に低くなった場合には、流れによっ
て生じる表面波によるドップラーシフトが発生し、これ
は、本来の平均流速を示すピーク山ｆD よりも低速側に
発生し、その周波数ｆｗは零から平均流速のピーク山ｆ
D の間で不定であり、パワーのピーク値も平均流速のピ
ーク山ｆD のピーク値（ピークパワー値）を越えること
がある。
【００１８】その為、表面波のパワースペクトルのピー
ク山ｆｗを平均流速のピーク山ｆDとして検出してしま
うという誤動作をおこすことがあるという問題点もあっ
た（図６参照）。
ク山ｆｗを平均流速のピーク山ｆDとして検出してしま
うという誤動作をおこすことがあるという問題点もあっ
た（図６参照）。
【００１９】従来は、受信信号の周波数の変移（シフ
ト）をＦＦＴ変換して得られたスペクトルを低周波域か
ら高周波域に向けて単純に検索し、最大パワー値をもつ
周波数をもって流速に換算するものであったため、上述
のように正確な流速を計測できないと言う問題点が発生
していた。
ト）をＦＦＴ変換して得られたスペクトルを低周波域か
ら高周波域に向けて単純に検索し、最大パワー値をもつ
周波数をもって流速に換算するものであったため、上述
のように正確な流速を計測できないと言う問題点が発生
していた。
【００２０】そこで、本発明はかかる問題点を解消でき
る超音波ドップラー流速・流量計を提供することを目的
とする。
る超音波ドップラー流速・流量計を提供することを目的
とする。
【００２１】
【課題を解決するための手段】前記目的を達成するため
に、請求項１の発明は、液体が流れる水路（１）と、該
水路（１）の底面、中央またはごみや泥を避けるために
底面、中央より側面方向に寄ったところに設置された超
音波の送信子と受信子とからなるセンサ（２）と、送信
子に高周波を供給する駆動回路（３）と、受信信号の増
幅部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての偏移周
波数信号を取り出すフィルタ部（６）と、偏移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、ＦＦＴ部（７）でフ
ーリエ変換した結果を単純加算平均したスペクトルのパ
ワー値をｆｆｔ（Ｎ）とする周波数ｆｎのある点Ｎの前
後に幅ＭｎをもってＮ−ＭｎからＮ＋Ｍｎの間で周波数
のパワーの平均値Ｍｅａｎ・ｆｆｔ（Ｎ）を取り、その
点Ｎの元来のパワー値ｆｆｔ（Ｎ）に置き換える。
に、請求項１の発明は、液体が流れる水路（１）と、該
水路（１）の底面、中央またはごみや泥を避けるために
底面、中央より側面方向に寄ったところに設置された超
音波の送信子と受信子とからなるセンサ（２）と、送信
子に高周波を供給する駆動回路（３）と、受信信号の増
幅部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての偏移周
波数信号を取り出すフィルタ部（６）と、偏移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、ＦＦＴ部（７）でフ
ーリエ変換した結果を単純加算平均したスペクトルのパ
ワー値をｆｆｔ（Ｎ）とする周波数ｆｎのある点Ｎの前
後に幅ＭｎをもってＮ−ＭｎからＮ＋Ｍｎの間で周波数
のパワーの平均値Ｍｅａｎ・ｆｆｔ（Ｎ）を取り、その
点Ｎの元来のパワー値ｆｆｔ（Ｎ）に置き換える。
【００２２】Ｎを零と最大値Ｎｍａｘとの間で同様にし
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルのピーク山に基づいて流速・換算
することを特徴とする超音波流速・流量計である。
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルのピーク山に基づいて流速・換算
することを特徴とする超音波流速・流量計である。
【００２３】この発明では、Ｎの零から最大値Ｎｍａｘ
までのパワー値を平均値Ｍｅａｎ・ｆｆｔ（Ｎ）に置き
換えること（移動平均）により、変移周波数の分布を滑
らかにして突発的なパワー値をもつ周波数を除去し、ス
ペクトルのピークを容易に検出できるようにしている
（図１（ｂ）参照）。
までのパワー値を平均値Ｍｅａｎ・ｆｆｔ（Ｎ）に置き
換えること（移動平均）により、変移周波数の分布を滑
らかにして突発的なパワー値をもつ周波数を除去し、ス
ペクトルのピークを容易に検出できるようにしている
（図１（ｂ）参照）。
【００２４】請求項２の発明は、請求項１の超音波ドッ
プラー流速・流量計において、ピーク山の検出のための
比較を周波数の最大値Ｎｍａｘより始め、第１のピーク
山Ｎｐ1 を検出した後、パワー値Ｍｅａｎ・ｆｆｔ
（Ｎ）がピーク山Ｎｐ1 のパワー値Ｍｅａｎ・ｆｆｔ
（Ｎｐ1 ）に対する比率で一定値以下になったところで
比較を中止し、第１のピーク山Ｎｐ1 の周波数を流速に
換算することを特徴とするものである。
プラー流速・流量計において、ピーク山の検出のための
比較を周波数の最大値Ｎｍａｘより始め、第１のピーク
山Ｎｐ1 を検出した後、パワー値Ｍｅａｎ・ｆｆｔ
（Ｎ）がピーク山Ｎｐ1 のパワー値Ｍｅａｎ・ｆｆｔ
（Ｎｐ1 ）に対する比率で一定値以下になったところで
比較を中止し、第１のピーク山Ｎｐ1 の周波数を流速に
換算することを特徴とするものである。
【００２５】第２の従来技術では、スペクトルのパワー
ピーク値の検出はＮ＝０の低周波域より始めてＮｍａｘ
まで各点のパワー値を１つずつ比較し、単に最大パワー
値をもつ点Ｎｐを探していたが、ごみ、泥、表面波によ
るスペクトルを誤って検出してしまっていた。
ピーク値の検出はＮ＝０の低周波域より始めてＮｍａｘ
まで各点のパワー値を１つずつ比較し、単に最大パワー
値をもつ点Ｎｐを探していたが、ごみ、泥、表面波によ
るスペクトルを誤って検出してしまっていた。
【００２６】水路内では流水より早く動くものが無いた
め、この請求項２の発明では、ピークの検出動作をＮｍ
ａｘから始めて周波数が小さくなる方向に行い、第１の
ピークＮｐ1 を検出したあと、パワー値Ｍｅａｎ・ｆｆ
ｔ（Ｎ）とＭｅａｎ・ｆｆｔ（Ｎｐ1 ）の比率が一定値
以下になったところで比較を中止し、このときのピーク
山Ｎｐ1 を流速によるスペクトルのピーク山と判断して
そのときの周波数ｆNP 1 を流速に換算する（図１（ｃ）
参照）。
め、この請求項２の発明では、ピークの検出動作をＮｍ
ａｘから始めて周波数が小さくなる方向に行い、第１の
ピークＮｐ1 を検出したあと、パワー値Ｍｅａｎ・ｆｆ
ｔ（Ｎ）とＭｅａｎ・ｆｆｔ（Ｎｐ1 ）の比率が一定値
以下になったところで比較を中止し、このときのピーク
山Ｎｐ1 を流速によるスペクトルのピーク山と判断して
そのときの周波数ｆNP 1 を流速に換算する（図１（ｃ）
参照）。
【００２７】こうすることで、表面波スペクトルのピー
クＮｐ2 、ごみや泥によるスペクトルのピークＮｐ3 等
を誤って検出しないようにした。なお、ピーク検出動作
としての比較の中止を判定するための比率Ｍｅａｎ・ｆ
ｆｔ（Ｎ）／Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）は測定場所の
条件に応じて値を変えても良い。
クＮｐ2 、ごみや泥によるスペクトルのピークＮｐ3 等
を誤って検出しないようにした。なお、ピーク検出動作
としての比較の中止を判定するための比率Ｍｅａｎ・ｆ
ｆｔ（Ｎ）／Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）は測定場所の
条件に応じて値を変えても良い。
【００２８】請求項３の発明は、請求項２の超音波ドッ
プラー流速・流量計において、前記一定値を１／４〜３
／４の間に定めたことを特徴とするものである。そし
て、請求項４の発明は、液体が流れる水路（１）と、該
水路（１）の底面、中央またはごみや泥を避けるために
底面、中央より側面方向に寄ったところに設置された超
音波の送信子と受信子とからなるセンサ（２）と、送信
子に高周波を供給する駆動回路（３）と、受信信号の増
幅部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての変移周
波数信号を取り出すフィルタ部（６）と、変移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、ＦＦＴ部（７）でフ
ーリエ変換した結果を単純加算平均したスペクトルのパ
ワー値をｆｆｔ（Ｎ）とする周波数ｆｎのある点Ｎの前
後に幅ＭｎをもってＮ−ＭｎからＮ＋Ｍｎの間で周波数
のパワーの平均値Ｍｅａｎ・ｆｆｔ（Ｎ）を取り、その
点Ｎの元来のパワー値ｆｆｔ（Ｎ）に置き換える。
プラー流速・流量計において、前記一定値を１／４〜３
／４の間に定めたことを特徴とするものである。そし
て、請求項４の発明は、液体が流れる水路（１）と、該
水路（１）の底面、中央またはごみや泥を避けるために
底面、中央より側面方向に寄ったところに設置された超
音波の送信子と受信子とからなるセンサ（２）と、送信
子に高周波を供給する駆動回路（３）と、受信信号の増
幅部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての変移周
波数信号を取り出すフィルタ部（６）と、変移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、ＦＦＴ部（７）でフ
ーリエ変換した結果を単純加算平均したスペクトルのパ
ワー値をｆｆｔ（Ｎ）とする周波数ｆｎのある点Ｎの前
後に幅ＭｎをもってＮ−ＭｎからＮ＋Ｍｎの間で周波数
のパワーの平均値Ｍｅａｎ・ｆｆｔ（Ｎ）を取り、その
点Ｎの元来のパワー値ｆｆｔ（Ｎ）に置き換える。
【００２９】Ｎを零と最大値Ｎｍａｘとの間で同様にし
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルの分布曲線を微分または差分し、
曲線が上に凸の変曲点のうち、Ｎが最大のＮｐ1 を流速
スペクトルのピーク山と判断して、該ピーク山の周波数
を流速に換算することを特徴とするものである。
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルの分布曲線を微分または差分し、
曲線が上に凸の変曲点のうち、Ｎが最大のＮｐ1 を流速
スペクトルのピーク山と判断して、該ピーク山の周波数
を流速に換算することを特徴とするものである。
【００３０】差分をＮｍａｘとＮｍａｘ−１より始め
て、次々にＮ＝１とＮ＝０まで行うと図１（ｂ）の実線
の曲線から図２に示す微分曲線が得られる。すなわち、
差分をＮｍａｘより行い、始めてプラスから零をクロス
する点Ｎｐ1 を流速によるスペクトルのピーク山と判断
してこのときの変移周波数ΔｆNP1 を流速に換算する。
て、次々にＮ＝１とＮ＝０まで行うと図１（ｂ）の実線
の曲線から図２に示す微分曲線が得られる。すなわち、
差分をＮｍａｘより行い、始めてプラスから零をクロス
する点Ｎｐ1 を流速によるスペクトルのピーク山と判断
してこのときの変移周波数ΔｆNP1 を流速に換算する。
【００３１】こうすることで、表面波によるピーク山Ｎ
ｐ2 とか、ごみ、泥によるＮｐ3 を誤って流速によるピ
ーク山と認識することを避けるようにした。請求項１〜
４の本願発明では、受信信号（ドップラー信号）の変移
周波数の周波数スペクトルから検出したピーク山の周波
数を数式（２）のΔｆに代入して、第２の従来技術の場
合と同様に流速Ｖを求める。つまり流速Ｖに換算する。
ｐ2 とか、ごみ、泥によるＮｐ3 を誤って流速によるピ
ーク山と認識することを避けるようにした。請求項１〜
４の本願発明では、受信信号（ドップラー信号）の変移
周波数の周波数スペクトルから検出したピーク山の周波
数を数式（２）のΔｆに代入して、第２の従来技術の場
合と同様に流速Ｖを求める。つまり流速Ｖに換算する。
【００３２】こうして求めた流速Ｖは平均流速に相当す
るので、これとマニングの式などの平均流速公式から水
位が定まり、流量を求める。なお、上記の説明は送信周
波数と同じかまたは低い周波数でヘテロダイン検波を行
った場合であり、高い周波数で行った場合は流速が大き
くなれば偏移周波数は低くなり、ピーク周波数検出は低
周波側から高周波側に向けて行うこととなり逆になる。
るので、これとマニングの式などの平均流速公式から水
位が定まり、流量を求める。なお、上記の説明は送信周
波数と同じかまたは低い周波数でヘテロダイン検波を行
った場合であり、高い周波数で行った場合は流速が大き
くなれば偏移周波数は低くなり、ピーク周波数検出は低
周波側から高周波側に向けて行うこととなり逆になる。
【００３３】
【発明の実施の形態】次に本発明の好ましい実施の形態
を図面の実施例に基づいて説明する。 〔実施例１〕図１（ａ）（ｂ）（ｃ）と図３に示す実施
例１で、水路１を流れる液体中へセンサ２の超音波送信
子から仰角θで超音波信号を送信する。センサ２には前
記送信子に隣接してほぼ同じ方向を向けた受信子が配設
され、液体中のごみや気泡などからの超音波の返射信号
を受信して電気信号に変換する。駆動回路３は送信子に
高周波を供給する。
を図面の実施例に基づいて説明する。 〔実施例１〕図１（ａ）（ｂ）（ｃ）と図３に示す実施
例１で、水路１を流れる液体中へセンサ２の超音波送信
子から仰角θで超音波信号を送信する。センサ２には前
記送信子に隣接してほぼ同じ方向を向けた受信子が配設
され、液体中のごみや気泡などからの超音波の返射信号
を受信して電気信号に変換する。駆動回路３は送信子に
高周波を供給する。
【００３４】受信子の電気信号（受信信号）は増幅部４
で増幅され、混合部１１へ供給されて、前記駆動回路３
からの送信信号と混合された後、検波回路５で検波さ
れ、更にフィルタ部６で受信信号と送信信号の差信号と
しての変移周波数（ドップラーシフト周波数）信号を取
り出す。
で増幅され、混合部１１へ供給されて、前記駆動回路３
からの送信信号と混合された後、検波回路５で検波さ
れ、更にフィルタ部６で受信信号と送信信号の差信号と
しての変移周波数（ドップラーシフト周波数）信号を取
り出す。
【００３５】ＦＦＴ部（高速フーリエ変換部）７は前記
変移周波数信号を高速フーリエ変換してスペクトル変換
する。この実施例ではフーリエ変換の結果出力されるデ
ータのポイント数は１０２４であり、周波数範囲を５ｋ
Ｈｚとすると、周波数の分解能は４．８８Ｈｚとなる。
変移周波数信号を高速フーリエ変換してスペクトル変換
する。この実施例ではフーリエ変換の結果出力されるデ
ータのポイント数は１０２４であり、周波数範囲を５ｋ
Ｈｚとすると、周波数の分解能は４．８８Ｈｚとなる。
【００３６】１）．スペクトルの移動平均 ＦＦＴ部７で得たスペクトルのデータはＭＰＵからなる
ＣＰＵを備えた演算部８で流速換算されて、表示部９に
測定結果が表示される。
ＣＰＵを備えた演算部８で流速換算されて、表示部９に
測定結果が表示される。
【００３７】１回のＦＦＴのスペクトル値は１／２５さ
れて演算部９の第１レジスタ（図３参照、以後第２，第
３レジスタについても図３参照）に０から順にＮｍａｘ
まで格納され、同様に次のＦＦＴ結果を１／２５にして
第１レジスタに加算し、これを２５回まで行うと２５回
の単純加算平均となる。
れて演算部９の第１レジスタ（図３参照、以後第２，第
３レジスタについても図３参照）に０から順にＮｍａｘ
まで格納され、同様に次のＦＦＴ結果を１／２５にして
第１レジスタに加算し、これを２５回まで行うと２５回
の単純加算平均となる。
【００３８】更に、第１レジスタのＮｘ点の前後にそれ
ぞれ２０ポイントの幅Ｍｎの間で平均を取り、第２レジ
スタのＮｘ点に格納する。Ｎ＝０からＮ＝Ｎｍａｘまで
順に同じ様な方法で平均（移動平均）をして行く。但
し、Ｎ−Ｍｎ＜０の場合は０〜Ｎ＋Ｍｎの間で、Ｎ＋Ｍ
ｎ＞Ｎｍａｘの場合はＮ−Ｍｎ〜Ｎｍａｘの間で平均
（移動平均）を行う（図１（ｂ）参照）。
ぞれ２０ポイントの幅Ｍｎの間で平均を取り、第２レジ
スタのＮｘ点に格納する。Ｎ＝０からＮ＝Ｎｍａｘまで
順に同じ様な方法で平均（移動平均）をして行く。但
し、Ｎ−Ｍｎ＜０の場合は０〜Ｎ＋Ｍｎの間で、Ｎ＋Ｍ
ｎ＞Ｎｍａｘの場合はＮ−Ｍｎ〜Ｎｍａｘの間で平均
（移動平均）を行う（図１（ｂ）参照）。
【００３９】なお、上記具体的にあげられた値、ＦＦＴ
のデータ数、その加算平均回数、移動平均の個数の組み
合わせは、実際の測定場所の条件に応じて決められるも
のである。
のデータ数、その加算平均回数、移動平均の個数の組み
合わせは、実際の測定場所の条件に応じて決められるも
のである。
【００４０】２）．ピーク検出 次に演算部９で行うスペクトルのピーク検出（ピーク山
の検出）は、ごみ、泥、表面波によるスペクトルピーク
が流速に基づく流速スペクトルピークよりも低周波域に
発生するため、高周波域から行う。
の検出）は、ごみ、泥、表面波によるスペクトルピーク
が流速に基づく流速スペクトルピークよりも低周波域に
発生するため、高周波域から行う。
【００４１】そのピーク検出方法は、図３の第２レジス
タのＮｍａｘとＮｍａｘ−１の比較から始め、その大き
い方のＮ0 とＭｅａｎ・ｆｆｔ（Ｎ0 ）を第３レジスタ
に格納し、順にＮ＝０に向かって第３レジスタのＭｅａ
ｎ・ｆｆｔ（Ｎ0 ）とＭｅａｎ・ｆｆｔ（Ｎｘ）を比較
し、パワー値の大きい方を第３レジスタに書き替えて行
く。
タのＮｍａｘとＮｍａｘ−１の比較から始め、その大き
い方のＮ0 とＭｅａｎ・ｆｆｔ（Ｎ0 ）を第３レジスタ
に格納し、順にＮ＝０に向かって第３レジスタのＭｅａ
ｎ・ｆｆｔ（Ｎ0 ）とＭｅａｎ・ｆｆｔ（Ｎｘ）を比較
し、パワー値の大きい方を第３レジスタに書き替えて行
く。
【００４２】このようにして得られた最初の大きなパワ
ー値Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）がそれ以後のＮｘとの
比較を行う中でＭｅａｎ・ｆｆｔ（Ｎｘ）が例えば２倍
されてもＭｅａｎ・ｆｆｔ（Ｎｐ1 ）よりも小さくなっ
たとき、即ち、 Ｍｅａｎ・ｆｆｔ（Ｎｘ）＜（１／２）Ｍｅａｎ・ｆｆ
ｔ（Ｎｐ1 ） となったとき、Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）をスペクト
ルの最高周波数側にある極大値と見做し、点Ｎｐ1 にお
ける変移周波数ΔｆNP1 を流速スペクトルのピーク周波
数とする（図１（ｃ）参照）。
ー値Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）がそれ以後のＮｘとの
比較を行う中でＭｅａｎ・ｆｆｔ（Ｎｘ）が例えば２倍
されてもＭｅａｎ・ｆｆｔ（Ｎｐ1 ）よりも小さくなっ
たとき、即ち、 Ｍｅａｎ・ｆｆｔ（Ｎｘ）＜（１／２）Ｍｅａｎ・ｆｆ
ｔ（Ｎｐ1 ） となったとき、Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）をスペクト
ルの最高周波数側にある極大値と見做し、点Ｎｐ1 にお
ける変移周波数ΔｆNP1 を流速スペクトルのピーク周波
数とする（図１（ｃ）参照）。
【００４３】そして、ΔｆNP1 を流速に換算し、流量を
演算する。なお、上記ピーク検出を中止するときの比率 Ｍｅａｎ・ｆｆｔ（Ｎｘ）／Ｍｅａｎ・ｆｆｔ（Ｎ
ｐ1 ） の値は、実際の測定場所の条件に応じて１／４〜３／４
の間で変化させる。
演算する。なお、上記ピーク検出を中止するときの比率 Ｍｅａｎ・ｆｆｔ（Ｎｘ）／Ｍｅａｎ・ｆｆｔ（Ｎ
ｐ1 ） の値は、実際の測定場所の条件に応じて１／４〜３／４
の間で変化させる。
【００４４】〔実施例２〕差分をＮｍａｘとＮｍａｘ−
１から始めて、次々にＮ＝１とＮ＝０まで行うと図１
（ｂ）の実線の曲線から図２に示す微分曲線が得られ
る。
１から始めて、次々にＮ＝１とＮ＝０まで行うと図１
（ｂ）の実線の曲線から図２に示す微分曲線が得られ
る。
【００４５】そこで、差分をＮｍａｘから行い、始めて
プラスから零をクロスする点Ｎｐ1を流速スペクトルの
ピークと見做して、このときの変移周波数ΔｆNP1 を流
速スペクトルのピーク周波数とする。
プラスから零をクロスする点Ｎｐ1を流速スペクトルの
ピークと見做して、このときの変移周波数ΔｆNP1 を流
速スペクトルのピーク周波数とする。
【００４６】そしてΔｆNP1 を流速に換算し、流量を演
算する。こうすることで、表面波によるピーク山Ｎｐ2
とか、ごみ、泥によるＮｐ3 を誤って流速によるピーク
山と認識することを避ける。
算する。こうすることで、表面波によるピーク山Ｎｐ2
とか、ごみ、泥によるＮｐ3 を誤って流速によるピーク
山と認識することを避ける。
【００４７】
【発明の効果】本発明は上述のように構成されているの
で、ＦＦＴによって得られたスペクトルを移動平均して
滑らかな周波数分布曲線にすることで、水路の流速分布
の乱れや泡によって生じる突発的なパワー値をもつ周波
数のデータを除去して、正確な平均流速を求めることが
できる。
で、ＦＦＴによって得られたスペクトルを移動平均して
滑らかな周波数分布曲線にすることで、水路の流速分布
の乱れや泡によって生じる突発的なパワー値をもつ周波
数のデータを除去して、正確な平均流速を求めることが
できる。
【００４８】また、ＦＦＴ結果を単純に加算平均する場
合に比較して、移動平均を行うことで短時間に周波数分
布曲線を滑らかにでき、計測時間を減らすことができ
る。更に、請求項２〜４の発明では、水路の表面波によ
るスペクトルピークを検出しないようにして、表面波に
よる誤動作を防止し、平均流速を正確に求められ、計測
精度が向上する。そのうえ、センサの前面に下水路では
避けて通ることのできないごみや泥が堆積しても、短期
間のうちに計測が中断する虞れがなく、長期間安定に測
定ができ、その面からも流速・流量計の信頼性が向上す
る。
合に比較して、移動平均を行うことで短時間に周波数分
布曲線を滑らかにでき、計測時間を減らすことができ
る。更に、請求項２〜４の発明では、水路の表面波によ
るスペクトルピークを検出しないようにして、表面波に
よる誤動作を防止し、平均流速を正確に求められ、計測
精度が向上する。そのうえ、センサの前面に下水路では
避けて通ることのできないごみや泥が堆積しても、短期
間のうちに計測が中断する虞れがなく、長期間安定に測
定ができ、その面からも流速・流量計の信頼性が向上す
る。
【図１】本発明の実施例で、（ａ）はブロック図、
（ｂ）はＦＦＴデータのパワー値を示す図、（ｃ）は流
速スペクトルのピークを検出する方法を説明するための
周波数分布曲線である。
（ｂ）はＦＦＴデータのパワー値を示す図、（ｃ）は流
速スペクトルのピークを検出する方法を説明するための
周波数分布曲線である。
【図２】図１（ｃ）の周波数分布曲線を差分した微分曲
線である。
線である。
【図３】図１の演算部８に含まれるレジスタを示す図で
ある。
ある。
【図４】従来技術のＦＦＴ結果のパワー値を示す図であ
る。
る。
【図５】超音波振動子の前面に付着したごみや泥による
周波数分布曲線の変化を説明する線図である。
周波数分布曲線の変化を説明する線図である。
【図６】流れによって生じる表面波によるパワースペク
トルのピーク山を説明する線図である。
トルのピーク山を説明する線図である。
１ 水路 ２ センサ ３ 駆動回路 ４ 増幅部 ５ 検波回路 ６ フィルタ部 ７ ＦＦＴ部 ８ 演算部 ９ 表示部
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成１０年３月３日
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】請求項１
【補正方法】変更
【補正内容】
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００２２
【補正方法】変更
【補正内容】
【００２２】Ｎを零と最大値Ｎｍａｘとの間で同様にし
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルのピーク山に基づいて流速換算す
ることを特徴とする超音波流速・流量計である。
て平均値に置き換えて移動平均すると共に、こうして移
動平均したスペクトルのピーク山に基づいて流速換算す
ることを特徴とする超音波流速・流量計である。
Claims (4)
- 【請求項１】 液体が流れる水路（１）と、該水路
（１）の底面、中央またはごみや泥を避けるために底
面、中央より側面方向に寄ったところに設置された超音
波の送信子と受信子とからなるセンサ（２）と、送信子
に高周波を供給する駆動回路（３）と、受信信号の増幅
部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての偏移周
波数信号を取り出すフィルタ部（６）と、偏移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、 ＦＦＴ部（７）でフーリエ変換した結果を単純加算平均
したスペクトルのパワー値をｆｆｔ（Ｎ）とする周波数
ｆｎのある点Ｎの前後に幅ＭｎをもってＮ−ＭｎからＮ
＋Ｍｎの間で周波数のパワーの平均値Ｍｅａｎ・ｆｆｔ
（Ｎ）を取り、その点Ｎの元来のパワー値ｆｆｔ（Ｎ）
に置き換える。Ｎを零と最大値Ｎｍａｘとの間で同様に
して平均値に置き換えて移動平均すると共に、こうして
移動平均したスペクトルのピーク山に基づいて流速・換
算することを特徴とする超音波流速・流量計。 - 【請求項２】 ピーク山の検出のための比較を周波数の
最大値Ｎｍａｘより始め、第１のピーク山Ｎｐ1 を検出
した後、パワー値Ｍｅａｎ・ｆｆｔ（Ｎ）がピーク山Ｎ
ｐ1 のパワー値Ｍｅａｎ・ｆｆｔ（Ｎｐ1 ）に対する比
率で一定値以下になったところで比較を中止し、第１の
ピーク山Ｎｐ1 の周波数を流速に換算することを特徴と
する請求項１記載の超音波ドップラー流速・流量計。 - 【請求項３】 前記一定値を１／４〜３／４の間に定め
たことを特徴とする請求項２記載の超音波ドップラー流
速・流量計。 - 【請求項４】 液体が流れる水路（１）と、該水路
（１）の底面、中央またはごみや泥を避けるために底
面、中央より側面方向に寄ったところに設置された超音
波の送信子と受信子とからなるセンサ（２）と、送信子
に高周波を供給する駆動回路（３）と、受信信号の増幅
部（４）と、受信信号と送信信号を検波する検波回路
（５）と、受信信号と送信信号の差信号としての変移周
波数信号を取り出すフィルタ部（６）と、変移周波数信
号をスペクトル変換するスペクトル変換部としてのＦＦ
Ｔ部（７）と、ＦＦＴ部（７）で得たスペクトルから流
速換算する演算部（８）と、必要に応じて設けられる測
定結果の表示部（９）とを具備し、 ＦＦＴ部（７）でフーリエ変換した結果を単純加算平均
したスペクトルのパワー値をｆｆｔ（Ｎ）とする周波数
ｆｎのある点Ｎの前後に幅ＭｎをもってＮ−ＭｎからＮ
＋Ｍｎの間で周波数のパワーの平均値Ｍｅａｎ・ｆｆｔ
（Ｎ）を取り、その点Ｎの元来のパワー値ｆｆｔ（Ｎ）
に置き換える。Ｎを零と最大値Ｎｍａｘとの間で同様に
して平均値に置き換えて移動平均すると共に、こうして
移動平均したスペクトルの分布曲線を微分または差分
し、曲線が上に凸の変曲点のうち、Ｎが最大のＮｐ1 を
流速スペクトルのピーク山と判断して、該ピーク山の周
波数を流速に換算することを特徴とする超音波ドップラ
ー流速・流量計。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP02466297A JP3752340B2 (ja) | 1997-02-07 | 1997-02-07 | 超音波ドップラー流速・流量計 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP02466297A JP3752340B2 (ja) | 1997-02-07 | 1997-02-07 | 超音波ドップラー流速・流量計 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10221141A true JPH10221141A (ja) | 1998-08-21 |
| JP3752340B2 JP3752340B2 (ja) | 2006-03-08 |
Family
ID=12144367
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP02466297A Expired - Fee Related JP3752340B2 (ja) | 1997-02-07 | 1997-02-07 | 超音波ドップラー流速・流量計 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3752340B2 (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2005083372A1 (ja) * | 2004-02-27 | 2007-11-22 | 富士電機システムズ株式会社 | パルスドップラ方式と伝搬時間差方式の両方式対応型超音波流量計、同流量計において測定方式を自動選択する方法およびプログラム、同流量計用の電子装置 |
| JP2008512653A (ja) * | 2004-09-03 | 2008-04-24 | ネフロス・インコーポレーテッド | ドップラー方式流速測定装置 |
| CN112639412A (zh) * | 2019-03-14 | 2021-04-09 | 欧姆龙株式会社 | 流量测量装置 |
| CN120405177A (zh) * | 2025-07-03 | 2025-08-01 | 杭州开闳流体科技有限公司 | 基于多指标融合的抑制气泡干扰声学多普勒测流方法、装置及其可读存储介质 |
-
1997
- 1997-02-07 JP JP02466297A patent/JP3752340B2/ja not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2005083372A1 (ja) * | 2004-02-27 | 2007-11-22 | 富士電機システムズ株式会社 | パルスドップラ方式と伝搬時間差方式の両方式対応型超音波流量計、同流量計において測定方式を自動選択する方法およびプログラム、同流量計用の電子装置 |
| JP2008512653A (ja) * | 2004-09-03 | 2008-04-24 | ネフロス・インコーポレーテッド | ドップラー方式流速測定装置 |
| CN112639412A (zh) * | 2019-03-14 | 2021-04-09 | 欧姆龙株式会社 | 流量测量装置 |
| CN120405177A (zh) * | 2025-07-03 | 2025-08-01 | 杭州开闳流体科技有限公司 | 基于多指标融合的抑制气泡干扰声学多普勒测流方法、装置及其可读存储介质 |
| CN120405177B (zh) * | 2025-07-03 | 2025-09-23 | 杭州开闳流体科技有限公司 | 基于多指标融合的抑制气泡干扰声学多普勒测流方法、装置及其可读存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3752340B2 (ja) | 2006-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6487916B1 (en) | Ultrasonic flow metering system | |
| US8452551B2 (en) | Method and apparatus for monitoring multiphase fluid flow | |
| US8151653B2 (en) | Coriolis flowmeter | |
| WO2005083372A1 (ja) | パルスドップラ方式と伝搬時間差方式の両方式対応型超音波流量計、同流量計において測定方式を自動選択する方法およびプログラム、同流量計用の電子装置 | |
| JP2006078362A (ja) | 同一軸型ドップラー超音波流速計 | |
| JPH039405B2 (ja) | ||
| US6595070B1 (en) | Acoustic flow meters | |
| JPH10221141A (ja) | 超音波ドップラー流速・流量計 | |
| DK0762086T3 (da) | Fremgangsmåde til ultralyd-måling af gennemstrømningsmængder af strømmende fluider | |
| JP2890188B2 (ja) | 流体速度演算方法およびその装置 | |
| JP2910815B2 (ja) | 超音波流速計を用いた管内流速の測定方法 | |
| JPH07260532A (ja) | 超音波流量計 | |
| JP2723291B2 (ja) | 超音波センサ | |
| JPH11183215A (ja) | 超音波ドップラー流速計 | |
| JPH09229734A (ja) | 超音波ドップラー流量計 | |
| JP2000180230A (ja) | 流量計及び流量測定方法 | |
| JP2000009508A (ja) | 超音波ドップラー流速計 | |
| JPS58151564A (ja) | 超音波流速計 | |
| JPH0587608A (ja) | 流体速度測定方法およびその装置 | |
| JP5043514B2 (ja) | 超音波式渦流量計 | |
| KR100559139B1 (ko) | 상호상관을 이용한 초음파 유량계 | |
| JP3036800B2 (ja) | 渦流量計 | |
| JP3116189B2 (ja) | 超音波式渦流量計 | |
| JPH0324607B2 (ja) | ||
| JPH04249716A (ja) | ポンプ吐き出し流量測定装置および測定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050822 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050830 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051027 |
|
| 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: 20051129 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20051212 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081216 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091216 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091216 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101216 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101216 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111216 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121216 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |