JP2016004007A - 流量算出方法 - Google Patents

Abstract

【課題】開放空間内であっても、壁面の穴から流出または流入する流れの流量を算出できる流量算出方法を提供することを目的とする。
【解決手段】本発明にかかる流量算出方法の代表的な構成は、開放空間の壁面１３４の穴１３６から流出または流入する流れの流量を算出する方法であって、穴の位置にレーザーシート１３０を形成し、流れの中にあるトレーサー１３２にレーザーシートの光が反射した輝点Ｄを連続的に撮影することにより流れの速度分布を算出し、穴を囲むように境界線１５２を設定し、境界線上のプロットの速度成分に単位面積をかけて単位流量を算出し、単位流量を積分して流れの流量を算出する。
【選択図】図４

Description

本発明は、ＰＩＶ計測システムを用いた流量算出方法に関するものである。

原子炉圧力容器内は、冷却材としての水が注入されている。そこで原子炉圧力容器内の状態を知るために、原子炉圧力容器内の水の流れを知りたいという要請がある。特許文献１には、原子炉圧力容器内の水の流れを知る方法として、放射線環境下でのＰＩＶ（Particle Imaging Velocity：粒子画像流速測定法）計測システムが提案されている。ＰＩＶ計測とは、流れの中にレーザーシートを形成してこれを撮影し、粒子（トレーサー）に反射した輝点の動きを観察することによって流れを把握する方法である。なお、このように流れを可視化する技術を総称してＦＶ（Flow Visualization）と呼ばれている。

特許文献２には、ＰＩＶ計測システムを用いた粒子画像流速測定装置が記載されている。この装置では、トレーサーが浮遊する空間に相互に平行かつ波長が異なる二つの２次元のレーザーシートを照射し、レーザーシートに対して直交する方向に配置した２台のカメラで各々のレーザーシート内に存在するトレーサーを撮像する。つぎに装置では、得られた二つの画像のトレーサーの輝度パターンを比較することで、レーザーシートの面内のＸ方向およびＹ方向の速度成分（面内速度成分）と、レーザーシートに直交するＺ方向の速度成分（面外速度成分）とを含む３速度成分を測定している。

特開２００６−０１７６１６号公報 特開平１０−１９９１９号公報

ところで、原子力発電所が過酷事故（シビアアクシデント）を生じている場合、原子炉圧力容器の内部のみならず、原子炉格納容器の内部にも水が注入されている場合がある。そしてこれらの場合には、通常状態よりもさらに内部の様子、例えば原子炉格納容器や原子炉圧力容器などの壁面の穴から流出または流入する流体の流れの様子を知りたいという要請が高い。なお原子炉格納容器や原子炉圧力容器のようなある程度の大きさを有する密閉空間では、その内部の流体の流れが管内のように一方向に生じるものとは異なっている。以下では、このような空間を開放空間と称する。

特許文献１は、耐放射線性能を備える伝送ケーブルを用いるなど、ＰＩＶ計測システムを放射線環境下でも使用可能とする技術である。しかしながら、特許文献１には、開放空間内で壁面の穴から流出または流入する流れの流量を算出する手法について具体的には開示されていない。

特許文献２は、半導体製造工程で使用されるクリーンルーム内で雰囲気中に浮遊する微粒子の発生源を、ＰＩＶ計測システムによって特定するものに過ぎない。すなわち特許文献２に記載の技術では、原子力発電所の過酷事故の状況下で、開放空間内の流れの流量を算出することはできない。

本発明は、このような課題に鑑み、開放空間内であっても、壁面の穴から流出または流入する流れの流量を算出できる流量算出方法を提供することを目的としている。

上記課題を解決するために、本発明にかかる流量算出方法の代表的な構成は、開放空間の壁面の穴から流出または流入する流れの流量を算出する方法であって、穴の位置にレーザーシートを形成し、流れの中にある粒子にレーザーシートの光が反射した反射光を連続的に撮影することにより流れの速度分布を算出し、穴を囲むように境界線を設定し、境界線上のプロットの速度成分に単位面積をかけて単位流量を算出し、単位流量を積分して流れの流量を算出することを特徴とする。

