JP2820448B2 - Volume measurement method and device - Google Patents

Volume measurement method and device

Info

Publication number
JP2820448B2
JP2820448B2 JP22900689A JP22900689A JP2820448B2 JP 2820448 B2 JP2820448 B2 JP 2820448B2 JP 22900689 A JP22900689 A JP 22900689A JP 22900689 A JP22900689 A JP 22900689A JP 2820448 B2 JP2820448 B2 JP 2820448B2
Authority
JP
Japan
Prior art keywords
volume
tank
main tank
signal
pressure
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 - Lifetime
Application number
JP22900689A
Other languages
Japanese (ja)
Other versions
JPH0392725A (en
Inventor
潤 西野
Original Assignee
株式会社カンセイ
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 株式会社カンセイ filed Critical 株式会社カンセイ
Priority to JP22900689A priority Critical patent/JP2820448B2/en
Publication of JPH0392725A publication Critical patent/JPH0392725A/en
Application granted granted Critical
Publication of JP2820448B2 publication Critical patent/JP2820448B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION 【産業上の利用分野】[Industrial applications]

この発明は、タンク内に収納された被測定物である液
体等の体積(容積)を測定する体積測定方法及びその装
置に関するものである。
The present invention relates to a volume measuring method for measuring a volume (volume) of a liquid or the like, which is an object to be measured, stored in a tank and an apparatus therefor.

【従来の技術】[Prior art]

従来のこの種の体積測定方法及びその装置としては、
先行技術として例えば特願平1−27808号に示されるも
のがあるので、それを第3図乃至第7図に基づいて説明
する。 まず第3図において、原理説明を行なうと、30は例え
ば液体,粉体,異形状物体等の被測定物を収納する異形
状のメインタンクであって、このメインタンク30には連
結パイプ32を介して該メインタンク30よりも容積の小さ
い補正タンク31が連結されている。また上記異形状メイ
ンタンク30の上部には小径の通気孔35が穿設されている
(なお、該通気孔35は設けられてなくてもよい)。上記
補正タンク31の上部には例えばピストン,スピーカ形状
(以下これをスピーカと云う)等の体積変化手段(機
構)33が設けられていて、この体積変化手段33の動作に
よって補正タンク31内の体積を変化させることができる
ようになっている。以上が本実施例の構成であって、次
にその構成による測定原理を説明する。 測定原理 (1) まず、小体積V1の補正タンク31と、大体積V2
メインタンク30とが連結される第3図に示すような連結
タンクシステムを考える。これはメインタンク30と補正
タンク31が流量抵抗r1のパイプ32によって接続されたも
ので、タンク30の通気孔35は流量抵抗r2である。双方の
タンク30,31内の気体の比熱比をγ、気体定数をR、タ
ンク31の熱時定数をτとする。タンク31にピストン,ダ
イヤフラム,ベローズ,スピーカ形状等を用いた体積変
化手段33を取り付け、この体積変化手段33によって実際
に発生する体積変化量をv(t)とする。 タンク30,31が剛体の場合、タンク30,31内気圧の加圧
又は減圧時にタンク30,31が歪まないので、ピストン,
ダイヤフラム,ベローズ等による体積変化量V0(t)と
実際に発生する体積変化量v(t)は等しい。もしタン
ク30がフレキシブルな場合、タンク30,31内気圧の加圧
又は減圧時にタンク30が歪むので、その拡張又は収縮に
よる体積変化量に応じた分だけ実際の体積変化量v
(t)は変化する。 v(t)=0のとき、タンク31の中の気体の絶対圧
力,温度,モル数をそれぞれP0,T1,n1,またタンク30の
中の気体の絶対圧力,温度,モル数をそれぞれP0,T2,n2
とする。測定環境が著しく変化しない場合は、通気孔35
を介してタンク30,31内外に気体が流通するのでタンク3
1,30内の絶対圧力P0は外気圧と等しく、その変化は非常
に遅く、外気圧と等しく変化する。 v(t)≠0のとき、圧力,温度,モル数も体積変化
手段33の状況に応じて変化し、 タンク31において、 圧力はP0+ΔP1(t), 温度はT1+ΔT1(t) モル数はn1−Δn12(t)と変化する。 タンク30においては、 圧力はP0+ΔP2(t), 温度はT2+ΔT2(t), モル数はn2+Δn12(t)−Δn2(t)と変化する。Δn
12(t)はタンク31からタンク30に流れた空気のモル
数、Δn2(t)はタンク30から通気孔35を介して外部に
漏れた空気のモル数である。 ここでこのシステムについて次の仮定を設定する。 1)タンク30,31内気体は理想気体である。 2)v(t)≪{V1,V2} 3)タンク30の熱容量は大きく、圧力変化ΔP2(t)に
伴うタンク内温度変化は体積変化量v(t)の変化の速
さに比べ非常に遅く無視できる。 4)体積変化量v(t)の変化の速さは、それに伴い変
化する圧力がタンク30,31内の至るところで等しい程度
とする。 5)タンク30内に被測定物を入れても、この物体により
タンク30内に2つ以上の閉じたガス空間、すなわち空洞
部分が構成されない。 以上の仮定は、大きな制約となるものではない。体積
変化量v(t)に対してΔP1(t),ΔP2(t),ΔT1
(t),ΔT2(t),Δn12(t)、Δn2(t)の変化
は本来、非線形方程式で表されるが、仮定2)よりその
大きさはP0,T1,T2,n1,n2に対して非常に小さく、そのた
め線形方程式で近似できる。静的状態においてタンク3
0,31内の気体の圧力,温度,モル数の関係は次の代数方
程式で表される。 P0V1=n1RT1,P0V2=n2RT2 (1a) また仮定1),3),4),5)より、動的状態においては
タンク30,31内の気体の圧力,温度,モル数の関係は、
次の線形常微分方程式で表される。 流量抵抗r、上式においてはr1,r2はパイプ32の長さ
lと直径dより次の実験式のように求められる。 この式は長さlが50〜650[mm]、直径dが2.0〜9.0
[mm]のアルミ製パイプを用いて実験的に求めたもので
ある。 また、体積変化量v(t)は次のように表される。 v(t)=v0(t)−Δv(t) (1h) ΔVはタンクの材質,形状,容積等から定まるタンク
固有の定数であり、Δv(t)は体積変化手段33の体積
変化量v0(t)の変化に伴うタンク30の拡張または収縮
による体積変化量である。 式(1a)〜(1i)にラプラス変換を施し、入力v
(t)から出力ΔP1(t)までの伝達関数を求めると次
のようになる。 ここで 係数r2V2/RT2,r1V2/RT2,r1V1/RT1,r1V2/RT1は、圧力
変化の時定数である。例えばr2V2/RT2はタンク30におけ
る空洞部分の絶対温度T2の気体が流量抵抗r2の通気孔35
を介してタンク30外に流れるときの圧力減衰の時定数で
ある。補正係数k2(s,r1,r2,V1,V2)は、メインタンク3
0の容積V2により変化するが、第4図に示されるk2(s,r
1,r2,V1,V2)の周波数特性において、適当な周波数、例
えば区間Aの4×10-4〜10-3Hzの周波数を選ぶことによ
り近似的に定数と見做せる。 r1<<<r2(r2は空気等の通気孔35の流量抵抗)で熱
時定数τとr2{V1+Min V2}/RT2が同程度の値なら次の
ような角周波数ωが存在する。 この条件において補正係数k2(s,r1,r2,V1,V2)は次
のように近似される。 |K2(iω,r1,r2,V1,V2)|≒1 ∠K2(iω,r1,r2,V1,V2)=0 (2d) ゆえに、式(2c)の条件を満たす場合、入力v(t)
から出力ΔP1(t)までの伝達関数は、γP0/(V1+V2
+ΔV)となる。 なお、体積変化手段33が角周波数ωで正弦波状に駆
動される場合、 が満たされていれば、パイプ32が閉塞されている状態と
同等と考えられる。 すなわち、 であればよい。 (2) 次に第5図のような単一タンクシステムを考え
る。これは第3図のタンク31,30の間を結合するパイプ3
2の断面積を非常に大きくしたもので、これによりパイ
プ32の流量抵抗r1の値が非常に小さくなる場合に相当す
る。これより、第5図のv(t)からΔP2(t)までの
伝達関数は式(2a)においてr1→0,T2=T1,ΔP2=ΔP1,
V2′=V1+V2としたものであり、次のようになる。 ここで、 となる。次のような角周波数ωを考える。 例えば第6図の周波数特性においてAに示す10-3Hz以
上の周波数である。このような周波数に設定することに
より補正係数K1(iω,r2,V3′)は次のように近似され
る。 |k1(iω,r2,V3′)|≒1, k1(iω,r2,V3′)≒0 (2i) このとき、伝達関数はγP0/(V2′+ΔV)となる。 次に上記原理を第7図に基づいて説明する。 第7図において33はスピーカ(体積変化手段)であ
り、この体積変化手段33によりメインタンク30と補正タ
ンク31が仕切られるように構成されている。さらにその
メインタンク30側にはコンデンサマイク(第1の圧力検
出手段)34bが設けられ、また補正タンク31側には上記
第1のマイクロホンより34b感度が低いダイナミックマ
イクロホン(第2の圧力検出手段)34aが設けられてい
る。また、前記スピーカ33は、所定の角周波数ωで駆
動する。この場合、第5図に示す単一タンクシステムの
理論が適用される。 ここで補正タンク31の容積をV1、メインタンク30内の
気体の体積をV2、メインタンク30内の液体の体積をVL
補正タンク31とメインタンク30の容積の和をVTとする。 補正タンク31の圧力変化ΔP1(t)は式(2h)を満た
す角周波数ωを用いると、 となる。また、メインタンク30を剛体、すなわちΔV=
0とし、角周波数ωが式(2h)を満たすとき、ΔP
1(t)は次のようになる。 またv(t)を角周波数ωで駆動したときの補正タ
ンク31側のΔP1′(t)の振幅を測定すると、 ここで体積変化手段33を、圧力変化が正弦波になるよ
うに駆動した場合の補正タンク31の圧力変化の振幅値を
A1、メインタンク30の圧力変化の振幅値をA2とすると、
これらの振幅値の比は式(3a),(3b)より となる。 これより となり、メインタンク30内に収納された液体の体積V
Lは、 となる。 なお、|k1/k2|,ΔVは予め実験により求めておくもの
とする。 次に上記原理に基づく具体的先行技術を説明する。 第7図は、補正タンク31とメインタンク30との間にス
ピーカ(体積変化手段)33を置いたもので、また、メイ
ンタンク30と補正タンク31の静圧を等しくするため、両
タンク30,31間がオリフィス47aを有する細いパイプ47で
接続され、メインタンク30と補正タンク31とが大気圧の
影響を受けないように、略完全に閉じられた系となって
いるため、スピーカ33の駆動角周波数ωがゆっくりし
た大気圧の変化に比べて非常に大きく、また通気孔35に
おける流体抵抗r2が非常に大きいことにより、スピーカ
33の駆動中は通気孔35が塞がれているように作用するの
で、タンク30,31内の気圧はタンク30,31内外の気圧差に
よる影響を全く受けない。ここでこのパイプ47の圧力伝
達の時定数は、スピーカ33による補正タンク31内の圧力
変化の時定数よりも十分大きく、またタンク30,31外の
大気圧、すなわち絶対圧力の圧力変化の時定数より十分
小さいものとする前提のもので、スピーカ33を駆動し、
v0 sin ω0tの関数で示されるような体積変化を補正タ
ンク31及びメインタンク30に与えると、すなわちメイン
タンク30および補正タンク31の双方に単一タンクシステ
ムの原理が適応される。そこで、 V(t)=v0 sin ω0t (4a) とすると、補正タンク31の圧力変化ΔP1(t)は、 となり、上式において角周波数ωが(2h)を満たす場
合、 となる。ΔP2(t)は次のようになる。 ここで上式において、メインタンク30が剛体で補正タ
ンク31と同様に角周波数ωが式(2h)を満たす場合、 となる。 すなわち、スピーカ33によって、メインタンク30と補
正タンク31のそれぞれの体積がV0 sinω0tだけ角周波数
ωで規則的に変動させられると、メインタンク30内と
補正タンク31内のそれぞれの圧力変動はそれぞれのタン
ク30,31に取付けられたダイナミックマイクロホン34a,
コンデンサマイクロホン34bによって検出され、メイン
タンク30内の圧力変動を検出したコンデンサマイクロホ
ン34bの出力は、ゲイン1、中心角周波数ωのバンド
パスフィルタ37bによって角周波数ωの信号成分が抽
出され、その後、第2の振幅検出器39bに供給され、γ
P0v0/V2が検出される出力される。また補正タンク31の
圧力変動を検出したダイナミックマイクロホン34aの出
力は、ゲインV1、中心角周波数ωのバンドパスフィル
タ37aによって角周波数ωの信号成分のみがV1倍され
て抽出され、その後第1の振幅検出器39aに供給され、
γP0v0が検出され出力される。その後、第1の振幅検出
器39aからの出力γP0v0は、第2の振幅検出器39bからの
出力γ P0v0/V2で割算器40によって除算され、メインタ
ンク30の空洞部分の体積V2が算出され、その演算結果
は、引算器41に供給され、設定されたメインタンク30の
全容積VTから引算され、その結果、メインタンク30内に
収納された液体の体積VLが算出される。
Conventional volume measurement methods and devices of this type include:
As a prior art, for example, there is one disclosed in Japanese Patent Application No. 1-28088, which will be described with reference to FIGS. First, referring to FIG. 3, the principle will be described. Reference numeral 30 denotes an irregularly shaped main tank for storing an object to be measured such as a liquid, powder, or an irregularly shaped object. A correction tank 31 having a smaller volume than the main tank 30 is connected via the main tank 30. Further, a small-diameter vent hole 35 is formed in the upper part of the irregularly shaped main tank 30 (the vent hole 35 may not be provided). A volume changing means (mechanism) 33 such as a piston and a speaker (hereinafter referred to as a speaker) is provided on the upper part of the correction tank 31, and the volume in the correction tank 31 is controlled by the operation of the volume changing means 33. Can be changed. The above is the configuration of the present embodiment. Next, the principle of measurement by the configuration will be described. Measurement principle (1) First, consider a correction tank 31 of small volume V 1, the connection tank system as shown in FIG. 3 that the main tank 30 of the large volume V 2 is connected. This is what the correct tank 31 main tank 30 is connected by a pipe 32 of the flow resistance r 1, the vent hole 35 of the tank 30 is a flow resistance r 2. The specific heat ratio of the gas in both tanks 30 and 31 is γ, the gas constant is R, and the thermal time constant of the tank 31 is τ. A volume changing means 33 using a piston, a diaphragm, a bellows, a speaker shape or the like is attached to the tank 31, and the volume change actually generated by the volume changing means 33 is defined as v (t). When the tanks 30 and 31 are rigid bodies, the tanks 30 and 31 are not distorted when the internal pressure of the tanks 30 and 31 is increased or decreased.
The volume change amount V 0 (t) due to the diaphragm, the bellows and the like is equal to the volume change amount v (t) actually generated. If the tank 30 is flexible, the tank 30 is distorted when the internal pressure of the tanks 30 and 31 is increased or decreased, so that the actual volume change amount v corresponds to the volume change amount due to expansion or contraction.
(T) changes. When v (t) = 0, the absolute pressure, temperature, and mole number of the gas in the tank 31 are P 0 , T 1 , n 1 , and the absolute pressure, temperature, and mole number of the gas in the tank 30 are, respectively. P 0 , T 2 , n 2 respectively
And If the measurement environment does not change significantly,
Gas flows into and out of tanks 30 and 31 via
The absolute pressure P 0 within 1,30 is equal to the outside pressure, its change is very slow and changes equal to the outside pressure. When v (t) ≠ 0, the pressure, temperature, and number of moles also change according to the state of the volume changing means 33. In the tank 31, the pressure is P 0 + ΔP 1 (t), and the temperature is T 1 + ΔT 1 (t The number of moles changes as n 1 −Δn 12 (t). In the tank 30, the pressure changes as P 0 + ΔP 2 (t), the temperature changes as T 2 + ΔT 2 (t), and the number of moles changes as n 2 + Δn 12 (t) −Δn 2 (t). Δn
12 (t) is the number of moles of air flowing from the tank 31 to the tank 30, and Δn 2 (t) is the number of moles of air leaking from the tank 30 to the outside through the vent hole 35. Here, the following assumptions are made for this system. 1) The gas in the tanks 30, 31 is an ideal gas. 2) v (t) ≪ {V 1 , V 2 } 3) The heat capacity of the tank 30 is large, and the temperature change in the tank due to the pressure change ΔP 2 (t) depends on the speed of change of the volume change amount v (t). Very slow to ignore. 4) The rate of change of the volume change amount v (t) is such that the pressure that changes accordingly is equal throughout the tanks 30 and 31. 5) Even if an object to be measured is placed in the tank 30, this object does not form two or more closed gas spaces, that is, hollow portions in the tank 30. The above assumption is not a major constraint. ΔP 1 (t), ΔP 2 (t), ΔT 1 with respect to the volume change amount v (t)
The changes in (t), ΔT 2 (t), Δn 12 (t), and Δn 2 (t) are originally represented by nonlinear equations, but from Assumption 2), the magnitudes are P 0 , T 1 , T 2 , n 1 , n 2 are very small and can therefore be approximated by linear equations. Tank 3 in static state
The relationship between the pressure, temperature, and number of moles of the gas in 0,31 is expressed by the following algebraic equation. P 0 V 1 = n 1 RT 1 , P 0 V 2 = n 2 RT 2 (1a) Assumptions 1), 3), 4) and 5) indicate that the gas in tanks 30 and 31 in the dynamic state The relationship between pressure, temperature and number of moles is
It is expressed by the following linear ordinary differential equation. The flow resistance r, in the above equation, r 1 and r 2 are determined from the length l and the diameter d of the pipe 32 as in the following experimental equation. In this formula, the length 1 is 50 to 650 [mm], and the diameter d is 2.0 to 9.0.
It was obtained experimentally using an aluminum pipe of [mm]. The volume change amount v (t) is expressed as follows. v (t) = v 0 (t) −Δv (t) (1h) ΔV is a constant peculiar to the tank determined by the material, shape, volume, etc. of the tank, and Δv (t) is a volume change due to expansion or contraction of the tank 30 due to a change in the volume change amount v 0 (t) of the volume changing means 33. Quantity. Expressions (1a) to (1i) are subjected to Laplace transform, and input v
The transfer function from (t) to the output ΔP 1 (t) is obtained as follows. here The coefficients r 2 V 2 / RT 2 , r 1 V 2 / RT 2 , r 1 V 1 / RT 1 and r 1 V 2 / RT 1 are time constants of the pressure change. For example, r 2 V 2 / RT 2 is a gas having an absolute temperature T 2 in a hollow portion of the tank 30 and a gas flow hole 35 having a flow resistance r 2.
Is the time constant of the pressure decay when flowing out of the tank 30 via. The correction coefficient k 2 (s, r 1 , r 2 , V 1 , V 2 ) is
Although it changes with the volume V 2 of 0, k 2 (s, r
1, the frequency characteristics of the r 2, V 1, V 2 ), the appropriate frequency, approximately constant and considered to by selecting the 4 × 10 -4 ~10 -3 Hz frequency, for example sections A. If the thermal time constant τ and r 2流量 V 1 + Min V 2 } / RT 2 in r 1 << r 2 (r 2 is the flow resistance of the air hole 35 such as air) are approximately the same, the following angle There is a frequency ω. Under this condition, the correction coefficient k 2 (s, r 1 , r 2 , V 1 , V 2 ) is approximated as follows. | K 2 (iω, r 1 , r 2 , V 1 , V 2 ) | ≒ 1∠K 2 (iω, r 1 , r 2 , V 1 , V 2 ) = 0 (2d) Therefore, equation (2c) If the condition is satisfied, the input v (t)
From the output ΔP 1 (t) is γP 0 / (V 1 + V 2
+ ΔV). Incidentally, if the volume change unit 33 is driven at an angular frequency omega 0 sinusoidally, Is satisfied, it is considered equivalent to a state where the pipe 32 is closed. That is, Should be fine. (2) Next, consider a single tank system as shown in FIG. This is the pipe 3 connecting the tanks 31 and 30 in FIG.
2 of cross-sectional area which was very large, thereby corresponds to the case where the value of the flow rate resistance r 1 of the pipe 32 becomes very small. Thus, the transfer function from v (t) to ΔP 2 (t) in FIG. 5 is r 1 → 0, T 2 = T 1 , ΔP 2 = ΔP 1 , in equation (2a).
V 2 ′ = V 1 + V 2 and is as follows. here, Becomes Consider the following angular frequency ω. For example, the frequency is 10 -3 Hz or more indicated by A in the frequency characteristic of FIG. By setting such a frequency, the correction coefficient K 1 (iω, r 2 , V 3 ′) is approximated as follows. | k 1 (iω, r 2 , V 3 ′) | ≒ 1, k 1 (iω, r 2 , V 3 ′) ≒ 0 (2i) At this time, the transfer function is γP 0 / (V 2 ′ + ΔV). Become. Next, the above principle will be described with reference to FIG. In FIG. 7, reference numeral 33 denotes a speaker (volume changing means), which is configured so that the main tank 30 and the correction tank 31 are separated by the volume changing means 33. Further, a condenser microphone (first pressure detecting means) 34b is provided on the main tank 30 side, and a dynamic microphone (second pressure detecting means) having lower sensitivity 34b than the first microphone is provided on the correction tank 31 side. 34a is provided. Further, the loudspeaker 33 is driven at a predetermined angular frequency omega 0. In this case, the theory of the single tank system shown in FIG. 5 applies. Here, the volume of the correction tank 31 is V 1 , the volume of the gas in the main tank 30 is V 2 , the volume of the liquid in the main tank 30 is V L ,
The sum of the volume of the compensation tank 31 and the main tank 30 and V T. Pressure change [Delta] P 1 of the correction tank 31 (t) is the use of the angular frequency omega H satisfying the equation (2h), Becomes Further, the main tank 30 is made rigid, that is, ΔV =
0 and the angular frequency ω H satisfies the equation (2h), ΔP
1 (t) is as follows. Also when measuring the amplitude of [Delta] P 1 of the correction tank 31 side '(t) when driven v a (t) at an angular frequency omega H, Here, the amplitude value of the pressure change in the correction tank 31 when the volume changing means 33 is driven so that the pressure change becomes a sine wave is
A 1 , if the amplitude value of the pressure change of the main tank 30 is A 2 ,
The ratio of these amplitude values is obtained from equations (3a) and (3b). Becomes Than this And the volume V of the liquid stored in the main tank 30
L is Becomes Note that | k 1 / k 2 |, ΔV is determined in advance by an experiment. Next, a specific prior art based on the above principle will be described. FIG. 7 shows a speaker (volume changing means) 33 placed between the correction tank 31 and the main tank 30. In order to equalize the static pressures of the main tank 30 and the correction tank 31, both tanks 30, 31 is connected by a thin pipe 47 having an orifice 47a, and the main tank 30 and the correction tank 31 are almost completely closed so as not to be affected by the atmospheric pressure. Since the angular frequency ω 0 is very large compared to the slowly changing atmospheric pressure, and the fluid resistance r 2 in the vent hole 35 is very large, the speaker
During the operation of 33, the air holes 35 act as if they are closed, so that the air pressure in the tanks 30, 31 is not affected at all by the pressure difference between the inside and outside of the tanks 30, 31. Here, the time constant of the pressure transmission of the pipe 47 is sufficiently larger than the time constant of the pressure change in the correction tank 31 by the speaker 33 and the time constant of the atmospheric pressure outside the tanks 30 and 31, ie, the pressure change of the absolute pressure. Assuming that it is much smaller, drive speaker 33,
When a volume change as shown by a function of v 0 sin ω 0 t is given to the correction tank 31 and the main tank 30, that is, the principle of the single tank system is applied to both the main tank 30 and the correction tank 31. Then, assuming that V (t) = v 0 sin ω 0 t (4a), the pressure change ΔP 1 (t) of the correction tank 31 becomes And when the angular frequency ω 0 satisfies (2h) in the above equation, Becomes ΔP 2 (t) is as follows. Here, in the above equation, when the main tank 30 is a rigid body and the angular frequency ω 0 satisfies the equation (2h) as in the case of the correction tank 31, Becomes That is, the speaker 33, the respective volumes of the main tank 30 correction tank 31 is regularly varied in angular frequency omega 0 only V 0 sinω 0 t, respectively pressure within the correction tank 31 main tank 30 Fluctuations are caused by the dynamic microphones 34a,
Is detected by the capacitor microphone 34b, the output of the condenser microphone 34b that detects the pressure fluctuations in the main tank 30, a gain 1, the signal component of the center angular frequency omega angular frequency omega 0 by a band-pass filter 37b of 0 is extracted, then , Supplied to the second amplitude detector 39b,
P 0 v 0 / V 2 is output is detected. The output of the dynamic microphone 34a detects a pressure variation of the correction tank 31, only the signal component of the gain V 1, the angular frequency omega 0 by a band-pass filter 37a of the center angular frequency omega 0 is extracted is 1 × V, then Supplied to the first amplitude detector 39a,
γP 0 v 0 is detected and output. After that, the output γP 0 v 0 from the first amplitude detector 39a is divided by the output γ P 0 v 0 / V 2 from the second amplitude detector 39b by the divider 40, and the cavity moiety is the volume V 2 of the calculation, the calculation result is supplied to the subtractor 41 is subtracted from the total volume V T of the main tank 30 that is set, as a result, the liquid stored in the main tank 30 volume V L is calculated for.

【発明が解決しようとする課題】[Problems to be solved by the invention]

しかしながら、このような従来の体積測定方法及びそ
の装置にあっては、スピーカ33よりの電気信号のうちか
らノイズ分を除去するバンドパスフィルタ37a,37bの出
力振幅値を検出する振幅検出器39a,39bなどによりマイ
クロホン34a,34bが検出した信号をアナログ処理する構
成であり、また信号に低周波数範囲で大きなノイズが重
畳されるので、計測誤差が発生する確率が大きいという
問題点があった。 この発明は上記のような問題点を解消するためになさ
れたもので、重畳ノイズを大幅に軽減し測定精度を向上
させた体積測定方法及びその装置を提供することを目的
とする。
However, in such a conventional volume measurement method and its device, the amplitude detector 39a, which detects the output amplitude value of the band-pass filters 37a, 37b for removing noise from the electric signal from the speaker 33, Since the signal detected by the microphones 34a and 34b by 39b or the like is processed in an analog manner, and a large noise is superimposed on the signal in a low frequency range, there is a problem that a measurement error is likely to occur. SUMMARY OF THE INVENTION The present invention has been made to solve the above problems, and an object of the present invention is to provide a volume measuring method and an apparatus therefor which greatly reduce superimposed noise and improve measurement accuracy.

【課題を解決するための手段】[Means for Solving the Problems]

第1請求項に係る体積測定方法は連通して設けられた
メインタンクと補正タンクとのそれぞれの内圧を位相が
一定区間毎に反転する所定周波数の信号により変化させ
て得られる検出出力と体積変化手段を駆動する信号とで
自己相関をとり、その相関値に基づいて上記メインタン
ク内に収納されている被測定物の体積を測定するように
したものである。 第2請求項に係る体積測定装置は所定周波数で発振す
る第1の発振器と、この第1の発振器の発振周波数の1/
nに同期して発振する第2の発振器と、これら第1およ
び第2の発振器の発振出力を乗算する乗算器と、被測定
物を収納するメインタンクと、このメインタンクに連通
して設けられた補正タンクと、これら両タンクのそれぞ
れの内圧を上記乗算器の乗算信号で変化させる体積変化
手段と、この体積変化手段による圧力変化を検出した検
出出力と上記体積変化手段を駆動する乗算信号との間で
相関を演算する自己相関手段と、この自己相関手段で得
られた相関値に基づいて上記メインタンク内に収納され
た被測定物の体積を測定する被測定物算出手段とで構成
されたものである。
The volume measurement method according to the first aspect is a detection output and a volume change obtained by changing respective internal pressures of a main tank and a correction tank provided in communication by a signal of a predetermined frequency whose phase is inverted every fixed section. An autocorrelation is obtained with a signal for driving the means, and the volume of the object to be measured stored in the main tank is measured based on the correlation value. A volume measuring device according to a second aspect includes a first oscillator that oscillates at a predetermined frequency, and 1/1 of an oscillation frequency of the first oscillator.
a second oscillator that oscillates in synchronization with n, a multiplier that multiplies the oscillation outputs of the first and second oscillators, a main tank that stores the device under test, and a communication unit that is provided in communication with the main tank. Correction tank, volume change means for changing the internal pressure of each of these tanks by a multiplication signal of the multiplier, a detection output of detecting a pressure change by the volume change means, and a multiplication signal for driving the volume change means. And an object-to-be-measured calculating means for measuring the volume of the object to be measured housed in the main tank based on the correlation value obtained by the auto-correlation means. It is a thing.

【作用】[Action]

第1請求項における体積測定方法は連通して設けられ
たメインタンクと補正タンクとのそれぞれの内圧を体積
変化手段で変化させ、それに伴う圧力変化を検出した検
出出力と上記体積変化手段を駆動する位相が一定区間毎
に反転する所定周波数の駆動信号とで相関を計算するよ
うにしたものである。 第2請求項における体積測定装置は連通して設けられ
たメインタンクと補正タンクとのそれぞれの内圧を変化
させて両タンクの圧力変化を検出しこの検出出力と体積
変化手段を駆動する駆動信号との相関を自己相関手段で
求めるようにしたものである。
In the volume measuring method according to the first aspect, the internal pressures of the main tank and the correction tank provided in communication with each other are changed by the volume changing means, and a detection output detecting the pressure change accompanying the driving and the volume changing means are driven. The correlation is calculated with a drive signal of a predetermined frequency whose phase is inverted every fixed section. The volume measuring device according to the second aspect changes the internal pressure of each of the main tank and the correction tank provided in communication with each other to detect a change in the pressure of both tanks, and outputs the detected output and a drive signal for driving the volume changing means. Are obtained by autocorrelation means.

【実施例】【Example】

以下、この発明を図面に基づいて詳細に説明する。第
1図はこの発明の第1実施例を示す回路図で、第1図に
おいて第7図と同一または均等な構成部分には同一符号
を付して重複説明を省略する。図において、50は、信号
V0sinωtを出力する第1の発振器、51は信号sin n ω
tを出力する第2の発振器で、この第2の発振器51の発
振周波数は第1の発振器50の発振周波数の1/nで同期し
て発振する(ただしnは偶数)。52は第1および第2の
発振器50,51の発振信号を掛け合せる第1の乗算器で、
この乗算信号により体積変化手段であるスピーカ33を駆
動する。53は位相差検出回路で、第1の乗算器52からの
出力信号V0sinωt・sin n ωtと第1のマイクロホン3
4aから出力される検出信号A1sin(ωt+φ)・sin(n
ωt+φ)との間の位相差φを検出する。54は第1の位
相合せ回路で、第1の乗算器52の出力信号V0sinωt+s
in n ωtを位相差検出回路53よりの検出出力の位相差
φにより第1のマイクロホン34aの出力の位相に合せ
る。55は第2の乗算器で、第1のマイクロホン34aの検
出出力A1sin(ωt+φ)・sin n ωtと第1の位相合
せ回路54の出力V0sin(ωt+φ)・sin(nωt+φ)
とを掛け合せる。56は第1の位相合せ回路54と同じ機能
を有する第2の位相合せ回路で、第2の発振器51の出力
信号sin n ωtを信号sin(nωtφ)として出力す
る。57はゼロクロス検出回路で、第2の位相合せ回路56
の出力sin(nωtφ)のある時点でのゼロクロスで積
分開始信号を出力し、それからm回目の所定時間(t)
後のゼロクロスで積分終了信号を出力する。58は第2の
乗算器55からの出力をゼロクロス検出回路57よりの積分
開始信号と積分終了信号とに基づいて積分してγP0v0
出力する第1の積分器で、この第1の積分器58は積分終
了信号が来るまでその積分出力はでない。59は第1の位
相合せ回路54よりの出力を反転させる反転回路、60は第
3の乗算器で、第2の乗算器55と同一機能を有する。61
は第2の積分器で、第1の積分器58と同一機能を有し、
ゼロクロス検出回路57から積分終了信号が供給されると
積分を停止してγ P0v0/V2を出力する。 次に動作について説明する。 第1の発振器50は第2図(a)に示されるV0sinωt
の出力信号を、また第2の発振器51は第2図(b)に示
されるsinωtの出力信号をそれぞれ第1の乗算器52に
出力し、第2図(c)に示される乗算信号V0sinωt・s
in nωtをスピーカ33に供給してスピーカ33を駆動す
る。第1の乗算器52の出力V0sinωt・sin nωtを基準
として第1のマイクロホン34aよりの検出出力A1sin(ω
t+φ)(例えば第2図(d))との間の位相差φを位
相差検出回路53により検出し、その検出された位相差φ
に基づいて第1の乗算器51からの出力V0sinωt・sin n
ωtの位相を修正してV0sin(ωt+φ)・sin(nωt
+φ)を出力する(第2図(e))。この修正された信
号V0sin(ωt+φ)・sin(nωt+φ)は第2の乗算
器55で第1のマイクロホン34aからの検出信号A1sin(ω
t+φ)と掛け合わされる。一方、第2の発振器51の出
力信号sin nωtは第2の位相合せ回路56により第2図
(f)のようにsin(nωt+φ)に変換され、ゼロク
ロス検出回路57に出力する。そして、第2の乗算器55よ
り第1の積分器58に供給される乗算信号をゼロクロス検
出回路57よりの積分開始信号に基づいて積分を開始し、
それから所定時間(t)後のゼロクロス検出回路57より
の積分終了信号に基づいて積分を終了する。この積分開
始から積分終了までの区間で得られるはずの信号成分を
Sとし、これに定常的なノイズNが重畳されている場
合、第2図(f)に示す正区間では が算出され、トータルの積分結果は によりノイズNがカットされ、信号成分Sだけが算出さ
れる。すなわちγP0v0を得る。 一方、第1の位相合せ回路54の出力V0sin(ωt+
φ)・sin(nωt+φ)は反転回路59で反転されてそ
の出力−V0sin・(ωt+φ)・sin(nωt+φ)が第
2のマイクロホン34bよりの検出信号−A2sin(ωt+
θ)・sin(nωt+φ)と第3の乗算器60で掛け合わ
され、ゼロクロス検出回路57よりの積分信号に基づいて
第2の積分器61において第1の積分器58と同じ働きのγ
P0v0/V2を得る。そののち、割算器40により第1の積分
器58の積分値γP0v0を第2の積分器61の積分値γ P0v0/
V2により割算し、メインタンク30の気体の成績V2を求
め、次いで減算器41によりメインタンク30の全体積VT
ら気体の体積V2を減算することによりメインタンク30内
の被計測物体の体積VLを求めることができる。
Hereinafter, the present invention will be described in detail with reference to the drawings. FIG. 1 is a circuit diagram showing a first embodiment of the present invention. In FIG. 1, the same or equivalent components as those in FIG. In the figure, 50 is a signal
A first oscillator that outputs V 0 sinωt, 51 is a signal sin n ω
The second oscillator 51 outputs t and oscillates in synchronization with the oscillation frequency of the second oscillator 51 at 1 / n of the oscillation frequency of the first oscillator 50 (where n is an even number). 52 is a first multiplier for multiplying the oscillation signals of the first and second oscillators 50 and 51,
The speaker 33 which is a volume changing means is driven by the multiplied signal. Reference numeral 53 denotes a phase difference detection circuit, which outputs an output signal V 0 sinωt · sin n ωt from the first multiplier 52 and the first microphone 3
4a, the detection signal A 1 sin (ωt + φ) · sin (n
ωt + φ) is detected. Reference numeral 54 denotes a first phase matching circuit, which is an output signal V 0 sinωt + s of the first multiplier 52.
in n ωt is matched with the phase of the output of the first microphone 34a by the phase difference φ of the detection output from the phase difference detection circuit 53. 55 is a second multiplier, which is the detection output A 1 sin (ωt + φ) · sin n ωt of the first microphone 34a and the output V 0 sin (ωt + φ) · sin (nωt + φ) of the first phase matching circuit 54.
Multiply with Reference numeral 56 denotes a second phase matching circuit having the same function as the first phase matching circuit 54, and outputs the output signal sin n ωt of the second oscillator 51 as a signal sin (nωtφ). 57 is a zero cross detection circuit, and a second phase matching circuit 56
Outputs an integration start signal at a zero crossing at a certain point in time of the output sin (nωtφ), and then performs an m-th predetermined time (t)
An integration end signal is output at the later zero cross. 58 is a first integrator that integrates the output from the second multiplier 55 based on the integration start signal and the integration end signal from the zero cross detection circuit 57 and outputs γP 0 v 0 . The integrator 58 does not output its integration until the integration end signal comes. 59 is an inverting circuit for inverting the output from the first phase matching circuit 54, and 60 is a third multiplier having the same function as the second multiplier 55. 61
Is a second integrator having the same function as the first integrator 58,
When the integration end signal is supplied from the zero cross detection circuit 57, the integration is stopped and γ P 0 v 0 / V 2 is output. Next, the operation will be described. The first oscillator 50 has V 0 sinωt shown in FIG.
And the second oscillator 51 outputs the output signal of sinωt shown in FIG. 2 (b) to the first multiplier 52, and the multiplied signal V 0 shown in FIG. 2 (c). sinωt ・ s
in nωt is supplied to the speaker 33 to drive the speaker 33. Based on the output V 0 sinωt · sin nωt of the first multiplier 52, the detection output A 1 sin (ω
t + φ) (for example, FIG. 2 (d)) is detected by the phase difference detection circuit 53, and the detected phase difference φ
Output V 0 sinωt · sin n from the first multiplier 51 based on
By correcting the phase of ωt, V 0 sin (ωt + φ) · sin (nωt
+ Φ) (FIG. 2 (e)). The corrected signal V 0 sin (ωt + φ) · sin (nωt + φ) is converted by a second multiplier 55 into a detection signal A 1 sin (ω) from the first microphone 34a.
t + φ). On the other hand, the output signal sin nωt of the second oscillator 51 is converted by the second phase matching circuit 56 into sin (nωt + φ) as shown in FIG. Then, the multiplication signal supplied from the second multiplier 55 to the first integrator 58 is integrated based on the integration start signal from the zero cross detection circuit 57, and
Then, the integration is terminated based on the integration termination signal from the zero-cross detection circuit 57 after a predetermined time (t). Let S be a signal component that should be obtained in the section from the start of integration to the end of integration, and when a steady noise N is superimposed on this, in the positive section shown in FIG. Is calculated, and the total integration result is , The noise N is cut off, and only the signal component S is calculated. That is, γP 0 v 0 is obtained. On the other hand, the output V 0 sin of the first phase matching circuit 54 (ωt +
φ) · sin (nωt + φ) is inverted by the inverting circuit 59, and the output −V 0 sin · (ωt + φ) · sin (nωt + φ) is the detection signal −A 2 sin (ωt +) from the second microphone 34b.
θ) · sin (nωt + φ) is multiplied by the third multiplier 60, and based on the integration signal from the zero-crossing detection circuit 57, the second integrator 61 performs the same operation as the first integrator 58 γ
Obtain P 0 v 0 / V 2. After that, the divider 40 integrates the integrated value γP 0 v 0 of the first integrator 58 with the integrated value γ P 0 v 0 /
Divided by V 2, obtains the results V 2 of the gas in the main tank 30, and then the measurement of the main tank 30 by subtracting the volume V 2 of gas from the total volume V T of the main tank 30 by the subtractor 41 The volume VL of the object can be determined.

【発明の効果】 以上説明したように第1請求項によればその構成を連
通して設けられたメインタンクと補正タンクとのそれぞ
れの内圧を位相が一定区間毎に反転する所定周波数の信
号で駆動される体積変化手段を用いて変化させ、この圧
力変化を圧力センサで検出し、この圧力センサの検出出
力と上記体積変化手段を駆動する上記信号とで相関をと
り、次いでその相関値に基づいて上記メインタンク内に
収納された被測定物の体積を測定する体積測定方法とし
たため、重畳ノイズを大幅に軽減し、測定精度を上げる
ことができるなどの効果が得られる。 また第2請求項によればその構成を所定周波数で発振
する第1の発振器と、この第1の発振器の発振周波数の
1/nに同期して発振する第2の発振器と、これら第1お
よび第2の発振器の発振出力を乗算する乗算器と、被測
定物を収納するメインタンクと、このメインタンクに連
通して設けられた補正タンクと、これら両タンクのそれ
ぞれの内圧を上記乗算器の乗算信号で変化させる体積変
化手段と、この体積変化手段による圧力変化を検出して
得た検出出力と上記体積変化手段を駆動する乗算信号と
の間で相関を演算する自己相関手段と、この自己相関手
段で得られた相関値に基づいて上記メインタンク内に収
納されている上記被測定物の体積を測定する被測定物体
積算出手段とを備えた体積測定装置としたためランダム
ノイズに対するS/Nを落すことなく定常ノイズに対するS
/Nを向上できるので、測定精度を上げることができるな
どの効果が得られる。
As described above, according to the first aspect, the internal pressure of each of the main tank and the correction tank provided in communication with the configuration is a signal of a predetermined frequency in which the phase is inverted every fixed section. The pressure change is detected by a pressure sensor, and the detected output of the pressure sensor is correlated with the signal for driving the volume changing means, and then based on the correlation value. Since the volume measurement method is used to measure the volume of the object stored in the main tank, effects such as greatly reducing superimposed noise and increasing measurement accuracy can be obtained. Further, according to the second aspect, a first oscillator oscillating the configuration at a predetermined frequency, and an oscillation frequency of the first oscillator.
A second oscillator that oscillates in synchronism with 1 / n, a multiplier that multiplies the oscillation outputs of the first and second oscillators, a main tank that stores the device under test, and communicates with the main tank. The correction tank provided, volume changing means for changing the internal pressure of each of these tanks with the multiplication signal of the multiplier, a detection output obtained by detecting a pressure change by the volume changing means, and the volume changing means. Autocorrelation means for calculating a correlation between the multiplication signal to be driven, and a measured object for measuring a volume of the object to be measured housed in the main tank based on the correlation value obtained by the autocorrelation means. The volume measurement device was equipped with an object volume calculation means, so the S / N for stationary noise was maintained without reducing the S / N for random noise.
Since / N can be improved, effects such as improvement in measurement accuracy can be obtained.

【図面の簡単な説明】[Brief description of the drawings]

第1図はこの発明に係る体積測定方法及びその装置の第
1実施例を示す回路ブロック説明図、第2図は第1図を
説明するためのタイムチャート図、第3図は先行技術の
原理説明図、第4図は第3図において体積変化手段の駆
動周波数とタンク内の被測定物体積とを変化させたとき
の伝達関数の係数K2の変化状態を示す特性図、第5図は
先行技術を説明するための原理説明図、第6図は第5図
において体積変化手段の駆動周波数とメインタンク内の
被測定物体積とを変化させたときの伝達関数の係数K1
変化状態を示す特性図、第7図は先行技術の具体的シス
テム説明図である。 30:メインタンク、31:補正タンク 33:スピーカ(体積変化手段) 50,51:発振器、52:第1の乗算器 53:位相差検出回路 54,56:位相合せ回路、55:第2の乗算器 57:ゼロクロス検出回路、58:第1の積分器 60:第3の乗算器、61:第2の積分器
FIG. 1 is a circuit block diagram showing a first embodiment of a volume measuring method and apparatus according to the present invention, FIG. 2 is a time chart for explaining FIG. 1, and FIG. 3 is a principle of the prior art. illustration, Figure 4 is a characteristic diagram showing a change in state of the coefficient K 2 of the transfer function when changing the object to be measured volume of the driving frequency and the tank of volume change means in Figure 3, Figure 5 is principle explanatory diagram for explaining a prior art, FIG. 6 is changing state of the coefficients K 1 of the transfer function when changing the object to be measured volume of the driving frequency and the main tank volume change means in Figure 5 FIG. 7 is an explanatory diagram of a specific system according to the prior art. 30: Main tank, 31: Correction tank 33: Speaker (volume changing means) 50, 51: Oscillator, 52: First multiplier 53: Phase difference detection circuit 54, 56: Phase matching circuit, 55: Second multiplication 57: Zero cross detection circuit, 58: First integrator 60: Third multiplier, 61: Second integrator

