(57)【特許請求の範囲】
【請求項１】
離散型に測定される瞬時流量を測定する瞬時流量測定手段を備えると共に、所定の微小流量域における流量出力を零流量出力とする流量の不感帯を設けたガスメーターであって、
前記瞬時流量のバラツキを表す統計関数、連続する前記瞬時流量の差分からなる数列の統計関数、もしくは単位時間あたりの瞬時流量の変化量からなる数列の統計関数のいずれか１種以上に基づいて、前記不感帯の上限値を設定する不感帯上限設定手段を備えたガスメーター。
【請求項２】
使用期間内においてメーター下流側に供給された流量である使用流量を、前記使用期間中に測定される前記瞬時流量から求めるに、
前記使用期間内に測定される前記瞬時流量に対して前記不感帯の上限値を適応して、前記使用流量を求める請求項１記載のガスメーター。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、離散型に測定される瞬時流量を測定する瞬時流量測定手段を備えたガスメーターであって、測定される瞬時流量から利用流量等を導出するガスメーターに関する。
【０００２】
【従来の技術】
従来、ガスメーターに使用される流量計としては、膜式のものが主流であるが、その利便性等との関係から、現今、超音波を利用したメーター、熱伝導式流量センサーを利用したメーターが提案されている。
このようなメーターにあっては、その測定原理から、取り込まれる流量は、時間的に離散して得られる瞬時流量であり、そのデーター列は離散型である。
【０００３】
一方、メーター下流側に供給されるガス流量を検出するためのガスメーターにあっては、流量の不感帯が設けられる。
この不感帯は、原則的には、測定原理等との関係から、測定の信頼性が確保できない場合に、この不安定な領域を零出力とするものである。
通常、この不感帯の上限は、上記のようにして設定され、ガスメーターの機種、号数（メーターにより測定される標準流量によって設定される）によって、予め設定される一定値とされる。
更に、従来、この不感帯は常時設定されており、上記の一定値より低い流量が検知された場合は、検知流量がないものとされる。
【０００４】
【発明が解決しようとする課題】
さて、この不感帯に関しては、以下のような事象が発生する。
（１）  不感帯の上限の設定値が低すぎる場合にあっては、メーターを取り付けた各客先にあって、本来の使用量より大きなガス流量を測定する、所謂、オーバーカウント状態を起こす。
（２）  不感帯の上限の設定値が高すぎる場合にあっては、ガス供給が行われているにもかかわらず、ガス供給が無いと測定する不測定ガス量が発生する。
（３）  更に、本来、メーターは、季節、経時変化等によって、その測定精度等が変化し、一般的な意味での誤差項を発生することとなるが、従来、不感帯は、一定に設定されていたために、このような変化に追随せず、不感帯が広い側に設定される。
（４）  流量の測定原理が、離散的に瞬時流量を求め、測定される多数の瞬時流量から利用流量等を求める構成のものにあっては、流れに脈動がある場合、この脈動と離散化された測定タイミングとの関係に起因する誤差が、流量の微小領域の流量測定結果に影響を与える場合があり、このような誤差の影響をも加味して、本来、不感帯は設定されるべきであるが、このような要因は、これまで、何ら顧みられていない。
【０００５】
本発明の目的は、離散型に測定される瞬時流量を測定する瞬時流量測定手段を備えたガスメーターにおいて、上記のような問題に対応できる合理的な不感帯を備えたガスメーターを得ることにある。
【０００６】
【課題を解決するための手段】
この目的を達成するための本発明による流量導出方法の特徴手段は、請求項１に記載されている様に、瞬時流量のバラツキを表す統計関数、連続する前記瞬時流量の差分からなる数列の統計関数、もしくは単位時間あたりの瞬時流量の変化量からなる数列の統計関数のいずれか１種以上に基づいて、前記不感帯の上限値を設定する不感帯上限設定手段を備えたことにある。
不感帯上限設定手段を備えたものにあっては、不感帯の上限が、前記統計関数のいずれかに基づいて設定される。このような設定にあたっては、不感帯の基底となる一定値の上に、統計関数分が乗るものとしても良いし、統計関数自体から不感帯の幅が決まるものとしてもよい。
さて、本願において使用する統計関数の性質であるが、基本的には、離散化されて測定される瞬時流量のバラツキ、ゆらぎ等を代表できる物理量であり、瞬時流量を離散的に測定する場合に発生しうる誤差を含んだものと見なすことができる。例えば、流量に脈動成分が含まれており、瞬時流量の測定タイミングを乱数化して分散設定する場合に、瞬時流量の測定所要時間は分散処理の対象とできないため、この分の誤差要因が残る。例えば、一定の積算時間Ｔ内で、ｎ回の測定を行う場合、一回の測定に要する測定所要時間をΔｔとすると、Ｔ−ｎ×Δｔの時間帯は分散可能で誤差要因は低いが、ｎ×Δｔの時間は分散の対象とできない。よって、この時間分の誤差がでる可能性が高い。
このような誤差要因は、本来メーターが測定すべき流量に影響すると共に、同時に、不感帯の信頼性にも影響を与える。従って、不感帯の上限（下限は当然０である）の設定にあたって、これを考慮するものとする。このようにすることで、例えば、測定の誤差要因を除去した状態での流量測定が可能となり、結果的に、不測定流量が過大でない状態でオーバーカウントのない流量測定を精度よく行うことが可能となる。ここで、瞬時流量のバラツキを表す統計関数、連続する前記瞬時流量の差分からなる数列の統計関数、もしくは単位時間あたりの瞬時流量の変化量からなる数列の統計関数のいずれかは、上記のようなゆらぎ等の影響を見るうえで好ましい物理量であるため、好適に採用できる。
【０００７】
上記の請求項１記載の導出方法において、請求項２に記載されているように、使用期間内においてメーター下流側に供給された流量である使用流量を、前記使用期間中に測定される前記瞬時流量から求めるに、
前記使用期間内に測定される前記瞬時流量に対して前記不感帯の上限値を適応して、前記使用流量を求めることが好ましい。
この構成を採用することで、利用流量として、過不足の無い流量を得ることができる。
【０００８】
さて、上記構成において、前記瞬時流量を測定するための平均測定タイミング間隔が予め設定されており、前記瞬時流量を測定した瞬時流量測定完了時点と後続する瞬時流量開始時点との間の時間である離散間隔に関して、多数の離散間隔の平均値と、前記瞬時流量の測定に要する測定所要時間との和が、前記平均測定タイミング間隔となるように、前記離散間隔をバラツキ処理する離散間隔分散手段を備え、前記統計関数から前記不感帯の上限値を設定するに、前記平均測定タイミング間隔及び前記測定所要時間に基づいて前記上限値を設定することが好ましい。
【０００９】
先にも説明したように、離散型の測定構造を有する場合、測定タイミングを乱数化処理等に基づいて分散させたとしても、乱数化できない測定に要する時間分（測定所要時間）に起因する誤差が発生する。そして、測定対象に脈動等があり、この脈動の周期が、平均測定タイミング間隔と一致していると、この脈動分が誤差要因となる前述の統計関数にかなり寄与することとなる。
従って、この寄与分を見越して、誤差としての性質を有する統計関数に基づいて、不感帯の上限を求める場合に、この影響を加味する意味から、平均測定タイミング間隔及び測定所要時間に基づいて上限値を設定すると、完全な分散処理が行えない場合の誤差要因を不感帯の上限値の設定で良好に反映させることができる。
【００１２】
【発明の実施の形態】
以下、本願の実施の形態を、図１に基づいて説明する。
ガスメーター１は、超音波を使用して所定の流れに関する情報を得るための測定管２と、この測定管２からの情報に基づいて、本願にいう瞬時流量を求める系及び得られた瞬時流量から本願にいう使用流量、安全機能用流量を求める系を備えて構成されている。
即ち、ガスメーター１は、前述の測定管２と、この測定管２の出力端２ａが接続されるメーター本体７とを備えており、コンピュータとして構成されるメーター本体７内には、前述の両者の処理系が、様々な手段として格納されている。
【００１３】
測定管２において、測定対象のガスｆは、導入部３から測定管２内に流入し、導出部４より排出される。同図において、測定管２内でのガスｆの流れ方向は、大矢印で示すように、図上、左から右である。左側が上流側、右側が下流側である。
測定管２には、この測定管２に対して所定の方向に位置づけて設けられる一対の送受波器５が備えられている。これら送受波器５の離間距離Ｌは一定とされる。
この送受波器５は、それぞれ、超音波の送波器５ａと受波器５ｂとを備えており、本体側からの指令に従って、一方側の送波器５ａでは、超音波を一方側から他方側へ送信可能に構成されていると共に、他方側の受波器５ｂでは、一方側から所定のタイミングで送信された超音波を受信することにより、その超音波の伝播時間（一方側の送波器５ａにおける送信タイミングと他方側の受波器５ｂにおける受信タイミングの時間差）を割り出すことが可能に構成されている。
この操作は、流れ方向に対して逆転した位置関係にある送波器５ａと受波器５ｂに関しても同様に実行可能とされ、流れ方向に沿った順方向（図上、左から右）と、逆方向（図上右から左）との両方向で、超音波が伝播され、各超音波の伝播時間を求めることができる構成が採用されている。
【００１４】
上記構造において、伝播時間の測定は、超音波が、音の伝播方向が流れに沿った方向である順方向に伝播されて、その順方向で超音波が送受波器間を伝播する順方向伝播時間ｔ１０を求める動作と、流れに逆らった方向である逆方向に伝播させ、その逆方向で超音波が送受波器間を伝播する逆方向伝播時間ｔ２０とを求める動作とを一対として、瞬時流量測定手段６によって実行される。
【００１５】
瞬時流量測定手段６は、出力されてくる順方向伝播時間ｔ１０と逆方向伝播時間ｔ２０とに基づいて、ガスの流速Ｖｘ（これは超音波の伝播方向に沿った流速成分）を求めると共に、測定部の断面積Ｓ１と、測定部の流量係数βを考慮して、単一回の流量を求める。
流速Ｖｘ、瞬時流量Ｑは、下記式で与えられる。
【００１６】
【数１】
Ｖｘ＝Ｌ×（１／ｔ１０−１／ｔ２０）／２
Ｑ＝Ｖｘ×Ｓ１×β×ｃｏｓ（θ）
【００１７】
この瞬時流量Ｑは、一回の測定により得られる値であり、離散的に測定される。瞬時流量の測定タイミングは、個々に、クロック８ａ及び分散手段８ｂを備えた測定時間制御手段８によって与えられる。
この測定時間制御手段８の働きを、図２、図３を参考にしながら説明する。図において、横軸上に離散点として個々の測定タイミングを示した。
個々の測定タイミングの指令を発するのが、この手段８の役割であるが、相互の測定タイミングの設定も、この手段８によって制御される。この制御は、分散手段８ｂにより以下のようなシーケンスを採用して実行される。
積算時間Ｔは、ガスメーターの用途等を考慮して、数十分〜数時間等、適切な一定値として予め設定される。また、この積算時間Ｔ内における測定タイミングの数（サンプリング数ｎ）も、使用目的等との関係から１００〜３００００に設定されている。従って、基準となるサンプリング間隔は、１秒程度である。
さて、各測定タイミングは、積算時間をＴ、超音波による順方向および逆方向を一対とした１回の測定に必要な時間（この測定時間は純粋に測定に必要な時間だけでなく、ＣＰＵ内の他の演算との兼ね合い等で発生する磁界測定に入れない禁止時間を含む）をΔｔ、測定タイミングの設定毎に毎回発生される一桁の乱数をＮとして以下の式に基づいて設定される。
この場合、ある一回の測定タイミングから次の測定タイミングまでの時間は次のようになる。
【００１８】
【数２】
Δｔ＋（Ｔ／ｎ−Δｔ）×Ｎ／｛（１＋Ｎｍａｘ）／２｝
ここで、乱数Ｎは１からＮｍａｘまでの整数の範囲で発生する乱数で、各数の発生する確率が同一のものとする。
このシーケンスを、順次、繰り返す。
【００１９】
即ち、分散手段８ｂにおける処理にあっては、瞬時流量を測定するための平均測定タイミング間隔Ｔ／ｎが予め設定されており、瞬時流量を測定した瞬時流量測定完了時点と後続する瞬時流量開始時点との間の時間である測定タイミング間隔（これを離散間隔と呼び、上記の場合はＴ／ｎ−Δｔである）に関して、多数の離散間隔の平均値と、瞬時流量の測定に要する測定所要時間Δｔとの和が、前記平均測定タイミング間隔Ｔ／ｎとなるように、離散間隔をバラツカせることとなる。
従って、この構成を採用することにより、積算時間Ｔ内でのサンプリング数ｎを統計的に確保できると共に、タイミングをバラツカせることで、論理的共鳴による誤差の発生確率を大幅に下げることができる。
図２、３にあっては、実測流量と記載したアナログ表示の流量線上に各点として瞬時流量の測定値を示している。
【００２０】
このようにして所定の積算時間Ｔ内における複数の瞬時流量Ｑが得られることとなる。
【００２１】
先ず、先に説明した瞬時流量Ｑから安全機能用流量Ｑｓａを求める系に関して説明する。この系で求まる誤差関数が、後に説明するように、本願の特徴である不感帯の設定に使用される。
【００２２】
安全機能用流量Ｑｓａを導出する部位には、積算時間Ｔ内に測定される複数の瞬時流量Ｑから、複数の瞬時流量Ｑのバラツキを表す第一統計関数を求める第一統計関数導出手段９と、積算時間Ｔ内に測定される複数の瞬時流量Ｑに関して、連続する瞬時流量の差分の統計関数である第二統計関数を求める第二統計関数導出手段１０とを備えている。
そして、第一統計関数導出手段９によって求められる第一統計関数と、第二統計関数導出手段１０によって求められる第二統計関数との、一方あるいは双方の関数である誤差関数を求める誤差関数導出手段１１が備えられている。
更に、誤差関数導出手段１１によって求められる誤差関数と、積算流量導出手段１２によって求められる積算流量とに基づいて、安全機能用の安全機能用流量Ｑｓａを求める安全機能用流量導出手段１３を備えている。
この安全機能用流量導出手段１３により得られた安全機能用流量Ｑｓａは、ガスメーター１が備えられている管路（図外）の遮断等の判断用に使用され、必要な場合は、管路の遮断が行われる。
【００２３】
以下、上記、図１に示す個々の手段に関して説明する。
イ  積算流量導出手段１２
瞬時流量測定手段６から、積算時間内に入力されてくる瞬時流量Ｑの平均と積算時間との積としての積算流量Ｑｓｕｍを導出する。この手段において、後述するように、本願の特徴構成である不感帯の処理が行われるが、この処理に関しては、再度、説明する。この積算流量Ｑｓｕｍは、安全機能用流量導出手段１３、不感帯処理、使用流量導出手段１６に送られる。
ロ  第一統計量導出手段９
瞬時流量測定手段６から、積算時間内に入力されてくる瞬時流量Ｑの標準偏差σを導出する。この標準偏差σは、誤差関数導出手段１１に送られる。
ハ  第二統計量導出手段１０
瞬時流量測定手段６から、積算時間内に入力されてくる連続する瞬時流量Ｑの差分の絶対値の平均ｄｍを導出する。この差分の平均ｄｍは、誤差関数導出手段１１に送られる。
ニ  安全機能用流量導出手段１３
第一統計関数としての標準偏差σ、第二統計関数としての差分の平均ｄｍを重み付けをして得られる誤差関数を、誤差関数導出手段１１から受け、入力されてくる積算流量Ｑｓｕｍに対して加減算を行って、安全機能用流量Ｑｓａを求める。
この安全機能用流量は、下記のように表現できる。
【００２４】
【数３】
Ｑｓａ＝Ｑｓｕｍ＋ａ１×σ＋ａ２×ｄｍ
ここで、ａ１、ａ２は、積算時間Ｔ等とバラツキ指標の種類に基づいて予め設定している定数である。
【００２５】
例えば、（ａ１、ａ２）として（−２、０）を採用する場合は、図２、図３における低位側の安全機能用流量（長破線で示す）を与え、（２、０）とする場合は、高位側の安全機能用流量（短破線で示す）を与える。
このような低位側の安全機能用流量を採用すると、脈動による誤作動の約９０％を防止することができる。
更に、（ａ１、ａ２）として、（０、−ｄｍ／ｓｑｒｔ（ｎ））としても良い。この場合も、ほぼ、上記と同様な効果を得ることができる。また、この手法にあっては、判断時間を短縮することができる。
【００２６】
（ａ１、ａ２）として、互いに有効な数字の組み合わせを採用することも可能であるし、例えば、積算時間内をさらに細かな時間に分割し、その分割された時間内の標準偏差と差分の平均との関係を判断して、脈動と機器のＯＮ／ＯＦＦによるものとを判別し、上記の組み合わせを変更するものとしても良い。
更に、脈動周期と平均測定タイミング間隔との関係において説明すると、脈動周期が平均測定タイミング間隔より小さい図２に示す状況では、瞬時流量の測定は、脈動の位相に対してほぼランダムに行われることになるから、差分の絶対値の平均は脈動振幅に比例するため、この平均値の倍数をとる方法は妥当性を持つ。
一方、脈動周期が平均測定タイミング間隔より大きい図３のような状況では、バラツキの程度を示す標準偏差と前記平均との両者間は有意差は存在しない。
従って、差分の絶対値の平均の倍数は、本願構造において容易且つ有用に選択・使用できる一つの関数である。
【００２７】
以上が、安全機能用流量を求めるための部分であるが、この部分とは並列に、所定の利用期間内で、メーターより下流側に供給されるガス流量を求めるための系が備えられている。
この系は、先の安全用の系と同様に、積算流量Ｑｓｕｍ、誤差関数を基本に、利用流量ΣＱｓｕｍを求める系である。
【００２８】
不感帯処理は、先の誤差関数を受けて、積算流量導出手段１２における積算流量の導出時に、不感帯処理をおこなうかどうかを設定するための情報を出力するようになっている。
即ち、図１に不感帯処理として示すように、前述の誤差関数を受けて不感帯の上限を設定する不感帯上限設定手段１４が設けられていると共に、この不感帯上限設定手段１４による設定を適応するかどうか（即ち、装置は、不感帯を設ける不感帯設定状態と、不感帯を設けない不感帯非設定状態とに切り換え自在に構成されている）の判断を行う不感帯設定切り換え手段１５を備えている。
このような不感帯に関する情報は、積算流量導出手段１２に送られ、この積算流量導出手段１２における、前述の積算流量Ｑｓｕｍの導出時に使用される。
以下、各手段を順に説明する。
【００２９】
ホ  不感帯上限設定手段１４
この手段は、連続する前記瞬時流量の差分からなる数列の統計関数（具体的には先の平均ｄｍ）に基づいて不感帯の上限を設定する。
この不感帯は、流量０からその上限値ＵＬｍａｘまでの領域であり、この上限は以下のように設定される。
【００３０】
【数４】
ＵＬｍａｘ＝ＵＬｍｉｎ＋ΔＱ
ΔＱ＝ｋ１×ｋ２×ｄｍ
【００３１】
ここで、ＵＬｍｉｎは、不感帯の最低値（０もしくは有限値）、
ｋ１は乱数発生を伴った分散処理における不完全性に依存する１以上の定数、
ｋ２は、積算時間と同じ周期を有するｃｏｓ（φ）で代表する脈動が入ってくる場合における超音波測定に必要となる時間Δｔの影響を代表する定数であり、乱数化されていない位相領域を−θから＋θとし、脈動の周期を−πから＋πとした場合に、以下のように定義される。
【００３２】
【数５】
ｋ２＝ｋ２１／ｋ２２
ｋ２１は、積分領域を−θから＋θとする∫ｃｏｓ（φ）ｄφ
ｋ２２は、積分領域を−πから＋πとする∫ｃｏｓ（φ）ｄφ
【００３３】
従って、誤差関数である統計関数から不感帯の上限値を設定するに、平均測定タイミング間隔（これは−πから＋πに包含される）及び測定所要時間（これは−θから＋θに包含される）に基づいて、その寄与分を見込んでいることとなる。
【００３４】
この構成を採用することにより、不感帯の上限ＵＬｍａｘは、瞬時流量の状態によって可変設定される。
この不感帯上限情報は、積算流量導出手段１２に送られる。
【００３５】
へ  不感帯設定切り換え手段１５
この手段は、積算流量導出手段１２において瞬時流量Ｑから積算流量Ｑｓｕｍを導出する場合に、不感帯を設定するかどうかに関する切り換えを行うための手段である。先に示したように、メーターは、不感帯を設ける不感帯設定状態と、不感帯を設けない不感帯非設定状態との両状態で動作可能であるが、この切り換え設定をこの手段が受け持ち、切り換え情報を、積算流量導出手段１２に送る。さらに、具体的には、この手段は、使用期間の始点にあって、不感帯設定状態とする切り換え情報を出力し、使用期間内に一定流量を超える流量の測定を行った場合に、不感帯非設定状態とする切り換え情報を出力する。
従って、もし、使用期間内に一定流量以下しかガス流量の測定が行われなければ、不感帯が使用期間中常に適応され、実質上、オーバーカウントを防止できる。
一方、一定量以上の使用があった場合は、不感帯の設定が行われなくなり、従来測定していなかった微小流量域まで、有意な流量とする。結果、不測定流量を発生しないようにできる。
【００３６】
イ  積算流量導出手段１２
さて、先に説明したように、この積算流量導出手段１２は、瞬時流量測定手段６から、積算時間内に入力されてくる瞬時流量Ｑの平均と積算時間との積としての積算流量Ｑｓｕｍを導出し、この積算流量Ｑｓｕｍを、所定の部位に送るのであるが、この積算流量Ｑｓｕｍの算定時に、前記不感帯の上限情報と、切り換え情報が使用される。
この処理は以下のように行われる。
イー１  不感帯非設定状態
不感帯の設定を一切行わず、積算時間内に入力されてくる瞬時流量を全て積算し、平均瞬時流量を求め、この結果に基づいて積算流量Ｑｓｕｍを求める。
イー２  不感帯設定状態
不感帯上限情報に基づいて、この上限より低い瞬時流量をすべて零瞬時流量出力として、平均瞬時流量を求め、この結果に基づいて積算流量Ｑｓｕｍを求める。この場合、不感帯の設定は行われるのであるが、数４において、ＵＬｍｉｎが有限の値の場合は、本願が対象とするΔＱの適応を行う状態と、行わない状態との両状態を選択可能としても良い。
このようにして求められた積算流量Ｑｓｕｍは逐次、使用流量導出手段１６へ送られる。
【００３７】
ト  使用流量導出手段１６
積算流量導出手段１２より逐次送られてくる積算流量Ｑｓｕｍを、所定の使用流量期間内で加算し、使用流量ΣＱｓｕｍを求め、出力する。
以上のように構成することで、不感帯の上限を測定状況によって可変することが可能となると共に、この不感帯を適応する場合と適応しない場合とを選択できる。
【００３８】
〔別実施の形態〕
（１）上記の実施の形態にあっては、超音波による流量測定をおこなうものを例に挙げたが、離散的に測定される瞬時流量から積算流量を求めるものとしては、熱伝導式流量センサーを流量測定の主要部とするものの他、電磁式等もあり、本願は、これら全てに対して適応可能である。
（２）上記の実施の形態にあっては、誤差関数は積算流量に対して定義されるものであるが、積算時間内が経過した後、得られている複数の瞬時流量から、本願の誤差関数に相当する関数を個々に瞬時流量に対して求め、これに基づいて瞬時流量自体を個々に誤差補正処理を施すことも考えられるが、この手法を採る場合は、瞬時流量の補正および、その補正後の値からの流量の導出という迂回手順をとることとなり、最終的には積算流量に対する安全機能用流量を、求めることが必要となるため、実質的に等価なものであると考えられる。
（３）上記の実施の形態にあっては、第一統計関数として標準偏差を採用しているが、積算流量のバラツキを表す統計関数としては、分散等を採用しても良い。（４）上記の実施の形態にあっては、第二統計関数として、連続する瞬時流量の差分の絶対値の平均値を採用したが、連続して測定される瞬時流量の差分からなる数列の平均値等の統計関数を第二統計関数、単位時間当たりの瞬時流量の変化量からなる数列の平均値等の統計関数を第二統計関数としても良い。更に、流量の一次微分に相当する前記関数の他、二次微分に相当する関数等も採用可能である。
（５）上記の実施の形態にあっては、積算流量と誤差関数とを線形に接続することで、安全機能用流量を得るものとしたが、誤差関数の値に対応して、積算流量に対して一定の割合となる流量を安全機能用流量としても良い。即ち、安全機能用流量は、積算流量と誤差関数から決めることができる。
（６）上記の実施の形態にあっては、不感帯の上限の設定にあたり、差分の平均ｄｍを使用したが、先に説明した標準偏差σ等の瞬時流量のバラツキを表す統計関数を使用しても良いし、連続する瞬時流量の差分からなる数列の統計関数、もしくは単位時間あたりの瞬時流量の変化量からなる数列の統計関数のいずれか１種以上、いずれを使用しても良い。
さらに  上記の実施例においては、数４に示す手順により不感帯の上限を設定するものとしたが、以下のような誤差要因を考慮してΔＱを設定するものとしても良い。
【００３９】
考慮する誤差
１  測定タイミング間隔の乱数化の不完全性に起因する誤差
ΔＱ１
この最大値は、平均流量測定間隔に同期する脈動の振幅をＡとして、
Ａ×ｓｉｎ（π×Δｔ×ｎ／Ｔ）／πとなる。
２  積算時間内におけるサンプル数が十分に大きくないことによる誤差
ΔＱ２＝Ａ×／ｓｑｒｔ（ｎ）
３  測定部を流れるガスの流速分布が定常流と異なることによる誤差
ΔＱ３
これは、ΔＱ１、ΔＱ２に対して充分に小さい。
ここで、Ａは、これまで説明してきた差分の二乗平均とすることができ、結果的に、この二乗平均を利用して、本願の誤差関数であるΔＱ１＋ΔＱ２を求め、これから不感帯の上限を設定することも可能である。
このようにして設定する場合の具体的な値は以下のようになる。
Ｔ／ｎ＝１ｓｅｃ，Δｔ＝０．１ｓｅｃ，ｎ＝４００，差分の二乗平均を１（ｍ3 ／ｈ）2 とすると、Ａ＝１ｍ3 ／ｈより、ＵＬｍｉｎ＝０として、
ΔＱ１＝０．０２８ｍ3 ／ｈ、ΔＱ２＝０．０５ｍ3 ／ｈなり、
ΔＱ＝０．０７８ｍ3 ／ｈすることができる。この値を不感帯の上限としても良い。
この誤差関数は、実際に現行のメーターに使用される不感帯の値に近い値を与え、本願手法が有効であることが判る。
【図面の簡単な説明】
【図１】本願のガスメーターのシステム構成を示す図
【図２】脈動周期が短い場合の状況を示す説明図
【図３】脈動周期が長い場合の状況を示す説明図
【符号の説明】
１      ガスメーター
２      測定管
３      導入部
４      導出部
５      送受波器
５ａ    送波器
５ｂ    受波器
６      瞬時流量測定手段
７      メーター本体
８      測定時間制御手段
８ａ    クロック
８ｂ    分散手段
９      第一統計関数導出手段
１０    第二統計関数導出手段
１１    誤差関数導出手段
１２    積算流量導出手段
１３    安全機能用流量導出手段
１４    不感帯上限設定手段
１５    不感帯設定切り換え手段
１６    使用流量導出手段
【図１】
0004165985.tif 000002
【図２】
0004165985.tif 000003
【図３】
0004165985.tif 000004