上記構成によれば、開放空間（例えば、原子炉格納容器など）の壁面の穴を囲むように境界線を設定することで、境界線上のプロットを通過する流れの速度成分を算出できる。そして、このプロットの速度成分に単位面積をかけて単位流量を算出し、この単位流量を穴の中心軸周りに積分することで、流れの流量を算出できる。なお単位面積は、境界線で囲まれた領域を境界線上のプロットで分割することで得られる。したがって、上記方法によれば、開放空間内であっても、壁面の穴から流出または流入する流れの流量を算出できる。

上記の速度分布は、マトリクス状のプロットについて算出されていて、境界線は、壁面に直交し穴の中心位置を通る中心軸を基準として互いに線対称となる２つの側線と、壁面に対向し側線の端同士をつなぐ上線とを含む矩形であるとよい。

このように、プロットがマトリクス状に配列されていて、境界線が矩形（コの字型）であるので、境界線の側線上および上線上にプロットを設定し易い。このため、側線上のプロットの速度成分（側線に対する垂直成分）を用いて、側線での流れの流量を算出できる。また、上線上のプロットの速度成分（上線に対する垂直成分）を用いて、上線での流れの流量を算出できる。さらに、側線および上線での流れの流量を加算することで、壁面の穴から流出または流入する流れの流量を算出できる。

上記の境界線は、その内側に最も速度が大きいプロットを含むように設定するとよい。レーザーシートは空間内の１つの断面しか照射しないので、境界線は基本的にはなるべく穴に近いのがよい。しかしながら、あまり穴に近い位置のプロットでは、流入する場合には流れに乱れが大きく、流出する場合にはトレーサーが穴に吸い込まれてしまうために、適切な速度を取れないおそれがある。そこで境界線の内側に最も速度が大きいプロットを含むように設定することにより、穴から適切に距離をおいて境界線を設定することができる。これにより壁面の穴から流出または流入する流れを最もよく捉えることができ、高い精度で流量を算出することが可能となる。

上記の流れの中にある粒子と理想的なトレーサーとを比較して補正係数を予め算出し、境界線上のプロットの速度に補正係数をかけて、速度を補正するとよい。ここで、例えば原子炉格納容器や原子炉圧力容器の内部には、流れの中に浮遊する粒子（水酸化鉄など）が存在している。このような粒子は、流れに追従する理想的なトレーサーに比べて追従性が劣る可能性が高い。そこで、上記構成では、流れの中にある粒子をトレーサーとして利用する場合、理想的なトレーサーとの追従性の差を考慮した補正係数を予め算出する。そして、この補正係数によって、境界線上のプロットの速度を補正することにより、流れの流量を正確に算出できる。

本発明によれば、開放空間内であっても、壁面の穴から流出または流入する流れの流量を算出できる流量算出方法を提供することができる。

本実施形態における流量算出方法が適用される原子炉の概略構成を説明する図である。 本実施形態におけるＰＩＶ計測の手法を概略的に示す図である。 撮像部で輝点を撮影する様子を例示する図である。 本実施形態における流量算出方法の各処理を示すフローチャートである。 穴から流出する流れを速度ベクトルで示す図である。 図５の速度ベクトルを解析する手法を示す図である。

以下に添付図面を参照しながら、本発明の好適な実施形態について詳細に説明する。かかる実施形態に示す寸法、材料、その他具体的な数値などは、発明の理解を容易とするための例示に過ぎず、特に断る場合を除き、本発明を限定するものではない。なお、本明細書及び図面において、実質的に同一の機能、構成を有する要素については、同一の符号を付することにより重複説明を省略し、また本発明に直接関係のない要素は図示を省略する。

図１は、本実施形態における流量算出方法が適用される原子炉の概略構成を説明する図である。以下では、理解を容易にするために沸騰水型原子炉（ＢＷＲ）を例示して説明するが、改良型沸騰水型原子炉（ＡＢＷＲ）についても本実施形態を適用可能である。

原子炉建屋１００には原子炉格納容器１０２が設置され、その中に原子炉圧力容器１０４が設置されている。原子炉格納容器１０２の下方にはサプレッションプール１０６（圧力抑制室ともいう）が設けられていて、過剰な蒸気をサプレッションプール１０６内の水に注入することによって圧力を下げることができる。

原子炉格納容器１０２および原子炉圧力容器１０４は、いずれも密閉容器であり、その内部に密閉空間を有している。なお、原子炉圧力容器１０４は高さ２０ｍ程度、直径６ｍ程度であり、また、原子炉格納容器１０２は高さ４０ｍ程度、直径１８ｍ程度の大きな設備である。このような大きな密閉空間を有する原子炉格納容器１０２や原子炉圧力容器１０４は、その内部の流体の流れが管内のように一方向に生じるものとは異なる、いわば開放空間を成している。

原子炉圧力容器１０４は、ウラン等からなる燃料棒を束ねた燃料集合体１０８を収容する。原子炉圧力容器１０４に収容された燃料集合体１０８は、臨界に達して、容器内部に充填された水（炉水）を加熱し蒸気を発生させる。

図１は事故時の状態を示している。図中では原子炉格納容器１０２の底部、原子炉圧力容器１０４の内部およびサプレッションプール１０６の部屋に水が溜まっているように描いている。なお通常の運転中は、原子炉格納容器１０２は窒素充填され、サプレッションプール１０６内部には適量の水が保有されている。過酷事故の状況下では、原子炉格納容器１０２および原子炉圧力容器１０４の内部には、放射線の遮蔽と炉心の冷却を図るために冷却水が注入されることになり、図示のように冷却水１１０、１１２が滞留する。

過酷事故によって燃料（燃料ペレット）が過熱すると、燃料被覆管や燃料集合体１０８のケース（チャネルボックス）が融解する（炉心溶融：メルトダウン）。すると燃料および炉の構造物の溶融物が混在した、いわゆるデブリが発生する。デブリは原子炉圧力容器１０４の下部に堆積すると考えられるが、程度によっては原子炉圧力容器の底をも融解させて原子炉格納容器１０２まで落下する場合もある（メルトスルー）。したがって図１に示すように、デブリ１１４が原子炉圧力容器１０４や原子炉格納容器１０２の底部に堆積すると考えられる。

また過酷事故の状況下では、原子炉格納容器１０２や原子炉圧力容器１０４などの壁面の穴（排水口や亀裂、スリット状の穴など）から流出または流入する冷却水１１０、１１２などの流体の流れの様子を知りたいという要請が高い。そこで、本実施形態では、ＰＩＶ計測システムを用いた流量算出方法によって、一例として冷却水１１０の流量を算出する。

過酷事故の状況下では、原子炉格納容器１０２の蓋を開くことができないため、原子炉格納容器１０２の壁面に設置されている検査孔１１６から、流量測定装置１１８の一部であるレーザーシート用の光学ユニットおよび撮影ユニットを挿入する必要がある。流量測定装置１１８は、検査孔１１６に挿入される挿入管１２０と、挿入管１２０の先端に接続されたヘッド１２２とを備える。

挿入管１２０は、可撓性を有するフレキシブルチューブであって、流量測定装置１１８の本体とヘッド１２２との間を接続する光ケーブルや信号ケーブルを挿通している。流量測定装置１１８は、ヘッド１２２から連続的に画像を取得して、画像間（フレーム間）で粒子の輝点の位置を比較することにより、輝点の移動方向と移動量を速度ベクトルで示し、これを解析することで流れの流量を算出する（詳細は後述）。

図２は、本実施形態におけるＰＩＶ計測の手法を概略的に示す図である。図２（ａ）は、流量測定装置１１８によるＰＩＶ計測の概略構成を示す図である。図２（ｂ）は、流量計測装置１１８の機能を示すブロック図である。

流量計測装置１１８のヘッド１２２は、図２（ａ）に示すように、治具１２４に固定された光学ユニット１２６と撮影ユニット（撮像部１２８）とを備える。光学ユニット１２６は、流体である冷却水１１０の流れの中にレーザーシート１３０を形成する。光学ユニット１２６は、挿入管１２０を通る光ケーブル介して流量計測装置１１８から送られてきたレーザー光を照射させてレーザーシート１３０を形成する。レーザーシート１３０が冷却水１１０の中に浮遊する水酸化鉄などの粒子（トレーサー１３２）に反射すると輝点（図３参照）を生じる。

