JP2873443B2 - 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法 - Google Patents

開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法

Info

Publication number
JP2873443B2
JP2873443B2 JP7342767A JP34276795A JP2873443B2 JP 2873443 B2 JP2873443 B2 JP 2873443B2 JP 7342767 A JP7342767 A JP 7342767A JP 34276795 A JP34276795 A JP 34276795A JP 2873443 B2 JP2873443 B2 JP 2873443B2
Authority
JP
Japan
Prior art keywords
converter
ultrasonic
flow
measuring
open channel
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.)
Expired - Fee Related
Application number
JP7342767A
Other languages
English (en)
Other versions
JPH09189590A (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.)
CHANMIN CO Ltd
Original Assignee
CHANMIN CO Ltd
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
Priority to US08/573,637 priority Critical patent/US5780747A/en
Priority to DE19548433A priority patent/DE19548433C2/de
Application filed by CHANMIN CO Ltd filed Critical CHANMIN CO Ltd
Priority to JP7342767A priority patent/JP2873443B2/ja
Publication of JPH09189590A publication Critical patent/JPH09189590A/ja
Application granted granted Critical
Publication of JP2873443B2 publication Critical patent/JP2873443B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related 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/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
    • 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/002Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow wherein the flow is in an open channel
    • 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/662Constructional details
    • 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/704Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • G01F1/708Measuring the time taken to traverse a fixed distance
    • G01F1/7082Measuring the time taken to traverse a fixed distance using acoustic detecting arrangements
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Acoustics & Sound (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measuring Volume Flow (AREA)

Description

【発明の詳細な説明】
【0001】
【発明が属する技術分野】本発明は、開水路用超音波流
速、流量測定技術と関連されたもので、大河川、大人工
開水路の流量を測定するのに主に活用される開水路用超
音波流量測定装置 、開水路用超音波流量測定方法、及び
開水路用超音波流量測定装置を修正しテストする方法に
関するものである。
【0002】
【従来の技術】開水路の流量を測定するために、超音波
にて水流れの水平平均流速を測定して、流量を算出する
方法が広く知られている。このような超音波流量測定装
置は、水流れの水平平均流速測定回線(channel)が一
つにて成されたものもあり、多回線(multi-channel)
にて成されたものもある。
【0003】例えば、Ultraflux社において、
開水路用多回線超音波流量計UF−210000が開発
され販売されている。英国ロンドン市を流れるテームス
江下流にも超音波流量測定所がある。
【0004】図1と図2に、広く知られている前記した
開水路超音波流量測定の基本原理図を図示した。図1に
示したところのように、河川の側面(麓)に沿って一定
なる間隔hにて、超音波発信受信子(以後転換器(Tran
sducer)と略称する)(1A、1B)、(2A、2B)、
(3B、3B)・・・を水流れ方向に対して、一定なる角
度ψi を成すべく設置されている。
【0005】図2に従って、水平線上の平均流速測定方
法を説明すれば次の通りである。転換器(1A)におい
て、転換器(1B)を向いて超音波パルスを発射する瞬
間から、超音波パルスが転換器(1B)に到達する時ま
での時間tABを測定し、同時に、転換器(1B)におい
て(1A)を向いて、超音波パルスを発射して、前記し
た時間tBAを測定する。
【0006】超音波パルスが(1A)から(1B)に、更
に、(1B)から(1A)に伝播する時間tAB、tBAは、
次と同様であると見る。
【数12】 (1) ここで、Lは転換器(1A)から(1B)間の間隔、Cは
川水における超音波流速、
【数13】 水平線上における平均流速、υ超音波パルスが伝播
する軌道線上における流速成分である。
【0007】このように仮定し、流速
【数14】 を次の式によって計算する。(Ultraflux社U
F2100流量計)。
【数15】 (2)
【0008】このような流速計算式にて流速を測定する
方法を、周波数差方法と言う。式(2)においてd/L
=cosψである。
【0009】このように流れの多くの水深において、水
平平均流速を測定し、これに基づいて河川の流量を次式
にて計算する。
【数16】 (3) ここで、Bは測定する水平線における河川水流れの幅で
ある。普通、n=10であれば、十分な正確度で河川流
量を測定するものと見ている。流量Qの計算式は多様で
あり、式(3)よりもっと正確な式もある。いかなる流
量計算式を使用するとしても、基本測定は、多くの水深
における水平平均流速の測定であり、流量測定誤差の大
小は流速測定誤差が大きく左右するものである。
【0010】前記した周波数差流速測定式(2)の外
に、時間差方法式も広く利用されている。時間差式は次
の通りである。
【数17】 超音波流速Cの代わりに、tAB・tBA=L2 2 、cos
ψ=d/Lの関係式を参酌したら次の様になる。
【0011】
【数18】 (4) 式(2)と式(4)は全く同一である。式4で見る様に
2 の項があるのでL測定誤差δL が流速測定誤差に2
δLが加算するのである。
【0012】このため、Lを正確に測定しなければなら
ないが、転換器が開水路の水の中に配置されるので、直
接にはどの様な精密距離測定器で測定すること出来
ない。Lの測定誤差のため水平平均流速測定誤差が大き
くなる場合がたびたびある。
【0013】本発明は、開水路河川の水中に配置されて
いる転換器間の距離Lを超音波で流量計自身が正確に測
できるようにした開水路用超音波流量測定方法と装
置、及び開水路用超音波流量測定装置の修正とテスト方
法を提案するものである。
【0014】本発明を詳細に説明する前に超音波の水流
での伝播軌道に対して説明しなければならない。図3
に、転換器1A(又は1B)から転換器B(又は1A)に
到達する超音波パルス軌道を図示した。該図面に図示し
たところのように、L線方向へ発射された超音波ビーム
は、点線にて示したところのように、転換器(1A)又
転換器(1B)に到達せずに、地点a又はbに到達す
る。
【0015】一方、転換器(1B)に転換器(1A)から
発射された超音波パルスは軌道L1を沿って到達し、又
転換器(1B)から発射した超音波パルスは、軌道L2
沿って転換器(1A)に到達するものである。軌道L1
2は同一でなくL2>L1である。
【0016】河川超音波流量計を設置する時、若しも、
転換器の指向性(directional patーtern)が狭い場合、
図2に示したところのように、転換器をL線と一致され
るべく、即ち、転換器(1 A )に対し転換器(1 B )を水
流の方向と角ψをなすように設置すれば、受信信号が弱
くなり、角ψよりもっと大きくなるように設置すれば、
(図4参照)受信信号が大きくなるということは経験的
に知られている。このような現象が、まさに超音波パル
ス伝播軌道が、図3に示したところのようになるとのこ
とを証明するのである。
【0017】実際の超音波伝播時間は次の様である。
【数19】 ここでL1、L2は幾何学的直線Lでなく伝播軌道の長
さである。又、L2>L1になるのでtBA>tABになるの
である。
【0018】又、L1とL2
【数20】 になる。ここでBは水路の幅、
【数21】 でもっとも速く到達する超音波の発射角である。
【0019】また、次のような関係式が成立することも
証明されている。
【数22】 (5) ここで
【数23】 は線上の流速成分である。式(5)から式(4)が導か
たのである。
【0020】式(5)からLは次のようになる。
【数24】 (6)
【0021】本発明では幾何学的距離を測定するのでは
なく超音波伝播距離L1又はL2を超音波で測定するので
ある。
【0022】そうなので、若しυが大きい時、υ/Cを
無視して、L 1 =Lとおいたならば、少ないけれども
差が発生する。
【0023】例えば、υ=1m/s の時とυ=6m/s の時
(流速が高い開水路である)を対比して見たら次のよう
になる
【数25】 超音波流速C=1425m/sとした時
【0024】従って、υ=1.0m/s の時は、Lの測定
誤差が0.07%になるので、これを無視することは出
来るが、υ=6m/sの時は、0.42%となり、流速測
定誤差が2×0.42=0.84%に増加することにな
り、これを無視することはできない。本発明では流速が
高い時もLを正確に測定できる方法も提供するのであ
る。
【0025】
【発明が解決しようとする課題】従って、本発明は、開
水路用超音波流量測定において、転換器間の距離を正確
に測定して所定式にて水平平均流速を算出することによ
り、流量測定誤差を顕著に減少せしめる、開水路用超音
流量測定装置及びその方法を提供することに、その主
な目的がある。
【0026】本発明の他の目的は、超音波パルスの伝播
軌道の距離を測定して、所定式にて流速を測定するよう
なす、開水路用超音波流量測定装置の修訂とテストの
方法を提供することにある
【0027】
【課題を解決するための手段】本発明による装置は、一
定の少ない間隔に配置された2個の第1、第2の転換器
が一体になっている複合転換器と、これから所定距離L
に配置された第3の転換器で構成した多数の水平平均流
速測定回線と;多数の前記第1、第2及び第3の転換器
が所定間隔にて設置されるべく成した設置台と; 前記第1、第2及び第3の転換器の超音波パルスを発生
せしめる、超音波パルス発振機と; 前記超音波パルスを、前記第1、第2及び第3の転換器
に印加し、それらより、超音波伝播パルスを順次に受信
して、出力すべくなすスイッチング機能を有する交換機
と; 前記交換機より信号を受信して波形整形する正形波成型
部と; 前記第1、第2及び第3の転換器の中のいずれか一つの
転換器と、他の転換器との間で超音波パルスが伝播する
時間間隔を測定する時間間隔測定器と: 所定周期のクラック信号を発生せしめる、同期信号発生
機を設け、前記時間間隔測定器よりの信号を受信して、
断面平均流速を算出して流量を測定する流速流量演算制
御装置にて構成される。
【0028】ここで、複合転換器は、外部のパルス発振
源に連結された中心部が空洞となっている環状形(Annu
lar ring)の第1のケース内に収納された第1の転換器
と、第2のケース内に収納された第2の転換器と、所定
の長さを有するパイプとを有し、前記パイプの末端が前
記第1ケースの中心部に固定され、前記パイプの先端に
第2のケースが所定間隔lに固定されている。
【0029】発明の開水路用超音波流量測定方法は、
第1の転換器(1 Ai )と、この第1の転換器(1 Ai )と
所定間隔を置いている第2の転換器(1 Ci )とを有する
複合転換器を、水路の一辺の側面(麓)に沿って多数設
置すると共に、水路の他の辺の側面(麓)に沿って第3
の超音波転換器(1 Bi )を多数設置することで、多数の
水平流速測定回線を構成するステップと、 前記の第1の
転換器(1 Ai )または第2の(1 Ci )と、第3の転換器
(1 B1 )との間で、その受信信号の振幅が最大となるよ
うにそれらの位置(発射角)を調節するステップと、
記第1の転換器(1 Ai )と第2の転換器(1 Ci )とから
前記第3の転換器(1 Bi )に向かって超音波パルスを発
射して、それぞれから伝播される超音波の 伝播時間t
AiBi 、t CiBi を測定して前記第1の転換器(1 Ai )と第
3の転換器(1 Bi )間の距離Liを測定するステップ
と、前記第3の転換器(1 Bi )から超音波パルスを発射
して前記第1の転換器(1 Ai )まで伝播する伝播時間t
BiAi を測定し、距離Li線上で水平平均流速
【数26】 測定するステップと、前記のように夫々の深さで水平
平均流速
【数27】 を測定し、その結果から水平平均流速分布曲線を作成し
て総平均流速を求め、これに流れ断面Sを乗じて流量を
求めるステップとを含んでいる。
【0030】また、本発明が提案するところの、開水路
用超音波流量測定装置の修正とテストの方法は、所定間
隔lをおいて離れた複合転換器の第1及び第2の転換器
を、水槽の中で各々の地点A及びCに置き、また第3の
超音波転換器を地点Bに置いて、前記第1及び第2超音
波転換器より第3の転換器へ超音波パルスを発射して、
その超音波パルスの伝播時間t AB 及びt CB を測定するス
テップと、地点Bにある第3の転換器を所定距離だけ移
動させて地点B′に置き、前述の如く、前記第1及び第
2超音波転換器より第3の転換器へ超音波パルスを発射
して、その超音波パルスの伝播時間t AB' 及びt CB' を測
定して、間隔lを測定すると同時に、時間間隔測定装置
の時間測定誤差を修正しテストするステップとを含んで
いる。
【0031】従って、本発明によれば、超音波を用いて
河川、開水路又は、規模が大きい閉水路の流量計を構成
するときに、直接測定する水流れの多くの断面線上の平
均流速測定の誤差を除去して、流速測定誤差を顕著に縮
小せしめ、流量測定誤差も減少せしめることができる。
【0032】特に、本発明によれば、河川又は他の開水
路、閉水路両面に配置されている転換器間の間隔の長さ
を測定する必要がなく、流速測定毎に流量計自体が測定
するようになされているので甚だ便利である。
【0033】
【発明の実施の形態】本発明を添付図面に基づいて詳細
に記述すれば、次の通りである。図4に示すように、開
水路の一定なる深さに該当される、開水路両面に、
1、第3の転換器1 A と1 B が配置され、又、超音波パル
スが発射される方向に沿って、前記第1の転換器1 A
距離lをおいて第2の転換器(1 C )を配置する。流速
方向へ、第1の転換器(1 A )が発射した超音波パルス
が第3の転換器(1 B )に到達する時間t AB 、又、第2
の転換器(1 c )が発射して第3の転換器(1 B )に到達
する時間t CB 、さらに、第3の転換器(1 B )が発射し
た超音波パルスが第1の転換器(1 C )と第2の転換器
(1 A )に到達する時間t BC とt BA を測定する。前記
の時間は次式に従い求められる。
【0034】
【数28】 (7)
【0035】前記式(7)において、tCB/tAB を求
めれば、
【数29】 なる。従って、
【0036】
【数30】 (8) これと同様にtCB/tAB によって、
【数31】 (9) を求めることができる。
【0037】流速が少ない時はL1=L又はL2=Lと見
ても良いが流速が高い場合は次の様にもっとも正確なL
値を求める。
【0038】先に式(8)で求められたL1を流速測定
式(4)に代入して流速成分υ′を計算する。式(4)
において
【数32】 にする。この様に求められたυ′は真の流速成分υとは
若干相違している。このυ′を式(6)に代入して正確
なLを求める。
【0039】
【数33】 (10) ここで、C’=L1/(tAB・tBA1/2 である。若
し、υ′/C′値が真のυ/Cに比べて5%の誤差が
り、又L 1 とLに対して1%の誤差をもっているとして
も、式(10)によって求めたL値の誤差は0.05%
となるのでとても正確である。
【0040】この様に測定された転換器間の距離Li
流速流量演算制御装置に入力記憶させたり、又は流速を
測定する時にLi を演算して流速測定式(4)に代入し
ても良い。
【0041】図5及び図6において見るところのように
複合転換器の構造は、環状リング形の圧電セラミックよ
りなる流速測定用転換器(1 A )と、円柱形の圧電セラ
ミックよりなり、前記転換器(1 A )より所定の長さl
の間隔を離れてパイプ(8)の端部近くに設置されてい
る超音波伝播軌道測定用転換器(1 C )とからなる。
【0042】パイプ管(8)は転換器のケース(7A)
を通過して転換器(1C)のケース(7c)と連結され
ている。ケース(7c)はパイプ管(8)に、例えばネ
ジで連結されていて、このパイプ管(8)に沿って移動
させることができ、転換器(1 A )と、転換器(1 C )と
の間隔lを調節することができる。
【0043】超音波発射効率を高めるために、転換器
(1A)と(1C)後面に緩衝材(demper)(6)が充満
されている。例えば、緩衝材として酸化タングステン粉
末を使用しても構わない。
【0044】図5及び図6において(9)は、転換器
(1A)と(1C)を連結する導線であり、(10)は、
複合転換器を河川の側面(麓)に設置した転換器固定台
に固着せしめる固着環である。転換器(1A)の内径D
A は、転換器(1C)のケース外径DC に比して、やや
大きくなければならない。
【0045】河川の幅が数十、数百メートルにもなる
合、超音波パルスの振動周波数は、普通200−500
kHz を使用するところ、水における波長は、わずか7〜
3mm程度になる。従って、転換器(1A )の前に3〜
7mm以上になる障害物があってはならない。しかし、
図5に示したところのように、2個の超音波転換器より
なる複合構造とすれば、転換器(1A)と(1C)が超音
波パルスを発射、受信するのに支障がない。
【0046】転換器(1A)と(1C)間の間隔lの選
は流量測定装置に内蔵されている時間間隔測定装置の誤
差によって決定される。式(4)で見るようにLの測定
誤差の如何によって流速測定誤差が左右される。間隔l
の選択関係式は次の通りである。
【数34】 (11) ここで、τは時間間隔測定絶対誤差、 max 河川水に
おいて期待することができる最大の超音波流速、δ L
L間隔測定許容相対誤差である。
【0047】例えば、τ=10-7 S、Cmax=1460m/
s,δl=0.001(0.1%)とすれば、これらの数
値を式(11)に代入すると、
【数35】 になる。従って、l=0.3mにて選択すれば十分であ
る。
【0048】図7には、一つの実例にて、河川用超音波
流量測定装置の構成を、概略的なブロック線図で図示し
た。流速流量演算制御装置(16)は、クロックパルス
を発生せしめる同期発振機(syncqenerator)(12)
と、それによって作動制御される無線送受信機17を設
け、以後、詳細に記述されるところのように、河川水の
流速及び流量を演算する。
【0049】超音波パルス発振機(13)は超音波パル
ス信号を発生せしめ、交換機(11)に印加する。交換
機(11)は超音波パルス信号を河川の一辺の側面
(麓)に設置した流速測定装置にて構成した多数の複合
転換器(1A1、1C1)、(1A2、1C2)...(1An、1
Cn )に超音波パルス信号を印加する。
【0050】パルス整形機(14)は以後記述されるよ
うに、交換機(11)からの受信信号を整形波に成形す
る。時間間隔測定装置(15)は、パルス信号発振機
(13)と、パルス整形機(14)からの信号間の時間
間隔、即ち、超音波伝播時間を測定して、流速流量演算
制御装置(16)に入力する。
【0051】同様に、他側河辺の側面(麓)にも、転換
器(1B1),(1B2)..(1Bn )が設置され、これと関
連して、交換機(11)、同期発生機(12)、パルス
発振機(13)、波形整形機(14)、と無線送受信機
(18)が設置される。
【0052】って、該多回線超音波流量測定装置は、
所定のアルコリズムとプログラムに従って、交換機(1
1)が複合転換器(1A1、1C1)、(1A2、1C2)・・
・(1An、1Cn )に超音波パルス信号を印加して、超音
波パルスが転換器(1B1)、(1B2)・・(1Bn )に伝
播するようにし、流速流量演算制御装置(16)によっ
て作動する無線送受信機(17)は、これを無線送受信
機(18)に通知することにより、交換機(11)は転
換器(1B1)、(1B2)・・(1Bn )が伝播信号を受信
するように成し、該信号を交換機(11)が波形整形機
(14)等を経由して、無線送受信機(18)が無線送
受信機(17)に伝達すべく成す。
【0053】この時、無線送受信機(17)は、流速流
演算制御装置(16)の制御によって、無線送受信機
(18)、交換機(11)及び超音波パルス発振機(1
3)を作動せしめ、転換器(1B1)、(1B2)・・(1
Bn より超音波パルスが転換器(1A1、1C1)、
(1A2、1C2)・・・(1An、1Cn に超音波パルスを
発射する。
【0054】って、転換器(1A1、1C1)、(1A2
1C2)・・・(1An、1Cn )によって受信された信号
交換機(11)を通じて波形整形機(14)におい
て整形波になり、時間間隔測定装置(15)に入力され
る。時間間隔測定装置(15)は、転換器1A1、1B1
C1に該当する伝播時間tAB1、tBA1、tCB1、tBC1
コード信号を発生させ、流速流量演算制御装置(16)
に入力する。流速流量演算制御装置(16)は、第1回
線水平平均流速1を演算して記憶する。同様に、該過
程が第2回線、第3回線、第n回線とに繰り返され、
【数36】 流速が測定記憶される。(19)は転換器(1Ai)、
(1Bi)、(1Ci)を設置固定せしめる設置台である。
【0055】流量Qは、従前に使用していた式(3)に
よって計算され得ることができるが、もっと正確な方法
は、
【数37】 測定データを利用して水平流速分布曲線を作成して総
平均流速 s を求め、これを次式のように、
【数38】 (12)近似積分し、 予め求めておいた水深と水流れ断面積Sと
の関係函数S(H)によって求められる断面積Sに乗じ
流量Qを求める方法がもっとも正確である。即ち、
【数39】 (13) である。図7に、水深Hを測定するレベルメートルは表
示しなかった。
【0056】以上におけるように、広い河川用の超音波
流量測定装置を構成するにおいて、両辺に配置された転
換器を高周波ケーブルにて第2次測定器(11−16)
と連結するのが非合理的である時、無線送受信機(1
7)、(18)を使用するのが好ましい。勿論、河川両
辺の側面(麓)に、転換器を配置せずに一方の側面
(麓)に、所定された間隔dにて転換器(1A1)、(1
C1)・・・、(1An)、(1Cn、さらに(1 B1 )、
(1 B2 )、・・・(1 Bn )を設置し、反の辺の側面
(麓)に超音波反射装置を装置するのも甚だ効果的であ
る。このような場合には、無線送受信器も必要無く同期
発振機(12)、超音波パルス発振機(13)、波形成
形機(14)、交換機(11)を二重に両方の辺の側面
(麓)に設置する必要も無い。
【0057】又、流量を測定する毎に転換器(1Ci)を
作動させてL1を測定しなくても良い。周期的にL1を測
定して流速流量演算制御装置(16)に入力記憶してお
いても良い。
【0058】又、多数の複合転換器を設置する代わり
に、開水路の一辺の側面(麓)に沿って転換器
(1 B1 )、(1 B2 )…(1 Bn )を設置した後、他の辺の
側面(麓)に沿って転換器(1 A1 )、(1 A2 )……(1
An )を設置する位置点で、複合転換器(1 A )、
(1 C )を移動させながら、(L 1 )、(L 2 )…(L
n )を測定し、その測定値を流量計の流速流量演算制御
装置に入力記憶させ、転換器(1 A1 )、(1 A2 )…(1
An )を設置することもできる。これはL1 が絶対に変わ
らない時に適用する。
【0059】また本発明は、複合転換器を構成する2個
の転換器間の間隔を修正しテストする方法と超音波伝播
時間測定装置と流速流量演算装置を修正しテストする方
法を提案している。
【0060】式(8)又は(9)を見ると複合転換器の
中に入っている転換器間の間隔lを正確に測定しなけれ
ばL測定誤差が大きくなって流速測定誤差も大きくな
る。従って、lを正確に測定する方法と又時間間隔測定
装置の誤差を修正しテストする方法が必要である。
【0061】勿論、精密な誤差が±10-6m程度の長さ
定器はあるが、このような長さ測定器にて式(8)に
代入される長さlを正確に測定することは難しい。
【0062】その理由は、長さlとは超音波パルスが伝
播する長さであるが、いくら超音波パルス送受信回路と
時間間隔測定回路を精密に構成したとしても、色々な原
因にて遅延時間が含まれる。更に、いくら遅延時間が小
さく無視することができる程度であったとしても、図5
に示した複合転換器、すなわち転換器(1Ai)と
(1Ci)の超音波パルス発射受信有効面間の間隔測定を
正確に行うことは難しい。
【0063】前記したことを参酌して、本発明に伴うl
測定の修正方法は次の通りである。図8に示すように、
水槽を周囲の温度が変わらない室内に設置し、この水槽
に水を充満させる。次に流量測定装置に連結された
(1A)と(1Cとからなる複合転換器を水中に設置し
固定し、転換器(1A)よりL1になる間隔に転換器(1
B)を設置する。この時、L1 間隔を正確に測定する必
要が無い。又、予め、長さ測定計にてl’を大略測定し
て流量計に入力せしめる。
【0064】流量計に内蔵された時間間隔測定装置の時
間間隔測定誤差をテストするために、精度(accuracy)
が高い基準時間間隔測定装置を並列に連結する。
【0065】このように事前準備を行い、流量計を作動
させ、1 間隔を水槽において超音波パルスが伝播する
時間tAB1、tCB1を繰り返し測定しながら、流量計が算
出する間隔L’1を記録する。即ち、
【数40】 (14) ここで、l’は真の長さlと異なるために、式(14)
によって計算されたL’1 は、真の長さL1 と異なるも
のである。
【0066】次に超音波転換器(1B)を正確に測定さ
れたL0間隔に移動する。移動した間隔L0 の正確な測
定は難しいものでは無い。転換器(1A)と(1B)間の
間隔は、Li 0 であり、転換器(1C)と(1B)間
の間隔はL 1 +L 0 −lとなる。
【0067】このように超音波転換器(1B)の位置を
変更し、超音波伝播時間t AB2 、t CB2 を測定して、超音
波転換器(1 A )と(1 B )間の間隔
【数41】 (15) を求める。式(15)から式(14)を減じれば、転換
器(1 C )と(1 B )間の間隔
【0068】
【数42】 (16)を求めることができる。
【0069】時間間隔測定誤差は甚だ小さいので、無視
すべき程度であれば、真の長さL0は、
【数43】 (17)となる。 ここで、lは求めようとする真の長さである。
従って、
【数44】 になり、
【数45】 (18) になるであろう。
【0070】このように、式(18)に従って、正確な
さlを求めることができる。このように求められたl
の値をl’代わりに流量測定装置の演算制御部に入力
、式(17)によってL0 を繰り返し測定する。繰り
返して測定した結果を分析すれば、流量測定装置が測定
した長さと、精密な長さ計測機にて測定した長さL0
の偏差が、偶然偏差であるか、あるいは固定偏差が含ま
れているのかが、容易に知り得る。固定(sysytemati
c)偏差が含まれているならば、前記した方法を繰り返
、更に式(18)によってもっと正確なlの値を求
め、流量計に入力、前記した測定を繰り返す。
【0071】このような過程を固定偏差が無くなるま
で、繰り返して最終的にl値を求め、流量測定装置の演
算制御装置に入力する
【0072】前記した長さlの測定方法の長点は、流量
測定装置と超音波転換器、更に両者間を連結するケーブ
ル線のすべてのものを連結して見る時、発生する遅延時
間等を補償するlの値を求めることが重要である。
【0073】前記した方法にて求められたlの値を流量
測定装置に入力、超音波転換器(1B)の位置を変更
させながら、流量測定装置が測定した変更された長さの
値と、精密長さ計測機にて測定した値と対比して、超
音波パルス伝播の長さ測定誤差をテストするようにな
る。同時に流量測定装置の時間間隔測定装置の時間間隔
測定誤差も、基準時間間隔測定装置と対比してテスト
れる。入力せしめたlの誤差、時間間隔測定誤差等を、
前記した水槽においてテストできるようになり、一定な
るL測定誤差、tAB、tCB、tBA、tBC測定誤差がある
場合の流速測定式(4)により測定される流速の誤差が
δviを評価するようになる。
【0074】ところが流速測定誤差を評価する時、次
状況に注目しなければならない。水槽で流量計が測定し
たL0の絶対誤差ΔL0で計算した相対誤差ΔL0/L0
δL0 現地開水路で測定する誤差ではない。L0はほぼ
1〜3mであるが現地では数十米になるのでLi測定相
対誤差はδL0より少なく、δ Li=ΔL0Liとなる。
【0075】超音波伝播時間測定誤差も水槽で求めた絶
対誤差Δtを現地で期待する伝播時間で除けば良い。δ
t=Δt/tこの様に計算したδLとδtを使って流速測定
誤差を評価するのである。
【0076】次は、流量計算を行うための積分式(1
2)を遂行する、演算制御部の誤差修正は、次のように
簡単に行うことができる。水平平均流速
【数46】 の多様な分布関数を作成して総平均流速
【数47】 を予め正確に積分計算しておく。次に流速分布曲線にお
いて、転換器が配置される地点における流速
【数48】 に該当する流速と、それらの間の間隔を流量測定装置の
流量演算部に入力し、そ の演算結果
【数49】 を正確に計算された
【数50】 と対比して、
【数51】 測定誤差を修正する。
【0077】超音波流量測定装置の最大流量測定相対誤
差は、
【数52】 (23) になり、
【数53】 が主に偶然誤差にて構成されているとすれば、次のよう
に評価される。
【数54】 (24) ここでδSは水流れの断面積誤差である。このように、
本発明の方法によれば、流速を測定する条件において、
河川用、開水路用、又は、規模が大きい閉水路用超音波
多回線流量計の流量測定誤差を間接的に非流入方法で、
高精度で修正することができるものである。
【0078】産業用流量計等は、標準流量測定装置にて
直接修正、テストすることができるけれども、開水路の
大流量を測定する流量計は間接的な修正テストにより測
定誤差を修正するより外に方法がないところ、正確な超
音波流速測定方法を使用することにより間接的な修正テ
ストの精度を高めることができる。
【図面の簡単な説明】
【図1】従来の多回線超音波河川流量測定の原理を説明
した断面図。
【図2】従前の水平平均流速測定の原理を説明した模式
図。
【図3】超音波が伝播する軌道の様子を示した模式図。
【図4】本発明に伴う水平平均流速測定方法を説明する
模式図。
【図5】本発明に伴う水平平均流速測定用複合超音波転
換器の断面図。
【図6】超音波伝播軌道測定用転換器(IC)を示す斜
視図。
【図7】本発明に伴う流量測定装置の一実施例を示す断
面図。
【図8】本発明の修正テスト方法を説明する模式図。
【符号の説明】
A ・・・第1転換器 1C ・・・第2転換器 1B ・・・第3転換器

Claims (5)

    (57)【特許請求の範囲】
  1. 【請求項1】 人工開水路又は河川の両側面(麓)の線
    上に沿って2地点間の水平平均流速を測定する多数の超
    音波平均流速測定回線により構成されている開水路用超
    音波流量測定装置において、 開水路の一辺の側面(麓)の線上に配置する超音波転換
    器(以下超音波は略する)は、所定間隔lをもって配置
    された2個の転換器(Transducer)が一体になっている
    複合転換器となっており、この複合転換器は、中心部が
    空洞となっている環状形(Annular ring)の第1のケー
    ス内に収納された第1の転換器(1 A )と、 第2のケース内に収納された第2の転換器(1 C )と、
    所定の長さを有するパイプとを有し、 前記パイプの末端が前記第1ケースの中心部に固定さ
    れ、前記パイプの先端に第2のケースが所定間隔lを微
    調整可能に固定されており、前記所定間隔lは、次の式
    により選択されることを特徴とする開水路用超音波流量
    測定装置 。 【数1】 ここでτ−超音波多回線流量計の超音波伝播時間測定装
    置の時間測定絶対誤差、Cmax −開水路水の中で期待す
    る最大の超音波、δL −水平平均流速測定用の超音波転
    換器間の距離測定許容相対誤差である。
  2. 【請求項2】 開水路用超音波流量測定方法において、第1の転換器(1 Ai )とこの第1の転換器(1 Ai )と所
    定間隔lをもって配置された第2の転換器(1 Ci )とが
    一体になっている複合転換器を開水路の一辺の側面
    (麓)に沿って設置すると共に、他の辺の側面(麓)に
    沿って、前記複合転 換器(1 Ai )(1 Ci )と結ぶ線が水
    平となり、かつ水流の方向と角Ψをなすように、第3の
    転換器(1 Bi )を設置することで、水平平均流速測定回
    線を構成するステップと、 前記第1の転換器(1 Ai )または第2の転換器(1 Ci
    と、第3の転換器(1 Bi )との間で、超音波が伝達され
    るとき、その受信信号の振幅が最大になる様に前記第
    1、第2及び第3の転換器(1 Ai )(1 Ci )(1 Bi )の
    発射角を調節するステップと、 前記第1の転換器(1 Ai )と第2の転換器(1 Ci )とか
    ら前記第3の転換器(1 Bi )に向かって超音波パルスを
    発射して、それぞれから伝播される超音波の伝播時間t
    AiBi 、t CiBi を測定して前記第1の転換器(1 Ai )と第
    3の転換器(1 Bi )間の距離Liを測定するステップ
    と、 前記第3の転換器(1 Bi )から超音波パルスを発射して
    前記第1の転換器(1 Ai )まで伝播する伝播時間t BiAi
    を測定し、距離Li線上で水平平均流速 【数2】 を測定するステップと、 前記の様に夫々の深さで水平平均流速 【数3】 を測定し、その結果から水平平均流速分布曲線を作成し
    て総平均流速を求め、これに流れ断面を乗じて流量を求
    めるステップとを含むことを特徴とする開水路用超音波
    流量測定方法尚、前記L i は、流速が高く無い時は次
    の式によって測定する; 【数4】 (a) (ここでlは第1と第2の転換器間の間隔)又流速が高
    い(υ> 3 m/s)ところでは次の式で測定する; 【数5】 (b)ここで、Cは、下記式で計算した超音波流速であり、 【数6】 υ′は、次式で計算されるL i 線上の流速成分であり、 【数7】 (c) この様に測定したLi を下記の流速測定式に代入して流
    速を演算する; 【数8】 (d) ここで di=Licosψ.
  3. 【請求項3】 距離L i と流量を毎回測定せず、一度L i
    を測定して流量計の流速流量演算制御装置に入力記憶し
    ておき、流速を測定する時は、超音波伝播時間t AiBi
    BiAi のみ測定し、L i は一定の周期で測定し、前に記
    憶させたL i 値と比較し修正することを特徴とする請求
    項2記載の開水路用超音波流量測定方法
  4. 【請求項4】 第1の転換器(1 Ai )と第3の転換器
    (1 Bi )との間の距離Liが変わらないことを条件とし
    て、開水路の一辺の側面(麓)に沿って転換器
    (1 B1 )、(1 B2 )…(1 Bn )を設置した後、他の辺の
    側面(麓)に沿って転換器(1 A1 )、(1 A2 )……(1
    An )を設置する位置点で、複合転換器(1 A )、
    (1 C )を移動させながら、(L 1 )、(L 2 )…(L
    n )を測定し、その測定値を流量計の流速流量演算制御
    装置に入力記憶させ、転換器(1 A1 )、(1 A2 )…(1
    An )を設置することを特徴とする請求項2記載の開水路
    用超音波流量測定方法
  5. 【請求項5】 水平平均流速測定回線を構成する、第1
    の転換器(1 A )と第2の転換器(1 C )でなっている複
    合転換器と、第3の転換器(1 B )とを、水槽に一定な
    間隔L A に配置し、前記転換器(1 A )と(1 C )間の間
    隔を長さ測定器で測定し、その値l′を流量計の流速流
    量演算制御装置に入力するステップと、 前記流量計の超音波伝播時間測定装置と並列して精密な
    基準時間間隔測定装置を連結するステップと 前記流量計を始動させて前記転換器(1 A )、(1 C )か
    ら(1 B )、又(1 B )から( 1 C )と(1 A )に超音波
    が伝播する時間t AB 、t CB 、t BA 、t BC を測定して間隔
    L′ A を次式で算出するステップと、 【数9】 (a) 前記転換器(1 B )を精密に測定したL 0 の間隔に移動し
    固定して前記の伝播時間を測定してL’ A +L’ 0 を次式
    により算出し、 【数10】 (b) 前記L’ A +L’ 0 からL’ A を引いてL’ 0 を得ること
    で、正確なl値を次式に従って算出するステップと、 【数11】 (c) 前記流速流量演算制御装置に入力されたl′に代えて、
    修正された正確なl値を前記流速流量演算制御装置に再
    入力し、前記の方法でL 0 を測定し、測定したL 0 の値
    L″ 0 が正確に測定された値と等しいかどうかをテスト
    し、 もしL″ 0 の値が許容範囲の誤差を越えるなら、前記式
    (C)にL″ 0 値を代入してもっと正確なlの値を計算
    するステップと、 前記の方法を繰り返して最終的なl値を求め、前記流速
    流量演算制御装置に入力するステップとを含むことを特
    徴とする開水路用超音波流量測定装置を修正しテストす
    る方法。
JP7342767A 1995-12-18 1995-12-28 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法 Expired - Fee Related JP2873443B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US08/573,637 US5780747A (en) 1995-12-18 1995-12-18 Open channel multichannel ultrasonic flowrate measurement apparatus and method
DE19548433A DE19548433C2 (de) 1995-12-18 1995-12-22 Multi-Kanalultraschalldurchflussmessvorrichtung und Verfahren zur Messung der Durchflussgeschwindigkeit eines Flusses oder eines offenen Kanals
JP7342767A JP2873443B2 (ja) 1995-12-18 1995-12-28 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US08/573,637 US5780747A (en) 1995-12-18 1995-12-18 Open channel multichannel ultrasonic flowrate measurement apparatus and method
DE19548433A DE19548433C2 (de) 1995-12-18 1995-12-22 Multi-Kanalultraschalldurchflussmessvorrichtung und Verfahren zur Messung der Durchflussgeschwindigkeit eines Flusses oder eines offenen Kanals
JP7342767A JP2873443B2 (ja) 1995-12-18 1995-12-28 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法

Publications (2)

Publication Number Publication Date
JPH09189590A JPH09189590A (ja) 1997-07-22
JP2873443B2 true JP2873443B2 (ja) 1999-03-24

Family

ID=27215773

Family Applications (1)

Application Number Title Priority Date Filing Date
JP7342767A Expired - Fee Related JP2873443B2 (ja) 1995-12-18 1995-12-28 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法

Country Status (3)

Country Link
US (1) US5780747A (ja)
JP (1) JP2873443B2 (ja)
DE (1) DE19548433C2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012098093A (ja) * 2010-10-29 2012-05-24 Hiroshima Univ 音響トモグラフィー計測システム及び音響トモグラフィー計測方法

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100293842B1 (ko) * 1998-06-10 2001-09-17 남상용 하천수평평균유속측정방법
JP2957576B1 (ja) 1998-09-16 1999-10-04 チャン ミン テック カンパニー リミテッド 水深・流速・水温測定装置
US6477900B2 (en) * 2001-01-08 2002-11-12 Jet Sensor Ltd. Fluidic gas metering system
US7117738B2 (en) * 2003-10-02 2006-10-10 Denso Corporation Liquid level detecting apparatus
DE102004031274B4 (de) * 2004-06-28 2007-07-12 Flexim Flexible Industriemesstechnik Gmbh Verfahren zur Kalibrierung von Ultraschall-Clamp-on-Durchflussmessgeräten
CN1828232B (zh) * 2006-04-13 2010-05-12 李佳宾 二线能坡流量计量法
US8397919B2 (en) * 2009-03-16 2013-03-19 Bruce H. Kittrick Continuous gravity assisted ultrasonic coal cleaner
IN2014DN10384A (ja) 2012-05-30 2015-08-14 Rubicon Res Pty Ltd
CN105051503B (zh) * 2012-10-11 2018-11-16 鲁比康研究有限公司 流量测量
US10161770B2 (en) 2016-06-30 2018-12-25 Ott Hydromet Gmbh Flow meter with adaptable beam characteristics
US10408648B2 (en) 2016-06-30 2019-09-10 Hach Company Flow meter with adaptable beam characteristics
US10295385B2 (en) 2016-06-30 2019-05-21 Hach Company Flow meter with adaptable beam characteristics
WO2018081503A1 (en) 2016-10-27 2018-05-03 Ysi, Inc. Ability to measure distance between stations using dgps/rtk gps in the velocity-area method (stationary) using acoustic doppler current profile
WO2019083662A1 (en) * 2017-10-27 2019-05-02 Daniel Measurement And Control, Inc. ADJUSTABLE TRANSDUCER ASSEMBLIES
CN108195436A (zh) * 2018-03-28 2018-06-22 上海中核维思仪器仪表有限公司 带自校准功能的***式气体超声波流量计测量装置及方法
CN109975578B (zh) * 2019-04-09 2021-12-31 武汉新烽光电股份有限公司 径流流速实时测量方法、***及控制装置
CN112903043B (zh) * 2019-12-03 2022-06-14 成都秦川物联网科技股份有限公司 一种多声道超声波流量计***
CN113917185A (zh) * 2021-10-22 2022-01-11 浙江天禹信息科技有限公司 一种声学全断面测流仪及其工作方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1181907A (en) * 1966-11-25 1970-02-18 Westinghouse Electric Corp Fluid Flow Measurement System and Method
US3546935A (en) * 1969-11-24 1970-12-15 Westinghouse Electric Corp Fluid measurement system and method
US4028938A (en) * 1976-01-26 1977-06-14 Ocean Research Equipment, Inc. Acoustical flow meter
DE2703439C3 (de) * 1977-01-28 1979-08-09 Danfoss A/S, Nordborg (Daenemark) Vorrichtung zum Messen von physikalischen Größen einer Flüssigkeit mit zwei Ultraschallwandlern
DE4016529C1 (en) * 1990-05-22 1991-11-07 Turbo-Werk Messtechnik Gmbh, 5000 Koeln, De Flowmeter for open duct - has flow sensor movable across channel to supply computer with flow speeds measured at several points
US5315880A (en) * 1992-08-13 1994-05-31 Henry Filters, Inc. Method for measuring fluid velocity by measuring the Doppler frequency shift or microwave signals
DE4320295A1 (de) * 1993-01-25 1994-07-28 Christian Dipl Ing Koelling Verfahren und Vorrichtung zum Messen des Durchflusses in teilgefüllten oder vollgefüllten Rohrleitungen und offenen Gerinnen
KR960013251B1 (ko) * 1993-08-25 1996-10-02 주식회사 창민물산 초음파 유량측정 방법과 장치
KR960003645B1 (ko) * 1993-08-25 1996-03-21 주식회사창민테크놀러지 하천 국부(局部) 유속측정방법 및 장치
CA2107043A1 (en) * 1993-09-27 1995-03-28 Wolfgang Stedtnitz Method and apparatus for measuring water velocities and discharges in rivers and the like
KR960009391B1 (en) * 1993-12-24 1996-07-18 Korea Inst Sci & Tech Apparatus and method of measuring flow velocity

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012098093A (ja) * 2010-10-29 2012-05-24 Hiroshima Univ 音響トモグラフィー計測システム及び音響トモグラフィー計測方法

Also Published As

Publication number Publication date
US5780747A (en) 1998-07-14
DE19548433A1 (de) 1997-06-26
JPH09189590A (ja) 1997-07-22
DE19548433C2 (de) 2003-07-24

Similar Documents

Publication Publication Date Title
JP2873443B2 (ja) 開水路用超音波流量測定装置、開水路用超音波流量測定方法、及び開水路用超音波流量測定装置を修正しテストする方法
US11333676B2 (en) Beam shaping acoustic signal travel time flow meter
JP4979838B2 (ja) 超音波流量測定ユニットの測定変換器を校正するための方法及び装置
EP1736741B1 (en) Ultrasonic method for measuring a flow rate of liquid and/or gaseous media and device for carrying out said method
US20110271769A1 (en) Flow measuring apparatus
JPH10104039A (ja) 多回線流量測定装置
JPS5948327B2 (ja) 流量計
CN101680797B (zh) 距离测量设备及其相关方法
JP2005241546A (ja) ドップラー式超音波流量計、その演算処理装置、プログラム
CA1131342A (en) Acoustic flowmeter with reynolds number compensation
GB2479115A (en) Ultrasonic flow measurement based on propagation time of Lamb waves
EP1726920B1 (en) Method for ultrasonic Doppler fluid flow measurement
CN112964319B (zh) 一种多频阵列超声波流量计
AU2010335057B2 (en) Measuring apparatus
KR0169091B1 (ko) 개수로용 초음파 유량측정 장치 및 방법
JP4827007B2 (ja) 超音波流量計
JPH11351928A (ja) 流量計および流量計測方法
JP2007322186A (ja) 超音波流量計
JP3438371B2 (ja) 流量計測装置
JPS6040916A (ja) 超音波流速・流量計の温度変化誤差の補正法
RU2447406C1 (ru) Ультразвуковой способ измерения расхода жидких и газообразных сред и устройство для его осуществления
JPS638515A (ja) 超音波流量測定装置
JP6354327B2 (ja) 超音波流量計
JP2005241407A (ja) 超音波流量計の校正方法、超音波流速分布計の校正方法、超音波流量計、超音波流速分布計、検出器
RU2139502C1 (ru) Ультразвуковой способ измерения расхода воды в открытом русле и устройство для его осуществления

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees