JP3685256B2

JP3685256B2 - Fluid flow measurement system

Info

Publication number: JP3685256B2
Application number: JP2001280508A
Authority: JP
Inventors: 治嗣森; 英昭手塚; 公一菱田; 昌信前田
Original assignee: Tokyo Electric Power Co Inc; Keio University
Current assignee: Tokyo Electric Power Co Inc; Keio University
Priority date: 2001-09-14
Filing date: 2001-09-14
Publication date: 2005-08-17
Anticipated expiration: 2021-09-14
Also published as: JP2003084005A

Description

【０００１】
【発明の属する技術分野】
本発明は、複雑な流動場における流体の流動を高精度かつ精密に測定する流体の流動計測技術に係り、特に閉空間内を流動する流体の流速および流れ方向を計測する流体の流動計測システムに関する。
【０００２】
【従来の技術】
複雑な流動場における流体の流動を高精度かつ精密に測定することは、原子力発電プラントや火力発電プラント、化学プラント等におけるプラント挙動予測の高精度化、プラント診断、プラント性能評価や機器改良にとって極めて重要である。実機レベルのプラント、特に外界と環境を異にする厳しい環境下での流体の流動計測技術は、複雑な流れをもつ流体流動計測条件のために、流動計測が極めて困難な状況にある。
【０００３】
中でも、原子炉圧力容器や火力発電プラントの熱交換器内のように、厳しい環境下にある熱流動場における流体の流速や流れ方向等の流体流動測定は、アクセス性を含め複雑な流体流量計測条件のために極めて困難である。
【０００４】
しかしながら、近年熱流動場における可視化技術の研究が進み、複雑な流動場における流体の流動を高精度かつ精密に測定可能な粒子画像流速計（ＰａｒｔｉｃｌｅＩｍａｇｅＶｅｌｏｃｉｍｅｔｒｙ：以下、ＰＩＶという。）が開発されつつある。このＰＩＶと光ファイバを束ねたイメージガイドを併用することにより、原子炉圧力容器内のように、外界と環境を異にする閉空間における流体の流動測定が期待できる。
【０００５】
【発明が解決しようとする課題】
しかし、イメージガイドを通してＣＣＤカメラによって撮像された粒子画像には、粒子の輝度パターンだけでなく、イメージガイドの光ファイバの配列（俵積みの構造）の輝度パターンも存在する。光ファイバの配列輝度パターンが相互相関法による解析の際に弊害となり、粒子の輝度パターンに対して誤対応（過誤ベクトル）を招く。
【０００６】
また、最小の空間分解能を与える光ファイバの単繊維一本よりも粒子像径を大きくすると、参照窓画像内に最低でも粒子は６〜８個あることが望ましいとされるため、参照窓画像を大きく取らなくてはならなくなり、計算量が莫大に増加してアルゴリズムとしては現実的ではないものとなる（探索する領域が増えるとその４乗に比例して計算量は増加する）。
【０００７】
本発明は、上述した事情を考慮してなされたもので、イメージガイドを利用した粒子軌跡追跡法（Ｐａｒｔｉｃｌｅｓｔｒｅａｋｔｒａｃｋｉｎｇｖｅｌｏｃｉｍｅｔｒｙ）により画像処理を行って流体の流動状態を検出するアルゴリズムの最適化を行い、流体の流動分布を精度よく正確に行うことができる流体の流動計測システムを提供することにある。
【０００８】
【課題を解決するための手段】
本発明に係る流体の流動計測システムは、上述した課題を解決するために、請求項１に記載したように、レーザ光を発振させるレーザ発振装置と、発振されたレーザ光を流体の流動場内にシート状に投入させるレーザシート形成用走査光学系と、この走査光学系からのレーザシート上の２次元粒子軌跡画像を撮像する画像撮像手段と、前記レーザ発振装置と画像撮像手段をタイミングをとって同期化駆動させるタイミングコントロール手段と、前記画像撮像手段で撮像された２時刻の粒子軌跡画像の輝度パターンを比較・解析し、個々の粒子の移動方向および移動量を計測する画像処理手段を備え、この画像処理手段による個々の粒子の移動量から流体の流動・流速分布を計測するように構成し、前記画像撮像手段は、走査光学系からのレーザシート上の２次元粒子軌跡画像を光伝送する画像伝送手段と、光伝送された２次元粒子軌跡画像を撮像する撮像手段とを備えたものであり、さらに、前記画像伝送手段は、多数本の光ファイバを束ねて一体化させ、両端面が平面加工された可撓性のイメージガイドで構成され、前記イメージガイドは、対物レンズにより一方のファイバ端面に結像された画像を光ファイバで各画素に分解して撮像手段側の他端面まで同一画像を伝送するように構成し、前記レーザ発振装置の発光時間を調整して粒子軌跡を前記イメージガイドの単繊維径よりも大きくしたものである。
【０００９】
また、上述した課題を解決するために、本発明の流体の流動計測システムにおいては、請求項２に記載したように、前記走査光学系と画像撮像手段は、一体的に組み合わせて構成され、流体の流動場で１次元あるいは２次元移動可能に構成されたものである。
【００１０】
さらに、上述した課題を解決するために、本発明の流体の流動計測システムにおいては、請求項３に記載したように、前記走査光学系は、直交するレーザシートを流体の流動場に形成し、この直交するレーザシート上の２次元粒子軌跡画像を光伝送する画像伝送手段をそれぞれ設けたものである。
【００１１】
【発明の実施の形態】
本発明に係る流体の流動計測システムの一実施形態について添付図面を参照して説明する。
【００１２】
図１は、本発明に係る流体の流動計測システムの構成を示す原理図である。
【００１３】
本発明の流体の流動計測システム１０には、高出力のレーザ光を発振させるレーザ発振装置１１が光源として備えられる。レーザ発振装置１１には小型で高出力レーザ光を発振する半導体レーザやパルスレーザ等が用いられる。例えば、出力１６Ｗの連続光で中心発振波長８０８ｎｍの半導体レーザが光源として用いられる。
【００１４】
レーザ発振装置１１から発振されたレーザ光は、送光用ファイバ１２を通ってレーザシート形成用走査光学系１３に案内される。走査光学系１３はレーザ発振装置１１から発振されたレーザ光を流体の流動場１４にシート状に投射させるようになっている。
【００１５】
流体の流動場１４は、外界と環境を異にする閉空間内に形成される。また、場合に応じて、外界に開放された形態で流体の流動場１４が形成されていてもよい。流動場１４は、気体、液体のうちの少なくともいずれかの流体からなる。具体的には、例えば、原子炉圧力容器のダウンカマ部、炉心シュラウド、火力発電プラントの熱交換器や蒸気発生器内等の閉空間における流体の速度場を対象とする。また、船舶、航空機あるいは自動車等のような、水流や気流と接触する構造部位や機材まわりにおける流体の流動場も対象とすることができる。なお、図１に示す原理図では、流体の流動場１４として水の流動場の実験例を示している。容器１５あるいは配管内にノズル１６より水を噴出させて流体としての水の流動場１４を形成している。
【００１６】
走査光学系１３は、流体の流動場１４にレーザ発振装置１１からの発振レーザ光をシート状に照射して流体流動場１４にレーザシート１７を形成している。このレーザシート１７により流体の流動場１４を可視化している。流動計測システム１０では、流体流動場１４の奥行き方向の空間分解能を確保するために、発振レーザ光をシート状に照射している。
【００１７】
走査光学系１３は、図２に示すようにレーザシート１７の作成のためのレーザ光照射光学系である。レーザ発振装置１１として例えば半導体レーザを用いると、半導体レーザからの発振レーザ光は、送光用ファイバ１２から大きな広がり角を持って出射される。出射されたレーザビームは軸対称の光となっておらず、非点較差があり、レーザビームの形状変換は複雑となる。
【００１８】
このため、走査光学系１３には、送光用ファイバ１２から出射されたレーザビームからレーザシート１７を形成するためのレーザシート形成用光学手段としてのレンズ群１８が用いられる。このレンズ群１８は、出射レーザビームを並行レーザビームにコリメートさせる２枚のコリメートレンズ１９ａ、１９ｂと、コリメートレンズ１９ａ、１９ｂによって並行光となったレーザビームをレーザシート１７の厚み方向で約１〜２ｍｍ程度にウエストさせる長焦点レンズ２０と、レーザシート１７を形成するシリンドリカルレンズ２１とを組み合わせてレーザシート形成用光学手段が構成される。コリメートレンズ１９ａ、１９ｂには、例えば、平凸レンズが用いられる。
【００１９】
走査光学系１３からシート状に照射されるレーザ光により、流体の流動場１４にレーザシート１７が形成され、流体の流動場１４がシート状に可視化される。そして、レーザシート１７の測定範囲１７ａに対向して画像伝送手段２２が設けられる。図１に示したように、画像伝送手段２２に近い容器１５壁には、レーザシート１７上の画像が画像撮像手段２９で撮像可能となるように、透過用窓ガラス２３が設置されている。
【００２０】
画像撮像手段２９は、レーザシート１７上の画像を伝送する画像伝送手段２２と、２次元粒子軌跡画像を撮像する撮像手段としてのＣＣＤカメラ２４とを有する。画像伝送手段２２はレーザシート１７上の２次元粒子軌跡画像を撮像手段としてのＣＣＤカメラ２４に伝送させるようになっている。画像伝送手段２２は、図３に示すように、数千本から数万本の光ファイバ２５を束ねて一体化（溶融）させた可撓性のイメージガイド２６が用いられる。イメージガイド２６は、図３に示すように、束ねられた光ファイバ２５の両端面の各位置が正確に対応するように並設されて束ねられる。イメージガイド２６の両端面は平面に仕上げられ、対物レンズ２７により一方のファイバ端面に結像された画像を各画素ファイバとしての光ファイバ２５に分解し、同一画像をＣＣＤカメラ２４側のファイバ端面まで伝送している。ファイバ端面からカメラレンズ２８を経てＣＣＤカメラに記録される。画像撮像手段２９はイメージガイドを備えずに、レーザシート１７上の２次元粒子軌跡画像を直接撮像できるようにしてもよい。
【００２１】
イメージガイド２４の画像の限界（解像度）は、単一の光ファイバ２５の直径と配列の仕方によって決まり、像の見え方は、イメージガイド２６に束ねられる光ファイバ２５の総本数と並べ方によって決まる。光ファイバ２５のコア径より小さな光量の光分布は、光ファイバ２５を伝わる間に平均化され、各コア断面内では明るさは一様となる。したがって、イメージガイド２６は光ファイバ２５のコア径が画素（ピクセル）としての空間分解能の最小単位となる。
【００２２】
実際のイメージガイド２６は、多数本の光ファイバ２５が俵積み（六方稠密配列）構造で整然とかつ稠密に配列されている。このイメージガイド２６で光ファイバ間隔とほぼ等しい間隔の格子状物体を観察すると、物体の像の強度分布と無関係なモアレ干渉縞が生じる。
【００２３】
また、光ファイバ２５間にクラッド部分が融着されているイメージガイド２６では、クラッドの厚さが薄くなると、１つの光ファイバ２５のコアに入った光が隣のコアに漏れてしまう場合が生じる。このため、像の境界が滲み、像のコントラストが低下する。このため、イメージガイド２６は、六方稠密配列構造を採用したとき、クラッドの厚さが薄過ぎないことが条件となる。
【００２４】
一方、流体の流動計測システム１０には、タイミングコントロール手段３０が備えられている。タイミングコントロール手段３０はタイミングスケジューラ３１とシンクロナイザ３２とから構成される。タイミングスケジューラ３１からのタイミングでレーザ発振装置１１が発振せしめられる一方、この発振タイミングとシンクロナイザ３２で同期をとってＣＣＤカメラ２４が駆動せしめられる。すなわち、タイミングコントロール手段３０は、レーザ発振装置１１と操作手段としてのＣＣＤカメラ２２との同期をとって駆動させるようになっている。
【００２５】
ＣＣＤカメラには、例えば、近赤外線領域に感度を持つ解像度ＶＧＡ（６４０×４８０ｐｉｘｅｌｓ）、フレームレート３０Ｈｚ、８ビット（モノクロ２５６階調）のものが用いられる。ＣＣＤカメラ２４は、レーザシート１７上の２次元粒子軌跡画像をレーザ発振装置１１からの発振タイミングと同期をとって撮像しており、ＣＣＤカメラ２４で撮像されたアナログの画像信号は画像処理手段３５に送られて画像処理され、後述する粒子軌跡追跡法により流体の流動場１４における流体の流速分布や流れ方向が計測される。ＣＣＤカメラ２４には、例えば、逐次読出方式（プログレッシブスキャン）のカメラが用いられ、配列されたＣＣＤの奇数ライン、偶数ラインの電荷蓄積時間は同時に行われる。このＣＣＤカメラ２４は通常のカメラに較べ垂直方向の解像度に優れている。
【００２６】
また、画像処理手段３５は、ＣＣＤカメラ２４からのアナログ画像信号をデジタル化するＡ／Ｄ変換器としてのフレームグラバボード３６と、このフレームグラバボード３６からデジタル画像信号がＰＣＩバス等の信号線３７を介して送られ、デジタル画像処理するコンピュータ３８とを有する。
【００２７】
この流体の流動計測システム１０は、ＣＣＤカメラ２４にてアナログ撮影し、画像処理手段３５によりデジタル処理し、デジタル録画するアナログ撮影・デジタル録画方式を採用している。
【００２８】
図１に示す流体の流動計測システム１０は実験室対応の例を示した。実用化に際しては、図４および図５に示すように、レーザシート用走査光学系１３と画像撮像手段２９は一体化される。図４では、走査光学系１３と一体化した画像撮像手段２９の画像伝送手段２２は側方視可能な可視化用プローブ４０として構成されている。プローブ４０先端に設置された光学窓（透過用窓ガラス）４１から４５°直角プリズム４３により光路を曲げ画像伝送手段を構成したファイバースコープ４４にて画像伝送される。図４（Ａ）は、この一体化した流体の流動計測システムを図１の容器１５内における流体としての水の流動場に用いた例を示す。図５では、一体化した走査光学系１３と画像撮像手段２９の画像伝送手段２２とを本体ケーシング４５に格納している。
【００２９】
次に、流体の流動計測方法を説明する。
【００３０】
図１に示された流体の流動計測システム１０では、照明に用いる半導体レーザの発光時間を比較的長くとり、ＣＣＤカメラ２４に粒子の流跡（軌跡）を撮像する。そして、この撮像された時間的に微小時間異なる２時刻の粒子軌跡画像を輝度パターンとみなして画像処理手段３５のコンピュータ３８により比較・解析の画像処理を行い、個々の粒子の移動量を推定する。この流体の流動計測装置１０においては、この手法を粒子軌跡追跡法（Ｐａｒｔｉｃｌｅｓｔｒｅａｋｔｒａｃｋｉｎｇｖｅｌｏｃｉｍｅｔｒｙ）と呼ぶ。この粒子軌跡追跡法は、粒子の流跡を記録することにより、必要以上に粒子像径を大きくすることなく、かつ最小の空間分解能となる光ファイバ単繊維に対して十分な粒子像の情報を与えることができる。
【００３１】
粒子軌跡画像は、レーザシート１７上の水（流体）中に浮遊し、拡散分布する個々の粒子（浮遊物等の異物）のある時刻での流跡の画像を表しており、デジタル化された輝度を持つピクセル（画素）によって構成される。画像処理手段３５による画像処理の過程では、各粒子軌跡画像のうち限定された矩形領域１７ａ（図１参照）を取り出して画像処理される。デジタル化された粒子軌跡画像中のある１点の値を「輝度値」と定義し、この輝度値が矩形領域１７ａの領域に分布しているものを「輝度パターン」と称する。
【００３２】
粒子軌跡画像から水（流体）中の個々の粒子の移動量を推定するアルゴリズムは、得られた粒子軌跡画像を輝度分布の変化であると判断し、レーザシート１７上の粒子の輝度パターンがある時間内に移動する量を個々の粒子について定量的に解析するものである。この粒子の移動量の解析には相互相関法を用いる。
【００３３】
そこで、本発明に係る流体の流動計測システム１０においては、粒子軌跡追跡法による流体の流動測定の最適化を図るために、相互相関計算の妨げとなる光ファイバ２５の背景処理を行う。
【００３４】
特願２０００−２０７４７６号公報に係る流体の流動計測システムおよびその計測方法では、粒子画像流速計測（ＰＩＶ）法を適用したものであり、このＰＩＶ法では、凍結画像が用いられる。通常のＰＩＶ法では、粒子像径がイメージガイド２６の光ファイバ２５単繊維一本より小さいと、もはや失われた粒子の情報は回復できるものではないので、算出される相関係数は俵積みの輝度パターンに強く影響を受け、速度検出に非常に大きな弊害となり、誤対応を招く。
【００３５】
凍結画像を用いるＰＩＶ法では、要求されるパラメータ（粒子像径、数密度）は非常に制限を受ける。そこで、本発明で用いられる粒子軌跡追跡法では、照明に用いる半導体レーザの発光時間を比較的長くとり、粒子の流跡を記録、粒子群ではなく個々に粒子の移動量を見積もるアルゴリズムを構築し、ソフトウェアに実装した。個々の粒子の流跡を記録することにより、必要以上に粒子像径を大きくとることなく、かつ最小の空間分解能となる光ファイバ２５の単繊維に対して十分な粒子像の情報を与えることができる。この粒子軌跡追跡法（Ｐａｒｔｉｃｌｅｓｔｒｅａｋｔｒａｃｋｉｎｇｖｅｌｏｃｉｍｅｔｒｙ）の手法は、レーザの発光時間を容易に制御することができる半導体レーザを用いることによって可能となる。
【００３６】
本発明で用いる粒子軌跡追跡法では、この背景処理も含めて、主に以下の６段階の手順を経て粒子の移動量の算出が行われる。
【００３７】
（１）ノイズの除去（背景処理）
イメージガイド２６を通してＣＣＤカメラ２４で撮像された粒子軌跡画像を処理し、時系列に得られた粒子軌跡画像を各ピクセル毎に輝度値の加算平均値を求め、各々の粒子軌跡画像の輝度値から加算平均値を引くことで、光ファイバ２５のファイバ配列の輝度パターン除去を行う。
【００３８】
光ファイバ配列の輝度パターン除去の画像処理アルゴリズムを以下に説明する。
【００３９】
測定データの粒子軌跡画像をｘとするとき、粒子軌跡画像ｘは、
【数１】
ｘ＝ｓ＋ｎ ……（１）
但し、ｓは信号成分（光ファイバ配列の輝度パターン）
ｎは雑音成分（粒子群の輝度パターン）
で表される。
【００４０】
すなわち、測定データの粒子軌跡画像ｘには、信号（光ファイバ配列の輝度パターン）に雑音（粒子の輝度パターン）が重なった画像を測定していると考える。加算平均を取ることによる光ファイバ配列の輝度パターン除去処理は、ＳＮ比が小さく、信号と雑音の構成周波数にも大差がない場合や同じ条件で何度も測定を繰返すことが可能な場合に有効な手法である。
【００４１】
式１を時系列要素を負荷して表示すると、式２のように表される。
【００４２】
【数２】
ｘ^ｉ（ｋ）＝ｓ^ｉ（ｋ）＋ｎ^ｉ（ｋ） ……（２）
ｉ＝ｉ枚目の粒子画像、
ｋ＝ｉ枚目の粒子画像内のｋ番目のピクセル輝度値、
をそれぞれ示す。
【００４３】
粒子軌跡画像の撮影において、取り込んだ粒子軌跡画像の枚数をＭ枚とすれば、Ｍ枚の粒子軌跡画像に対する加算平均ｘ（ｋ）は次のように表せる。
【００４４】
【数３】