撮像部１２８は、撮影方向Ａに沿ってこれらの輝点を撮影する。撮像部１２８は、カメラ（ＣＣＤやＣＭＯＳなどの撮像素子）を有し、挿入管１２０を通る信号ケーブルによって画像データが流量測定装置１１８の本体に伝達される。あるいは撮像部１２８にレーザーシート１３０を撮影するための対物レンズなどを設置し、カメラは流量測定装置１１８の本体に備えるようにしてもよい。この場合には、挿入管１２０を通る光ケーブルによって画像（光）が流量測定装置１１８の本体に伝達される。なお図２（ａ）中の矢印Ｂは、開放空間である原子炉格納容器１０２の壁面１３４の穴１３６から流出する流体の流れの方向を示している。また領域Ｃは、撮像部１２８による撮影範囲を示している。

流量測定装置１１８は、本実施形態にかかる流量算出方法を実現するものであり、図２（ｂ）に示すように、上記撮像部１２８と、流量算出部１３８と、出力部１４０とを備える。流量算出部１３８は、画像処理部１４２と、ベクトル解析部１４４と、補正係数格納部１４６とを備える。

図３は、撮像部１２８で輝点Ｄを撮影する様子を例示する図である。撮像部１２８は、図示のように、壁面１３４の穴１３６の位置に形成されたレーザーシート１３０が撮影範囲（領域Ｃ）内に位置するように撮影方向Ａを適宜調整する。すなわち、撮像部１２８の撮影範囲には、図示のようにレーザーシート１３０のうち穴１３６を横切る照射ライン１４８が含まれている。また輝点Ｄは、冷却水１１０の中の上記トレーサー１３２にレーザーシート１３０の光が反射した反射光である。

図４は、本実施形態における流量算出方法の各処理を示すフローチャートである。まず、流量測定装置１１８では、光学ユニット１２６を用いて穴１３６の位置にレーザーシート１３０を形成する（ステップＳ１００）。つぎに、撮像部１２８は、図３に示す輝点Ｄを連続的に撮影する（ステップＳ１０２）。

流量算出部１３８の画像処理部１４２は、撮像部１２８から画像データを例えば１／３０秒間隔で５００枚分取得し、取得したデータをベクトル解析部１４４に出力する。続いて、ベクトル解析部１４４は、画像処理部１４２から得た連続する各時刻の画像データの相関に基づいて、輝点Ｄの移動方向と移動量を示す速度ベクトルを求め、これを解析し流れの速度分布を算出する（ステップＳ１０４）。ステップＳ１０４の後、ベクトル解析部１４４は、穴１３６を囲むように境界線（図５参照）を設定する（ステップＳ１０６）。以下、ステップＳ１０４、Ｓ１０６について説明する。

図５は、穴１３６から流出する流れを速度ベクトルＶで示す図である。なお図中の領域Ｅは、いわゆる斜角補正が施された状態で出力部１４０に表示されたものである。斜角補正では、ベクトル解析部１４４が補正係数格納部１４６から斜角補正係数を読み出し、これを用いて撮像部１２８からの距離が遠いほど速度ベクトルを伸ばし、距離が近いほど速度ベクトルを縮める補正を行う。

領域Ｅは、マトリクス状のプロット１５０と、穴１３６を囲むように設定された矩形（コの字型）の境界線１５２とを含む。このように、プロット１５０がマトリクス状に配列されていて、境界線１５２が矩形であるので、境界線１５２上にプロットが設定し易くなる。ここで、ステップＳ１０４で算出される速度分布とは、図中の速度ベクトルＶのうち境界線１５２上のプロットを通過する流れの速度成分である（図６参照）。

境界線１５２は、図６（ａ）に示すように、その内側に最も速度が大きいプロットｆを含むように設定する。なおプロットｆを通過する流れは、速度ベクトルＶｆとして示している。ここで、レーザーシート１３０は空間内の１つの断面しか照射しないので、境界線１５２は基本的にはなるべく穴１３６に近いほうが好ましい。しかしながら、穴１３６にあまりに近い位置のプロットでは、穴１３６に流体が流入する場合には流れに乱れが大きく、穴１３６から流体が流出する場合にはトレーサー１３２が穴１３６に吸い込まれてしまうために、適切な速度を得られない可能性がある。

そこで、境界線１５２の内側に最も速度が大きいプロット（ここでは、プロットｆ）を含むように設定することにより、穴１３６から適切に距離をおいて境界線１５２を設定できる。これにより、壁面１３４の穴１３６から流出または流入する流れを最もよく捉えることができる。

