(57)【要約】      （修正有）【課題】ガウス求積法を用いて、少ない計測点でも精度よく流量を計測する方法及び装置を提供する。【解決手段】流路断面内の推定流速における断面法線方向成分ｖz0（ｘ，ｙ）の、ｙ座標に関する平均値をＶz0（ｘ）とし、Ｌ（ｘ）＝ｙmax（ｘ）−ｙmin（ｘ）として、ガウス求積法の重み関数をＶz0（ｘ）Ｌ（ｘ）とする。さらに、実際の流速の断面法線方向成分ｖz（ｘ，ｙ）の、ｙに関する平均値をＶz（ｘ）とし、被積分関数をＶz（ｘ）／Ｖz0（ｘ）とする。これらを用いて、ガウス求積法により流量の近似値を決定する。また、仮想的に流速分布のゼロ点をずらして逆流等の問題を見かけ上回避するための方法。【選択図】図５
2015078962.tif 000002
【特許請求の範囲】
【請求項１】
流路内を流れる流体の流量として、以下の式：
【数１】
2015078962.tif 000063
ただし、
2015078962.tif 000064
2015078962.tif 000065
（Ｓは流路断面を表わし、該流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、前記流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における該流体の流速のうち、該位置（ｘ，ｙ）における該法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ該流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。）
によって表わされる流量Ｑを計測する方法であって、
前記流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の前記法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、以下の式：
【数２】
2015078962.tif 000066
ただし、
2015078962.tif 000067
によって表わされるＶz0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する段階と、
各々の分点ｘiについて、Ｖz（ｘi）を決定するとともにＶz0（ｘi）を決定する段階と、
前記流量Ｑの近似値として、以下の式：
【数３】
2015078962.tif 000068
を用いて表わされる、近似値ＱNを決定する段階と
を備えた方法。
【請求項２】
流路内を流れる流体の流量として、以下の式：
【数４】
2015078962.tif 000069
ただし、
2015078962.tif 000070
2015078962.tif 000071
（Ｓは流路断面を表わし、該流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、前記流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における該流体の流速のうち、該位置（ｘ，ｙ）における該法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ該流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。）
によって表わされる流量Ｑを計測する方法であって、
前記流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の前記法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、オフセット関数Ｖz1（ｘ）を用いて以下の式：
【数５】
2015078962.tif 000072
ただし、
2015078962.tif 000073
によって表わされるＶ’z0（ｘ）を用いて、以下の式：
【数６】
2015078962.tif 000074
によって表わされるＶ’z0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する段階と、
各々の分点ｘiについて、Ｖz（ｘi）を決定するとともに、Ｖ’z0（ｘi），Ｖz1（ｘi），及び、以下の式：
【数７】
2015078962.tif 000075
2015078962.tif 000076
で表わされるＱz0，Ｑz1を決定する段階と、
以下の式：
【数８】
2015078962.tif 000077
を用いて表わされる、近似値ＱNを決定する段階と
を備えた方法。
【請求項３】
流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以上となるよう、又は流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以下となるよう、前記Ｖz1（ｘ）が選定される、請求項２に記載の方法。
【請求項４】
請求項１乃至３のいずれか１項に記載の方法により前記近似値ＱNを決定する段階と、
前記近似値ＱNよりも信頼性の低い、第２の近似値ＱSUBを用いて、以下の式：
【数９】
2015078962.tif 000078
に基づき、前記近似値ＱNの誤差範囲を推定する段階と
を備えた流量計測方法。
【請求項５】
流路内を流れる流体の流量として、以下の式：
【数１０】
2015078962.tif 000079
ただし、
2015078962.tif 000080
2015078962.tif 000081
（Ｓは流路断面を表わし、該流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、前記流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における該流体の流速のうち、該位置（ｘ，ｙ）における該法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ該流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。）
によって表わされる流量Ｑを計測する装置であって、
前記流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の前記法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、以下の式：
【数１１】
2015078962.tif 000082
ただし、
2015078962.tif 000083
によって表わされるＶz0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する、ガウス求積演算部と、
前記ガウス求積演算部により決定された各々の分点ｘiについてＶz（ｘi）を決定する、Ｖz（ｘi）決定部と、
前記ガウス求積演算部により決定された各々の分点ｘiについてＶz0（ｘi）を決定する、Ｖz0（ｘi）決定部と、
前記流量Ｑの近似値として、以下の式：
【数１２】
2015078962.tif 000084
を用いて表わされる近似値ＱNを決定する、ＱN決定部と
を備えた装置。
【請求項６】
流路内を流れる流体の流量として、以下の式：
【数１３】
2015078962.tif 000085
ただし、
2015078962.tif 000086
2015078962.tif 000087
（Ｓは流路断面を表わし、該流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、前記流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における該流体の流速のうち、該位置（ｘ，ｙ）における該法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ該流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。）
によって表わされる流量Ｑを計測する装置であって、
前記流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の前記法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、オフセット関数Ｖz1（ｘ）を用いて以下の式：
【数１４】
2015078962.tif 000088
ただし、
2015078962.tif 000089
によって表わされるＶ’z0（ｘ）を用いて、以下の式：
【数１５】
2015078962.tif 000090
によって表わされるＶ’z0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する、ガウス求積演算部と、
前記ガウス求積演算部により決定された各々の分点ｘiについて、Ｖz（ｘi）を決定する、Ｖz（ｘi）決定部と、
前記ガウス求積演算部により決定された各々の分点ｘiについてのＶ’z0（ｘi），Ｖz1（ｘi），及び、以下の式：
【数１６】
2015078962.tif 000091
2015078962.tif 000092
で表わされるＱz0，Ｑz1を決定する、Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1決定部と、
以下の式：
【数１７】
2015078962.tif 000093
を用いて表わされる近似値ＱNを決定する、ＱN決定部と
を備えた装置。
【請求項７】
前記ガウス求積演算部は、流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以上となるよう、又は流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以下となるよう、選定された前記Ｖz1（ｘ）を用いる、請求項６に記載の装置。
【請求項８】
前記Ｖz（ｘi）決定部が、
前記流路を規定する流路配管上で、前記各々の分点ｘiに対応する位置に設置される超音波素子のペアと、
前記超音波素子のペアに対して超音波の送受信を指示する、超音波送受信回路と、
前記超音波の送受信を示す信号を記録する、超音波送受信信号記録回路と、
前記超音波送受信信号記録回路に記録された、前記超音波の送受信を示す信号を用い、前記各々の分点ｘiについてＶz（ｘi）を計算する、流速計算回路と
を備えることを特徴とする、請求項５乃至７のいずれか１項に記載の装置。
【請求項９】
前記近似値ＱNよりも信頼性の低い、第２の近似値ＱSUBを用いて、以下の式：
【数１８】
2015078962.tif 000094
に基づき、前記近似値ＱNの誤差範囲を推定する誤差範囲推定部
を更に備えた、請求項５乃至８のいずれか１項に記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、典型例として、液体ロケットエンジンやその地上試験設備等、十分に整流することが難しい流路において、高精度で流量を計測するための方法及び装置に関する。
【背景技術】
【０００２】
一般に体積流量は、流路断面を通って単位時間内に流れる流体の体積と定義することができる。一例として、図１に示す流路配管１内を流れる流体の流量Ｑは、流路断面上の各点における流速のうち、流路断面の法線方向成分の積分値として以下のとおり表される。
【数１】
2015078962.tif 000003
…（１）
ただし、
2015078962.tif 000004
…（２）
2015078962.tif 000005
…（３）
【０００３】
上記式（１）〜（３）において、Ｓは流路断面を表わし、流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における流体の流速のうち、位置（ｘ，ｙ）における法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。また、図１において流路断面Ｓ1は流路配管１の軸に直交する平面であるが、流路断面が軸に対して傾いていたり、曲面であったりしても、上記式（１）〜（３）により流量を表わすことができる。
【０００４】
上記式（２）から明らかなとおり、ｘ座標としてあるｘiを定めたとき、Ｖz（ｘi）は、流路断面Ｓ内においてｘ座標がｘiに等しい領域、即ち線分ＡＢ（図２）上の各点（ｘ，ｙ）についてのｖz（ｘ，ｙ）の平均値を表しており、例えば超音波方式によって計測することが可能である。
【０００５】
一例においては、図３に示すとおり、点Ａ，Ｂからそれぞれｚ方向に所定の大きさだけシフトした点をＡ＊，Ｂ＊として、超音波が点Ａ＊から点Ｂ＊に伝播する場合と点Ｂ＊から点Ａ＊に伝播する場合とで流速に応じた見かけ上の音速差が生じることを利用して、それぞれの方向について計測された伝播時間から、Ａ＊，Ｂ＊間におけるｖz（ｘ，ｙ）の平均値を求めることができる。通常、これを線分ＡＢ上の平均流速Ｖz（ｘi）とみなすことができる。
【０００６】
具体的には、図４に示すとおり、点Ａ＊から超音波を送信する超音波素子２ａと、点Ｂ＊から超音波を送信する超音波素子２ｂとを流路配管１上に配置し、超音波素子２ａから送信された超音波が超音波素子２ｂによって受信されるまでの伝播時間と、超音波素子２ｂから送信された超音波が超音波素子２ａによって受信されるまでの伝播時間とを計測する。単純な一例においては、特許文献１に記載されているとおり、超音波素子２ａから超音波素子２ｂに対して超音波を送信したとき、長さＬの経路を通って超音波素子２ｂに到達するまでの伝播時間をｔ1とし、超音波素子２ｂから超音波素子２ａに対して超音波を送信したとき超音波素子２ａに到達するまでの伝播時間をｔ2として、以下の式（４）により、ｖz（ｘ，ｙ）の平均値を決定することができる（超音波の送受信経路とｚ方向との間の角度をθとする。）。
【数２】
2015078962.tif 000006
…（４）
【０００７】
通常、この平均値をＶz（ｘi）として用いることができる。
【０００８】
実用上、計測できるＶz（ｘi）は離散的であることから、多くの場合には数値積分を用いてＱの近似値を求める必要がある。また、図４に示すとおり超音波素子のペアを用いて計測する場合には、計測点ｘiの数に応じた数だけ、超音波素子のペアを流路配管上に配置することが考えられるが、この場合、実装できる計測点には限りがあるため、少ない計測点でも精度のよい計算手法を用いることが必要である。
【０００９】
上記式（１）で表わされる流量Ｑの近似値を、離散的に計測した流速からガウス求積法（ガウス−ルジャンドル，チェビチェフ，ロバート）によって求める方法が、特許文献２に開示されている。特許文献２が開示する方法は、被積分関数をＶz（ｘ）Ｌ（ｘ）とし、重み関数を１（定数）としてガウス求積法を用いることにより、以下の近似式
【数３】
2015078962.tif 000007
…（５）
によって流量Ｑの近似値ＱNを決定するものである。ここで、ｘi，ｗiはそれぞれガウス求積法の分点、係数であり、各分点ｘ＝ｘiにおいて計測された離散的な流速Ｖz（ｘi）を用いて上記式（５）が計算される。
【００１０】
特許文献３，４においては、半径Ｒの円管流路を流れる流体の流量の近似値を離散的に計測した流速から求める方法が開示されている。半径Ｒの円管流路では、流路断面内でｘ座標を定めたとき、上記式（３）に示すＬ（ｘ）が
【数４】
2015078962.tif 000008
…（６）
と表わされる。これを利用し、特許文献３，４が開示する方法においては重み関数を
【数５】
2015078962.tif 000009
…（７）
とし、被積分関数を２ＲＶz（ｘ）として、ガウス求積法（ガウス−チェビチェフ）に従い分点ｘi，係数ｗiが計算される。各分点ｘ＝ｘiにおいて計測された離散的な流速Ｖz（ｘi）を用いて、以下の近似式
【数６】
2015078962.tif 000010
…（８）
により流量Ｑの近似値ＱNが決定される。
【００１１】
従来技術においては被積分関数に流速分布そのものが含まれているため、流速分布が大きく歪んでいる場合等では、精度のよい近似値を得るために、より高次の求積法（即ち、より多くの計測点）が必要となるという課題がある。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００１−２０８５８４号公報
【特許文献２】米国特許第３５６４９１２号明細書
【特許文献３】米国特許第３９４０９８５号明細書
【特許文献４】特開昭５１−１２９２６１号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
以上に鑑み、本発明は、少ない計測点でも精度よく流量を計測する方法及び装置を提供することを課題とする。
【課題を解決するための手段】
【００１４】
上記課題を解決するため、本発明は、流路内を流れる流体の流量として、以下の式：
【数７】
2015078962.tif 000011
…（９）
ただし、
2015078962.tif 000012
…（１０）
2015078962.tif 000013
…（１１）
（Ｓは流路断面を表わし、流路断面Ｓ内の位置を（ｘ，ｙ）で表わす。ｖzは、流体の流速のうち、流路断面Ｓの法線方向（ｚ方向）成分を表わし、ｖz（ｘ，ｙ）は、流路断面Ｓ内位置（ｘ，ｙ）における流体の流速のうち、位置（ｘ，ｙ）における法線方向の成分を表わす。ｙmin（ｘ），ｙmax（ｘ）は、それぞれ流路断面Ｓ内でｘ座標を定めたときのｙ座標の最小値、最大値を表わし、ｘmin，ｘmaxは、それぞれ流路断面Ｓ内でのｘ座標の最小値、最大値を表わす。後述の各式においても同様。）
によって表わされる流量Ｑを計測する方法であって、
流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、以下の式：
【数８】
2015078962.tif 000014
…（１２）
ただし、
2015078962.tif 000015
…（１３）
によって表わされるＶz0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する段階と、
各々の分点ｘiについて、Ｖz（ｘi）を決定するとともにＶz0（ｘi）を決定する段階と、
流量Ｑの近似値として、以下の式：
【数９】
2015078962.tif 000016
…（１４）
を用いて表わされる、近似値ＱNを決定する段階と、を備えた方法を提供する。
【００１５】
上記のとおり、本発明の方法は、流路断面内Ｓの推定流速分布（理論計算、シミュレーション、予備実験等により、予め求めておくことができる。）を用いて重み関数Ｖz0（ｘ）Ｌ（ｘ）を構成し、ガウス求積法により流量の近似値を決定することを特徴とする。このとき、被積分関数は
【数１０】
2015078962.tif 000017
…（１５）
となり、あるｘ座標についての、推定流速分布の平均値に対する実際の流速分布の平均値の比を表わしている。推定流速分布が確からしいほど被積分関数は定数関数に近くなるため、近似精度が向上する。
【００１６】
上記方法は、適当なオフセット関数Ｖz1（ｘ）によってＶz0（ｘ）を修正することで得られるＶ′z0（ｘ）を用いても実施可能である。
【００１７】
すなわち本発明は、流路内を流れる流体の流量として、上記式（９）〜（１１）によって表わされる流量Ｑを計測する方法であって、
流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、オフセット関数Ｖz1（ｘ）を用いて以下の式：
【数１１】
2015078962.tif 000018
…（１６）
ただし、
2015078962.tif 000019
…（１７）
で表わされるＶ’z0（ｘ）を用いて、以下の式：
【数１２】
2015078962.tif 000020
…（１８）
によって表わされるＶ’z0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する段階と、
各々の分点ｘiについて、Ｖz（ｘi）を決定するとともに、Ｖ’z0（ｘi），Ｖz1（ｘi），及び、以下の式：
【数１３】
2015078962.tif 000021
…（１９）
2015078962.tif 000022
…（２０）
で表わされるＱz0，Ｑz1を決定する段階と、
以下の式：
【数１４】
2015078962.tif 000023
…（２１）
を用いて表わされる、近似値ＱNを決定する段階とを備えた方法を提供する。
【００１８】
この方法は、オフセット関数Ｖz1（ｘ）を用いることにより、後述のとおりＶz0（ｘ），Ｖz（ｘ）のゼロ点をずらしてＶ′z0（ｘ），Ｖ’z（ｘ）を導入し、これらを用いて、ゼロ点がずらされた流量Ｑ’にガウス求積法を適用するという発想に基づく。後述のとおり、上記式（９）中、Ｖz（ｘ）の代わりにＶ′z（ｘ）を用いてゼロ点をずらした流量Ｑ′を導入し、この流量Ｑ’に、重み関数をＶ′z0（ｘ）Ｌ（ｘ）としたガウス求積法を適用して近似値Ｑ’Nを多項式近似の形で表し、更にＱ’NとＱNとの間の関係式を用いてＱNを多項式近似で表すという手順で、上記式（２１）を導出することが可能である。
【００１９】
上記のＶz1（ｘ）を用いる方法においては、流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以上となるよう、又は流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以下となるよう、前記Ｖz1（ｘ）を選定することが好ましい。
【００２０】
いわゆる逆流等によって重み関数の符号が一定とはならない場合、ガウス求積法の分点ｘiを必ずしも積分区間内（流路断面Ｓ内のｘ座標がとる範囲）に取ることができないため、Ｖz（ｘi）の計測点ｘiやＶ0（ｘi）の計算点ｘiを実装できない可能性がある。この場合であっても、重み関数Ｖ′z0（ｘ）Ｌ（ｘ）の符号が積分区間内で一定となるよう（例えば、Ｖz0（ｘ）の最小値の絶対値をＶz1（ｘ）としたり、Ｖz0（ｘ）の最大値の絶対値×（―１）をＶz1（ｘ）としたりすれば、重み関数Ｖ′z0（ｘ）Ｌ（ｘ）の符号を一定にできる。）Ｖz1（ｘ）を選択すれば、ガウス求積法の原理から分点ｘiも必ず積分区間内に実装できるため、見かけ上逆流を回避して高精度で流量を決定することが可能となる。
【００２１】
ただし、後述のとおり、重み関数の符号が積分区間内で変化する場合であっても、積分区間内に分点ｘiをとることができる場合もある。したがって、上記のとおり重み関数の符号が一定となるようＶz1（ｘ）を選定することは本発明にとって必須ではない。
【００２２】
また本発明は、上述の本発明が提供する方法により近似値ＱNを決定する段階と、近似値ＱNよりも信頼性の低い、第２の近似値ＱSUBを用いて、以下の式：
【数１５】
2015078962.tif 000024
…（２２）
に基づき、近似値ＱNの誤差範囲を推定する段階とを備えた流量計測方法を提供する。後述のとおり、このような方法で近似値ＱNの誤差を見積もることが可能である。ここにおいて、「近似値ＱNよりも信頼性の低い」第２の近似値ＱSUBとは、本発明の方法を用いる場合には計測点数Ｎを減らしたり、あるいは後述のとおり、ガウス求積法の被積分関数を各分点ｘiについて個別に計算するのではなく、その最大値又は最小値で一律に置き換えたりして決定された、近似値ＱNよりも誤差が大きいと推定される近似値である。
【００２３】
また本発明は、流路内を流れる流体の流量として、上記式（９）〜（１１）式によって表わされる流量Ｑを計測する装置であって、流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、上記式（１２），（１３）によって表わされるＶz0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する、ガウス求積演算部と、ガウス求積演算部により決定された各々の分点ｘiについてＶz（ｘi）を決定する、Ｖz（ｘi）決定部と、ガウス求積演算部により決定された各々の分点ｘiについてＶz0（ｘi）を決定する、Ｖz0（ｘi）決定部と、流量Ｑの近似値として、上記式（１４）を用いて表わされる近似値ＱNを決定する、ＱN決定部と、を備えた装置を提供する。
【００２４】
また本発明は、流路内を流れる流体の流量として、上記式（９）〜（１１）によって表わされる流量Ｑを計測する装置であって、流路断面Ｓ内の各位置（ｘ，ｙ）における、流速の法線方向成分の推定値をｖz0（ｘ，ｙ）としたとき、オフセット関数Ｖz1（ｘ）を用いて上記式（１６），（１７）によって表わされるＶ’z0（ｘ）を用いて、上記式（１８）によって表わされるＶ’z0（ｘ）Ｌ（ｘ）を重み関数として、１からＮ（Ｎは１以上の任意の自然数）までの自然数ｉについてガウス求積法の分点ｘi及び係数ｗiを決定する、ガウス求積演算部と、ガウス求積演算部により決定された各々の分点ｘiについて、Ｖz（ｘi）を決定する、Ｖz（ｘi）決定部と、ガウス求積演算部により決定された各々の分点ｘiについてのＶ’z0（ｘi），Ｖz1（ｘi），及び、上記式（１９），（２０）で表わされるＱz0，Ｑz1を決定する、Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1決定部と、上記式（２１）を用いて表わされる近似値ＱNを決定する、ＱN決定部とを備えた装置を提供する。
【００２５】
上記のＶz1（ｘ）を用いる装置においては、前記ガウス求積演算部が、流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以上となるよう、又は流路断面Ｓ内でＶ′z0（ｘ）が常にゼロ以下となるよう選定された前記Ｖz1（ｘ）を用いることが好ましい。ただし、既に述べたとおり重み関数の符号が積分区間内で変化する場合であっても、積分区間内に分点ｘiをとることができる場合もあるため、重み関数の符号が一定となるようＶz1（ｘ）を選定することは必須ではない。
【００２６】
上記本発明が提供する装置において、Ｖz（ｘi）決定部を、流路を規定する流路配管上で各々の分点ｘiに対応する位置に設置される超音波素子のペアと、超音波素子のペアに対して超音波の送受信を指示する超音波送受信回路と、超音波の送受信を示す信号を記録する超音波送受信信号記録回路と、超音波送受信信号記録回路に記録された、超音波の送受信を示す信号を用い、各々の分点ｘiについてＶz（ｘi）を計算する流速計算回路とを備えるよう構成することができる。
【００２７】
また、上記本発明が提供する装置に、近似値ＱNよりも信頼性の低い第２の近似値ＱSUBを用いて、上記式（２２）に基づき近似値ＱNの誤差範囲を推定する、誤差範囲推定部を更に備えることができる。
【発明の効果】
【００２８】
本発明によれば、任意の流速分布に対して最適なガウス求積法を構成することができるため、少ない計測点でも精度よく流量を計測することができる。本発明の方法、装置は、整流されていない歪んだ流速分布に対して特に効果的である。
【００２９】
本発明が用いるガウス求積法において、被積分関数は、流速分布そのものではなく、推定流速分布に対する実際の流速分布の比を表している。ある程度確からしい流速分布を予め推定することができれば、流速分布が大きく歪んでいる場合等でも被積分関数は定数関数に近くなるため、精度のよい流量の近似値を得ることができる。
【００３０】
さらに、本発明が教示する誤差範囲推定の手法によれば、実際に計測した流量の不確かさ（積分の近似誤差範囲に相当）について、簡易的に見当を付けることができる。これにより、必ずしも基準器との比較を行わずに、積分誤差に起因する流量の不確かさを評価することが可能となる。ただし、当然ながら、可能ならば基準器との比較検定を行う方が確実である。なお、この誤差範囲推定手法は従来技術にも適用可能である。
【図面の簡単な説明】
【００３１】
【図１】本発明による流量計測を実施可能な流路の一例を示す図である。
【図２】図１に示される流路の流路断面Ｓ1を示す図である。
【図３】図１に示される流路において、超音波が伝播する経路Ａ＊Ｂ＊を示す図である。
【図４】図１に示される流路において、超音波が伝播する経路Ａ＊Ｂ＊を、流路の軸に直交する方向から見たときの図である。
【図５】流路配管１に接続された、本発明に係る流量計測装置の一実施形態を示す構成図である。なお、超音波素子２ａ，２ｂのペアは、計測点ｘiの数に応じた数の組だけ流路配管１上に設置されるが、図５においては代表して１組だけ描かれている。
【図６】円管流路における流速分布の例を示すグラフである。
【図７】Ｎ＝３として本発明による流量計測を実施するときの、超音波の送受信地点Ａi＊，Ｂi＊を示す図である。ただし、計測点ｘ1〜ｘ3の位置は正確ではない。またＡi＊，Ｂi＊は、図３と同様に流路断面Ｓ1からｚ方向へ所定の大きさだけシフトされているが、ｚ座標は図７中に示されていない。
【図８】流路配管１に接続された、本発明に係る流量計測装置の一実施形態を示す構成図である。なお、超音波素子２ａ，２ｂのペアは、計測点ｘiの数に応じた数の組だけ流路配管１上に設置されるが、図８においては代表して１組だけ描かれている。
【図９Ａ】図８中の流路配管１の一例として、邪魔板が設けられた円管流路を示す図である。
【図９Ｂ】邪魔板をｘｙ平面内で見たときの平面図である。
【図１０】従来のガウス求積法（ガウス‐チェビチェフ）によって、図９Ａ中のＡ地点の断面で流量を計測した時の誤差と、本発明の方法によって同じくＡ地点の断面で流量を計測した時の誤差とを示すグラフである。
【図１１】流路断面の別の例として、流路の軸に対して角度φだけ傾いた面Ｓ2を示した図である。
【図１２】流路断面Ｓ2を、その法線方向から見たときの図である。
【図１３】流路断面Ｓ2を用いる場合に、流速ｖの法線方向成分ｖzを求めるための方法を示す図である。
【発明を実施するための形態】
【００３２】
これより図面を用いて、本発明の流量計測方法、流量計測装置の具体例を説明する。ただし、本発明の流量計測方法、流量計測装置は、以下に説明する具体的態様に限定されるわけではなく、本発明の範囲内で適宜変更可能である。例えば、実施例１においては半径１の円管流路を流れる流体の流量を計測するが、これ以外の任意の流路を流れる流体の流量であっても、本発明に従って計測することが可能である（ガウス求積法の積分区間は慣例的に［−１，１］とされているが、流路断面のｘ区間がこれとは異なる場合であっても、適切な変数変換により積分区間を［−１，１］とすることが可能である。）。また、以下の各実施例では、流速の流路断面法線方向成分の平均値を計測する手段として、流路配管上に１計測点ごとに２つの超音波素子を設置し、超音波素子間での双方向における超音波伝播時間の差から流量を決定する方式（伝播時間差方式）を用いるが、この方式に限らず、不透明な液体に適したパルスドップラー方式など、任意の方式を用いてよい。超音波を用いる方式以外の方式として、公知のあらゆる方式、及び将来的に実施可能となるあらゆる方式を用いることも可能である。さらに、以下の実施例では流量Ｑを体積流量として計測しているが、計測された体積流量に密度を乗じるなどして質量流量を求めることも可能である。そのような質量流量の計測は必然的に体積流量の計測を伴うものであり、即ち本発明の技術的範囲内にある。また、図５，図８で示される各回路は、物理的に隔離された別個の回路であってもよいし、１つの回路中で各機能を担う概念的機能部であってもよい。各回路は複数の回路から構成されてもよく、例えば図８中のＶ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路は、Ｖ’z0（ｘi）を計算するための回路と、Ｖz1（ｘi）を計算するための回路と、Ｑz0を計算するための回路と、Ｑz1を計算するための回路とからなる回路群であってもよい。本発明を実施するために行われるそれぞれの処理を、１以上の任意の要素に対して任意に分配することが可能である。
【００３３】
本明細書における、ガウス求積法
本明細書、及び特許請求の範囲中で用いられる「ガウス求積法」という用語は、重み関数を用いて有限個の分点及び係数を決定し、これらを用いた多項式近似として実行される数値積分法一般を指す総称である。分点及び係数の具体的な値は、重み関数等の条件に依存する。以下、このような一般的意味におけるガウス求積法を説明する。
【００３４】
（１）重み関数をｗ（ｘ），積分区間をＸとし、
【数１６】
2015078962.tif 000025
…（２３）
とおく。
【００３５】
（２）以下のｐk：
【数１７】
2015078962.tif 000026
…（２４）
2015078962.tif 000027
…（２５）
を係数とする直交多項式：
【数１８】
2015078962.tif 000028
…（２６）
に対して、ＰN（ｘ）＝０となる根をガウス求積法の分点ｘi（ｉ＝１，２，…，Ｎ）とする。
【００３６】
（３）被積分関数がＮ次未満の多項式である場合に正しい積分値が得られるよう、以下の式（２７）によってガウス求積法の係数ｗiを求める。
【数１９】
2015078962.tif 000029
…（２７）
（４）一般の被積分関数ｆ（ｘ）に対して、以下のとおり有限和で積分値を近似する。
【数２０】
2015078962.tif 000030
…（２８）
誤差項ＲNは２Ｎ次導関数を含む形で表されるため、被積分関数ｆ（ｘ）が２Ｎ次未満の多項式である場合、近似誤差は０となる。
【実施例１】
【００３７】
本発明による流量計測方法
本発明の一実施形態による流量計測方法によれば、上述のガウス求積法において重み関数ｗ（ｘ）をＶz0（ｘ）Ｌ（ｘ）とし（式（１２））、式（２３）〜（２７）を用いて分点ｘi及び係数ｗiを決定し、被積分関数ｆ（ｘ）をＶz（ｘ）／Ｖz0（ｘ）（式（１５））として式（２８）で演算することにより、上記式（１４）に示される近似値ＱNを得ることができる。
【００３８】
本発明による誤差範囲推定
また本発明の一実施形態によれば、以下のとおり誤差Ｑ−ＱNの範囲を推定することが可能である。
【００３９】
上記本発明の一実施形態による流量計測方法に従い、最も確からしいと考えられる流量Ｑの近似値ＱNを求める。併せて、ＱNよりも信頼性の低い近似値ＱSUB、すなわち以下の式：
【数２１】
2015078962.tif 000031
…（２９）
が成り立つと推定される近似値ＱSUBを、利用する計測点数を減らした上で本発明の流量計測方法を用いたり、あるいは別の任意の方法を用いたりして求める。
【００４０】
ここで、まずＱSUB＜ＱNである場合について考えると、
（１）ＱN≦Ｑのとき、
【数２２】
2015078962.tif 000032
…（３０）
より、以下の式が得られる。
【数２３】
2015078962.tif 000033
…（３１）
（２）ＱN＞Ｑのとき、
まず、仮にＱSUB≧Ｑとすると、ＱSUB＜ＱNからＱSUB−Ｑ＜ＱN−Ｑとなり、ＱSUBがＱNよりも信頼性が低いという仮定に矛盾するので、ＱSUB＜Ｑとする。このとき、当該仮定から
【数２４】
2015078962.tif 000034
…（３２）
なので、上記式（３２）の両辺にＱNを加えて整理することにより上記式（３１）が得られる。
【００４１】
次に、ＱSUB＞ＱNである場合について考えると、
（１）ＱN≦Ｑのとき、
まず、仮にＱSUB≦Ｑとすると、ＱSUB＞ＱNからＱ−ＱSUB＜Ｑ−ＱNとなり、ＱSUBがＱNよりも信頼性が低いという仮定に矛盾するので、ＱSUB＞Ｑとする。このとき、当該仮定から
【数２５】
2015078962.tif 000035
…（３３）
なので、上記式（３３）の両辺からＱNを引いて整理することにより、以下の式：
【数２６】
2015078962.tif 000036
…（３４）
が得られる。
（２）ＱN＞Ｑのとき、
【数２７】
2015078962.tif 000037
…（３５）
より、上記式（３４）が得られる。
【００４２】
以上をまとめると、
【数２８】
2015078962.tif 000038
…（３６）
として、ＱNの誤差範囲を推定することができる。
【００４３】
以下、上記式（３６）を用いた誤差範囲推定の一例を示す。本実施例１の流量計測方法に従い、Ｎ個の計測点で計測したＶz（ｘ）を用いて得られる流量の近似値ＱNよりも信頼性の低い近似値としては、Ｎ個の計測点ｘiについてそれぞれ求められる、Ｎ個の被積分関数値Ｖz（ｘi）／Ｖz0（ｘi）の最小値、最大値をそれぞれｆmin，ｆmaxとして、
【数２９】
2015078962.tif 000039
2015078962.tif 000040
…（３７）
で表されるＱmin，Ｑmaxを考えることができる。Ｑmin＜ＱN＜Ｑmaxであり、また、実用上十分な確からしさで、ＱNはＱmin，Ｑmaxよりも精度が高いものと考えられる。よって、上記（３６）を用いれば、以下の式：
【数３０】
2015078962.tif 000041
…（３８）
によりＱNの誤差範囲を概算することができる。
【００４４】
本発明による流量計測装置
本発明の流量計測方法を実施することができる装置の一例を、図５に示す。図５の流量計測装置は、ガウス求積演算回路、流路配管１に設置された少なくとも１組の超音波素子２ａ，２ｂのペア、超音波送受信回路、超音波送受信信号記録回路、流速計算回路、Ｖz0（ｘi）計算回路、ＱN決定回路、及び任意で誤差範囲推定回路を備える。これらの要素のうち、少なくとも１組の超音波素子２ａ，２ｂのペア、超音波送受信回路、超音波送受信信号記録回路、流速計算回路によってＶz（ｘi）決定部が構成される。
【００４５】
ガウス求積演算回路には、推定流速分布として、流路断面内の各位置（ｘ，ｙ）における流速の、法線方向成分の推定値ｖz0（ｘ，ｙ）が入力される。図１に示されるとおり、流路配管１の軸に対して直交する流路断面Ｓ1を用いる場合には、法線方向（ｚ方向）は軸方向に一致するのであり、推定流速分布としては、流路断面Ｓ1内の各位置（ｘ，ｙ）における軸方向の流速の推定値ｖz0（ｘ，ｙ）が入力される。ガウス求積演算回路は、入力されたｖz0（ｘ，ｙ）から重み関数ｗ（ｘ）＝Ｖz0（ｘ）Ｌ（ｘ）を構成し、上述のガウス求積法により分点ｘi及び係数ｗiを決定する。
【００４６】
超音波素子２ａ，２ｂは、図４を用いて既に説明したとおり、流路配管１上に傾斜しつつ対向して配置された超音波素子のペアである。具体的には、ガウス求積演算回路により決定された、１からＮまでのｉに対する各々の分点ｘiについて、図２に示すとおりｙ座標の最大値、最小値に対応する流路断面上の点Ａ，Ｂが決定され、図３に示すとおり、点Ａ，Ｂからそれぞれｚ方向に所定の大きさだけシフトした点をＡ＊，Ｂ＊として、Ａ＊，Ｂ＊間で超音波の送受信を行うことができる位置に超音波素子２ａ，２ｂのペアが設置される。すなわち、超音波素子２ａ，２ｂのペアはＮ組設置されることになる。あるいは、予め多数組の超音波素子２ａ，２ｂペアを流路配管１上に設置しておき、ガウス求積演算回路により決定された分点ｘiに対応する位置にある超音波素子２ａ，２ｂペアのみを超音波送受信回路により動作させてもよい。
【００４７】
超音波送受信回路は、超音波素子２ａ，２ｂに超音波の送受信を指示する信号を送信し、超音波素子２ａ，２ｂが受信した超音波信号を受信するとともに、超音波送受信信号記録回路に対しては超音波の送受信を示す信号を送信する回路である。
【００４８】
超音波送受信信号記録回路は、超音波送受信回路から送信された超音波の送受信を示す信号を記録する回路である。一例においては、上記式（４）中のｔ1，ｔ2を示す信号が超音波送受信回路から送信されて、超音波送受信信号記録回路に記録される。
【００４９】
流速計算回路は、超音波送受信信号記録回路から超音波の送受信を示す信号を受信し、これを演算して、各々の分点ｘiに対するＶz（ｘi）を計算する回路である。一例において、流速計算回路は、超音波送受信信号記録回路から上記式（４）中のｔ1，ｔ2を示す信号を受信して、上記式（４）に基づいてｖz（ｘ，ｙ）の平均値を計算により決定し、この平均値をＶz（ｘi）として出力する。なお、上記式（４）中のＬやθは、例えばガウス求積演算回路により分点ｘi、係数ｗiとともに幾何学計算されて、流速計算回路に入力されるものとする。
【００５０】
Ｖz0（ｘi）計算回路は、ガウス求積演算部により決定された各々の分点ｘiについて、推定流速分布を用いてＶz0（ｘi）を計算により決定する。各々のｘ座標に対応するｙmin（ｘ），ｙmax（ｘ）は流路断面形状により決定されるため、Ｖz0（ｘi）計算回路に予め入力しておくことが可能である。分点ｘiはガウス求積演算部から入力される。またｖz0（ｘ，ｙ）は、あらかじめ入力されていても、ガウス求積演算部から入力されてもよい。
【００５１】
ＱN決定回路は、ガウス求積演算回路により決定された係数ｗi、流速計算回路により決定されたＶz（ｘi）、及びＶz0（ｘi）計算回路によって決定されたＶz0（ｘi）の、各々の演算回路、計算回路からの入力を受け、これらを用いて、式（１４）により流量の近似値ＱNを計算により決定する回路である。
【００５２】
なお、ガウス求積法を初めとして各回路により行われる演算には解析解が存在しないことが多いため、典型的には計算機等による数値的な演算が行われる。
【００５３】
誤差範囲推定回路は、ＱN決定回路からＱNの入力を受けるとともに、ＱNよりも信頼性の低い近似値ＱSUBの入力をＱN決定回路又は任意の外部回路等（不図示）から受けて、上記式（３６）により誤差範囲を推定する回路である。
【００５４】
実施例１の具体例
次に、半径１の円管流路において、実際の流速分布が以下の式（３９），（４０）：
（層流）
【数３１】
2015078962.tif 000042
…（３９）
（指数則）
2015078962.tif 000043
ただし、ｎ＝５，７，９
…（４０）
で表わされるそれぞれの場合（図６）に、Ｎ＝３として本発明の方法により流量の近似値Ｑ3を得るための、流量計測装置の動作を説明する。
【００５５】
本実施例１においては、流路断面として図１等に示される流路断面Ｓ1を用い、流速の法線方向成分の推定値ｖz0（ｘ，ｙ）が、
【数３２】
2015078962.tif 000044
…（４１）
として、予め与えられているとする。この推定値ｖz0（ｘ，ｙ）がガウス求積演算回路に入力される。ガウス求積演算回路は、入力されたｖz0（ｘ，ｙ）から重み関数：
【数３３】
2015078962.tif 000045
…（４２）
を構成し、上述のガウス求積法により分点ｘi及び係数ｗiを決定する。本具体例においてはＮ＝３であるため、ｘ1〜ｘ3及びｗ1〜ｗ3が決定される。
【００５６】
上記決定された各々の分点ｘ1〜ｘ3に対応して、超音波の送受信点Ａ1＊〜Ａ3＊，Ｂ1＊〜Ｂ3＊が決定される（図７）。具体的には、各々の分点ｘ1〜ｘ3について図２に示す要領で点Ａ1〜Ａ3，Ｂ1〜Ｂ3を決定し（図２のｘiを分点ｘ1として得られる、図２中のＡ，ＢをＡ1，Ｂ1とし、図２のｘiを分点ｘ2として得られる、図２中のＡ，ＢをＡ2，Ｂ2とし、図２のｘiを分点ｘ3として得られる、図２中のＡ，ＢをＡ3，Ｂ3とする。）、これらの点をｚ方向に所定の大きさだけシフトさせることにより（図３参照。点Ａiと点Ｂiを逆方向にシフトさせる。）、点Ａ1＊〜Ａ3＊，Ｂ1＊〜Ｂ3＊が決定される。Ａ1＊，Ｂ1＊間、Ａ2＊，Ｂ2＊間、及びＡ3＊，Ｂ3＊間のそれぞれで超音波の送受信を行うことができる位置に、３組の超音波素子２ａ，２ｂのペアが配置される。超音波素子の配置は、ガウス求積演算回路に接続された任意のアクチュエータ等（不図示）によって行うことができる。既に述べたとおり、予め多数組の超音波素子２ａ，２ｂペアを流路配管１上に設置しておき、ガウス求積演算回路により決定された分点ｘiに対応する位置にある超音波素子２ａ，２ｂペアのみを超音波送受信回路により動作させるならば、超音波素子の配置操作は不要である。
【００５７】
次に、超音波送受信回路が、上記３組の超音波素子２ａ，２ｂのペアに対して、超音波の送受信を指示する信号を送信する。３組の超音波素子２ａ，２ｂは、それぞれＡ1＊，Ｂ1＊間、Ａ2＊，Ｂ2＊間、及びＡ3＊，Ｂ3＊間で双方向に超音波の送受信を行い、各超音波素子が受信した超音波信号を超音波送受信回路へと送信する。超音波送受信回路は、受信した超音波信号から、各送受信経路における、それぞれの方向についての伝播時間ｔ1，ｔ2を計測し、超音波の送受信を示す信号として超音波送受信信号記録回路に送信する。
【００５８】
超音波送受信信号記録回路は、上記各送受信経路における、それぞれの方向についての伝播時間ｔ1，ｔ2を記録する。流速計算回路は、超音波送受信信号記録回路からこれら伝播時間ｔ1，ｔ2を受信して、上記式（４）に従い、分点ｘ1〜ｘ3に対するＶz（ｘ1）〜Ｖz（ｘ3）を計算する。
【００５９】
また、Ｖz0（ｘi）計算回路は、分点ｘ1〜ｘ3について、上記式（４１）の推定流速分布を用いてＶz0（ｘ1）〜Ｖz0（ｘ3）を計算により決定する（この場合、ｙmax（ｘ）＝√（１−ｘ2）であり、ｙmin（ｘ）＝−√（１−ｘ2）である。）。ＱN決定回路は、ガウス求積演算回路により決定された係数ｗ1〜ｗ3、流速計算回路により決定されたＶz（ｘ1）〜Ｖz（ｘ3）、及びＶz0（ｘi）計算回路によって決定されたＶz0（ｘ1）〜Ｖz0（ｘ3）を用いて、式（１４）により流量の近似値Ｑ3を計算により決定する。
【００６０】
誤差範囲推定回路は、ＱN決定回路からＱ3の入力を受けるとともに、ＱNよりも信頼性の低い近似値ＱSUBの入力をＱN決定回路又は任意の外部回路等から受けて、上記式（３６）により誤差範囲を推定する。近似値ＱSUBは、上記各回路を用いて３以外のＮについて決定された流量の近似値であってもよいし、本発明とは異なる任意の方法により決定された流量の近似値であってもよい。
【００６１】
本発明による流量計測方法を用いた実験
以上の動作のうち、誤差範囲推定以外の動作について、計算機シミュレーションによる模擬実験を行った。ただし、超音波素子２ａ，２ｂを用いた伝播時間ｔ1，ｔ2の計測は行わず、上記（３９），（４０）式のｖz（ｘ，ｙ）を用いてＶz（ｘi）を直接計算した。併せて、従来のガウス求積法により同様に流量の近似値を計算した。以下、それらの計算結果を比較する。
【００６２】
表１に、従来のさまざまなガウス求積法、及び本発明の方法において、Ｎ＝３について計算された分点ｘi及び係数ｗiを示す。なお、本発明の求積法では、上記式（４１）のｖz0（ｘ，ｙ）を用いて、上記式（１２），（１３）に従い重み関数Ｖz0（ｘ）Ｌ（ｘ）を構成し、上記式（２３）〜（２７）を用いて分点ｘi及び係数ｗiを決定した。
【表１】
2015078962.tif 000046
【００６３】
表１中、Ｎ＝３のときのｘ2，ｘ3は正負が逆の関係にあり（例えばガウスールジャンドルの従来例ではｘ1＝０．００００，ｘ2，ｘ3＝±０．７７４６０）、またｗ2＝ｗ3である（ガウス−ルジャンドルの従来例では、ｗ1＝０．８８８８９，ｗ2＝ｗ3＝０．５５５５６。他の従来例、本発明についても同様）。なお、ロバートの求積法における両端の固定分点は、被積分関数をゼロとみなすことができるため表１中に含まれていない。
【００６４】
表１の分点ｘi及び係数ｗiを用いて、各従来例、及び本発明それぞれの方法を用いて流量の近似値Ｑ3を決定した。さらに、実際の流量Ｑを上記式（３９），（４０）から計算して、それぞれの方法で決定された近似値Ｑ3の実際の誤差を計算した。これらの結果を表２に示す。
【表２】
2015078962.tif 000047
【００６５】
従来技術と本発明の技術に対する積分誤差の理論値を比較すると、ほとんどの場合に本発明の技術が優位である。特に、流速分布がｖz0（ｘ，ｙ）に近い場合、その効果が大きい。
【実施例２】
【００６６】
既に述べたとおり、逆流等によって重み関数の符号が一定とはならない場合、ガウス求積法の分点ｘiを必ずしも積分区間内にとれない場合がある。このときに、仮想的に流速分布のゼロ点をずらして逆流を見かけ上回避するための方法を、以下に説明する。
【００６７】
この方法は、適当な関数Ｖz1（ｘ）を用いて、
【数３４】
2015078962.tif 000048
…（４３）
2015078962.tif 000049
…（４４）
と、Ｖz0（ｘ），Ｖz（ｘ）のゼロ点をずらし、特にＶ′z0（ｘ）が流路断面Ｓ内で常にゼロ以上となるように、又は流路断面Ｓ内で常にゼロ以下となるようにＶz1（ｘ）を選定することで重み関数の符号を一定に保つ、という発想に基づく。ここで、
【数３５】
2015078962.tif 000050
2015078962.tif 000051
…（４５）
とする。
【００６８】
これらゼロ点をずらしたＶ′z0（ｘ），Ｖ′z（ｘ）を用いて、ゼロ点をずらした流量Q’を
【数３６】
2015078962.tif 000052
…（４６）
と構成し、更に流量Ｑ’に対してガウス求積法を適用するための被積分関数を
【数３７】
2015078962.tif 000053
…（４７）
と、重み関数を
【数３８】
2015078962.tif 000054
…（４８）
と構成することができる。上記式（４８）の重み関数を用いてガウス求積法の分点ｘi及び係数ｗiを決定すれば、これらを用いて流量Ｑ’の近似値Ｑ’Nを
【数３９】
2015078962.tif 000055
…（４９）
と表わすことができる。
【００６９】
また、上記式（４４）〜（４６）より
【数４０】
2015078962.tif 000056
…（５０）
が得られるので、Ｑ，Ｑ’をそれぞれ近似値ＱN，Ｑ’Nで置き換えることにより、
【数４１】
2015078962.tif 000057
…（５１）
が得られる。さらに、上記式（４４），（４９），（５１）から
【数４２】
2015078962.tif 000058
…（５２）
が得られる。
【００７０】
上記式（５２）によって、逆流等がある場合でも流量の近似値ＱNを決定することが可能となる。具体的には、ゼロ点がずらされた重み関数Ｖ’z0（ｘ）Ｌ（ｘ）を用いてガウス求積法の分点ｘi，係数ｗiを決定し、決定された計測点ｘiについてＶz（ｘi）を、上述のとおり超音波計測等により決定するとともに、Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1を、それぞれ定義式に従い計算により決定して上記式（５２）に代入することで近似値ＱNを決定できる。
【００７１】
ここで、既に述べたとおり、重み関数の符号が積分区間内で変化する場合であっても、積分区間内に分点ｘiをとることができる場合もある。一例として、推定流速分布が
【数４３】
2015078962.tif 000059
…（５３）
で与えられる場合は、−１≦ｘ，ｙ≦１の半径１の円管流路において、Ｖz0（ｘ）の符号が［−１，１］区間で変化するが、重み関数をＶz0（ｘ）Ｌ（ｘ）として上記式（２３）〜（２７）のガウス求積法を適用すると（Ｎ＝３とする。）、分点ｘ1，ｘ2，ｘ3＝０．００，０．７５，−０．７５、係数ｗ1，ｗ2，ｗ3＝０．０５，−０．３６，−０．３６となり、積分区間内に分点ｘiをとることができる。なお、このときｘ1〜ｘ3に対応するＶz0（ｘi）＝０．１９，−０．８１，−０．８１であり、重み関数の符号は確かに変化している。またオフセット関数Ｖz1（ｘ）を導入する場合も、例えば推定流速分布を
【数４４】
2015078962.tif 000060
…（５４）
とし、オフセット関数Ｖz1（ｘ）＝１とすれば、このＶz1（ｘ）及び上記式（５４）のｖz0（ｘ，ｙ）を用いて得られるＶ’z0（ｘ）は上記式（５３）を用いて上記式（１３）に従い得られる（オフセットを導入しない時の）Ｖz0（ｘ）と等しくなるため、Ｖ’z0（ｘ）Ｌ（ｘ）を重み関数としてガウス求積法を実行したときには積分区間内に分点ｘi（ｘ1，ｘ2，ｘ3＝０．００，０．７５，−０．７５）が得られる。したがって、実施例１，２の両方において、重み関数の符号が積分区間内で変化する場合であっても、理論的には、ガウス求積法により流量の近似値を得ることができる場合がある。
【００７２】
なお、上記式（５２）を用いて得られた近似値ＱNについても、当該ＱNよりも信頼性の低い近似値ＱSUBを用いて上記式（３６）により誤差範囲を推定することが可能である。
【００７３】
本発明による流量計測装置
本実施例２に示す、本発明の流量計測方法を実施することができる装置の一例を、図８に示す。図８の装置においては、図５の装置中のＶz0（ｘi）計算回路に代わってＶ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路が用いられる。以下、図８の装置の構成、動作を説明するが、図５の装置と同様の部分については説明を省略する。
【００７４】
図８の流量計測装置は、ガウス求積演算回路、流路配管１に設置された少なくとも１組の超音波素子２ａ，２ｂのペア、超音波送受信回路、超音波送受信信号記録回路、流速計算回路、Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路、ＱN決定回路、及び任意で誤差範囲推定回路を備える。これらの要素のうち、少なくとも１組の超音波素子２ａ，２ｂのペア、超音波送受信回路、超音波送受信信号記録回路、流速計算回路によってＶz（ｘi）決定部が構成される。
【００７５】
ガウス求積演算回路には、推定流速分布として、流路断面内の各位置（ｘ，ｙ）における流速の、法線方向成分の推定値ｖz0（ｘ，ｙ）が入力されるとともに、好ましくは（ただし必須ではない）上記式（１６）で表されるＶ’z0（ｘ）が流路断面Ｓ内で常にゼロ以上となるよう、又は流路断面Ｓ内で常にゼロ以下となるよう選定されたオフセット関数Ｖz1（ｘ）が入力される。ガウス求積演算回路は、入力されたｖz0（ｘ，ｙ），Ｖz1（ｘ）から重み関数ｗ（ｘ）＝Ｖ’z0（ｘ）Ｌ（ｘ）を構成し、上述のガウス求積法により分点ｘi及び係数ｗiを決定する。
【００７６】
超音波素子２ａ，２ｂは、図５の装置に関して説明したものと同様の素子であってよい。これらを配置すべき計測点ｘiは、修正された重み関数ｗ（ｘ）＝Ｖ’z0（ｘ）Ｌ（ｘ）を用いて決定される。超音波送受信回路、超音波送受信信号記録回路、流速計算回路も、それぞれ図５の装置に関して説明したものと同様の回路であってよい。
【００７７】
Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路は、ガウス求積演算部により決定された各々の分点ｘiについて、Ｖ’z0（ｘi），Ｖz1（ｘi）を計算するとともに、Ｑz0，Ｑz1を計算により決定する回路である。分点ｘiはガウス求積演算部から入力される。またｖz0（ｘ，ｙ）は、あらかじめ入力されていても、ガウス求積演算部から入力されてもよい。
【００７８】
ＱN決定回路は、ガウス求積演算回路により決定された係数ｗi、流速計算回路により決定されたＶz（ｘi）、及びＶ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路によって決定されたＶ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1の、各々の演算回路、計算回路からの入力を受け、これら用いて、式（２１）により流量の近似値ＱNを計算により決定する回路である。実施例１と同様に、各回路により行われる演算は典型的に計算機等を用いて行われる。
【００７９】
誤差範囲推定回路は、ＱN決定回路からＱNの入力を受けるとともに、ＱNよりも信頼性の低い近似値ＱSUBの入力をＱN決定回路又は任意の外部回路等から受けて、上記式（３６）により誤差範囲を推定する回路である。
【００８０】
実施例２における、流量計測装置の動作
次に、流速分布がｖz（ｘ，ｙ）で表わされる場合に、Ｎ＝３として本発明の方法により流量の近似値Ｑ3を得るための、流量計測装置の動作を説明する。ただし、実施例１で説明した動作と同様の部分については説明を省略する。
【００８１】
本実施例２においては、流路断面として図１等に示される流路断面Ｓ1を用い、流速の法線方向成分の推定値がｖz0（ｘ，ｙ）であるとする。この推定値ｖz0（ｘ，ｙ）がガウス求積演算回路に入力される。さらに、Ｖz1（ｘ）として、好ましくは（ただし必須ではない）流路断面Ｓ1内でＶ’z0（ｘ）＝Ｖz0（ｘ）＋Ｖz1（ｘ）が常にゼロ以上となるよう、又は常にゼロ以下となるよう選定されたＶz1（ｘ）も、ガウス求積演算回路に入力される。
【００８２】
ガウス求積演算回路は、入力されたｖz0（ｘ，ｙ），Ｖz1（ｘ）から重み関数ｗ（ｘ）＝Ｖ’z0（ｘ）Ｌ（ｘ）を構成し、上述のガウス求積法により分点ｘi及び係数ｗiを決定する。本具体例においてはＮ＝３であるため、ｘ1〜ｘ3及びｗ1〜ｗ3が決定される。
【００８３】
上記決定された各々の分点ｘ1〜ｘ3に対応して、超音波の送受信点Ａ1＊〜Ａ3＊，Ｂ1＊〜Ｂ3＊が決定される（図７。ただしｘ1〜ｘ3の実際の値は実施例１の場合と異なる。）。すなわち、各々の分点ｘ1〜ｘ3について図２に示す要領で点Ａ1〜Ａ3，Ｂ1〜Ｂ3を決定し（図２のｘiを分点ｘ1として得られる、図２中のＡ，ＢをＡ1，Ｂ1とし、図２のｘiを分点ｘ2として得られる、図２中のＡ，ＢをＡ2，Ｂ2とし、図２のｘiを分点ｘ3として得られる、図２中のＡ，ＢをＡ3，Ｂ3とする。）、これらの点をｚ方向に所定の大きさだけシフトさせることにより（図３参照。点Ａiと点Ｂiを逆方向にシフトさせる。）、点Ａ1＊〜Ａ3＊，Ｂ1＊〜Ｂ3＊が決定される。Ａ1＊，Ｂ1＊間、Ａ2＊，Ｂ2＊間、及びＡ3＊，Ｂ3＊間のそれぞれで超音波の送受信を行うことができる位置に、３組の超音波素子２ａ，２ｂのペアが配置される。
【００８４】
次に、超音波送受信回路が、上記３組の超音波素子２ａ，２ｂのペアに対して、超音波の送受信を指示する信号を送信する。３組の超音波素子２ａ，２ｂは、それぞれＡ1＊，Ｂ1＊間、Ａ2＊，Ｂ2＊間、及びＡ3＊，Ｂ3＊間で双方向に超音波の送受信を行い、各超音波素子が受信した超音波信号を超音波送受信回路へと送信する。超音波送受信回路は、受信した超音波信号から、各送受信経路における、それぞれの方向についての伝播時間ｔ1，ｔ2を計測し、超音波の送受信を示す信号として超音波送受信信号記録回路に送信する。
【００８５】
超音波送受信信号記録回路は、上記各送受信経路における、それぞれの方向についての伝播時間ｔ1，ｔ2を記録する。流速計算回路は、超音波送受信信号記録回路からこれら伝播時間ｔ1，ｔ2を受信して、上記式（４）に従い、分点ｘ1〜ｘ3に対するＶz（ｘ1）〜Ｖz（ｘ3）を計算する。
【００８６】
また、Ｖ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路は、分点ｘ1〜ｘ3について、Ｖ’z0（ｘi），Ｖz1（ｘi）を計算するとともに、Ｑz0，Ｑz1を計算により決定する。ＱN決定回路は、ガウス求積演算回路により決定された係数ｗ1〜ｗ3、流速計算回路により決定されたＶz（ｘ1）〜Ｖz（ｘ3）、及びＶ’z0（ｘi），Ｖz1（ｘi），Ｑz0，Ｑz1計算回路によって決定されたＶ’z0（ｘ1）〜Ｖ’z0（ｘ3），Ｖz1（ｘ1）〜Ｖz1（ｘ3），Ｑz0，Ｑz1を用いて、式（２１）により流量の近似値Ｑ3を計算により決定する。
【００８７】
誤差範囲推定回路は、ＱN決定回路からＱ3の入力を受けるとともに、ＱNよりも信頼性の低い近似値ＱSUBの入力をＱN決定回路又は任意の外部回路等から受けて、上記式（３６）により誤差範囲を推定する。近似値ＱSUBは、上記各回路を用いて３以外のＮについて決定された流量の近似値であってもよいし、本発明とは異なる任意の方法により決定された流量の近似値であってもよい。
【００８８】
本発明による流量計測方法を用いた実験
実用的な流体設備では、流路中の種々の構成機器、構造物によって流速分布が乱される場合があり、流量計測に誤差を生じる要因となる。そのような整流されていない流路における流量計測精度について、本発明の流量計測方法を用いた計測結果と従来技術による計測結果とを比較することにより、本発明の有効性を例示する。
【００８９】
（実施方法）
整流されていない流路の典型例として、曲管（肘部）及び邪魔板（曲管下流の流路断面を部分的に封止）を配置した円管流路を考える（図９Ａ，図９Ｂ）。邪魔板は図９Ａ中の位置Ｘに配置されており、その断面は、図９Ｂに示すとおり半円形状である。なお、流路断面は直径約８３ｍｍの円形である。
【００９０】
図９Ａの円管流路には、５（リットル／秒）の一定流量で試験流体（水）が流入し続けている。実施例２に従い、Ａ断面（曲管及び邪魔板の直下。位置Ｘから軸方向に約２００ｍｍの位置。）において流量を計測する。ただし、推定流速分布は事前のシミュレーションによって求める。邪魔板下流に逆流が発生するため、仮想的に流速分布のゼロ点をずらしてガウス求積法を構成する必要がある。本件においては、Ｖz1（ｘ）を定数関数（推定流速の流路断面平均の１倍）と取ることにより、上記式（１６）〜（１８）に基づき重み関数を構成し、実施例２の方法に従い流量を計測した。
【００９１】
また比較対象として、特許文献３，４の従来技術（ガウス‐チェビチェフ求積法）に従い、Ａ断面において流量を計測した。さらに、流量の基準値として、汎用の流量計測装置を用いてＢ断面（曲管及び邪魔板の十分下流。位置Ｘから軸方向に３０００ｍｍ以上の位置。）においても流量を計測した。なお、本発明の実施例２の方法により流量を計測したときの、ガウス求積法の分点ｘi及び係数ｗi、Ｖz0（ｘi），Ｑz0，Ｖz1（ｘi），Ｑz1の値はそれぞれ以下のとおりである。ただし流路半径を１とし、推定流速は、流路断面平均が１となるよう正規化している。なお、ガウス−チェビチェフ求積法（Ｎ＝３）により流量を計測したときの分点ｘiはｘ1，ｘ2，ｘ3＝０．００，０．７１，−０．７１であり、係数ｗｉはｗ１，ｗ２，ｗ３＝０．７９，０．３９，０．３９であった。
【表３】
2015078962.tif 000061
2015078962.tif 000062
【００９２】
上記Ｂ断面で計測した基準流量値と比較することにより、Ａ断面において実施例２の方法と従来技術とでそれぞれ計測した流量の誤差を評価した。
【００９３】
図１０に、Ａ断面での流量計測方法として従来のガウス‐チェビチェフ求積法を用いた場合と実施例２の方法を用いた場合とのそれぞれにおける、「Ｂ断面において計測した流量」と「Ａ断面において計測した流量」との差の絶対値を示す。本発明により、従来技術と比較して、流量計測誤差を１／１０程度に低減することができることが示されている。
【００９４】
なお、既に述べたとおり、流路断面が傾いた面や曲面であったとしても、本発明によって同様の原理で流量を計測することが可能である。例えば、図１１に示すとおり、流路の軸に対して角度φだけ傾いた流路断面Ｓ2（流路断面Ｓ1を、図１のｘ軸の廻りにφ回転させることで得られる。）を用いる場合には、図１２に示すとおり、この流路断面Ｓ2上でｘ座標、ｙ座標を規定し、既に説明した方法と同様の方法で流量を計測することができる。このとき、図１１に示すとおり流路配管１の軸方向に流体が流れているならば、流速のうち流路断面Ｓ2の法線方向（ｚ方向）成分は、流速の大きさをｖ（ｘ，ｙ）とすればｖz（ｘ，ｙ）＝ｖ（ｘ，ｙ）・ｃｏｓφとなる（図１３）。また、超音波の送受信を行う位置Ａ＊，Ｂ＊は、図１２に示すＡ，Ｂを、それぞれ流路配管１の軸方向逆向きに所定の大きさだけシフトした位置とすればよい。この場合、上述の伝播時間差方式においては、上記式（４）を用いて得られた平均流速にｃｏｓφを乗じればＶz（ｘi）を得られる。流路断面が曲面である場合にも、流路断面内の各位置（ｘ，ｙ）に応じた法線方向の成分を用いる（各位置（ｘ，ｙ）で異なるφを用いる）ことにより、同様にｖz（ｘ，ｙ）を決定して流量計測を行うことができる。あるいは、最小二乗法等によって当該曲面を最もよく近似する平面を決定し、当該決定された平面を仮想的な流路断面として用いる等して、（４）式を用いて（傾斜角の代表値をφとして、ｃｏｓφを乗じる。）平均流速を計測し、流量を決定することも可能である。
【産業上の利用可能性】
【００９５】
本発明の計測方法、計測装置を、高精度の流量計測一般に利用できる可能性がある。例えば、実際のプラント等で利用する場合、従来とは異なり整流されていない流速分布下で正確な計測が可能となるため、整流化するために延長していた配管等の圧縮、流量計取り付け位置の自由度の向上等が期待できる。
【符号の説明】
【００９６】
１                流路配管
２ａ，２ｂ        超音波素子
【図１】
2015078962.tif 000095
【図２】
2015078962.tif 000096
【図３】
2015078962.tif 000097
【図４】
2015078962.tif 000098
【図５】
2015078962.tif 000099
【図６】
2015078962.tif 000100
【図７】
2015078962.tif 000101
【図８】
2015078962.tif 000102
【図９Ａ】
2015078962.tif 000103
【図９Ｂ】
2015078962.tif 000104
【図１０】
2015078962.tif 000105
【図１１】
2015078962.tif 000106
【図１２】
2015078962.tif 000107
【図１３】
2015078962.tif 000108