(57)【要約】      （修正有）【課題】流速の評価に要する時間を短縮化する。【解決手段】流速計測システム１００は、加熱部（熱交換器）２と、プレヒーター２ｂと、温度計測部３と、制御ユニット（流速算出部）４と、流速計５と、を含む。配管１０は、蒸気製造装置２０（ボイラーなど）と負荷設備３０との間に配設されている。加熱部２は、例えば、リング状のヒーターから構成されており、配管１０の表面の所定位置において周方向に亘って配置し、当該配管１０の表面を均一に加熱する。加熱部２は、制御ユニット４に電気的に接続されており、その動作が制御される。温度計測部３は、複数の温度センサー群から構成され、各温度センサー群は配管１０の表面において、配管１０の管軸方向に沿って配置される。ＣＰＵは、温度計測デー（温度分布）およびメモリに記憶された情報に基づき、配管１０の内部を流れる蒸気の流速を算出する。【選択図】図１
2017072424.tif 000002
【特許請求の範囲】
【請求項１】
配管の内部を流れる流体の流速を評価する評価方法であって、
１点以上の流速について温度分布の特徴量の解析値に基づいて、流速を決定する、
流速の評価方法。
【請求項２】
温度分布の特徴量は、流速の増加に応じて減少する特徴量である、
請求項１に記載の流速の評価方法。
【請求項３】
温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である、
請求項１または請求項２のいずれか１項に記載の流速の評価方法。
【請求項４】
前記配管の表面の所定部分で熱交換を行う熱交換工程と、
前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算工程と、
前記計算工程により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定工程と、
前記流速設定工程により設定された流速を用いて、流速を決定する流速決定工程と、
を含む、請求項１から請求項３のいずれか１項に記載の流速の評価方法。
【請求項５】
前記流速設定工程は、想定された流速の範囲の中点の流速について前記計算工程により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する、
請求項４に記載の流速の評価方法。
【請求項６】
前記流速設定工程は、３点以上の流速を設定し、
前記流速決定工程は、前記流速設定工程により設定された前記３点以上の流速を用いて、流速を決定する、
請求項４または請求項５のいずれか１項に記載の流速の評価方法。
【請求項７】
前記流速決定工程は、前記流速設定工程で使用された流速の一部または全部を用いる、
請求項４から請求項６のいずれか１項に記載の流速の評価方法。
【請求項８】
前記流速決定工程は、
前記流速設定工程により設定された流速のそれぞれにおいて、温度分布の解析値と実測値との２乗誤差を計算する第１−１工程と、
前記流速と前記２乗誤差との関係を求める第１−２工程と、
流速の候補を決定する第１−３工程と、
決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程と、
前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程と、
前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程と、
を含む、
請求項４から請求項７のいずれか１項に記載の流速の評価方法。
【請求項９】
前記流速決定工程は、
前記流速設定工程により設定された流速のそれぞれにおいて、流速と温度分布の前記特徴量との関係の解析値を計算する第２−１工程と、
計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程と、
決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程と、
前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程と、
を含む、
請求項４から請求項７のいずれか１項に記載の流速の評価方法。
【請求項１０】
配管の内部を流れる流体の流速を評価する評価システムであって、
１点以上の流速について温度分布の特徴量の解析値に基づいて、流速を決定する、
流速の評価システム。
【請求項１１】
温度分布の特徴量は、流速の増加に応じて減少する特徴量である、
請求項１０に記載の流速の評価システム。
【請求項１２】
温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である、
請求項１０または請求項１１のいずれか１項に記載の流速の評価システム。
【請求項１３】
前記配管の表面の所定部分で熱交換を行う熱交換器と、
前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算部と、
前記計算部により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定部と、
前記流速設定部により設定された流速を用いて、流速を決定する流速決定部と、
を含む、請求項１０から請求項１２のいずれか１項に記載の流速の評価システム。
【請求項１４】
前記流速設定部は、想定された流速の範囲の中点の流速について前記計算部により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する、
請求項１３に記載の流速の評価システム。
【請求項１５】
前記流速設定部は、３点以上の流速を設定し、
前記流速決定部は、前記流速設定部により設定された前記３点以上の流速を用いて、流速を決定する、
請求項１３または請求項１４のいずれか１項に記載の流速の評価システム。
【請求項１６】
前記流速決定部は、前記流速設定部で使用された流速の一部または全部を用いる、
請求項１３から請求項１５のいずれか１項に記載の流速の評価システム。
【請求項１７】
前記流速決定部は、
前記流速設定部により設定された流速のそれぞれにおいて、温度分布の解析値と実測値との２乗誤差を計算する第１−１工程と、
前記流速と前記２乗誤差との関係を求める第１−２工程と、
流速の候補を決定する第１−３工程と、
決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程と、
前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程と、
前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程と、
を実行する、
請求項１３から請求項１５のいずれか１項に記載の流速の評価システム。
【請求項１８】
前記流速決定部は、
前記流速設定部により設定された流速のそれぞれにおいて、流速と温度分布の前記特徴量との関係の解析値を計算する第２−１工程と、
計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程と、
決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程と、
前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程と、
を実行する、
請求項１３から請求項１５のいずれか１項に記載の流速の評価システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、流速の評価方法および評価システムに関する。
【背景技術】
【０００２】
ヒーター法を用いて蒸気の流量を計測する方法が検討されている（特許文献１参照）。
特許文献１では、配管内を流れる流体の流速計測方法が開示されている。この流速計測方法では、配管の表面の所定部分で熱交換を行い、当該配管の管軸方向における温度分布を計測し、計測した温度分布に基づいて、配管の内部を流れる流体の流速を求める。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１５−１４８５０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
温度分布の解析では、例えば、有限要素法の処理に時間がかかる。一例として、１個の点（ポイント）について処理に２０分かかる場合、３個の点だと１時間かかり、さらに周辺の多くの点で処理を行うと総じて多大な時間がかかる。
このため、流速の評価に多大な時間を要する場合があった。
【０００５】
具体例として、ヒーター法では、温度分布の実測値に理論値をフィッティングして、流速を求める。理論値の温度分布は、管内の流速を入力値として、数値計算を用いて、解く必要があるため、１回の計算にもかなりの時間を要する。そして、フィッティングのためには、多くの流速に対応した温度分布を計算する必要があるため、短時間での解析が困難であった。
【０００６】
本発明は、流速の評価に要する時間を短縮化することが可能な流速の評価方法および評価システムを提供する。
【課題を解決するための手段】
【０００７】
本発明の一態様は、配管の内部を流れる流体の流速を評価する評価方法であって、１点以上の流速について温度分布の特徴量の解析値に基づいて、流速を決定する、流速の評価方法である。
【０００８】
本発明の一態様は、温度分布の特徴量は、流速の増加に応じて減少する特徴量である、流速の評価方法である。
【０００９】
本発明の一態様は、温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である、流速の評価方法である。
【００１０】
本発明の一態様は、前記配管の表面の所定部分で熱交換を行う熱交換工程と、前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算工程と、前記計算工程により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定工程と、前記流速設定工程により設定された流速を用いて、流速を決定する流速決定工程と、を含む、流速の評価方法である。
【００１１】
本発明の一態様は、前記流速設定工程は、想定された流速の範囲の中点の流速について前記計算工程により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する、流速の評価方法である。
【００１２】
本発明の一態様は、前記流速設定工程は、３点以上の流速を設定し、前記流速決定工程は、前記流速設定工程により設定された前記３点以上の流速を用いて、流速を決定する、流速の評価方法である。
【００１３】
本発明の一態様は、前記流速決定工程は、前記流速設定工程で使用された流速の一部または全部を用いる、流速の評価方法である。
【００１４】
本発明の一態様は、前記流速決定工程は、前記流速設定工程により設定された流速のそれぞれにおいて、温度分布の解析値と実測値との２乗誤差を計算する第１−１工程と、前記流速と前記２乗誤差との関係を求める第１−２工程と、流速の候補を決定する第１−３工程と、決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程と、前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程と、を含む、流速の評価方法である。
【００１５】
本発明の一態様は、前記流速決定工程は、前記流速設定工程により設定された流速のそれぞれにおいて、流速と温度分布の前記特徴量との関係の解析値を計算する第２−１工程と、計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程と、決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程と、を含む、流速の評価方法である。
【００１６】
本発明の一態様は、配管の内部を流れる流体の流速を評価する評価システムであって、１点以上の流速について温度分布の特徴量の解析値に基づいて、流速を決定する、流速の評価システムである。
【００１７】
本発明の一態様は、温度分布の特徴量は、流速の増加に応じて減少する特徴量である、流速の評価システムである。
【００１８】
本発明の一態様は、温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である、流速の評価システムである。
【００１９】
本発明の一態様は、前記配管の表面の所定部分で熱交換を行う熱交換器と、前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算部と、前記計算部により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定部と、前記流速設定部により設定された流速を用いて、流速を決定する流速決定部と、を含む、流速の評価システムである。
【００２０】
本発明の一態様は、前記流速設定部は、想定された流速の範囲の中点の流速について前記計算部により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する、流速の評価システムである。
【００２１】
本発明の一態様は、前記流速設定部は、３点以上の流速を設定し、前記流速決定部は、前記流速設定部により設定された前記３点以上の流速を用いて、流速を決定する、流速の評価システムである。
【００２２】
本発明の一態様は、前記流速決定部は、前記流速設定部で使用された流速の一部または全部を用いる、流速の評価システムである。
【００２３】
本発明の一態様は、前記流速決定部は、前記流速設定部により設定された流速のそれぞれにおいて、温度分布の解析値と実測値との２乗誤差を計算する第１−１工程と、前記流速と前記２乗誤差との関係を求める第１−２工程と、流速の候補を決定する第１−３工程と、決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程と、前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程と、を実行する、流速の評価システムである。
【００２４】
本発明の一態様は、前記流速決定部は、前記流速設定部により設定された流速のそれぞれにおいて、流速と温度分布の前記特徴量との関係の解析値を計算する第２−１工程と、
計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程と、決定された流速の候補において、温度分布の解析値と実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程と、を実行する、流速の評価システムである。
【発明の効果】
【００２５】
上記した流速の評価方法および評価システムによれば、流速の評価に要する時間を短縮化することが可能である。
【図面の簡単な説明】
【００２６】
【図１】本実施形態に係る流速計測システムの概略構成を示す図である。
【図２】流速計測システムの要部構成を示す図である。
【図３】制御ユニット４を示す模式図である。
【図４】配管１０の領域のセルへの分割例を示す説明図である。
【図５】配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速５ｍ／ｓのとき）の例を示す図である。
【図６】配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速２０ｍ／ｓのとき）の例を示す図である。
【図７】配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速３５ｍ／ｓのとき）の例を示す図である。
【図８】温度分布の特性における半値幅Ｐ１の例を示す図である。
【図９】温度分布の特性における最高温度Ｐ２の例を示す図である。
【図１０】温度分布の特性における温度分布波形の面積Ｐ３の例を示す図である。
【図１１】流速と半値幅との関係の例を示す図である。
【図１２】温度分布の実測値４０１と温度分布の解析値４１１との関係の例を示す図である。
【図１３】流速と２乗誤差（温度分布の実測値と温度分布の解析値との２乗誤差）との関係の例を示す図である。
【図１４】第１の処理群の一例を示す図である。
【図１５】第１の処理群における流速設定処理の一例を示す図である。
【図１６】第２の処理群の第１の例を示す図である。
【図１７】流速と２乗誤差（温度分布の実測値と温度分布の解析値との２乗誤差）との関係の例を示す図である。
【図１８】第２の処理群の第２の例を示す図である。
【発明を実施するための形態】
【００２７】
以下、本発明の実施形態を説明するが、以下の実施形態は特許請求の範囲にかかる発明を限定するものではない。また、実施形態の中で説明されている特徴の組み合わせのすべてが発明の解決手段に必須であるとは限らない。
本実施形態に係る流速計測システム（流速の評価システムの一例）は、例えば、ボイラーなどの蒸気製造装置と負荷設備との間に配設される配管内を流れる流体（例えば、蒸気）の流速を計測可能なシステムである。
【００２８】
図１は、本実施形態に係る流速計測システム１００の概略構成を示す図である。図２は、流速計測システム１００の要部構成を示す図である。
本実施形態に係る流速計測システム１００は、図１に示すように、加熱部（熱交換器）２と、プレヒーター２ｂと、温度計測部３と、制御ユニット（流速算出部）４と、流速計５と、を含む。図１において、配管１０は、蒸気製造装置２０（ボイラーなど）と負荷設備３０との間に配設されている。蒸気製造装置２０からの蒸気が配管１０を流れ、負荷設備３０に送られる。負荷設備３０において、蒸気または蒸気の熱が利用される。負荷設備３０から排出された蒸気はドレンとして回収され、還水槽（不図示）に集約された後、蒸気製造装置２０に再度給水される。また、配管１０の周りには断熱材が巻かれている。
【００２９】
ここで、本実施形態では、流速計測システム１００に流速計５を備える構成を示すが、他の構成例として、流速計５が備えられない構成が用いられてもよい。つまり、流速計５を用いて流速を計測する構成が用いられてもよく、または、他の情報から流速を計算する（推定的な計算でもよい）構成が用いられてもよい。
【００３０】
加熱部２は、配管１０の表面１０ａと熱交換することで所定部分を加熱するためのものである。本実施形態において、加熱部２は、例えば、リング状のヒーターから構成されており、図２に示すように、配管１０の表面１０ａの所定位置において周方向に亘って配置される。これにより、加熱部２は、配管１０の所定部分（当該加熱部２の設置部分１１）において、当該配管１０の表面１０ａを均一に加熱する。加熱部２は、制御ユニット４に電気的に接続されており、その動作が制御される。
なお、ここでいう熱交換器は、温度の高い物体から温度の低い物体へ熱を移動させるものである。ヒーターは、熱交換器の一例である。
なお、図２において、矢印Ｂ１１は、流体が配管１０内を流れる向きの例を示している。
【００３１】
プレヒーター２ｂは、配管１０内を流れる流体を加熱する。このプレヒーター２ｂは、配管１０内を流れる気体の流体の液化を防止するために設けられている。例えば、配管１０内を流れる流体が蒸気である場合に、プレヒーター２ｂが当該蒸気を加熱することで、当該蒸気が液化して配管１０内や負荷設備３０内に水滴が付着することを防止する。
【００３２】
温度センサー群３Ａは、配管１０における表面１０ａの加熱部２の設置部分１１の両側（上流側および下流側）に配置されている。各温度センサー群３Ａは、上記設置部分１１からの距離に応じて設置位置が決定される。例えば、設置部分１１の上流側を例に挙げると、各温度センサー群３Ａは、設置部分１１の端面からの距離が０ｍｍ（ミリメートル）、６ｍｍ、１４ｍｍ、２４ｍｍ、３６ｍｍ、５０ｍｍ、６６ｍｍ、８４ｍｍ、１０４ｍｍ、１２６ｍｍ、１５０ｍｍ、１７６ｍｍに設置されている。ここで、設置部分１１の端面からの距離が０ｍｍとは、温度センサー群３Ａが加熱部２の端面に沿って配置されることを意味する。なお、図２では、配管１０の断面構造として、設置部分１１の下流側端面の近傍（Ａ−Ａ矢視による断面）と、設置部分１１の下流側端面からの距離２４ｍｍ近傍（Ｂ−Ｂ矢視による断面）とを図示した。
【００３３】
図２に示されるように、温度センサー群３Ａは、設置部分１１から離間するに従って、隣接する温度センサー群３Ａ間の距離が２ｍｍずつ大きくなるように配置されている。したがって、温度センサー群３Ａは、設置部分１１（加熱部２）に近い程、センサーが密集して配置されたものとなっている。これにより、設置部分１１の近傍において配管１０の表面１０ａの温度を精度良く検出することが可能とされている。
【００３４】
温度計測部３は、複数（本実施形態では、例えば、１２個）の温度センサー群３Ａから構成される。各温度センサー群３Ａは、配管１０の表面１０ａにおいて、当該配管１０の管軸方向に沿って配置される。各温度センサー群３Ａは、それぞれ配管１０の表面１０ａの温度を計測する温度センサー３ａを複数含む。本実施形態において、各温度センサー群３Ａは４つの温度センサー３ａから構成される。４つの温度センサー３ａは、配管１０の表面１０ａにおいて、周方向に均等に配置されている。すなわち、４つの温度センサー３ａは、配管１０の周方向において、９０度ずつ位置を違えるように配置されている。各温度センサー群３Ａは、４つの温度センサー３ａが計測した値の平均値を計測値として出力する。このように温度センサー群３Ａは、配管１０の表面１０ａにおける複数個所を計測した値の平均を計測値とすることで信頼性の高い計測結果（温度）を出力することが可能である。
【００３５】
なお、各温度センサー群３Ａを構成する温度センサー３ａの数は、図２に示す４つに限られず１つ以上であればよい。例えば、各温度センサー群３Ａが、２つの温度センサー３ａから構成され、２つの温度センサー３ａが、配管１０の周方向に左右均等に（配管１０の周方向において１８０度位置を違えるように）配置されていてもよい。上記のように、各温度センサー群３Ａが複数の温度センサー３ａを備え、配管１０の表面１０ａにおける複数個所を計測した値の平均を計測値とすることで信頼性の高い計測結果（温度）を出力することが可能である。
【００３６】
上述した構成に基づき、温度計測部３は、各温度センサー群３Ａの計測結果から配管１０の管軸方向における表面１０ａの温度分布を計測することが可能である。温度計測部３が計測した温度分布は、制御ユニット４に送信される。
【００３７】
配管１０は、表面１０ａの少なくとも一部が保温材１２により覆われている。本実施形態において、保温材１２は、配管１０の表面に設けられた加熱部２および温度計測部３（各温度センサー３ａ）を覆うように管軸方向に亘って設置されている。
【００３８】
図３は、制御ユニット４を示す模式図である。図３において、計算装置５０は、例えばコンピュータシステムである。制御ユニット４は、計算装置５０に加え、入力装置６０、および表示装置（出力装置）６４を有する。計算装置５０は、Ａ／Ｄ変換器等の変換器６１、演算処理手段（ＣＰＵ：Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）６２、およびメモリ６３等を有する。流速計測システム１００の温度計測部３から送られる計測データ（温度分布）が、必要に応じて変換器６１等で変換され、ＣＰＵ６２に取り込まれる。また、初期設定値、および仮データなどが入力装置６０などを介して計算装置５０に取り込まれる。表示装置６４は、入力されたデータに関する情報、および計算に関する情報などを表示することができる。
【００３９】
ＣＰＵ６２は、計測データ、およびメモリ６３に記憶された情報に基づき、配管１０の内部を流れる蒸気の流速を算出することができる。ＣＰＵ６２は、例えば、温度計測部３の計測結果（配管１０の表面１０ａにおける温度分布）を用い、メモリ６３に記憶された情報から配管１０の内部を流れる蒸気の流速を算出する。制御ユニット４は、流速算出部の一例である。
具体的には、メモリ６３は、配管１０の表面１０ａにおける温度分布と、配管１０の内部を流れる流体の流速との関係を示す関係情報を予め（流速の計測よりも前に）記憶しておく。メモリ６３は、記憶部の一例である。メモリ６３が記憶する関係情報では、流体の流速ごとに、当該流速と温度分布とが対応付けられている。
そして、温度計測部３が配管１０の表面１０ａにおける温度分布を計測すると、ＣＰＵ６２は、関係情報を参照して、温度計測部３が計測した温度分布に最も近い（例えば、差の絶対値の合計が最も小さい）温度分布に対応付けられている流速を読み出して、流速の計測値とする。
【００４０】
このように、メモリ６３は、温度分布と流速とを対応付けた関係情報を予め記憶しておく。そして、加熱部２が、リングヒーター（リング状のヒーター）で所定箇所を加熱し、配管１０の温度分布が定常状態になった状態で、温度計測部３が、配管１０における温度分布を計測する。そして、制御ユニット４は、関係情報を参照して、計測で得られた温度分布に最も近い温度分布に対応付けられている流速を流速計測値とする。
【００４１】
次に、メモリ６３が予め記憶しておく関係情報の取得方法について説明する。
関係情報は、配管１０の温度定常状態における熱の伝わりを、有限要素法で解析して求めることができる。以下では、制御ユニット４が有限要素法の計算を行う場合を例に説明するが、他のコンピュータを用いて有限要素法の解析を行うようにしてもよい。
配管１０に有限要素法を適用するために、配管１０の領域をセル（部分領域）に分割する。
【００４２】
図４は、配管１０の領域のセルへの分割例を示す説明図である。
同図に示す領域Ａ１１は、管内の領域（流体が流れる領域）を示す。領域Ａ１２は、配管１０の領域（管壁の領域）を示す。領域Ａ１３は、配管１０の周りに巻かれた断熱材の領域を示す。領域Ａ１４は、断熱材の外側の空気の領域を示す。また、矢印Ｂ１１は、流体が配管１０内を流れる向きの例を示している。
【００４３】
有限要素法を適用するためのセル分割にて、図４に示すように、配管１０の管壁の領域Ａ１２を均等な厚みで３層に分割する。この３層と、流体の領域Ａ１１、断熱材の領域Ａ１３、および、外部の空気の領域Ａ１４とで、配管１０の半径方向に６層に分割されている。
また、配管１０の軸方向（長手方向）に関しては、例えば３ｍｍ幅など比較的小さい幅で均等に分割する。
また、加熱部２としてリング状のヒーターを用い、図４に示すように、配管１０の長手方向における加熱部２の厚みを無視する。特に、加熱部２が、配管１０の管壁の領域Ａ１２のセルのうち１つのみに入熱するものとして近似する。また、加熱部２自体の容量の影響は無視する（容量が十分に小さいものとする）。
【００４４】
隣接する接点との熱移動は、管内面と管内の流体との対流熱伝達、管壁内（図４の領域Ａ１２内）での熱伝導、管壁から断熱材への熱伝導および断熱材内での熱伝導、断熱材外表面と周囲空気との対流熱伝達による熱移動とする。また、加熱部２から十分離れた管端の部分のセルの外縁を断熱条件とする。
以下の熱バランスによる方程式（以下の式（１））をセルごとに設定しておき、制御ユニット４が、有限要素法を用いて解析することで、セル間の温度差を算出する。有限要素法の解法（連立方程式の解法）として、例えばNewton-Raphson法を用いることができる。但し、制御ユニット４が用いる解法は、Newton-Raphson法に限られず、連立方程式に適用することが可能ないろいろな解法を用いることができる。
【００４５】
ここで、定常状態では、隣接するセルからの入熱量の和が０になる。なお、熱量の放出は、入熱量マイナスとして表す。
座標（ｉ，ｊ）に位置するセルにおける熱バランスは、式（１）のように表される。
【００４６】
Ｑｉ，ｊ−１  ＋  Ｑｉ，ｊ＋１  ＋  Ｑｉ−１，ｊ  ＋  Ｑｉ＋１，ｊ  ＝  ０
・・（１）
【００４７】
ここで、座標（ｉ，ｊ）に位置するセルに隣接するセルの座標を、（ｉ，ｊ−１）、（ｉ，ｊ＋１）、（ｉ−１，ｊ）、（ｉ＋１，ｊ）とする。また、Ｑｉ，ｊ−１、Ｑｉ，ｊ＋１、Ｑｉ−１，ｊ、Ｑｉ＋１，ｊは、それぞれ、添え字で示す座標に位置するセルからの入熱量を示す。なお、隣接するセルが無い場合は、当該セルからの入熱量を０とする。
上記のように、セルごとに熱バランスによる方程式（式（１））を設定する。当該方程式の設定は、例えば流速計測システム１００のユーザーが行って、制御ユニット４のメモリ６３に記憶させる。そして、制御ユニット４は、流速の設定値ごとに有限要素法による解析を行って、温度分布と流速との関係を示す関係情報を取得する。
【００４８】
式（１）のＱｉ，ｊ−１、Ｑｉ，ｊ＋１、Ｑｉ−１，ｊ、Ｑｉ＋１，ｊには、例えば以下の式（２）、式（３）、式（４）のいずれかの右辺を適用する。
【００４９】
対流熱伝達：配管１０のセルのうち最内側のセルへの流体（配管１０内の流体）からの入熱量は、式（２）のように示される。
【００５０】
対流熱伝達による入熱量  ＝  α・Ａ・Δｔ    ・・（２）
【００５１】
ここで、α［Ｗ（ワット）／（ｍ２（平方メートル）・Ｋ（ケルビン））］は、対流熱伝達率を示す。
Ａ［ｍ２］は、伝熱面積を示す。ここでは、配管１０のセルが管内の流体に接する面積である。
Δｔ［Ｋ］は、隣接するセルとの温度差を示す。ここでは、配管１０と管内の流体との接触部分における温度差を示す。
【００５２】
半径方向熱伝導：配管１０内における配管１０の半径方向の熱伝導による入熱量は、式（３）のように示される。すなわち、管の半径方向に隣接する管のセル（管をメッシュに切ったセル）からの入熱量は、式（３）のように示される。
【００５３】
半径方向熱伝導による入熱量  ＝  ２π・λ・Ｌ・Δｔ／ｌｎ（ｒｏ／ｒｉ）
・・（３）
【００５４】
ここで、πは、円周率を示す。
λ［Ｗ／（ｍ（メートル）・Ｋ）］は、配管１０の素材（例えば鋼鉄）の熱伝導率を示す。
Ｌ［ｍ］は、管軸方向（配管１０の軸方向）におけるセルの長さを示す。
Δｔ［Ｋ］は、上記のように、隣接するセルとの温度差を示す。ここでは、配管１０の半径方向に隣接する配管１０のセル同士の接触部分における温度差を示す。
ｌｎは、自然対数を示す。
ｒｏ／ｒｉ［ｍ］は、半径方向におけるセル間の距離（例えば、セルの中心間の距離）を示す。
【００５５】
管軸方向熱伝導：配管１０内における配管１０の軸方向（長手方向）の熱伝導による入熱量は式（４）のように示される。すなわち、管軸方向に隣接する配管１０のセルからの入熱量は式（４）のように示される。
【００５６】
管軸方向熱伝達による入熱量  ＝  λＡΔｔ／Ｌ    ・・（４）
【００５７】
ここで、λ、Ｌは、上記のとおりである。
Ａ［ｍ２］は、上記のように、伝熱面積を示す。ここでは、配管１０のセル同士（配管１０の半径方向に隣接する配管１０のセル）が接する面積である。
Δｔ［Ｋ］は、上記のように、隣接するセルとの温度差を示す。ここでは、管軸方向に隣接する配管１０のセル同士の接触部分における温度差を示す。
【００５８】
また、対流熱伝達に関して、以下の式（５）を用いる。
【００５９】
Ｎｕｄ  ＝  ０．０２２・ＸＷＴ・Ｒｅｄ０．８・Ｐｒ０．４    ・・（５）
【００６０】
ここで、Ｎｕｄは、ヌセルト数（Nusselt Number）を示す。
Ｒｅｄは、レイノルズ数（Reynolds Number）を示す。
Ｐｒは、プラントル数（Prandtl Number）を示す。
ＸＷＴは、温度を修正する係数であり、加熱部２（ヒーター）からの距離に応じた値を取る。特に、ＸＷＴは、加熱部２の近傍について温度を高くする（すなわち、熱伝達を大きくする）。以下では、ＸＷＴを熱伝達率修正係数と称する。
【００６１】
なお、本実施形態では、熱伝達率修正係数ＸＷＴを使用することで精度を向上させる場合を示すが、他の構成例として、精度が十分であれば、熱伝達率修正係数ＸＷＴを使用しない構成（熱伝達率修正係数ＸＷＴ＝１とみなす構成）が用いられてもよい。
【００６２】
ヌセルト数Ｎｕｄは、式（６）のように示される。
【００６３】
Ｎｕ＝αｉ・ｄｉ／λ    ・・（６）
【００６４】
ここで、λは、上記のとおりである。
ｄｉ［ｍ］は、配管１０の内径を示す。
αｉ［Ｗ／（ｍ２・Ｋ）］は、配管１０の内部における熱伝導率を示す。
【００６５】
また、レイノルズ数Ｒｅｄは、式（７）のように示される。
【００６６】
Ｒｅ＝ｕ×ｄｉ／ν    ・・（７）
【００６７】
ここで、ｕ［ｍ／ｓ（秒）］は、配管１０内を流れる流体の流速を示す。
ν［ｍ２／ｓ］は、配管１０内を流れる流体の動粘性係数を示す。
ｄｉ［ｍ］は、配管１０の内径を示す。
【００６８】
また、プラントル数Ｐｒは、式（８）のように示される。
【００６９】
Ｐｒ＝ν×ρ×Ｃｐ／λ    ・・（８）
【００７０】
ν、λは、上記のとおりである。
ρ［ｋｇ（キログラム）／ｍ３（立方メートル）］は、流体の密度を示す。
Ｃｐ［ｋＪ（キロジュール）／（ｋｇ・Ｋ）］は、流体の比熱を示す。
【００７１】
式（５）のヌセルト数Ｎｕｄは、式（２）の対流熱伝達率αと比例し、ヌセルト数Ｎｕｄから対流熱伝達率αを求めることができる。式（５）のように熱伝達率修正係数ＸＷＴを導入してヌセルト数Ｎｕｄの値を高精度に求めることで、対流熱伝達率αの値を高精度に求めることができる。対流熱伝達率αの値を高精度に得られることで、有限要素法を用いての解析でセル間の温度差を高精度に求めることができ、これにより、配管１０における温度分布の解析値の精度を高めることができる。
【００７２】
また、対流熱伝達に関して、以下の式（９）も用いる。
【００７３】
Ｗ  ＝  ＹＨＴ・Ｇ・Ｃｐ・ΔＴ    ・・（９）
【００７４】
ここで、Ｃｐは、上記のとおりである。
Ｗ［Ｊ（ジュール）／ｓ］は、加熱部２による加熱量を示す。
Ｇ［ｋｇ／ｓ］は、配管１０を流れる流体の全流量を示す。
ΔＴ［℃（度）］は、加熱部２のヒーターの加熱による流体の上昇温度を示す。
ＹＨＴは、管断面を流れる流体の全量のうち、加熱部２のヒーターからの熱の伝達に寄与する流体の量（加熱部２のヒーターの加熱による熱を受けた流体の量）の割合を示す。具体的には、管断面における流体全体の面積をＳとし、管断面における流体のうちヒーターの加熱による熱を受けた部分の面積をＳ１として、ＹＨＴ＝Ｓ１／Ｓと表される。以下では、ＹＨＴを温度境界係数と称する。
【００７５】
なお、本実施形態では、温度境界係数ＹＨＴを使用することで精度を向上させる場合を示すが、他の構成例として、精度が十分であれば、温度境界係数ＹＨＴを使用しない構成（温度境界係数ＹＨＴ＝１とみなす構成）が用いられてもよい。
【００７６】
上述したように、配管１０内に蒸気を流して温度分布および流速を実測し、温度分布の解析値（一般的な熱伝達の式を適用して有限要素法を用いた解析で得られた値）と実測値とを比較したところ、加熱部２のヒーター設置位置の近傍で、温度実測値が解析値よりも高くなった。上述した加熱部２のヒーター設置位置から上流側での温度のずれに加えて、加熱部２のヒーター設置位置から下流側でも、温度実測値が解析値よりも高くなった。この温度のずれの一因として、配管１０内を流れる流体のうち、ヒーターからの熱を伝達するのは配管１０の内面に近い一部のみであることが考えられる。そこで、温度境界係数ＹＨＴを導入して再計算をおこなったところ、解析値と実測値とがよりよく一致した。特に、熱伝達率修正係数ＸＷＴ、温度境界係数ＹＨＴの両方を導入することで、解析値と実測値とがよりよく一致した。
【００７７】
温度境界係数ＹＨＴの値を１より小さくすることは、ヒーターからの熱の伝達に寄与する流体の量を少なくすることに相当する。ヒーターからの熱の伝達に寄与する流体の量が少ないと、式（２）の温度差Δｔの値（配管１０と配管内の流体との接触部分における温度差）が大きく算出される。温度境界係数ＹＨＴを導入して式（２）の温度差Δｔの値を高精度に算出することで、配管１０における温度分布の解析値の精度を高めることができる。
【００７８】
流速計５は、流速計測部の一例であり、配管１０内を流れる流体の流速を計測する。流速計５を備えることで、流速計測システム１００は、配管１０における温度分布と、配管１０内を流れる流体の流速とを計測する。これにより、流速計測システム１００は、配管１０の表面１０ａにおける温度分布と、配管１０の内部を流れる流体の流速との実測値における対応関係を取得する。
本実施形態では、流速計測システム１００の制御ユニット４は、上述した、配管１０の温度定常状態における熱の伝わりを有限要素法で解析する機能を有している。
【００７９】
ここで、配管１０内を流れる流体の様々な流速について温度分布および流速を計測すれば、メモリ６３に記憶させる関係情報を取得することができる。しかしながら、関係情報をすべて実測にて取得しようとすると計測回数が多くなり、流速計測システム１００のユーザーにとって、流速計測システム１００を設定する負担（例えば、流体の流速を調節する負担）が大きくなる。
そこで、制御ユニット４が、配管１０の温度定常状態における熱の伝わりを流体の流速ごとに有限要素法で解析して関係情報を取得する。これにより、ユーザーが流速計測システム１００を設定する負担を低減させることができる。制御ユニット４は、関係情報設定部の一例である。
【００８０】
制御ユニット４は、上述した式（１）〜式（９）に基づいて、配管１０の温度定常状態における熱の伝わりを解析する。その際、制御ユニット４は、熱伝達率修正係数ＸＷＴの値、および、温度境界係数ＹＨＴの値を予め（関係情報を取得するための解析を行う前に）設定しておく。制御ユニット４は、温度境界係数値取得部の一例である。
制御ユニット４は、例えば、熱伝達率修正係数ＸＷＴの値と温度境界係数ＹＨＴの値との組み合わせを複数用意しておく。そして、制御ユニット４は、流速計測システム１００が実測した流速（流速計５が計測した流速）について、熱伝達率修正係数ＸＷＴの値と温度境界係数ＹＨＴの値との組み合わせごとに、有限要素法による解析を行って配管１０における温度分布を算出する。そして、制御ユニット４は、配管１０における温度分布の実測値に最も近い解析値を得られた、熱伝達率修正係数ＸＷＴの値と温度境界係数ＹＨＴの値との組み合わせを採用する。
【００８１】
このように、本実施形態では、配管１０の温度分布から配管１０の内部を流れる流体の流速を求める方法として、温度分布と流速との関係を予め求めておき、得られた関係を用いて温度分布を流速に換算する方法が用いられ、この場合に、流速を精度良く求めるために、温度分布と流速との関係を精度良く求めておく。
【００８２】
ここで、本実施形態では、配管１０において、加熱部２および温度計測部３（各温度センサー３ａ）が保温材１２で覆われた構成を例に挙げたが、これに限定されることは無い。例えば、制御ユニット４が配管１０の表面１０ａからの放熱を考慮して温度計測部３から送られる測定データ（温度分布）を補正する態様であれば、配管１０の表面１０ａを保温材１２で被覆しなくてもよい。或いは、表面１０ａの一部（温度計測部３の設置部分）のみを保温材１２で被覆する構成であってもよい。
【００８３】
また、配管１０内を流れる流体は蒸気または空気に限られない。例えば、配管１０内を流れる熱水の流速を計測する場合にも、本実施形態に係る構成を適用することが可能である。また、配管１０内を流れる流体がフロン、アンモニア、ＬＮＧ（Liquefied Natural Gas）等であってもよく、これらの流体の流速を計測する場合にも本実施形態に係る構成を適用することが可能である。
【００８４】
また、本実施形態では、配管１０と熱交換を行う熱交換器として加熱部２を例示したが、これに限定されない。例えば、熱交換器としては、配管１０の表面１０ａを冷却する冷却器を用いてもよく、冷却することで配管１０の管軸方向に生じた温度分布に基づいて配管１０内を流れる蒸気の流速を計測しても良い。この場合において、蒸気が飽和蒸気或いはそれに近い過熱蒸気の時は、凝縮が生じる可能性が有ることから熱伝達率算出の際はそれを考慮する必要がある。
【００８５】
また、流速計測システム１００（図１）において、プレヒーター２ｂは必須の構成ではない。流速計測システム１００の構成を、プレヒーター２ｂを備えていない構成としてもよい。
【００８６】
［フィッティングを使用する解析処理の説明］
図５は、配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速５ｍ／ｓのとき）の例を示す図である。
図６は、配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速２０ｍ／ｓのとき）の例を示す図である。
図７は、配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係（流速３５ｍ／ｓのとき）の例を示す図である。
【００８７】
図５〜図７に示されるグラフにおいて、横軸は配管１０の軸方向の距離［ｍ］を表わしており、縦軸は配管１０の表面１０ａの温度［℃］を表わしている。図５〜図７のそれぞれにおいて、配管１０の軸方向の距離と配管１０の表面１０ａの温度との関係の特性（流速５ｍ／ｓのときの温度分布の特性２０１、流速２０ｍ／ｓのとき温度分布の特性２０２、流速３５ｍ／ｓのとき温度分布の特性２０３）の例を示してある。
各特性２０１〜２０３は、横軸に関して、左右対称の形状（または、左右対称に近い形状）を有しており、当該左右対称の横軸の中心の位置（または、横軸の中心に近い位置）が加熱部２の位置に対応する。
【００８８】
図８は、温度分布の特性２０１における半値幅Ｐ１の例を示す図である。
図９は、温度分布の特性２０１における最高温度Ｐ２の例を示す図である。最高温度Ｐ２は、例えば、ピークの温度を示す。
図１０は、温度分布の特性２０１における温度分布波形の面積Ｐ３の例を示す図である。
温度分布波形の面積Ｐ３は、本実施形態では、温度のベースライン２５１と温度分布波形（山形の線）とで囲まれる領域の面積としてある。この場合、面積Ｐ３は、ベースライン２５１の温度に左右されず、感度が良い。図１０の例において、温度のベースライン２５１としては、左のベースライン上の点（例えば、右端の方の点）と右のベースライン上の点（例えば、左端の方の点）とを結んだ線が用いられている。
なお、他の構成例として、温度分布波形の面積Ｐ３としては、他の面積が用いられてもよい。具体例として、温度分布波形の面積Ｐ３としては、横軸（縦軸の値が０）と温度分布波形（山形の線）とで囲まれる領域の面積が用いられてもよく、または、当該横軸の代わりに縦軸の値が０でない線（当該横軸に平行な線、または、非平行な線）が用いられてもよく、また、面積を求める領域を区切る横軸の値の最小値または最大値のうちの一方または両方が設定されてもよい。
ここで、図８〜図１０は、図５に示されるグラフに対応する例を示した図である。
通常、流速が大きくなると、半値幅Ｐ１は小さくなり、最高温度Ｐ２は小さくなり、温度分布波形の面積Ｐ３は小さくなる。
【００８９】
ここで、上記のように、解析に使用することが可能なパラメータ（温度分布の特徴量）として、半値幅、最高温度、または、温度分布波形の面積などがある。
本実施形態では、解析のパラメータとして、半値幅を使用する場合を示すが、他の構成例として、最高温度が使用されてもよく、または、温度分布波形の面積が使用されてもよく、または、半値幅、最高温度、温度分布波形の面積のうちの２以上を組み合わせたものが使用されてもよい。
【００９０】
図１１は、流速と半値幅との関係の例を示す図である。
図１１に示されるグラフにおいて、横軸は流速［ｍ／ｓ］を表わしており、縦軸は半値幅を表わしている。当該グラフに、流速と半値幅との関係に関し、実測値３０１の例と、解析値３１１の例を示してある。ここで、解析値３１１は、例えば、複数の実測値３０１に対してフィッティングの処理を行った結果として得られた特性（流速と半値幅との関係の特性）である。本例では、解析値３１１を求めるフィッティングにおいて、２次関数の解析値が用いられている。
【００９１】
図１２は、温度分布の実測値４０１と温度分布の解析値４１１との関係の例を示す図である。
図１２に示されるグラフにおいて、横軸は配管１０の軸方向の距離（流れ方向の距離）［ｍ］を表わしており、縦軸は温度［℃］を表わしている。
例えば、複数の実測値４０１のそれぞれの点（ポイント）において、実測値４０１と解析値４１１との差分が誤差となる。また、複数の実測値４０１のそれぞれの点における誤差（実測値４０１と解析値４１１との差分）を２乗した値を、当該複数の実測値４０１の点について総和すると、２乗誤差となる。本例では、解析値４１１を求めるフィッティングにおいて、２次関数の解析値が用いられている。
【００９２】
図１３は、流速と２乗誤差（温度分布の実測値と温度分布の解析値との２乗誤差）との関係の例を示す図である。
図１３に示されるグラフにおいて、横軸は流速［ｍ／ｓ］を表わしており、縦軸は２乗誤差を表わしている。当該グラフに、流速と２乗誤差との関係に関し、実測値５０１の例と、解析値５１１の例を示してある。ここで、解析値５１１は、例えば、複数の実測値５０１に対してフィッティングの処理を行った結果として得られた特性（流速と２乗誤差との関係の特性）である。本例では、解析値５１１を求めるフィッティングにおいて、２次関数の解析値が用いられている。
【００９３】
＜フィッティングを使用する解析処理の第１の例＞
図１４〜図１７を参照して、フィッティングを使用する解析処理の第１の例を説明する。
本実施形態では、説明の便宜上、全体の処理群を第１の処理群と第２の処理群に分けて説明するが、第１の処理群と第２の処理群とは一体化されてもよい。第１の処理群の処理が行われた後に、第２の処理群の処理が行われる。本実施形態では、制御ユニット４において、計算装置５０のＣＰＵ６２がメモリ６３等を用いて、全体の処理群に含まれる各処理を実行する。
【００９４】
図１４は、第１の処理群の一例を示す図である。
図１５は、第１の処理群における流速設定処理の一例を示す図である。
図１６は、第２の処理群の第１の例を示す図である。
図１７は、流速と２乗誤差（温度分布の実測値と温度分布の解析値との２乗誤差）との関係の例を示す図である。
【００９５】
図１４を参照して、第１の処理群を説明する。
（ステップＳ１）
制御ユニット４において、所定の測定値（計測値）を入力する。この測定値は、例えば、配管１０内の流体（例えば、蒸気等）の温度の測定値、および当該流体の圧力の測定値である。当該流体は、例えば、加熱部２で加熱される前のものである。
【００９６】
（ステップＳ２）
制御ユニット４において、所定のパラメータを入力する。このパラメータは、例えば、配管１０の内径、肉厚、材質、熱伝導率、保温材の厚さ等のパラメータである。ここで、本実施形態では、当該パラメータが入力装置６０から入力されるが、他の構成例として、当該パラメータのうちの一部または全部が、予めメモリ６３に記憶されていて使用されてもよい。
なお、（ステップＳ１）の処理と（ステップＳ２）の処理とは、逆の順序で行われてもよい。また、（ステップＳ１）の処理および（ステップＳ２）の処理の段階では、特定の流速を想定していない。
【００９７】
（ステップＳ３）
制御ユニット４において、３点の流速を設定する。
ここで、本実施形態では、３点の流速を設定するが、他の構成例として、４点以上の流速が設定されてもよい。つまり、３点以上の任意の数の流速が設定されてもよい。
【００９８】
図１５を参照して、図１４に示される第１の処理群における流速設定処理（ステップＳ３の処理）の一例を説明する。
（ステップＳ２１）
まず、制御ユニット４において、１点の流速を設定する。
【００９９】
（ステップＳ２２）
次に、制御ユニット４において、設定した流速で半値幅を計算する。具体的には、制御ユニット４において、設定した流速で有限要素法を用いて温度分布の解析値（理論値）を計算し、計算した温度分布の解析値に基づいて半値幅の解析値（理論値）を計算する。
【０１００】
（ステップＳ２３）
次に、制御ユニット４において、半値幅の実測値に基づいて、１点の流速（ステップＳ２２〜ステップＳ２４の最初のループでは、２点目の流速）を設定する。具体的には、制御ユニット４において、既に設定した流速における半値幅の解析値と、計測された半値幅の実測値に基づいて、新たに１点の流速を設定する。
【０１０１】
（ステップＳ２４）
次に、制御ユニット４において、流速の設定を終了するか否かを判定する。制御ユニット４において、流速の設定を終了すると判定した場合には、本処理を終了する。一方、制御ユニット４において、流速の設定を終了しないと判定した場合には、ステップＳ２２の処理へ移行する。
ここで、流速の設定を終了する条件としては、例えば、所定の数の流速が設定された場合に終了するという条件が用いられてもよい。所定の数としては、例えば、３以上の任意の数が用いられてもよい。本実施形態では、所定の数として、３が設定されている。
【０１０２】
図１５に示される流速設定処理（ステップＳ２１〜ステップＳ２４の処理）の具体例を示す。
図１５に関する具体例を示す。
例えば、制御ユニット４において、３点以上の流速を設定する場合、まず、１点以上の流速で温度分布の半値幅の解析値を計算する。次に、制御ユニット４において、計測された半値幅の実測値を取得する。そして、制御ユニット４において、当該解析値と当該実測値に基づいて、新たな流速（次の点の流速）を決定する。
なお、このように流速を決定する処理が、さらに続けて、１回以上繰り返して行われてもよい。
また、これら３点以上の流速は、例えば、以降の処理を行うために仮に設定された流速（ダミーの流速）と捉えることも可能である。
【０１０３】
図１５に関する詳細な具体例を示す。
まず、制御ユニット４において、１点目の流速として、想定される流速の範囲のうちの中点の流速を設定する。一例として、想定される流速の範囲が０〜４０［ｍ／ｓ］である場合、中点の流速は２０［ｍ／ｓ］である。次に、制御ユニット４において、設定した中点の流速における半値幅の解析値を計算し、また、計測された半値幅の実測値を取得する。
この結果、半値幅の実測値が解析値よりも小さい場合には、制御ユニット４において、２点目の流速として、中点の流速よりも大きい流速を決定して設定する。前記一例では、この流速としては、２０〜４０のうちの任意の流速が設定され、例えば、２０〜４０のうちの中点である３０の流速が設定されてもよい。
一方、半値幅の実測値が解析値よりも大きい場合には、制御ユニット４において、２点目の流速として、中点の流速よりも小さい流速を決定して設定する。前記一例では、この流速としては、０〜２０のうちの任意の流速が設定され、例えば、０〜２０のうちの中点である１０の流速が設定されてもよい。
【０１０４】
本実施形態では、説明を簡易化するために、一例として、想定される流速の範囲が０〜４０［ｍ／ｓ］であるとしたが、０［ｍ／ｓ］は現実的でなく、実際には、例えば、配管１０内の蒸気等の流速は２〜４０［ｍ／ｓ］程度であると考えられる。
【０１０５】
図１５に関する詳細な具体例の続きを示す。
さらに、制御ユニット４において、３点目の流速を設定する。
具体的には、制御ユニット４において、１点目の流速における半値幅の解析値と２点目の流速における半値幅の解析値との間に半値幅の実測値が位置する場合には、１点目の流速と２点目の流速との間に位置する流速を、３点目の流速として決定して設定する。１点目の流速と２点目の流速との間に位置する流速としては、一例として、１点目の流速と２点目の流速との中点の流速が用いられてもよい。
一方、制御ユニット４において、１点目の流速における半値幅の解析値と２点目の流速における半値幅の解析値との間に半値幅の実測値が位置しない場合には、１点目の流速から見て２点目の流速を超えた位置の流速を、３点目の流速として決定して設定する。１点目の流速から見て２点目の流速を超えた位置の流速としては、一例として、１点目の流速から見て、２点目の流速までの距離の１．５倍の位置にある流速が用いられてもよい。この場合に、例えば、流速の範囲（最大値と最小値）が想定されて設定されてもよい。そして、１点目の流速から見て２点目の流速の方が大きい場合、前記した１．５倍の位置の流速が想定範囲を超えているときには、当該想定範囲のうちで最大の流速が用いられてもよい。逆に、１点目の流速から見て２点目の流速の方が小さい場合、前記した１．５倍の位置の流速が想定範囲未満であるときには、当該想定範囲のうちで最小の流速が用いられてもよい。
なお、同様にして、制御ユニット４において、４点目以降の流速を設定してもよい。
【０１０６】
ここで、他の構成例として、図１４および図１５に示される処理において３点以上の流速を設定する方法として、他の方法が用いられてもよい。一例として、想定される流速の範囲（一例として、０〜４０［ｍ／ｓ］）で、任意の３点以上の流速を選択（決定）して設定する方法が用いられてもよい。具体的には、例えば、５［ｍ／ｓ］、２０［ｍ／ｓ］、３５［ｍ／ｓ］という流速の組み合わせのように、すべての流量の範囲にわたって偏りなく流量が分布するように３点以上の流量を選択（決定）して設定する方法が用いられてもよい。
【０１０７】
図１６を参照して、第２の処理群の第１の例を説明する。図１６に示される第２の処理群の処理は、図１４に示される第１の処理群の処理の後に行われる。
（ステップＳ４１）
まず、制御ユニット４において、設定した３点以上の流速のうちの各流速において、有限要素法を用いて得られる温度分布の解析値（理論値）と、計測された温度分布の実測値とについて、両者の２乗誤差を計算する。
【０１０８】
（ステップＳ４２）
次に、制御ユニット４において、流速と２乗誤差との関係（解析値）をフィッティングにより求める。フィッティングでは、例えば、２次関数が用いられる。
（ステップＳ４３）
次に、制御ユニット４において、フィッティングの結果（解析値）に基づいて、流速の候補を決定する。
一例として、制御ユニット４において、フィッティングの結果（解析値）において２乗誤差が最も小さくなる流速を、流速の候補として決定する。
ここで、図１７には、図１３と同様なグラフにおいて、２乗誤差が最小値の流速Ｑ１を示してある。図１７に示されるグラフにおいて、横軸は流速［ｍ／ｓ］を表わしており、縦軸は２乗誤差を表わしている。当該グラフに、流速と２乗誤差との関係に関し、実測値６０１の例と、解析値６１１の例を示してある。
【０１０９】
（ステップＳ４４）
次に、制御ユニット４において、決定した流速の候補において、有限要素法を用いて得られる温度分布の解析値（理論値）と、計測された温度分布の実測値とについて、両者の２乗誤差を計算する。そして、制御ユニット４において、計算された２乗誤差が所定の規定値以下であるか否かを判定する。規定値としては、任意の値が用いられてもよく、例えば、予め設定されてもよい。
【０１１０】
（ステップＳ４５）
制御ユニット４において、上記した（ステップＳ４４）の判定の結果、２乗誤差が所定の規定値以下であると判定した場合には、候補とした流速を、正式に決定する。そして、制御ユニット４において、本処理を終了する。
【０１１１】
（ステップＳ４６）
一方、制御ユニット４において、上記した（ステップＳ４４）の判定の結果、２乗誤差が所定の規定値を超えると判定した場合には、流速（３点以上の流速の組み合わせ）を変更して、上記した（ステップＳ４１）の処理へ移行する。
ここで、流速（３点以上の流速の組み合わせ）を変更する方法としては、一例として、もとの３点以上の流速のうちの１点の流速を除いて、その代わりに候補とした１点の流速を加える方法が用いられてもよい。除く対象とする１点の流速としては、一例として、これら複数点の流速のうちで最も２乗誤差が大きい点の流速が用いられてもよい。
また、（ステップＳ４１）〜（ステップＳ４４）、（ステップＳ４６）の処理のループは、何回行われてもよいが、通常、例えば、１〜５回程度で２乗誤差が規定値以下になると予想される。
【０１１２】
ここで、本実施形態では、図１５に示される処理において新たな点の流速を設定するために使用された流速（例えば、３点の流速を設定する場合には、１点目の流速および２点目の流速）は、図１６に示される（ステップＳ４１）の処理における３点以上の流速に含められて使用されている。これにより、図１５に示される処理において設定された流速が、図１６に示される（ステップＳ４１）の処理において有効に使用されることが実現される。
なお、他の構成例として、図１５に示される処理において新たな点の流速を設定するために使用された流速の一部（例えば、４点以上の流速を設定する場合には、１点目の流速など）が除かれて、図１６に示される（ステップＳ４１）の処理における３点以上の流速が設定されてもよい。例えば、図１５に示される処理において４点以上のＺ（Ｚは４以上の整数）点の流速が設定される場合には、（Ｚ−３）以下の数の点の流速が除かれて、図１６に示される（ステップＳ４１）の処理が行われてもよい。
【０１１３】
以上のように、図１４〜図１６に示されるフィッティングを使用する解析処理の第１の例では、流速の評価に要する時間を短縮化することが可能である。また、流速の評価の精度を向上させることが可能である。
例えば、図１３に示される流速と２乗誤差との関係のフィッティングに関し、精度の向上、および処理速度の向上を図ることが可能である。
なお、本実施形態では、２乗誤差が用いられたが、他の構成例として、他の誤差が用いられてもよい。
【０１１４】
＜フィッティングを使用する解析処理の第２の例＞
図１４、図１５、図１８を参照して、フィッティングを使用する解析処理の第２の例を説明する。
本実施形態では、説明の便宜上、全体の処理群を第１の処理群と第２の処理群に分けて説明するが、第１の処理群と第２の処理群とは一体化されてもよい。第１の処理群の処理が行われた後に、第２の処理群の処理が行われる。本実施形態では、制御ユニット４において、計算装置５０のＣＰＵ６２がメモリ６３等を用いて、全体の処理群に含まれる各処理を実行する。
【０１１５】
ここで、図１４および図１５に示される処理（第１の処理群の処理）については、上記した＜フィッティングを使用する解析処理の第１の例＞の場合と同様である。
まず、制御ユニット４において、図１４および図１５に示される処理（第１の処理群の処理）を行う。
【０１１６】
図１８は、第２の処理群の第２の例を示す図である。
図１８を参照して、第２の処理群の第２の例を説明する。図１８に示される第２の処理群の処理は、図１４に示される第１の処理群の処理の後に行われる。
（ステップＳ６１）
まず、制御ユニット４において、設定した３点以上の流速のうちの各流速において、有限要素法を用いて温度分布の解析値（理論値）を計算し、流速と温度分布の半値幅（解析値）との関係をフィッティングにより計算する。フィッティングでは、例えば、２次関数が用いられる。
【０１１７】
（ステップＳ６２）
次に、制御ユニット４において、流速と温度分布の半値幅との関係のフィッティングの結果（解析値）と、計測された温度分布の半値幅の実測値とに基づいて、フィッティングの結果において半値幅の実測値に対応する流速を求めて、当該流速を、流速の候補として決定する。
【０１１８】
（ステップＳ６３）
次に、制御ユニット４において、決定した流速の候補において、有限要素法を用いて得られる温度分布の解析値（理論値）と、計測された温度分布の実測値とについて、両者の２乗誤差を計算する。そして、制御ユニット４において、計算された２乗誤差が所定の規定値以下であるか否かを判定する。規定値としては、任意の値が用いられてもよく、例えば、予め設定されてもよい。
【０１１９】
（ステップＳ６４）
制御ユニット４において、上記した（ステップＳ６３）の判定の結果、２乗誤差が所定の規定値以下であると判定した場合には、候補とした流速を、正式に決定する。そして、制御ユニット４において、本処理を終了する。
一方、制御ユニット４において、上記した（ステップＳ６３）の判定の結果、２乗誤差が所定の規定値を超えると判定した場合には、上記した（ステップＳ６１）の処理へ移行する。このとき、制御ユニット４において、（ステップＳ６１）の処理を行う場合に、３点以上の流速の組み合わせを、既に使用した流速の組み合わせとは異なるものに変更する。このような流速の変更の処理は、任意の方法により行われてもよく、一例として、３点以上の流速のうちの所定の１点以上の流速に対して所定値（例えば、プラスの値、または、マイナスの値）を加えた結果の値を使用する方法が用いられてもよい。
【０１２０】
ここで、本実施形態では、図１５に示される処理において新たな点の流速を設定するために使用された流速（例えば、３点の流速を設定する場合には、１点目の流速および２点目の流速）は、図１８に示される（ステップＳ６１）の処理における３点以上の流速に含められて使用されている。これにより、図１５に示される処理において設定された流速が、図１８に示される（ステップＳ６１）の処理において有効に使用されることが実現される。
なお、他の構成例として、図１５に示される処理において新たな点の流速を設定するために使用された流速の一部（例えば、４点以上の流速を設定する場合には、１点目の流速など）が除かれて、図１８に示される（ステップＳ６１）の処理における３点以上の流速が設定されてもよい。例えば、図１５に示される処理において４点以上のＺ（Ｚは４以上の整数）点の流速が設定される場合には、（Ｚ−３）以下の数の点の流速が除かれて、図１８に示される（ステップＳ６１）の処理が行われてもよい。
【０１２１】
以上のように、図１４、図１５および図１８に示されるフィッティングを使用する解析処理の第２の例では、流速の評価に要する時間を短縮化することが可能である。また、流速の評価の精度を向上させることが可能である。
例えば、図１１に示される流速と半値幅との関係のフィッティングに関し、精度の向上、および処理速度の向上を図ることが可能である。
なお、本実施形態では、２乗誤差が用いられたが、他の構成例として、他の誤差が用いられてもよい。
【０１２２】
［以上の実施形態のまとめ］
一構成例として、配管１０の内部を流れる流体の流速を評価する評価方法であって、１点以上の流速について温度分布の特徴量の解析値（理論値）に基づいて、流速を決定する、流速の評価方法である。
一構成例として、流速の評価方法において、温度分布の特徴量は、流速の増加に応じて減少する特徴量である。
一構成例として、流速の評価方法において、温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である。
一構成例として、流速の評価方法において、前記配管１０の表面１０ａの所定部分で熱交換を行う熱交換工程と、前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算工程と、前記計算工程により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定工程と、前記流速設定工程により設定された流速を用いて、流速を決定する流速決定工程と、を含む。
一構成例として、流速の評価方法において、前記流速設定工程は、想定された流速の範囲の中点の流速について前記計算工程により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する（ステップＳ２３）。
一構成例として、流速の評価方法において、前記流速設定工程は、３点以上の流速を設定し、前記流速決定工程は、前記流速設定工程により設定された前記３点以上の流速を用いて、流速を決定する。
一構成例として、流速の評価方法において、前記流速決定工程は、前記流速設定工程で使用された流速の一部または全部を用いる。
一構成例として、流速の評価方法において、前記流速決定工程は、前記流速設定工程により設定された流速のそれぞれにおいて、（有限要素法を用いて得られた）温度分布の解析値と（計測された温度分布の）実測値との２乗誤差を計算する第１−１工程（ステップＳ４１）と、前記流速と前記２乗誤差との関係を（フィッティングにより）求める第１−２工程（ステップＳ４２）と、流速の候補を決定する第１−３工程（ステップＳ４３）と、決定された流速の候補において、温度分布の解析値と（計測された温度分布の）実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程（ステップＳ４４）と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程（ステップＳ４５）と、前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程（ステップＳ４６）と、を含む。
一構成例として、流速の評価方法において、前記流速決定工程は、前記流速設定工程により設定された流速のそれぞれにおいて、流速と（有限要素法を用いて得られた）温度分布の前記特徴量との関係の解析値を計算する第２−１工程（ステップＳ６１）と、計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程（ステップＳ６２）と、決定された流速の候補において、温度分布の解析値と（計測された温度分布の）実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程（ステップＳ６３）と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程（ステップＳ６３、ステップＳ６４）と、を含む。
【０１２３】
一構成例として、配管１０の内部を流れる流体の流速を評価する評価システムであって、１点以上の流速について温度分布の特徴量の解析値（理論値）に基づいて、流速を決定する、流速の評価システム（実施形態では、流速計測システム１００）である。
一構成例として、流速の評価システムにおいて、温度分布の特徴量は、流速の増加に応じて減少する特徴量である。
一構成例として、流速の評価システムにおいて、温度分布の特徴量は、半値幅、最高温度、または、温度分布波形の面積のうちの１以上である。
一構成例として、流速の評価システムにおいて、前記配管１０の表面１０ａの所定部分で熱交換を行う熱交換器（実施形態では、加熱部２）と、前記１点以上の流速について温度分布の特徴量の前記解析値を計算する計算部（実施形態では、制御ユニット４の機能）と、前記計算部により計算された前記解析値と、前記配管に関して計測された温度分布の特徴量の実測値に基づいて、流速を設定する流速設定部（実施形態では、制御ユニット４の機能）と、前記流速設定部により設定された流速を用いて、流速を決定する流速決定部（実施形態では、制御ユニット４の機能）と、を含む。
一構成例として、流速の評価システムにおいて、前記流速設定部は、想定された流速の範囲の中点の流速について前記計算部により計算された前記解析値と、計測された温度分布の特徴量の前記実測値との大小関係に応じて、新たな点の流速を設定する。
一構成例として、流速の評価システムにおいて、前記流速設定部は、３点以上の流速を設定し、前記流速決定部は、前記流速設定部により設定された前記３点以上の流速を用いて、流速を決定する。
一構成例として、流速の評価システムにおいて、前記流速決定部は、前記流速設定部で使用された流速の一部または全部を用いる。
一構成例として、流速の評価システムにおいて、前記流速決定部は、前記流速設定部により設定された流速のそれぞれにおいて、（有限要素法を用いて得られた）温度分布の解析値と（計測された温度分布の）実測値との２乗誤差を計算する第１−１工程と、前記流速と前記２乗誤差との関係を（フィッティングにより）求める第１−２工程と、流速の候補を決定する第１−３工程と、決定された流速の候補において、温度分布の解析値と（計測された温度分布の）実測値との２乗誤差が所定の規定値以下であるか否かを判定する第１−４工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する第１−５工程と、前記２乗誤差が前記所定の規定値を超えると判定した場合には、流速を変更して、再び前記第１−１工程へ移行する第１−６工程と、を実行する。
一構成例として、流速の評価システムにおいて、前記流速決定部は、前記流速設定部により設定された流速のそれぞれにおいて、流速と（有限要素法を用いて得られた）温度分布の前記特徴量との関係の解析値を計算する第２−１工程と、計算された前記解析値と、計測された温度分布の特徴量に基づいて、流速の候補を決定する第２−２工程と、決定された流速の候補において、温度分布の解析値と（計測された温度分布の）実測値との２乗誤差が所定の規定値以下であるか否かを判定する第２−３工程と、前記２乗誤差が前記所定の規定値以下であると判定した場合には、前記流速の候補を正式に決定する一方、他の場合には、流速を変更して、再び前記第２−１工程へ移行する第２−４工程と、を実行する。
【０１２４】
なお、本実施形態では、以下のような構成が用いられてもよい。
一構成例として、配管の内部を流れる流体の流速と前記配管の表面の温度分布との関係を示す関係情報を求める関係情報設定方法であって、前記配管の表面の所定部分で熱交換を行う熱交換工程と、前記配管を流れる流体のうち、前記配管の表面の所定部分で行われる熱交換による熱の伝達に寄与する流体の量的割合を示す温度境界係数の値を取得する温度境界係数値取得工程と、前記温度境界係数値取得工程で取得した温度境界係数の値に基づいて前記関係情報を求める関係情報設定工程と、を含む関係情報設定方法が用いられてもよい。
一構成例として、関係情報設定方法において、前記温度境界係数値取得工程では、前記流体の流速ごとに前記温度境界係数の値を取得する、構成が用いられてもよい。
一構成例として、関係情報設定方法において、前記所定部分で熱交換が行われた前記配管の管軸方向における前記配管の表面の温度分布を計測する温度計測工程と、前記流体の流速を計測する流速計測工程と、を含み、前記温度境界係数値取得工程では、前記温度計測工程で計測した前記温度分布、及び、前記流速計測工程で計測した流速に基づいて前記温度境界係数の値を設定する、構成が用いられてもよい。
一構成例として、配管の内部を流れる流体の流速が計測対象の流速となっている状態で、前記配管の表面の所定部分で熱交換を行う流速計測時熱交換工程と、前記流速計測時熱交換工程にて前記所定部分で熱交換が行われた前記配管の管軸方向における前記表面の温度分布を計測する流速計測時温度分布計測工程と、前記流速計測時温度分布計測工程で計測した前記温度分布、及び、関係情報設定方法にて得られた関係情報に基づいて、前記配管の内部を流れる前記流体の流速を求める流速取得工程と、を含む流速計測方法が用いられてもよい。
【０１２５】
一構成例として、配管の内部を流れる流体の流速と前記配管の表面の温度分布との関係を示す関係情報を求める関係情報設定システムであって、前記配管の表面の所定部分で熱交換を行う熱交換器と、前記配管を流れる流体のうち、前記配管の表面の所定部分で行われる熱交換による熱の伝達に寄与する流体の量的割合を示す温度境界係数の値を取得する温度境界係数値取得部と、前記温度境界係数値取得部が取得した温度境界係数の値に基づいて前記関係情報を求める関係情報設定部と、を備える関係情報設定システムが用いられてもよい。
一構成例として、関係情報設定システムにおいて、前記温度境界係数値取得部は、前記流体の流速ごとに前記温度境界係数の値を取得する、構成が用いられてもよい。
一構成例として、関係情報設定システムにおいて、前記所定部分で熱交換が行われた前記配管の管軸方向における前記配管の表面の温度分布を計測する温度計測部と、前記流体の流速を計測する流速計測部と、を備え、前記温度境界係数値取得部は、前記温度計測部が計測した前記温度分布、及び、前記流速計測部が計測した流速に基づいて前記温度境界係数の値を設定する、構成が用いられてもよい。
一構成例として、関係情報設定システムの前記関係情報設定部が設定した関係情報を記憶する記憶部と、配管の内部を流れる流体の流速が計測対象の流速となっている状態で、前記配管の表面の所定部分で熱交換を行う熱交換器と、前記所定部分で熱交換が行われた前記配管の管軸方向における前記配管の表面の温度分布を計測する温度計測部と、前記温度計測部が計測した温度分布、及び、前記記憶部が記憶している関係情報に基づいて、前記流体の流速を算出する流速算出部と、を備える流速計測システムが用いられてもよい。
【０１２６】
以上、本発明の実施形態について図面を参照して詳述したが、具体的な構成はこの実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【０１２７】
また、本実施形態のように、以上に示した実施形態に係る装置（例えば、計算装置５０）の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、処理を行ってもよい。
なお、ここでいう「コンピュータシステム」とは、オペレーティング・システム（ＯＳ：Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）や周辺機器等のハードウェアを含むものであってもよい。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、フラッシュメモリ等の書き込み可能な不揮発性メモリ、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｋ）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。
さらに、「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（例えばＤＲＡＭ（Ｄｙｎａｍｉｃ  Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ））のように、一定時間プログラムを保持しているものも含むものとする。
また、上記のプログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。
また、上記のプログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、上記のプログラムは、前述した機能をコンピュータシステムに既に記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【符号の説明】
【０１２８】
２…加熱部、２ｂ…プレヒーター、３…温度計測部、３ａ…温度センサー、３Ａ…温度センサー群、４…制御ユニット、５…流速計、１０…配管、１０ａ…表面、１１…設置部分、１２…保温材、２０…蒸気製造装置、３０…負荷設備、５０…計算装置、６０…入力装置、６１…変換器、６２…ＣＰＵ、６３…メモリ、６４…表示装置、１００…流速計測システム、２０１〜２０３…特性、２５１…ベースライン、３０１、４０１、５０１、６０１…実測値、３１１、４１１、５１１、６１１…解析値
【図１】
2017072424.tif 000003
【図２】
2017072424.tif 000004
【図３】
2017072424.tif 000005
【図４】
2017072424.tif 000006
【図５】
2017072424.tif 000007
【図６】
2017072424.tif 000008
【図７】
2017072424.tif 000009
【図８】
2017072424.tif 000010
【図９】
2017072424.tif 000011
【図１０】
2017072424.tif 000012
【図１１】
2017072424.tif 000013
【図１２】
2017072424.tif 000014
【図１３】
2017072424.tif 000015
【図１４】
2017072424.tif 000016
【図１５】
2017072424.tif 000017
【図１６】
2017072424.tif 000018
【図１７】
2017072424.tif 000019
【図１８】
2017072424.tif 000020