【００４５】
光ファイバ配列の輝度パターンｓは、同一の実験条件（送光系・受光系）であれば、連続する粒子軌跡画像においては同じパターンが出現するため、Ｍ回の加算によりＭ倍の輝度値となり、これを平均する（Ｍで除する）加算平均処理を行うと、もとの輝度値のままである。すなわち、加算平均処理により、粒子の存在しない光ファイバ配列のみ映った画像を得ることができ、光ファイバ配列の存在する粒子軌跡画像から各画素毎に光ファイバ配列の輝度値（輝度パターン）を引けば、光ファイバ配列の輝度パターンの除去が可能となる。
【００４６】
この背景処理によって、取得した画像の背景輝度は画面全体にわたってほぼ均一になるので、相関計算処理に対して有効だけでなく、二値化処理においても有効な手段となる。
【００４７】
（２）二値化・膨張化
画像処理粒子軌跡画像から粒子の識別を行うための前処理として、二値化および膨張化画像処理を行う。水（流体）中に浮遊する粒子の軌跡像は、光ファイバ２５のクラッドの影響を受け、連続的なものとはならず、断続的な画像として記録される（図６（ａ）参照）。この断続的粒子画像を二値化処理しても、クラッドによって情報の失われた部分は回復されない（図６（ｂ）参照）。これらの分断された粒子の軌跡像を同一の粒子の移動軌跡として認識させるために、膨張化画像処理が行われる。膨張化画像処理とは、注目の画素が存在すれば、その周囲のすべての画素を０から２５５に変換するものである（８連結の膨張処理）。本アルゴリズムでは、この膨張化処理を２回繰返すことで分断された断続的粒子画像は、連続性のある軌跡画像として回復される（図６（ｃ）参照）。
【００４８】
（３）ラベリングおよび粒子画像位置の検出
１枚の粒子軌跡画像には複数（多数）の粒子画像が存在するので、それぞれを個々の粒子として認識させるラベリングを行う。それぞれの連結した画素には同じ番号を、異なった連結成分には別の番号をつける。このような個々の領域に番号付けを行うことをラベリングと呼ぶ。ラベリングによって識別された粒子画像は、元画像に戻って中心位置が算出される。中心座標は、以下の重心法（Ｐａｒｔｉｃｌｅｉｍａｇｅｃｅｎｔｒｏｉｄ）によって求められる。
【００４９】
Ｘ_０、Ｙ_０を求める重心、ｉ、ｊを各点の座標、Ｉをその座標の輝度値とすると、重心Ｘ_０、Ｙ_０は、
【数４】