図６は、図５の速度ベクトルＶを解析する手法を示す図である。図６（ａ）は、境界線１５２上のプロットの速度成分を示す図である。図６（ｂ）は、図６（ａ）のプロットの速度成分を用いて流れの流量を算出する過程を示す概念図である。

矩形の境界線１５２は、図６（ａ）に示すように、２つの側線１５４ａ、１５４ｂと、上線１５６とを含む。側線１５４ａ、１５４ｂは、壁面１３４に直交し穴１３６の中心位置を通る中心軸Ｇを基準として互いに線対称となっている。上線１５６は、壁面１３４に対向し側線１５４ａ、１５４ｂの端同士をつないでいる。図中では、境界線１５４の側線１５４ａにプロットａ、ｂ、ｃが設定され、上線１５６にプロットｄ、ｅが設定された状態を例示している。

ここで側線１５４ａ上のプロットａ、ｂを通過する流れは、速度ベクトルＶａ、Ｖｂとして示している。また、プロットａ、ｂの速度成分とは、速度ベクトルＶａ、Ｖｂの側線１５４ａに対するＸ方向に沿った垂直成分であり、図６（ａ）ではＶａｘ、Ｖｂｘとして示している。上線１５６上のプロットｄ、ｅを通過する流れは、速度ベクトルＶｄ、Ｖｅとして示している。また、プロットｄ、ｅの速度成分とは、速度ベクトルＶｄ、Ｖｅの上線１５６に対するＹ方向に沿った垂直成分であり、図６（ａ）ではＶｄｙ、Ｖｅｙとして示している。

再び図３に戻り説明する。ベクトル解析部１４４は、境界線１５２上のプロットの速度成分に単位面積をかけて単位流量を算出し、単位流量を積分して流れの流量を算出する（ステップＳ１０８）。以下、ステップＳ１０８について説明する。

図６（ｂ）には、一例として、領域Ｈに速度成分Ｖｍ、Ｖｎおよび単位面積ΔＳ、ΔＲを示している。なお領域Ｈは、境界線１５２で囲まれた領域（境界線１５２を穴１３６の中心軸Ｇを中心として回転させたときに囲まれる領域）とされる。速度成分Ｖｍは、プロットａ、ｂの速度成分Ｖａｘ、Ｖｂｘと方向が同一で大きさが（Ｖａｘ＋Ｖｂｘ）／２となるベクトルである。速度成分Ｖｎは、プロットｄ、ｅの速度成分Ｖｄｙ、Ｖｅｙと方向が同一で大きさが（Ｖｄｙ＋Ｖｅｙ）／２となるベクトルである。

単位面積ΔＳ、ΔＲは、領域Ｈを境界線１５２上のプロット（ここでは、プロットａ、ｂ、ｃ、ｄ、ｅ）で分割することでそれぞれ得られる。例えば、単位面積ΔＳは、Ｙ方向に関し側線１５４ａ上のプロットａ、ｂ、ｃにより４つに分割されている。単位面積ΔＳは、円周方向に関しては、適宜の数で分割してよく、分割数が増えるほど流量の算出の精度を高めることができる。また、上線１５６での流れの流量が穴１３２の中心軸Ｇ周りに対称であると仮定することで、単位面積ΔＲは、中心軸Ｇを中心とするプロットｄ、ｅで囲まれた同心円の面積となる。

つまり、ステップＳ１１０では、領域Ｈに示すように、速度成分Ｖｍに単位面積ΔＳをかけることで単位流量を算出し、これを中心軸Ｇ周りに積分することで、領域Ｈの側面での流れの流量を算出できる。なお側線１５４ａ、１５４ｂは中心軸Ｇに対して線対称であるため、側線１５４ａ、１５４ｂでの流れの流量は同一と仮定される。また領域Ｈに示す速度成分Ｖｎに単位面積ΔＲをかけることで単位流量を算出し、これを積分することで、領域Ｈの上面での流れの流量を算出できる。さらに、ベクトル解析部１４４は、ステップＳ１０８において、側線１５４ａ、１５４ｂおよび上線１５６での流れの流量を加算することで、壁面１３４の穴１３６から流出または流入する流れの流量を算出できる。