Claims (2)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】連通して設けられたメインタンクと補正タ
ンクとのそれぞれの内圧を位相が一定区間毎に反転する
所定周波数の信号で駆動される体積変化手段を用いて変
化させ、この圧力変化を圧力センサで検出し、この圧力
センサの検出出力と上記体積変化手段を駆動する信号と
で相関をとり、次いでその相関値に基づいて上記メイン
タンク内に収納された被測定物の体積を測定する体積測
定方法。
An internal pressure of each of a main tank and a correction tank provided in communication with each other is changed using volume changing means driven by a signal of a predetermined frequency whose phase is inverted every predetermined section, and the pressure change is performed. Is detected by a pressure sensor, and the detected output of the pressure sensor is correlated with the signal for driving the volume changing means, and then the volume of the measured object stored in the main tank is measured based on the correlation value. Volume measurement method.
【請求項2】所定周波数で発振する第1の発振器と、こ
の第1の発振器の発振周波数の1/nに同期して発振する
第2の発振器と、これら第1および第2の発振器の発振
出力を乗算する乗算器と、被測定物を収納するメインタ
ンクと、このメインタンクに連通して設けられた補正タ
ンクと、これら両タンクのそれぞれの内圧を上記乗算器
の乗算信号で変化させる体積変化手段と、この体積変化
手段による圧力変化を検出して得た検出出力と上記体積
変化手段を駆動する乗算信号との間で相関を演算する自
己相関手段と、この自己相関手段で得られた相関値に基
づいて上記メインタンク内に収納されている上記被測定
物の体積を測定する被測定物体積算出手段とを備えた体
積測定装置。
2. A first oscillator that oscillates at a predetermined frequency, a second oscillator that oscillates in synchronization with 1 / n of an oscillation frequency of the first oscillator, and oscillations of the first and second oscillators. A multiplier for multiplying the output, a main tank for storing the device under test, a correction tank provided in communication with the main tank, and a volume for changing the internal pressure of each of these tanks by a multiplication signal of the multiplier Changing means, autocorrelation means for calculating a correlation between a detection output obtained by detecting a pressure change by the volume changing means and a multiplication signal for driving the volume changing means, and autocorrelation means. A volume measuring device comprising: a measuring object volume calculating means for measuring a volume of the measuring object stored in the main tank based on a correlation value.
JP22900689A 1989-09-04 1989-09-04 Volume measurement method and device Expired - Lifetime JP2820448B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP22900689A JP2820448B2 (en) 1989-09-04 1989-09-04 Volume measurement method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP22900689A JP2820448B2 (en) 1989-09-04 1989-09-04 Volume measurement method and device