で表される。Ｙ_０についても同様である。
【００５０】
（４）相関計算による整数移動量の算出
（１）〜（３）の粒子画像前処理によって、粒子軌跡画像の中心位置が検出されたので、個々の粒子の移動量を相互相関法により求める。
【００５１】
なお、流体の流動計測システム１０のＣＣＤカメラ２４で撮像された第１時刻における２次元粒子軌跡画像をデジタル処理した粒子軌跡画像を「参照画像」と、「参照画像」のうち限定された矩形領域を「参照窓画像」と、第１時刻とは微小時間異なる第２時刻における２次元粒子軌跡画像をデジタル処理した粒子軌跡画像を「探索画像」と、「探索画像」のうち限定された矩形領域を「探索窓画像」とそれぞれ定義する。
【００５２】
相互相関法による粒子の移動量推定は、「探索窓画像」と「参照窓画像」との間で次式で表される相関値Ｒをもつことが知られており、この相関値Ｒをもって個々の粒子の輝度パターンの類似度を求め、比較検討して２画像間の個々の粒子の移動量を求める手法である。
【００５３】
【数５】

【００５４】
【数６】

【００５５】
式５および式６に示すように、相関値Ｒは、探索窓画像および参照窓画像の輝度値の二乗平均で除する（正規化）ことによって、−１から１までの値をとる。相関値Ｒが１のときに２画像は完全に一致し、相関値Ｒの値が大きいほど窓画像間の類似度が大きいことを意味する。
【００５６】
式５および式６において、Ｉ_１およびＩ_２は参照窓画像と探索窓画像の各画素の輝度を表し、ξ、ηは探索窓画像と参照窓画像の相対的な位置で表す。相関値が最大となる位置ξ、ηが画像内での粒子の移動量ΔＸ’、ΔＹ’に相当する。参照窓画像の大きさｎ×ｍおよび参照窓画像との相対的な位置は、予測される最小および最大速度から決定される。
【００５７】
（５）サブピクセル移動量の算出
（４）の相関計算により、粒子画像の移動量の整数値が求まるので、相関係数のピークの近傍３点を用いた正規分布近似（Ｇａｕｓｓｉａｎｐｅａｋ−ｆｉｔ）により、サブピクセル移動量を見積もる。ｘ、ｙ方向についてそれぞれ独立に算出される。
【００５８】
（６）較正
以上の手順により、粒子軌跡画像の移動量がピクセル（ｐｉｘｅｌ）単位で求まるので、較正により実空間での移動量に変換される。
【００５９】
図７にトレーサー粒子の混入したグリセリン溶液を回転円盤上に設置し、擬似的な流れを撮影した粒子軌跡画像を、図８にこの粒子軌跡画像６０フレーム分を粒子軌跡追跡法により画像処理して個々の粒子の移動量を算出した結果を示す。
【００６０】
【実施例】
以下に、粒子軌跡追跡法による画像処理を行う本発明に係る流体の流動計測システムの性能の評価、および流動計測システムの実流動場への適用性について説明する。
【００６１】
【実施例１】
（微動台による移動量の検定）
平板上に描かれた粒子画像を微動台（メカニカルステージ）によって移動させ、その移動量を７０ｍｍ×７０ｍｍの較正板を使って測定し、実際に移動させた真の値と、本発明の流動計測システムによりＳｔｒｅａｋ法により計測した測定値とを比較し、誤差を実測してその精度の検定を行った。
【００６２】
較正板は、７０ｍｍ×７０ｍｍの矩形内に６４点の格子点を有するものを用いた。したがって、計測範囲は７０ｍｍ×７０ｍｍである。このときの撮像倍率は、５．０１１ｍｍ／ｐｉｘｅｌであった。
【００６３】
粒子は粒子径３、６、８ｐｉｘｅｌｓの３種類を用いた。粒子像は、流体中にシーディングされた粒子をレーザシート照明によって撮像し、これをプリンタによってプリントアウトして平板上に接着させたものである。
【００６４】
粒子像の撮像は、基準位置（較正における基準位置）に設置された平板上の粒子像と微動台による移動後の粒子像とをＣＣＤカメラにより撮像し、これを一つの画像対とした。画像対から移動量ベクトルの算出は相互相関法によるものであり、実流動場の計測に用いられる方法と同等のものである。一組の画像対から算出されるベクトル数は約１５０個であった。
【００６５】
粒子像の描かれた平板の移動は、ｘ方向に５０μｍづつ１０ｍｍまで行い、２００枚の画像を取得し、１００組の画像対を得た。移動量算出は、相互相関法および正規分布近似（Ｇａｕｓｓｉａｎｐｅａｋ−ｆｉｔ）により行った。
【００６６】
最小の空間分解能となるイメージガイド２６の光ファイバ単繊維径は、イメージガイド２６の周波数解析結果からｄｆ＝３．１２ｐｉｘｅｌｓであるので、３種類の粒子像径３、６、８ｐｉｘｅｌｓに対して、それぞれｄτ／ｄｆ＝０．９６、１．９２、２．５６倍となる。
【００６７】
図９、図１０および図１１にそれぞれ、Ｓｔｒｅａｋ法（粒子軌跡追跡法）による計測結果の偏り誤差（Ｂｉａｓｅｒｒｏｒ）、ばらつき誤差（ＲＭＳｅｒｒｏｒ）および確率密度分布（ＰＤＦ）を示す。各図において、ｄτ／ｄｆ＝０．９６の曲線をａ、ｄτ／ｄｆ＝１．９２の曲線をｂ、ｄτ／ｄｆ＝２．５６の曲線をｃとし、図中の横軸は真の移動量（Ａｃｔｕａｌｄｉｓｐｌａｃｅｍｅｎｔ）［ｐｉｘ．］を示す。図１１に示されたＰＤＦは、真の移動量がそれぞれ０．５ｐｉｘｅｌｓ、２．５ｐｉｘｅｌｓ、４．５ｐｉｘｅｌｓにおけるものである。
【００６８】
図９は、平板をｘ方向に移動させたときのｘ方向移動量の計測値の偏り誤差を示しており、真値の移動ピクセル数からどけだけはずれて測定したかをプロットしたもので、偏り誤差は計測範囲内の全ベクトルの平均値と真値との差を示している。
【００６９】
また、図１０は計測範囲内の計測値のルート二乗平均をとったものであり、図１１は、真の移動量が、左の山から０．５ピクセル、２．５ピクセル、４．５ピクセルの時、Ｓｔｒｅａｋ法（粒子追跡法：ＰＴＶ）による測定値が、全体の誤差数を１０（％）とした時、何個（％）、真値よりずれたところで実際計測されたか、示している。この山の幅が狭ければ狭いほど、誤差は少ないことを表している。例えば、真中の山の２．５ピクセル移動の時、２ピクセルと計測したのは、１００個のうち約９個ということを意味する。
【００７０】
図９ないし図１１に示されたように、粒子軌跡追跡法では、粒子径による計測値の相違はあまり現れず、偏り誤差βは約β＝±０．２ｐｉｘｅｌｓの範囲に収まった。
【００７１】
実験で用いた光ファイバ２５には光を通さない（透過率が比較的低い）クラッドがあるために、粒子像の欠落が原因となり、振動効果が見られる。図９に示された偏り誤差には、粒子像径に依る値の変化は見られず、これは粒子軌跡によって光ファイバ径に対して十分な情報が与えることができることによる。
【００７２】
図１０に示されたＳｔｒｅａｋ法によるばらつき誤差も粒子像径に依らずほぼ同一の値を示しているが、比較的大きな値となった。これは、確率密度分布を見ても分かるように、ヒストグラムの形状が正規分布とはならず、真値付近で水平に近い分布となることが原因している。これは粒子軌跡像の相関をとっているためであり、すなわち、粒子軌跡の輝度値は比較的なだらかであり、その画像の相関係数分布もまたなだらかなものとなるためである。粒子軌跡追跡法では、その粒子を通常の粒子画像に比べ大きくとっていることが、ばらつき誤差ε_ｋを大きくする（０．３〜０．５ｐｉｘｅｌｓ）。しかし、粒子像径によってＰＤＦが大きく割れてしまうような誤差の発生の仕方はしないことが示され、また、偏り誤差の精度は有効であることが確認された。
【００７３】
【実施例２】
（実流動場での計測）
図１に示すような実験容器１５内に水（流体）の流動場１４を設けて本発明の流動計測システム１０により軸対称噴流の計測を行い、実流動場計測における流体の流動計測システムの性能を検討した。
【００７４】
流体の流動計測システムの検討は、ドップラ式レーザ流速計（ＬＤＶ）による計測と比較して行った。
【００７５】
容器１５は、例えば縦２５０ｍｍ、横２５０ｍｍ、高さ２５０ｍｍの立方体容器であり、中央に矩形噴流出口（ノズル）１６を持たせた。作動流体には水道水を用い、これにトレーサー粒子として３〜３０μｍの球形ポリエチレン粒子を混入した。図１に示すように、実験設備では作動流体を、下部の貯溜タンク５０に一旦ためた後、ポンプ５１により上部のリザーバタンク５２へ供給した。リザーバタンク５２に供給された作動流体は一部を仕切量５２ａを介してオーバーフローさせることにより、水頭圧（ヘッド）が一定に保たれる。作動流体はリザーバタンク５２から流動調節バルブ５３、整流用のハニカムおよびメッシュ（図示せず）を経て整流化し、軸対称ノズル１６から設備容器１５内に噴出させて水（流体）の流動場を構成した。リザーバタンク５２内でオーバーフローした水は戻り管５４を経て貯溜タンク５０内に戻される。
【００７６】
一方、容器１５内のノズル１６は面積比が約１０：１に縮流され、安定した速度分布が得られ、ほぼ定常になる状態において流体（水）の流動計測実験を行った。水（流体）の流動場の座標系はノズル１６出口の中心を原点とし、噴流方向（高さ方向）をＺ、半径方向をｒ、周方向をθ、それぞれの方向の速度成分をｗ、ｕ、ｖとした。
【００７７】
ノズル１６からの作動流体の出口速度をｗ＝１６５ｍｍ／秒とし、このときの出口直径（Ｄ＝７ｍｍ）に対するレイノズル数はＲｅ＝１１５５であった。容器１５をオーバーフローする水も、受けタンク５６で一旦受けた後、戻り水配管５７を経て貯溜タンク５０に戻され、水（流体）の循環系５８を構成した。
【００７８】
走査光学系１３と画像伝送手段２２は、図４に示すような走査光学系１３と画像伝送手段２２とが一体化した図４に示すような可視化用プローブ４０を用いた。図４に示すように、可視化用プローブ４０を容器１５の上方から水（流体）の流動場１４に挿入し、レーザ光は上方から鉛直に投入され、対向するプローブ４０からファイバースコープ４４により粒子軌跡画像を撮影した。レーザシート１７と対向するプローブ４０の受光面との距離は約３０ｍｍである。計測範囲は較正板の格子点が存在する範囲である、例えば１０ｍｍ×１０ｍｍの矩形範囲であり、レーザシート１７の厚さは２ｍｍである。
【００７９】
撮影対象領域は、ノズル１６の直径Ｄに対して、水の噴流のポテンシャルコアが消滅し、発達領域といわれ、周囲流体との運動量交換が噴流軸中心にまでおよび、比較的なだらかな速度分布の得られるＺ＝６Ｄの領域に設定した。
【００８０】
撮像された粒子軌跡画像は画像処理手段３５で粒子軌跡追跡法（Ｓｔｒｅａｋ法）により画像処理を行った。
【００８１】
一方、ＬＤＶ（ドップラ式レーザ流速計）による計測もＺ＝６Ｄにおける噴流軸に垂直な断面において別個に行った。これらの計測は、流動場が定常であると仮定し流量を同一条件に調整して行った。
【００８２】
図１２にＬＤＶで測定された噴流出口条件（レイノルズ数：Ｒｅ＝１１５５）を示す。図１２には、平均速度（ＬＤＶ（ｍｅａｎ））と変動速度（ＬＤＶ（ｒｍｓ））の曲線ｄ，ｅを示す。また、図１３に粒子軌跡追跡法によるノズル口から高さ方向Ｚ＝６Ｄ付近の速度ベクトルマップを示す。
【００８３】
粒子軌跡画像では、前述の粒子軌跡追跡法（Ｓｔｒｅａｋ法）によって速度を算出した。図１４および図１５に、それぞれＺ＝６Ｄにおける軸方向平均速度分布ｆと、その変動速度分布ｇを示す。横軸は半径方向、縦軸は軸方向である。本システムによる計測とＬＤＶ計測は、得られた速度分布である軸方向平均速度分布ｈとその変動速度分布ｉよりこの領域を抽出し、統計処理を行って２つの計測法による結果の比較を行った。この流体の流動計測システムにおけるファイバースコープによる軌跡画像は、個々の粒子に対応して算出されること、また、比較的低数密度のシーディングであるため、一組の画像対から得られる速度データが少ない。このような統計量の少なさを補うため、時間方向に多量の画像を撮る必要があり、本実験では１６０００枚の画像から統計量を算出した。図１４および図１５に示す軸方向平均および変動速度分布ｆ，ｇは各点で統計量約１０００であった。ＬＤＶによる計測では、統計量３０００から平均・変動速度データを算出した。
【００８４】
図１５に示される軸方向変動速度分布ｇ，ｉを比較すると、Ｓｔｒｅａｋ法はＬＤＶと傾向が一致し、その大きさにはバイアスが見られた。このバイアス量に関しては、噴流中心（ｒ＝０）において、ＬＤＶに対しては、送り台（微動台）による実験で示されたものよりも大きい量が算出された。図１４に示された軸方向平均速度分布ｆ，ｋについては信頼しうる計測値を有することが確認された。
【００８５】
変動速度分布を参照すると、Ｓｔｒｅａｋ法はＬＤＶと傾向が一致するものの、ＬＤＶの変動速度分布ｉに比べＳｔｒｅａｋ法の変動速度分布ｇは大きな値を示した。レーザの照射時間から換算し、Ｓｔｒｅａｋ法はｒ＝−２．２ｍｍにおいて、ＬＤＶに対して０．６３ｐｉｘｅｌｓとなり、送り台検定結果による精密度に収まっており、比較的せん断の影響を受けにくいことが確かめられた。噴流軸付近では、変動が大きくなっており、速度勾配の小さい領域での微小な変動までは捉えられなかったものと考えられる。また、ｒ＝３．３ｍｍおよびｒ＝−３．３ｍｍ付近において変動が周囲に比べ大きな値を示しているが、平均流速が低速であり、相対的に変動が大きく計測されていると考えられる。
【００８６】
図１６および図１７に、Ｓｔｒｅａｋ法とＬＤＶによる噴流軸中心（ｒ＝０ｍｍ）およびせん断層（ｒ＝２ｍｍ）における測定値の軸方向速度確率密度分布（ＰＤＦ）ｊ，ｋおよびｌ，ｍをそれぞれ示す。噴流軸中心では、両計測法とも正常に計測が行われていることが分かる。また、せん断層（ｒ＝２ｍｍ）においてもほぼ計測は正常に行われていることが分かる。
【００８７】
粒子軌跡追跡法による画像処理を行う本発明の流動計測システムの性能と実流動場への適用性が確認され、この流動計測システムにより流体の流動分布（流速および流れ方向）の測定を、精度よく正確に行うことができる。これにより、原子炉圧力容器内のような外界と環境を異にする厳しい閉空間内の流動場の流体（冷却材）の流動状態を正確に測定することができる。
【００８８】
例えば、原子炉圧力容器のダウンカマ部や炉心シュラウド内の流体（冷却材）の流動状態（冷却材の流速や流れ方向、炉心流量、よどみ状態）の計測を行うことができる。原子炉圧力容器内や蒸気発生器内、火力発電プラントの熱交換器内等の熱流動状態、すなわち、流れが一様でない熱流動（流動分布）の測定を正確に精度よく行うことができるので、構成部品の熱疲労や腐食の進捗状況を正確に管理することができる。
【００８９】
具体的には、蒸気発生器や熱交換器等の内部の伝熱管を支える支持板と熱交換器等の内壁等との間隙に形成される流体のよどみ状態を計測することによって、構成部品の熱疲労や腐食の進捗状況を管理することができるので、より安全に、効率的に操業することができる。また、局部的な流体の流動状態が明らかにされるので熱交換器などの最適化の設計が可能となる。例えば、図５に示された流動計測システムによれば、複雑な間隙においても容易に計測することができる。
【００９０】
また、例えば、図５に示された流動計測システムにおいて、レーザシート形成用走査光学系１３からシート状のレーザ光を２方向に直交するように照射して直行する２方向のレーザシートを形成し、直交する２つのレーザシート上の粒子の流跡をそれぞれの画像伝送手段を介してＣＣＤカメラ２４で撮像することにより、粒子の動きを３次元的に測定でき、ひいては、流体の流動状態を３次元的に測定できる。例えば、原子炉圧力容器のダウンカマ部の冷却材流量を正確に測定するには、ダウンカマ部に本体ケーシング４５（図５参照）を半径方向および周方向の少なくとも一方に移動自在に設置し、本体ケーシング４５を半径方向または周方向、すなわち、Ｘ方向またはＹ方向に移動させつつ個々の粒子の移動量を立体的に正確に測定し、流体の流動状態を測定するようにすればよい。
【００９１】
図１に示す流動計測システム１０において、図４および図５に示すように走査光学系と画像伝送系とを一体化することにより実用機として応用することができ、さらに広範囲に適用することができる。例えば、原子炉内のインターナルポンプまわりのよどみ状態等を計測することによりインターナルポンプの流動較正の高精度化を図ることができるとともに、流体のよどみ等による構成部品の熱疲労や腐食の進捗状況を正確に管理することができる。また、ＡＢＷＲの改良に伴う炉内構造物の大型化、シュラウドレス化等に伴う流動試験測定に適用可能であり、さらなる情報の高精度化により設計の合理化、またプラントの効率向上に寄与することも可能である。
【００９２】
また、この流動計測システムは、船舶、航空機あるいは自動車等において、水流あるいは気流と接触する構造部位や機材まわりにおける流体の流動場の計測に適用することができる。この場合、外界に開放された形態で流動場が形成されていてもよい。
【００９３】
例えば、船舶の底部や壁面における水流の流動状態（流速や流れ方向、流量、よどみ状態）を計測することができる。流れが一様でない流動状態の測定をも正確に精度よく行うことができる。これにより船舶の底部や壁面等における流れのよどみ箇所等に発生する貝類等の付着汚染や機材の疲労等の進捗状況を正確に管理することができる。また、船体の形状、構造、部材等に伴う流動試験測定により、船舶設計の最適化、合理化を図ることも可能である。
【００９４】
航空機あるいは自動車においても、流動試験測定により高速化、走行安全性等の最適化設計などに寄与することができる。
【００９５】
【発明の効果】
本発明に係る流体の流動計測システムによれば、必要以上に粒子像径を大きくすることなく、かつ最小の空間分解能となる光ファイバ単繊維に対して十分な粒子像の情報を与えることができ、複雑な流動場における個々の粒子の流速や流れ方向を測定して流体の流速や流れ方向等の流動状態を精度よく正確に測定することができる。
【図面の簡単な説明】
【図１】本発明に係る流体の流動計測システムの一実施形態を示す概念図。
【図２】本発明に係る流体の流動計測システムに組み込まれるレーザシート用走査光学系の原理を示す図。
【図３】本発明に係る流体の流動計測システムに組み込まれる画像伝送手段の原理を示す図。
【図４】本発明に係る流体の流動計測システムを実用的な計測システムに適用した概念図。
【図５】本発明に係る流体の流動計測システムを実用的な計測システムに適用した概念図。
【図６】本発明に係る流体の流動計測システムで撮像された粒子軌跡画像の二値化処理および膨張化処理の結果を示す図。
【図７】本発明に係る流体の流動計測システムで撮像された粒子軌跡画像を示す図。
【図８】図７に示された粒子軌跡画像を画像処理した計算結果を示す図。
【図９】本発明に係る流体の流動計測システムによる計測結果の偏り誤差を示す図。
【図１０】本発明に係る流体の流動計測システムによる計測結果のばらつき誤差を示す図。
【図１１】本発明に係る流体の流動計測システムによる計測結果の確率密度分布を示す図。
【図１２】軸対称噴流の計測実験における噴流出口条件を示す図。
【図１３】本発明の流動計測システムにおいて、粒子軌跡追跡法で画像処理した計算結果を示す図。
【図１４】本発明の流動計測システムとＬＤＶによる計測で得られたＺ＝６Ｄにおける軸方向平均速度分布を示す図。
【図１５】本発明の流動計測システムとＬＤＶによる計測で得られたＺ＝６Ｄにおける軸方向変動速度分布を示す図。
【図１６】本発明の流動計測システムとＬＤＶによる計測で得られたｒ＝０ｍｍにおける軸方向速度確率密度分布を示す図。
【図１７】本発明の流動計測システムとＬＤＶによる計測で得られたｒ＝２ｍｍにおける軸方向速度確率密度分布を示す図。
【符号の説明】
１０流体の流動計測システム
１１レーザ発振装置（光源）
１２送光用ファイバ
１３走査光学系
１４流体の流動場
１５容器
１６ノズル
１７レーザシート
１８レンズ群（レンズ光学系）
１９ａ、１９ｂコリメートレンズ
２０長焦点レンズ
２１シリンドリカルレンズ
２２画像伝送手段
２４ＣＣＤカメラ（撮像手段）
２５光ファイバ
２６イメージガイド
２７対物レンズ
２９画像撮像手段
３０タイミングコントロール手段
３１タイミングスケジューラ
３２シンクロナイザ
３５画像処理手段
３６フレームグラバボード
３７ＰＣＩバス
３８コンピュータ
４０可視化用プローブ
４１透過用窓ガラス
４３４５°直角プリズム
４４ファイバースコープ（画像伝送手段）
４５本体ケーシング
５０貯溜タンク
５１ポンプ
５２リザーバタンク
５３バルブ
５８流体の循環系[0001]
BACKGROUND OF THE INVENTION
  The present invention relates to a fluid flow measurement technique for measuring a fluid flow in a complicated flow field with high accuracy and precision, and in particular, measures a flow velocity and a flow direction of a fluid flowing in a closed space.For fluid flow measurement systemRelated.
[0002]
[Prior art]
  Accurate and precise measurement of fluid flow in complex flow fields is extremely important for improving plant behavior prediction, plant diagnosis, plant performance evaluation and equipment improvement in nuclear power plants, thermal power plants, chemical plants, etc. is important. The fluid flow measurement technology in an actual plant level, particularly in a harsh environment where the environment is different from the outside, is in a situation where flow measurement is extremely difficult due to fluid flow measurement conditions having a complex flow.
[0003]
  Above all, fluid flow measurement such as fluid flow velocity and flow direction in heat flow field under severe environment, such as in reactor pressure vessel and thermal power plant heat exchanger, is complicated fluid flow measurement including accessibility Extremely difficult due to conditions.
[0004]
  However, in recent years, research on visualization technology in a thermal flow field has progressed, and a particle image velocimetry (hereinafter referred to as PIV) capable of measuring fluid flow in a complex flow field with high accuracy and precision is being developed. is there. By using this PIV and an image guide in which optical fibers are bundled together, it is possible to expect fluid flow measurement in a closed space where the environment is different from the environment, such as in a reactor pressure vessel.
[0005]
[Problems to be solved by the invention]
  However, the particle image captured by the CCD camera through the image guide includes not only the luminance pattern of the particles but also the luminance pattern of the optical fiber array (stacked structure) of the image guide. The array luminance pattern of the optical fiber is detrimental when analyzed by the cross-correlation method, and causes erroneous correspondence (error vector) to the particle luminance pattern.
[0006]
  In addition, if the particle image diameter is made larger than one single fiber of an optical fiber that gives the minimum spatial resolution, it is desirable that there are at least 6 to 8 particles in the reference window image. The amount of calculation must be large and the amount of calculation increases enormously, making the algorithm impractical (the amount of calculation increases in proportion to the fourth power when the area to be searched increases).
[0007]
  The present invention has been made in consideration of the above-mentioned circumstances, and performs an optimization of an algorithm for detecting a fluid flow state by performing image processing by a particle trajectory tracking method using an image guide (Particle stroke tracking velocimetry). An object of the present invention is to provide a fluid flow measurement system capable of accurately and accurately performing fluid flow distribution.
[0008]
[Means for Solving the Problems]
  In order to solve the above-described problems, a fluid flow measurement system according to the present invention includes a laser oscillation device that oscillates laser light and the oscillated laser light in a fluid flow field. A scanning optical system for forming a laser sheet to be thrown into a sheet, an image capturing unit that captures a two-dimensional particle trajectory image on the laser sheet from the scanning optical system, and the laser oscillation device and the image capturing unit are timed. Timing control means for synchronous driving, and image processing means for comparing and analyzing the luminance patterns of the two-time particle trajectory images imaged by the image imaging means, and measuring the movement direction and movement amount of each particle, The flow and flow velocity distribution of fluid is measured from the amount of movement of individual particles by this image processing means.The image imaging means comprises image transmission means for optically transmitting a two-dimensional particle trajectory image on the laser sheet from the scanning optical system, and imaging means for imaging the optically transmitted two-dimensional particle trajectory image. Further, the image transmission means is composed of a flexible image guide in which a large number of optical fibers are bundled and integrated, and both end faces are processed to be flat. The image formed on the fiber end face is decomposed into each pixel by an optical fiber and the same image is transmitted to the other end face on the imaging means side, and the particle locus is adjusted by adjusting the light emission time of the laser oscillation device. Larger than the single fiber diameter of the image guideIs.
[0009]
  Also,In order to solve the above-described problem, in the fluid flow measurement system of the present invention,Claim 2As described above, the scanning optical system and the image pickup means are configured to be integrally combined and configured to be movable one-dimensionally or two-dimensionally in a fluid flow field.
[0010]
  Furthermore, in order to solve the above-described problem, in the fluid flow measurement system of the present invention, a third aspect is provided.As described above, the scanning optical system includes an image transmission means for forming an orthogonal laser sheet in a fluid flow field and optically transmitting a two-dimensional particle trajectory image on the orthogonal laser sheet. is there.
[0011]
DETAILED DESCRIPTION OF THE INVENTION
  According to the present inventionOf fluid flow measurement systemAn embodiment will be described with reference to the accompanying drawings.
[0012]
  FIG. 1 is a principle diagram showing the configuration of a fluid flow measurement system according to the present invention.
[0013]
  The fluid flow measurement system 10 of the present invention includes a laser oscillation device 11 that oscillates a high-power laser beam as a light source. As the laser oscillation device 11, a small semiconductor laser, a pulse laser, or the like that oscillates high-power laser light is used. For example, a semiconductor laser having a central oscillation wavelength of 808 nm with continuous light of 16 W output is used as the light source.
[0014]
  The laser light oscillated from the laser oscillation device 11 is guided to the laser sheet forming scanning optical system 13 through the light transmission fiber 12. The scanning optical system 13 projects the laser beam oscillated from the laser oscillation device 11 onto a fluid flow field 14 in a sheet shape.
[0015]
  The fluid flow field 14 is formed in a closed space having a different environment from the outside. Further, according to circumstances, the fluid flow field 14 may be formed in a form open to the outside. The flow field 14 is made of at least one of gas and liquid. Specifically, for example, the velocity field of fluid in a closed space such as a downcomer portion of a reactor pressure vessel, a core shroud, a heat exchanger or a steam generator of a thermal power plant is targeted. In addition, a fluid flow field around a structural part or equipment in contact with a water flow or an air current, such as a ship, an aircraft, an automobile, or the like, can also be targeted. The principle diagram shown in FIG. 1 shows an experimental example of a water flow field as the fluid flow field 14. Water is ejected from a nozzle 16 into a container 15 or piping to form a water flow field 14 as a fluid.
[0016]
  The scanning optical system 13 forms a laser sheet 17 in the fluid flow field 14 by irradiating the fluid flow field 14 with the oscillation laser light from the laser oscillation device 11 in a sheet shape. The fluid flow field 14 is visualized by the laser sheet 17. In the flow measurement system 10, in order to ensure the spatial resolution in the depth direction of the fluid flow field 14, the oscillation laser light is irradiated in a sheet shape.
[0017]
  The scanning optical system 13 is a laser beam irradiation optical system for creating the laser sheet 17 as shown in FIG. When, for example, a semiconductor laser is used as the laser oscillation device 11, the oscillation laser light from the semiconductor laser is emitted from the light transmission fiber 12 with a large divergence angle. The emitted laser beam is not axisymmetric light, has an astigmatism difference, and the shape conversion of the laser beam becomes complicated.
[0018]
  Therefore, the scanning optical system 13 uses a lens group 18 as laser sheet forming optical means for forming the laser sheet 17 from the laser beam emitted from the light transmission fiber 12. The lens group 18 includes two

collimating lenses

19 a and 19 b that collimate the emitted laser beam into a parallel laser beam, and the laser beam that has become parallel light by the

collimating lenses

19 a and 19 b in the thickness direction of the laser sheet 17. The long focal lens 20 that is waisted to about 2 mm and the cylindrical lens 21 that forms the laser sheet 17 are combined to constitute a laser sheet forming optical means. For the

collimating lenses

19a and 19b, for example, plano-convex lenses are used.
[0019]
  A laser sheet 17 is formed in the fluid flow field 14 by the laser light irradiated in a sheet form from the scanning optical system 13, and the fluid flow field 14 is visualized in a sheet form. And the image transmission means 22 is provided facing the measurement range 17a of the laser sheet 17. As shown in FIG. 1, a transmissive window glass 23 is installed on the wall of the container 15 near the image transmission means 22 so that an image on the laser sheet 17 can be picked up by the image pickup means 29.
[0020]
  The image capturing unit 29 includes an image transmitting unit 22 that transmits an image on the laser sheet 17 and a CCD camera 24 as an image capturing unit that captures a two-dimensional particle trajectory image. The image transmission means 22 transmits the two-dimensional particle trajectory image on the laser sheet 17 to a CCD camera 24 as an imaging means. As shown in FIG. 3, the image transmission means 22 uses a flexible image guide 26 in which thousands to tens of thousands of optical fibers 25 are bundled and integrated (melted). As shown in FIG. 3, the image guides 26 are bundled side by side so that the positions of both end faces of the bundled optical fibers 25 correspond exactly. Both end faces of the image guide 26 are finished flat, and an image formed on one fiber end face by the objective lens 27 is decomposed into optical fibers 25 as respective pixel fibers, and the same image is extended to the fiber end face on the CCD camera 24 side. Is transmitting. Recording is performed on the CCD camera from the end face of the fiber via the camera lens 28. The image capturing unit 29 may be configured to directly capture a two-dimensional particle trajectory image on the laser sheet 17 without providing an image guide.
[0021]
  The limit (resolution) of the image of the image guide 24 is determined by the diameter of the single optical fiber 25 and the arrangement method, and the appearance of the image is determined by the total number of optical fibers 25 bundled in the image guide 26 and the arrangement. The light distribution having a light quantity smaller than the core diameter of the optical fiber 25 is averaged while traveling through the optical fiber 25, and the brightness is uniform within each core cross section. Accordingly, in the image guide 26, the core diameter of the optical fiber 25 is the minimum unit of spatial resolution as a pixel.
[0022]
  In the actual image guide 26, a large number of optical fibers 25 are arranged in an ordered and dense manner in a stacked (hexagonal dense arrangement) structure. When a lattice-like object having an interval substantially equal to the optical fiber interval is observed with the image guide 26, moire interference fringes unrelated to the intensity distribution of the image of the object are generated.
[0023]
  Further, in the image guide 26 in which the clad portion is fused between the optical fibers 25, when the thickness of the clad becomes thin, the light entering the core of one optical fiber 25 may leak to the adjacent core. . For this reason, the boundary of the image is blurred and the contrast of the image is lowered. For this reason, when the hexagonal close-packed structure is adopted for the image guide 26, it is a condition that the thickness of the clad is not too thin.
[0024]
  On the other hand, the fluid flow measurement system 10 is provided with a timing control means 30. The timing control means 30 includes a timing scheduler 31 and a synchronizer 32. While the laser oscillation device 11 is oscillated at the timing from the timing scheduler 31, the CCD camera 24 is driven in synchronization with this oscillation timing by the synchronizer 32. That is, the timing control means 30 is driven in synchronization with the laser oscillation device 11 and the CCD camera 22 as the operation means.
[0025]
  For example, a CCD camera having a resolution VGA (640 × 480 pixels) having sensitivity in the near infrared region, a frame rate of 30 Hz, and 8 bits (monochrome 256 gradations) is used. The CCD camera 24 images a two-dimensional particle trajectory image on the laser sheet 17 in synchronization with the oscillation timing from the laser oscillation device 11, and an analog image signal captured by the CCD camera 24 is image processing means 35. The flow velocity distribution and the flow direction of the fluid in the fluid flow field 14 are measured by a particle trajectory tracking method described later. For example, a sequential readout system (progressive scan) camera is used as the CCD camera 24, and the charge accumulation times of the odd and even lines of the arranged CCDs are simultaneously performed. The CCD camera 24 is superior in resolution in the vertical direction as compared with a normal camera.
[0026]
  Further, the image processing means 35 includes a frame grabber board 36 as an A / D converter for digitizing an analog image signal from the CCD camera 24, and a digital image signal from the frame grabber board 36 is a signal line 37 such as a PCI bus. And a computer 38 for digital image processing.
[0027]
  This fluid flow measurement system 10 employs an analog photographing / digital recording method in which analog photographing is performed by a CCD camera 24, digital processing is performed by an image processing means 35, and digital recording is performed.
[0028]
  The fluid flow measurement system 10 shown in FIG. 1 shows an example corresponding to a laboratory. In practical use, as shown in FIGS. 4 and 5, the laser sheet scanning optical system 13 and the image pickup means 29 are integrated. In FIG. 4, the image transmission means 22 of the image pickup means 29 integrated with the scanning optical system 13 is configured as a visualization probe 40 that can be laterally viewed. An optical window (transmission window glass) 41 installed at the tip of the probe 40 is bent at the optical path by a 45 ° right-angle prism 43, and an image is transmitted by a fiberscope 44 constituting an image transmission means. FIG. 4A shows an example in which this integrated fluid flow measurement system is used for the flow field of water as fluid in the container 15 of FIG. In FIG. 5, the integrated scanning optical system 13 and the image transmission means 22 of the image pickup means 29 are stored in a main body casing 45.
[0029]
  Next, a fluid flow measurement method will be described.
[0030]
  In the fluid flow measurement system 10 shown in FIG. 1, the emission time of a semiconductor laser used for illumination is set relatively long, and a flow track (trajectory) of particles is imaged on the CCD camera 24. The captured particle trajectory images at two times that are slightly different in time are regarded as luminance patterns, and comparison / analysis image processing is performed by the computer 38 of the image processing means 35 to estimate the movement amount of each particle. . In the fluid flow measuring apparatus 10, this method is called a particle track tracking velocimetry. This particle trajectory tracking method records sufficient particle image information for an optical fiber monofilament with minimum spatial resolution without recording the particle image diameter more than necessary by recording the particle trajectory. Can be given.
[0031]
  The particle trajectory image represents an image of a trajectory at a certain time of individual particles (foreign matter such as suspended matter) floating and dispersed in water (fluid) on the laser sheet 17 and digitized. It is composed of pixels having luminance. In the course of image processing by the image processing means 35, a limited rectangular region 17a (see FIG. 1) is extracted from each particle trajectory image and image processing is performed. A value of one point in the digitized particle trajectory image is defined as “brightness value”, and a case where this luminance value is distributed in the rectangular region 17a is referred to as a “brightness pattern”.
[0032]
  An algorithm for estimating the amount of movement of individual particles in water (fluid) from the particle trajectory image determines that the obtained particle trajectory image is a change in luminance distribution, and there is a luminance pattern of particles on the laser sheet 17. The amount of movement in time is quantitatively analyzed for each particle. The cross-correlation method is used to analyze the amount of particle movement.Is used.
[0033]
  there,In the fluid flow measurement system 10 according to the present invention, the optical fiber 25 hinders the cross-correlation calculation in order to optimize the fluid flow measurement by the particle trajectory tracking method.Background processingDo.
[0034]
  In the fluid flow measurement system and the measurement method according to Japanese Patent Application No. 2000-207476, a particle image flow velocity measurement (PIV) method is applied. In this PIV method, a frozen image is used. In the normal PIV method, if the particle image diameter is smaller than one optical fiber 25 of the image guide 26, the lost particle information can no longer be recovered. It is strongly affected by the luminance pattern, which is a very bad effect on speed detectionInvite.
[0035]
  Frozen imageIn the PIV method using, required parameters (particle image diameter, number density) are very limited. Therefore, in the particle trajectory tracking method used in the present invention, the emission time of the semiconductor laser used for illumination is made relatively long, the particle trajectory is recorded, and an algorithm for estimating the movement amount of the particles individually rather than the particle group is constructed. Implemented in software. By recording the trajectory of each particle, sufficient particle image information can be given to a single fiber of the optical fiber 25 with the smallest spatial resolution without taking a particle image diameter larger than necessary. it can. This method of particle track tracking velocimetry is possible by using a semiconductor laser that can easily control the laser emission time.
[0036]
  In the particle trajectory tracking method used in the present invention, the movement amount of particles is calculated mainly through the following six steps including this background processing.
[0037]
  (1) Noise removal (background processing)
The particle trajectory image picked up by the CCD camera 24 through the image guide 26 is processed, and the particle trajectory image obtained in time series is obtained for each pixel as an average value of luminance values, and from the luminance value of each particle trajectory image. The luminance pattern of the optical fiber 25 is removed by subtracting the average value.
[0038]
  An image processing algorithm for removing the luminance pattern of the optical fiber array will be described below.
[0039]
  When the particle trajectory image of the measurement data is x, the particle trajectory image x is
[Expression 1]
    x = s + n (1)
        Where s is signal component (luminance pattern of optical fiber array)
              n is a noise component (luminance pattern of particle group)
It is represented by
[0040]
  That is, it is considered that an image in which noise (particle luminance pattern) is superimposed on a signal (luminance pattern of the optical fiber array) is measured in the particle trajectory image x of the measurement data. Luminance pattern removal processing of the optical fiber array by taking the averaging is effective when the signal-to-noise ratio is small and there is no large difference in the frequency of signal and noise, or when measurement can be repeated many times under the same conditions. It is a technique.
[0041]
  When Expression 1 is displayed with time series elements loaded, it is expressed as Expression 2.
[0042]
[Expression 2]
    xⁱ(K) = sⁱ(K) + nⁱ(K) ...... (2)
        i = i-th particle image,
        k = k-th pixel luminance value in the i-th particle image,
Respectively.
[0043]
  When capturing the particle trajectory image, if the number of captured particle trajectory images is M, the addition average x (k) for the M particle trajectory images can be expressed as follows.
[0044]
[Equation 3]

[0045]
  The luminance pattern s of the optical fiber array is M times the luminance value by adding M times since the same pattern appears in the continuous particle trajectory images under the same experimental conditions (light transmission system / light reception system). When the addition averaging process of averaging (dividing by M) is performed, the original luminance value remains unchanged. In other words, an image showing only the optical fiber array without particles can be obtained by the averaging process, and the luminance value (luminance pattern) of the optical fiber array is subtracted for each pixel from the particle trajectory image with the optical fiber array. For example, the luminance pattern of the optical fiber array can be removed.
[0046]
  By this background processing, the background luminance of the acquired image becomes substantially uniform over the entire screen, so that it is effective not only for correlation calculation processing but also for binarization processing.
[0047]
  (2) Binarization and expansion
  As preprocessing for identifying particles from the image processing particle trajectory image, binarization and expansion image processing are performed. The trajectory image of the particles floating in the water (fluid) is affected by the cladding of the optical fiber 25 and is not continuous, but is recorded as an intermittent image (see FIG. 6A). Even if this intermittent particle image is binarized, a portion where information is lost by the clad is not recovered (see FIG. 6B). In order to recognize the trajectory image of these divided particles as the movement trajectory of the same particle, the expanded image processing is performed. In the dilated image processing, if a pixel of interest exists, all the surrounding pixels are converted from 0 to 255 (8-connected dilation processing). In this algorithm, the intermittent particle image divided by repeating this expansion process twice is recovered as a continuous trajectory image (see FIG. 6C).
[0048]
  (3) Labeling and particle image position detection
  Since one particle trajectory image includes a plurality (a large number) of particle images, labeling is performed so that each particle image is recognized as an individual particle. Each connected pixel is assigned the same number, and different connected components are assigned different numbers. This numbering of individual areas is called labeling. The particle image identified by labeling returns to the original image and the center position is calculated. The center coordinates are obtained by the following centroid method (Particle image centroid).
[0049]
  X₀, Y₀Centroid, i, j are the coordinates of each point, and I is the luminance value of the coordinates, the centroid X₀, Y₀Is
[Expression 4]

It is represented by Y₀The same applies to.
[0050]
  (4) Calculation of integer movement by correlation calculation
  Since the center position of the particle trajectory image is detected by the particle image preprocessing (1) to (3), the movement amount of each particle is obtained by the cross correlation method.
[0051]
  The particle trajectory image obtained by digitally processing the two-dimensional particle trajectory image at the first time imaged by the CCD camera 24 of the fluid flow measurement system 10 is defined as a “rectified image” and a limited rectangular area among the “reference image”. A “reference window image”, a particle trajectory image obtained by digitally processing a two-dimensional particle trajectory image at a second time slightly different from the first time as a “search image”, and a limited rectangular area among the “search image” Are defined as “search window images”.
[0052]
  It is known that the amount of movement of particles by the cross-correlation method has a correlation value R expressed by the following equation between the “search window image” and the “reference window image”. This is a method for obtaining the similarity of the luminance patterns of the particles and obtaining the amount of movement of each particle between the two images by comparison.
[0053]
[Equation 5]

[0054]
[Formula 6]

[0055]
As shown in Expression 5 and Expression 6, the correlation value R takes a value from −1 to 1 by dividing (normalizing) the luminance value of the search window image and the reference window image by the root mean square. When the correlation value R is 1, the two images are completely matched, and the larger the correlation value R is, the higher the similarity between the window images is.
[0056]
In

Equations

5 and 6, I₁And I₂Represents the luminance of each pixel of the reference window image and the search window image, and ξ and η represent the relative positions of the search window image and the reference window image. The positions ξ and η at which the correlation values are maximum correspond to the movement amounts ΔX ′ and ΔY ′ of the particles in the image. The size nxm of the reference window image and the relative position to the reference window image are determined from the predicted minimum and maximum velocities.
[0057]
  (5) Calculation of sub-pixel movement
  Since the integer value of the movement amount of the particle image is obtained by the correlation calculation of (4), the sub-pixel movement amount is estimated by normal distribution approximation (Gaussian peak-fit) using three points near the peak of the correlation coefficient. It is calculated independently for the x and y directions.
[0058]
  (6) Calibration
  By the above procedure, the movement amount of the particle trajectory image is obtained in units of pixels, and is converted into the movement amount in the real space by calibration.
[0059]
  Fig. 7 shows a particle trajectory image in which a glycerin solution mixed with tracer particles is placed on a rotating disk and images a pseudo flow, and Fig. 8 shows image processing of 60 frames of the particle trajectory image by the particle trajectory tracking method. The result of calculating the movement amount of each particle is shown.
[0060]
【Example】
  The performance evaluation of the fluid flow measurement system according to the present invention that performs image processing by the particle trajectory tracking method and the applicability of the flow measurement system to an actual flow field will be described below.
[0061]
[Example 1]
(Verification of the amount of movement using a fine motion table)
  The particle image drawn on the flat plate is moved by a fine movement table (mechanical stage), the movement amount is measured using a calibration plate of 70 mm × 70 mm, and the true value actually moved and the flow measurement of the present invention are measured. The system was compared with the measurement value measured by the Strak method, the error was measured, and the accuracy was tested.
[0062]
  A calibration plate having 64 grid points in a 70 mm × 70 mm rectangle was used. Therefore, the measurement range is 70 mm × 70 mm. The imaging magnification at this time was 5.011 mm / pixel.
[0063]
  Three types of particles with particle sizes of 3, 6, and 8 pixels were used. The particle image is obtained by capturing particles seeded in a fluid with laser sheet illumination, printing them out with a printer, and bonding them onto a flat plate.
[0064]
  The image of the particle image was obtained by capturing a particle image on a flat plate placed at a reference position (reference position in calibration) and a particle image after movement by a fine moving table with a CCD camera, and forming this as one image pair. The calculation of the movement amount vector from the image pair is based on the cross-correlation method, which is equivalent to the method used for measuring the actual flow field. The number of vectors calculated from one set of image pairs was about 150.
[0065]
  The movement of the flat plate on which the particle image was drawn was performed up to 10 mm in 50 μm increments in the x direction, 200 images were acquired, and 100 image pairs were obtained. The amount of movement was calculated by the cross-correlation method and normal distribution approximation (Gaussian peak-fit).
[0066]
  The optical fiber single fiber diameter of the image guide 26 that provides the minimum spatial resolution is df = 3.12 pixels based on the frequency analysis result of the image guide 26. Therefore, for the three types of

particle image diameters

3, 6, and 8 pixels, respectively. dτ / df = 0.96, 1.92, 2.56 times.
[0067]
  FIG. 9, FIG. 10, and FIG. 11 show the bias error (Bias error), variation error (RMS error), and probability density distribution (PDF) of the measurement results by the Strak method (particle trajectory tracking method), respectively. In each figure, the curve of dτ / df = 0.96 is a, the curve of dτ / df = 1.92 is b, the curve of dτ / df = 2.56 is c, and the horizontal axis in the figure is the true movement Amount (Actual displacement) [pix. ] Is shown. The PDF shown in FIG. 11 is for true movement amounts of 0.5 pixels, 2.5 pixels, and 4.5 pixels, respectively.
[0068]
  FIG. 9 shows the deviation error of the measured value of the amount of movement in the x direction when the flat plate is moved in the x direction, and is a plot of how far the measured value deviates from the true number of moving pixels. The error indicates the difference between the average value and the true value of all vectors within the measurement range.
[0069]
  FIG. 10 shows the root mean square of the measurement values within the measurement range, and FIG. 11 shows that the true movement amount is 0.5 pixel, 2.5 pixel, 4.5 pixel from the left mountain. In this case, it is shown how many (%) measured values by the Strak method (particle tracking method: PTV) were actually measured at a deviation from the true value when the total error number was 10 (%). . The narrower the mountain, the smaller the error. For example, when 2.5 pixels of the middle mountain are moved, measuring 2 pixels means about 9 out of 100 pixels.
[0070]
  As shown in FIGS. 9 to 11, in the particle trajectory tracking method, the difference in the measured value due to the particle diameter does not appear so much, and the deviation error β is within the range of about β = ± 0.2 pixels.
[0071]
  Since the optical fiber 25 used in the experiment has a cladding that does not transmit light (having a relatively low transmittance), a loss of particle image is caused and a vibration effect is observed. In the bias error shown in FIG. 9, no change in value depending on the particle image diameter is observed, which is because sufficient information can be given to the optical fiber diameter by the particle trajectory.
[0072]
  The variation error by the Strak method shown in FIG. 10 also shows almost the same value regardless of the particle image diameter, but it is a relatively large value. This is because, as can be seen from the probability density distribution, the shape of the histogram does not become a normal distribution, but a distribution near the horizontal near the true value. This is because the particle trajectory image is correlated, that is, the luminance value of the particle trajectory is comparatively gentle, and the correlation coefficient distribution of the image is also gentle. In the particle trajectory tracking method, it is assumed that the particle is larger than the normal particle image._k(0.3 to 0.5 pixels). However, it has been shown that there is no way of generating an error that causes the PDF to crack greatly depending on the particle image diameter, and it has been confirmed that the accuracy of the bias error is effective.
[0073]
[Example 2]
(Measurement in actual flow field)
  A flow field 14 of water (fluid) is provided in an experimental container 15 as shown in FIG. 1, and an axisymmetric jet is measured by the flow measurement system 10 of the present invention, and the performance of the fluid flow measurement system in actual flow field measurement is measured. It was investigated.
[0074]
  The fluid flow measurement system was examined in comparison with measurement using a Doppler laser velocimeter (LDV).
[0075]
  The container 15 is, for example, a cubic container having a length of 250 mm, a width of 250 mm, and a height of 250 mm, and has a rectangular jet outlet (nozzle) 16 at the center. Tap water was used as the working fluid, and 3-30 μm spherical polyethylene particles were mixed as tracer particles. As shown in FIG. 1, in the experimental facility, the working fluid was temporarily stored in the lower storage tank 50 and then supplied to the upper reservoir tank 52 by the pump 51. The hydraulic fluid supplied to the reservoir tank 52 partially overflows through the partition amount 52a, so that the hydraulic head pressure (head) is kept constant. The working fluid is rectified from the reservoir tank 52 through a flow regulating valve 53, a rectifying honeycomb and a mesh (not shown), and ejected from the axisymmetric nozzle 16 into the equipment container 15 to form a flow field of water (fluid). did. The water overflowed in the reservoir tank 52 is returned to the storage tank 50 through the return pipe 54.
[0076]
  On the other hand, the nozzle 16 in the container 15 was contracted to an area ratio of about 10: 1, a stable velocity distribution was obtained, and a fluid (water) flow measurement experiment was performed in a substantially steady state. The coordinate system of the flow field of water (fluid) has the center of the outlet of the nozzle 16 as the origin, the jet direction (height direction) is Z, the radial direction is r, the circumferential direction is θ, and velocity components in the respective directions are w and u. , V.
[0077]
  The outlet speed of the working fluid from the nozzle 16 was w = 165 mm / sec, and the number of lay nozzles with respect to the outlet diameter (D = 7 mm) at this time was Re = 1155. The water overflowing the container 15 was once received by the receiving tank 56 and then returned to the storage tank 50 through the return water pipe 57 to constitute a water (fluid) circulation system 58.
[0078]
  As the scanning optical system 13 and the image transmission means 22, a visualization probe 40 as shown in FIG. 4 in which the scanning optical system 13 and the image transmission means 22 are integrated as shown in FIG. 4 is used. As shown in FIG. 4, the visualization probe 40 is inserted into the water (fluid) flow field 14 from above the container 15, and the laser beam is injected vertically from above, and the particle trajectory from the opposing probe 40 by the fiber scope 44. An image was taken. The distance between the laser sheet 17 and the light receiving surface of the probe 40 facing the laser sheet 17 is about 30 mm. The measurement range is a range where the grid points of the calibration plate exist, for example, a rectangular range of 10 mm × 10 mm, and the thickness of the laser sheet 17 is 2 mm.
[0079]
  The region to be imaged is said to be a development region where the potential core of the water jet disappears with respect to the diameter D of the nozzle 16, and the momentum exchange with the surrounding fluid extends to the center of the jet axis and has a relatively gentle velocity distribution. It set to the area | region of Z = 6D obtained.
[0080]
  The captured particle trajectory image was subjected to image processing by the image processing means 35 by a particle trajectory tracking method (Strak method).
[0081]
  On the other hand, measurement by LDV (Doppler laser velocimeter) was also performed separately in a cross section perpendicular to the jet axis at Z = 6D. These measurements were made by assuming that the flow field is steady and adjusting the flow rate to the same conditions.
[0082]
  FIG. 12 shows jet outlet conditions (Reynolds number: Re = 1155) measured by LDV. FIG. 12 shows curves d and e of the average speed (LDV (mean)) and the fluctuation speed (LDV (rms)). FIG. 13 shows a velocity vector map in the vicinity of the height direction Z = 6D from the nozzle opening by the particle trajectory tracking method.
[0083]
  In the particle trajectory image, the velocity was calculated by the particle trajectory tracking method (Strak method) described above. FIG. 14 and FIG. 15 show the axial average speed distribution f and the fluctuation speed distribution g at Z = 6D, respectively. The horizontal axis is the radial direction, and the vertical axis is the axial direction. In this system measurement and LDV measurement, this area is extracted from the obtained average speed distribution h in the axial direction and its fluctuation speed distribution i, and statistical processing is performed to compare the results of the two measurement methods. It was. The trajectory image by the fiberscope in this fluid flow measurement system is calculated corresponding to each particle, and since it is a relatively low density seeding, velocity data obtained from a set of image pairs Less is. In order to compensate for such a small amount of statistics, it is necessary to take a large amount of images in the time direction. In this experiment, statistics were calculated from 16000 images. The axial average and fluctuation speed distributions f and g shown in FIGS. 14 and 15 were about 1000 statistics at each point. In the measurement by LDV, average / fluctuation speed data was calculated from the statistic 3000.