つぎに、ベクトル解析部１４４は、補正係数格納部１４６からトレーサー補正係数を読み出し、この補正係数を用いて流れの流量を補正する（ステップＳ１１０）。ここで、流れの中にあるトレーサー１３２は、流れに追従する理想的なトレーサーに比べて追従性が劣る可能性が高い。一例としてトレーサー１３２は、理想的なトレーサーに比べて鉛直方向の流れに対し、その速度が約２／３程度となる場合もある。

そこで、本実施形態では、流れの中にある水酸化鉄などをトレーサー１３２として利用する場合、理想的なトレーサーとの追従性の差を考慮したトレーサー補正係数を予め算出し、これを補正係数格納部１４６に格納している。そして、ステップＳ１１０において、ベクトル解析部１４４は、トレーサー補正係数による流量の補正を行う。なおステップＳ１１０では、ステップＳ１０８で算出された流れの流量を補正してもよく、あるいは、ステップＳ１０４で得られた境界線１５２上の各プロットの速度ベクトル自体を予め補正してもよい。

したがって、本実施形態によれば、原子炉格納容器１０２などの開放空間内であっても、壁面１３４の穴１３６を囲むように境界線１５２を設定し、ベクトル解析を行うことで、穴１３６から流出または流入する流れの流量を算出できる。さらに、トレーサー補正係数による流量の補正を行うことで、流れの流量をより正確に算出できる。

上記実施形態では、境界線１５２は、矩形としたがこれに限られず、境界線１５２上に各プロットが設定できるのであれば、多角形や円形などの適宜の形状であってもよい。また、上記実施形態では、原子炉格納容器１０２の穴１３６から流出する流れの流量を算出する場合を例示したが、これに限られず、穴１３６から原子炉格納容器１０２に流入する流れの流量を算出してもよい。なお上記実施形態においては原子炉格納容器１０２の内部を例に用いて説明したが、全く同様にして原子炉圧力容器１０４の内部の穴から流出または流入する流れの流量を算出することも可能である。さらに、本実施形態にかかる流量算出方法は、開放空間を有する原子炉以外の適宜の構造物に適用してもよく、それらの構造物の壁面の穴から流出または流入する流れの流量を算出してもよい。

以上、添付図面を参照しながら本発明の好適な実施形態について説明したが、本発明は係る例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。

本発明は、ＰＩＶ計測システムを用いた流量算出方法に利用することができる。

１００…原子炉建屋、１０２…原子炉格納容器、１０４…原子炉圧力容器、１０６…サプレッションプール、１０８…燃料集合体、１１０、１１２…冷却水、１１４…デブリ、１１６…検査孔、１１８…流量測定装置、１２０…挿入管、１２２…ヘッド、１２４…治具、１２６…光学ユニット、１２８…撮像部、１３０…レーザーシート、１３２…トレーサー、１３４…壁面、１３６…穴、１３８…流量算出部、１４０…出力部、１４２…画像処理部、１４４…ベクトル解析部、１４６…補正係数格納部、１４８…照射ライン、１５０…プロット、１５２…境界線、１５４ａ、１５４ｂ…側線、１５６…上線

Claims (4)

1. 開放空間の壁面の穴から流出または流入する流れの流量を算出する方法であって、
   前記穴の位置にレーザーシートを形成し、
   流れの中にある粒子に前記レーザーシートの光が反射した反射光を連続的に撮影することにより流れの速度分布を算出し、
   前記穴を囲むように境界線を設定し、
   前記境界線上のプロットの速度成分に単位面積をかけて単位流量を算出し、該単位流量を積分して前記流れの流量を算出することを特徴とする流量算出方法。
2. 前記速度分布は、マトリクス状のプロットについて算出されていて、
   前記境界線は、前記壁面に直交し前記穴の中心位置を通る中心軸を基準として互いに線対称となる２つの側線と、前記壁面に対向し該側線の端同士をつなぐ上線とを含む矩形であることを特徴とする請求項１に記載の流量算出方法。
3. 前記境界線は、その内側に最も速度が大きいプロットを含むように設定することを特徴とする請求項１または２に記載の流量算出方法。
4. 流れの中にある前記粒子と理想的なトレーサーとを比較して補正係数を予め算出し、
   前記境界線上のプロットの速度に前記補正係数をかけて、該速度を補正することを特徴とする請求項１から３のいずれか１項に記載の流量算出方法。

Images