Publications (2)

Publication Number Publication Date
JPH0392725A JPH0392725A (en) 1991-04-17
JP2820448B2 true JP2820448B2 (en) 1998-11-05

Family

ID=16885283

Family Applications (1)

Application Number Title Priority Date Filing Date
JP22900689A Expired - Lifetime JP2820448B2 (en) 1989-09-04 1989-09-04 Volume measurement method and device

Country Status (1)

Country Link
JP (1) JP2820448B2 (en)

Also Published As

Publication number Publication date
JPH0392725A (en) 1991-04-17

Similar Documents

Publication Publication Date Title
JPS5838738B2 (en) pressure gauge
JPS6239935B2 (en)
US5824892A (en) Acoustic device for measuring volume difference
JP2717011B2 (en) Circuit device for temperature compensation of capacitive and differential pressure sensors
Sell et al. Real-time monitoring of a high pressure reactor using a gas density sensor
JP2820448B2 (en) Volume measurement method and device
JP2856772B2 (en) Volume measurement method and device
JP3329480B2 (en) Differential pressure measuring device
CN114008414A (en) Method for determining a physical parameter of a charged liquid
JP2778971B2 (en) Volume measurement method and device
JPH09292290A (en) Sensor circuit
RU2627544C2 (en) Piezo resonant-viscous vacuum gauge
JP3790873B2 (en) Acoustic surface area meter
JP2907965B2 (en) Volume measurement device
JPH09159459A (en) Vibrating gyro
JP4379759B2 (en) Acoustic volume meter
JP3281102B2 (en) Volume measuring method and volume measuring device
JPH0635140Y2 (en) Volume measuring device
JPS63165719A (en) Volumeter
JP3455283B2 (en) Fluid volume measuring method and fluid volume measuring device
JP2000171282A (en) Acoustic volume meter for measuring volume difference
JP2897630B2 (en) Volumetric meter
JPH0447231A (en) Volume measuring apparatus
JPH07311071A (en) Gas type volumeter
JPH04296635A (en) Temperature correcting method of density sensor in vibratory gas density meter