[0084]
  When comparing the axial direction fluctuation speed distributions g and i shown in FIG. 15, the trend method agrees with the tendency of LDV, and a bias is seen in the magnitude. Regarding the bias amount, an amount larger than that shown in the experiment with the feed base (fine movement base) was calculated for the LDV at the jet center (r = 0). It was confirmed that the axial average velocity distributions f and k shown in FIG. 14 have reliable measurement values.
[0085]
  Referring to the fluctuation speed distribution, the Stroke method has the same tendency as that of LDV, but the fluctuation speed distribution g of the Stroke method shows a larger value than the fluctuation speed distribution i of LDV. Converted from the laser irradiation time, the Strak method is 0.63 pixels with respect to LDV at r = -2.2 mm, which is within the precision based on the test result of the feed table, and is relatively less susceptible to shear. It was confirmed. In the vicinity of the jet axis, the fluctuation is large, and it is considered that even minute fluctuations in the region where the velocity gradient is small were not captured. Moreover, although the fluctuation | variation has shown the large value compared with the circumference | surroundings in the vicinity of r = 3.3mm and r = -3.3mm, it is thought that the average flow velocity is low speed and the fluctuation | variation is measured comparatively large.
[0086]
  FIGS. 16 and 17 show axial velocity probability density distributions (PDF) j, k, and l, m of measured values at the stream axis center (r = 0 mm) and shear layer (r = 2 mm) by the Strak method and LDV, respectively. Show. At the center of the jet axis, it can be seen that both measurement methods are measuring normally. It can also be seen that the measurement is normally performed in the shear layer (r = 2 mm).
[0087]
  The flow measurement system of the present invention that performs image processing by the particle trajectory tracking method and the applicability to the actual flow field have been confirmed. With this flow measurement system, the flow distribution (flow velocity and flow direction) of the fluid can be accurately measured. Can be done accurately. Thereby, it is possible to accurately measure the flow state of the fluid (coolant) in the flow field in a strict closed space in which the environment is different from the outside, such as in the reactor pressure vessel.
[0088]
  For example, it is possible to measure the flow state (coolant flow velocity, flow direction, core flow rate, stagnation state) of the fluid (coolant) in the downcomer portion of the reactor pressure vessel or the core shroud. Because it is possible to accurately and accurately measure the heat flow state in reactor pressure vessels, steam generators, heat exchangers of thermal power plants, etc., that is, heat flow (flow distribution) with non-uniform flow. The progress of thermal fatigue and corrosion of components can be managed accurately.
[0089]
  Specifically, by measuring the stagnation state of the fluid formed in the gap between the support plate supporting the heat transfer tubes inside the steam generator, heat exchanger, etc. and the inner wall, etc. of the heat exchanger, the component parts Since the progress of thermal fatigue and corrosion can be managed, it is possible to operate more safely and efficiently. Further, since the local fluid flow state is clarified, it is possible to design an optimization of a heat exchanger or the like. For example, according to the flow measurement system shown in FIG. 5, it is possible to easily measure even in a complicated gap.
[0090]
  Further, for example, in the flow measurement system shown in FIG. 5, a laser sheet forming scanning optical system 13 irradiates a sheet-shaped laser beam so as to be orthogonal to the two directions, thereby forming a two-way laser sheet. In addition, the movement of the particles can be measured three-dimensionally by imaging the track of the particles on the two orthogonal laser sheets with the CCD camera 24 via the respective image transmission means, and thus the fluid flow state is 3 It can be measured dimensionally. For example, in order to accurately measure the coolant flow rate in the downcomer portion of the reactor pressure vessel, a main body casing 45 (see FIG. 5) is installed in the downcomer portion so as to be movable in at least one of the radial direction and the circumferential direction. The amount of movement of each particle may be accurately measured three-dimensionally while moving 45 in the radial direction or the circumferential direction, that is, the X direction or the Y direction, and the fluid flow state may be measured.
[0091]
  The flow measurement system 10 shown in FIG. 1 can be applied as a practical machine by integrating the scanning optical system and the image transmission system as shown in FIGS. 4 and 5, and can be applied in a wider range. . For example, it is possible to improve the accuracy of internal pump flow calibration by measuring the stagnation state around the internal pump in the reactor, and the progress of thermal fatigue and corrosion of components due to fluid stagnation, etc. The situation can be managed accurately. In addition, it can be applied to flow test measurement due to the increase in the size of the reactor internal structure due to the improvement of ABWR, shroudless, etc., and contribute to the rationalization of design and the improvement of plant efficiency by further improving the accuracy of information. Is also possible.
[0092]
  In addition, this flow measurement system can be applied to measurement of a fluid flow field around a structural part or equipment in contact with water flow or air flow in a ship, an aircraft, an automobile, or the like. In this case, the flow field may be formed in a form open to the outside.
[0093]
  For example, it is possible to measure the flow state (flow velocity, flow direction, flow rate, stagnation state) of the water flow at the bottom or wall surface of the ship. It is possible to accurately and accurately measure the flow state where the flow is not uniform. As a result, it is possible to accurately manage progress such as adhesion contamination of shellfish or the like generated at the stagnation point of the flow at the bottom or wall surface of the ship, and fatigue of the equipment. It is also possible to optimize and rationalize ship design by measuring flow tests associated with the shape, structure, members, etc. of the hull.
[0094]
  Even in an aircraft or automobile, flow test measurement can contribute to optimization design such as high speed and driving safety.
[0095]
【The invention's effect】
  According to the present inventionAccording to the fluid flow measurement system, sufficient particle image information can be given to an optical fiber monofilament with minimum spatial resolution without increasing the particle image diameter more than necessary,By measuring the flow velocity and flow direction of individual particles in a complex flow field, the flow state such as the flow velocity and flow direction of the fluid can be accurately and accurately measured.
[Brief description of the drawings]
FIG. 1 is a conceptual diagram showing an embodiment of a fluid flow measurement system according to the present invention.
FIG. 2 is a diagram showing the principle of a laser sheet scanning optical system incorporated in a fluid flow measurement system according to the present invention.
FIG. 3 is a diagram showing the principle of image transmission means incorporated in a fluid flow measurement system according to the present invention.
FIG. 4 is a conceptual diagram in which the fluid flow measurement system according to the present invention is applied to a practical measurement system.
FIG. 5 is a conceptual diagram in which the fluid flow measurement system according to the present invention is applied to a practical measurement system.
FIG. 6 is a diagram showing a result of binarization processing and expansion processing of a particle trajectory image captured by the fluid flow measurement system according to the present invention.
FIG. 7 is a diagram showing a particle trajectory image captured by the fluid flow measurement system according to the present invention.
8 is a diagram showing a calculation result obtained by performing image processing on the particle locus image shown in FIG. 7;
FIG. 9 is a diagram showing a bias error of a measurement result by the fluid flow measurement system according to the present invention.
FIG. 10 is a view showing a variation error of a measurement result by the fluid flow measurement system according to the present invention.
FIG. 11 is a diagram showing a probability density distribution of a measurement result obtained by the fluid flow measurement system according to the present invention.
FIG. 12 is a diagram showing jet outlet conditions in an axially symmetric jet measurement experiment.
FIG. 13 is a diagram showing a calculation result obtained by performing image processing by a particle trajectory tracking method in the flow measurement system of the present invention.
FIG. 14 is a diagram showing an axial average velocity distribution at Z = 6D obtained by measurement using the flow measurement system of the present invention and LDV.
FIG. 15 is a diagram showing an axial direction fluctuation velocity distribution at Z = 6D obtained by measurement using the flow measurement system of the present invention and LDV.
FIG. 16 is a diagram showing an axial velocity probability density distribution at r = 0 mm obtained by measurement using the flow measurement system of the present invention and LDV.
FIG. 17 is a diagram showing an axial velocity probability density distribution at r = 2 mm obtained by measurement using the flow measurement system of the present invention and LDV.
[Explanation of symbols]
10 Fluid flow measurement system
11 Laser oscillator (light source)
12 Optical fiber
13 Scanning optical system
14 Fluid flow field
15 containers
16 nozzles
17 Laser sheet
18 Lens group (lens optical system)
19a, 19b Collimating lens
20 Long focus lens
21 Cylindrical lens
22 Image transmission means
24 CCD camera (imaging means)
25 optical fiber
26 Image Guide
27 Objective lens
29 Image pickup means
30 Timing control means
31 Timing scheduler
32 Synchronizer
35 Image processing means
36 frame grabber board
37 PCI bus
38 computers
40 Visualization probe
41 Window glass for transmission
43 45 ° right angle prism
44 Fiberscope (image transmission means)
45 Body casing
50 storage tank
51 pump
52 Reservoir tank
53 Valve
58 Fluid Circulation System

Claims

A laser oscillation device that oscillates laser light, a scanning optical system for forming a laser sheet that oscillates the oscillated laser light into a fluid flow field, and a two-dimensional particle trajectory image on the laser sheet from the scanning optical system The image pickup means for picking up images, the timing control means for synchronously driving the laser oscillation device and the image pickup means at the timing, and the luminance pattern of the particle trajectory image at two times picked up by the image pickup means are compared. It comprises an image processing means for analyzing and measuring the moving direction and moving amount of each particle, and configured to measure the fluid flow and flow velocity distribution from the moving amount of each particle by this image processing means ,
The image imaging means includes an image transmission means for optically transmitting a two-dimensional particle locus image on a laser sheet from a scanning optical system, and an imaging means for imaging an optically transmitted two-dimensional particle locus image.
Further, the image transmission means is composed of a flexible image guide in which a large number of optical fibers are bundled and integrated, and both end surfaces are processed to be flat.
The image guide is configured to transmit the same image to the other end surface on the imaging means side by decomposing an image formed on one fiber end surface by an objective lens into each pixel by an optical fiber,
A fluid flow measurement system, wherein a light emission time of the laser oscillation device is adjusted so that a particle locus in the particle locus image at the two times is larger than a single fiber diameter of the image guide.

2. The fluid flow measurement system according to claim 1 , wherein the scanning optical system and the image capturing unit are configured to be integrated and configured to be movable in a one-dimensional or two-dimensional manner in a fluid flow field.

3. The fluid flow according to claim 2 , wherein the scanning optical system is provided with image transmission means for forming orthogonal laser sheets in a fluid flow field and optically transmitting a two-dimensional particle image on the orthogonal laser sheets. Measuring system.