JP2006258553A - 流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラム - Google Patents

Abstract

【課題】 粉体を含む様々な流動性媒体の流れを評価することが可能な、流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムを提供する。
【解決手段】 流動性媒体の流れ評価方法は、粉体を含む流動性媒体の流れを、該流動性媒体と異なる少なくとも一種類のトレーサー及び該流動性媒体を透過すると共に該少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価する。流動性媒体の流れ評価装置は、上記の流動性媒体の流れ評価方法を実行する手段を含む。流動性媒体の流れ評価プログラムは、上記の流動性媒体の流れ評価方法をコンピュータに実行させる。
【選択図】 図１

Description

本発明は、流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムに関する。

電子写真装置は、感光体に形成された静電潜像に対して、現像ユニットから有色のトナーを含む現像剤を供給し、静電潜像を帯電した有色の現像剤で現像することによって、可視画像を形成する装置である。また、現像剤は、現像ユニット内におけるスクリューで攪拌及び輸送され、現像剤の摩擦帯電によって所定の帯電量に帯電させられる。ここで、現像ユニット内における有色の現像剤の流れを評価又は計測することは、現像剤の帯電量を制御する上で重要であると考えられる。

従来、流体中に粒子状のトレーサーを混入させ、流体とトレーサーの流れが一致するという仮定に基づいて、流体の流れをトレーサーの流れを追跡する方法が開示されている。特に、可視の粒子状のトレーサーを用いる場合には、流体の流れを可視化することができる。

例えば、特許文献１には、微小なトレーサー粒子を流体中に混入し、レーザ光を流体に照射し、流体中に混入されているトレーサー粒子をテレビカメラで撮像し、得られた画像を処理することによって、流体の流れを可視化する、流体の流れの可視化方法が開示されている。

ここで、特許文献１に開示されるような先行技術においては、トレーサー粒子をテレビカメラで追跡する必要がある。ところが、電子写真装置で用いられる現像剤は、一般に有色である。特許文献１に開示されるように、有色である現像剤にトレーサー粒子を混入して可視のレーザー光をその現像剤に照射する場合には、可視のレーザー光は、表面における現像剤によって遮断されてしまい、現像剤の内部に混入したトレーサー粒子にまでそのレーザー光が到達することができない。このため、現像剤の内部に混入したトレーサー粒子をテレビカメラで撮像することができない。よって、現像剤の表面側における現像剤の流れをテレビカメラで撮像することはできるが、現像剤の内部における現像剤の流れを撮像することはできない。

従って、有色の現像剤の表面及び内部の両方における現像剤の流れを評価することができる方法が望まれる。

また、可視光ではなく、Ｘ線及びガンマ線のような高エネルギーの電磁波を用いた流体の流れに関する先行技術も開示されている。

例えば、特許文献２には、粒子状造影剤または分子量１０万以上の高分子量化造影剤を用いて、コンピュータ断層装置により、液体流通装置内の液体流動状態を検査する、液体流通装置内の液体流動状態の非破壊検査方法が、開示されている。

また、特許文献３には、流れに追従でき、陽電子を放射する同位体でラベル化された多量の粒子を流体中に射出し、陽電子が静止される際放射されるガンマ線を検知し、流体の流れの表示を生成させる、流体流れの表示方法も開示されている。

しかしながら、特許文献２及び３に開示されるような先行技術は、液体などの流体の流れを検査又は表示する技術であり、乾式トナーを含む乾式の現像剤のような粉体の流れを評価することを目的としていない。

さらに、特許文献４には、複数方向からのＸ線の平行光線の透過光による流体中の粉体の像を用いて、流体中の被検出物である粉体の３次元空間分布の立体像および当該分布の３次元移動速度分布の簡易再構築装置が開示されている。

しかしながら、特許文献４に開示される先行技術は、流体中における被検出物である粉体の３次元空間分布及び３次元移動速度分布を得る技術であり、被検出物である粉体（又は流体のいずれか）が、Ｘ線と相互作用する必要がある。よって、粉体（又は流体のいずれか）が、Ｘ線と乏しい相互作用を有する場合には、特許文献４に開示される先行技術を用いることができない。例えば、乾式トナーを含む乾式の現像剤は、通常、Ｘ線を透過する粉体のみからなるため、乾式トナーを含む乾式の現像剤に、特許文献４に開示される先行技術を用いることはできない。
特開平６−６６６７５号公報 特開２００４−５３５５２号公報 特開昭６１−８６７２号公報 特開平１１−３２６００８号公報

本発明は、粉体を含む様々な流動性媒体の流れを評価することが可能な、流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムを提供することを目的とする。

請求項１に記載の発明は、流動性媒体の流れ評価方法において、粉体を含む流動性媒体の流れを、該流動性媒体と異なる少なくとも一種類のトレーサー及び該流動性媒体を透過すると共に該少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価することを特徴とする。

請求項２に記載の発明は、請求項１に記載の流動性媒体の流れ評価方法において、前記少なくとも一種類のトレーサーを含む前記流動性媒体に前記電磁波を照射し、該電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、該画像を用いて前記粉体を含む流動性媒体の流れを評価することを含むことを特徴とする。

請求項３に記載の発明は、請求項１又は２に記載の流動性媒体の流れ評価方法において、前記流動性媒体は、粉体であることを特徴とする。

請求項４に記載の発明は、請求項１乃至３のいずれか一項に記載の流動性媒体の流れ評価方法において、前記粉体は、トナーを含む現像剤であることを特徴とする。

請求項５に記載の発明は、請求１乃至４のいずれか一項に記載の流動性媒体の流れ評価方法において、前記トレーサーは、粒子であり、前記流動性媒体の流れと共に移動する該粒子の位置の対応付けを行うことを含むことを特徴とする。

請求項６に記載の発明は、請求項１乃至４のいずれか一項に記載の流動性媒体の流れ評価方法において、前記トレーサーは、液体又は粒子集合体であり、前記流動性媒体の流れと共に移動する該液体又は粒子集合体の対応付けを行うことを含むことを特徴とする。

請求項７に記載の発明は、請求項１乃至６のいずれか一項に記載の流動性媒体の流れ評価方法において、前記少なくとも一種類のトレーサーは、複数種類のトレーサーであることを特徴とする。

請求項８に記載の発明は、請求項１乃至７のいずれか一項に記載の流動性媒体の流れ評価方法において、前記電磁波は、互いに異なる波長を備えた複数の電磁波を含むことを特徴とする。

請求項９に記載の発明は、請求項１乃至６のいずれか一項に記載の流動性媒体の流れ評価方法において、前記少なくとも一種類のトレーサーは、複数種類のトレーサーであり、前記電磁波は、互いに異なる波長を備えた複数の電磁波を含み、該互いに異なる波長を備えた複数の電磁波は、前記少なくとも一種類のトレーサーの一つの吸収端の波長よりも短い波長を備えた第一の電磁波及び該吸収端の波長よりも長い波長を備えた第二の電磁波を含み、前記第一の電磁波によって形成される前記複数種類のトレーサーを含む前記流動性媒体の第一の画像及び前記第二の電磁波によって形成される前記複数種類のトレーサーを含む前記流動性媒体の第二の画像に少なくとも減算処理を施して、前記複数種類のトレーサーを識別することを含むことを特徴とする。

請求項１０に記載の発明は、請求項８又は９に記載の流動性媒体の流れ評価方法において、前記互いに異なる波長を備えた複数の電磁波は、それぞれ、異なる時間に前記少なくとも一種類のトレーサーを含む前記流動性媒体に照射されることを特徴とする。

請求項１１に記載の発明は、請求項８乃至１０のいずれか一項に記載の流動性媒体の流れ評価方法において、前記互いに異なる波長を備えた複数の電磁波は、単一の電磁波の発生源から選択的に発生させられることを特徴とする。

請求項１２に記載の発明は、請求項８乃至１０のいずれか一項に記載の流動性媒体の流れ評価方法において、前記互いに異なる波長を備えた複数の電磁波は、複数の電磁波の発生源から発生させられることを特徴とする。

請求項１３に記載の発明は、請求項１乃至１２のいずれか一項に記載の流動性媒体の流れ評価方法において、前記電磁波は、エックス線であることを特徴とする。

請求項１４に記載の発明は、請求項１３に記載の流動性媒体の流れ評価方法において、前記エックス線の波長は、前記少なくとも一種類のトレーサーの少なくとも一つの吸収端の波長を含むことを特徴とする。

請求項１５に記載の発明は、請求項１乃至１４のいずれか一項に記載の流動性媒体の流れ評価方法において、前記電磁波を、前記少なくとも一種類のトレーサーを含む前記流動性媒体に対して複数の異なる方向から照射し、該複数の方向に対して、該電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、前記複数の方向に対して取得された該画像の対応付けを行うことを含むことを特徴とする。

請求項１６に記載の発明は、請求項１５に記載の流動性媒体の流れ評価方法において、前記複数の異なる方向から照射される前記電磁波は、複数の電磁波の発生源から発生させられることを特徴とする。

請求項１７に記載の発明は、請求項１５に記載の流動性媒体の流れ評価方法において、前記複数の異なる方向から照射される前記電磁波は、移動可能な電磁波の発生源から発生させられる。

請求項１８に記載の発明は、請求項１乃至１４のいずれか一項に記載の流動性媒体の流れ評価方法において、前記少なくとも一種類のトレーサーを含む前記流動性媒体の全体を移動させ、前記少なくとも一種類のトレーサーを含む前記流動性媒体の移動の前後に、前記電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、該移動の前後に取得された該画像の対応付けを行うことを含むことを特徴とする。

請求項１９に記載の発明は、流動性媒体の流れ評価装置において、請求項１乃至１８のいずれか一項に記載の流動性媒体の流れ評価方法を実行する手段を含むことを特徴とする。

請求項２０に記載の発明は、流動性媒体の流れ評価プログラムにおいて、請求項１乃至１８のいずれか一項に記載の流動性媒体の流れ評価方法をコンピュータに実行させることを特徴とする。

本発明によれば、粉体を含む様々な流動性媒体の流れを評価することが可能な、流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムを提供することができる。

次に、本発明の実施の形態を図面と共に説明する。

まず、本発明の第一の態様である流動性媒体の流れ評価方法は、粉体を含む流動性媒体の流れを、流動性媒体と異なる少なくとも一種類のトレーサー及び流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価する。本発明の第一の態様である流動性媒体の流れ評価方法によれば、粉体を含む様々な流動性媒体の流れを評価することが可能な流動性媒体の流れ評価方法を提供することができる。

また、本発明の第二の態様である流動性媒体の流れ評価装置は、上記の流動性媒体の流れ評価方法を実行する手段を含む。すなわち、流動性媒体の流れ評価装置は、粉体を含む流動性媒体の流れを、流動性媒体と異なる少なくとも一種類のトレーサー及び流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価する手段を含む。本発明の第二の態様である流動性媒体の流れ評価装置によれば、粉体を含む様々な流動性媒体の流れを評価することが可能な流動性媒体の流れ評価装置を提供することができる。

さらに、本発明の第三の態様である流動性媒体の流れ評価プログラムは、上記の流動性媒体の流れ評価方法をコンピュータに実行させる。すなわち、流動性媒体の流れ評価プログラムは、粉体を含む流動性媒体の流れを、流動性媒体と異なる少なくとも一種類のトレーサー及び流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価する手順をコンピュータに実行させる。本発明の第三の態様である流動性媒体の流れ評価プログラムによれば、粉体を含む様々な流動性媒体の流れを評価することが可能な流動性媒体の流れ評価プログラムを提供することができる。なお、本発明の第三の態様である流動性媒体の流れ評価プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することもできる。

ここで、粉体を含む流動性媒体は、流れの評価の対象であり、粉体のみからなる流動性媒体並びに粉体及び流体の両方からなる流動性媒体を含む。粉体（粉末）は、１ｍｍよりも小さい（平均）直径を有する複数の固体粒子の集合体又は凝集体である。流体は、液体及び気体並びに液体及び気体の混合物を含む。粉体及び流体を含む流動性媒体は、流動性媒体に照射される電磁波を透過する性質を有する材料からなる。ここで、流動性媒体が、粉体であるときには、粉体の流れを評価することができる。例えば、粉体が、トナーを含む現像剤である場合には、現像剤（を含む流動性媒体）の流れを評価することができる。

また、少なくとも一種類のトレーサーは、流動性媒体の流れの評価に使用するためのものであって、流動性媒体と異なる。少なくとも一種類のトレーサーは、流動性媒体に照射される電磁波の少なくとも一部を吸収する性質を有する材料を含む。

トレーサーとしては、例えば、流動性媒体に照射される電磁波が、エックス線であるときには、トレーサーとしてエックス線を吸収することが可能な金属を用いることができる。エックス線を吸収することが可能な金属としては、例えば、鉄（Ｆｅ）、銅（Ｃｕ）、銀（Ａｇ）、鉛（Ｐｂ）、タングステン（Ｗ）、水銀（Ｈｇ）などが挙げられる。なお、トレーサーの材料の吸収端が、流動性媒体用の容器の材料の吸収端と近似する場合には、トレーサーによるエックス線の吸収と流動性媒体用の容器の材料によるエックス線の吸収とを明確に区別するため（トレーサーと背景とのコントラストを向上させるため）、トレーサーの材料を、流動性媒体用の容器の材料の吸収端と異なる吸収端を備えた材料に変更することが好ましい。

さらに、電磁波は、流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波である。電磁波としては、例えば、紫外線（１ｎｍを超える４００ｎｍ以下の波長を備えた電磁波）、エックス線（０．１ｎｍ以上１ｎｍ以下の波長を備えた電磁波）、ガンマ線（０．１ｎｍ未満の波長を備えた電磁波）などが挙げられる。しかしながら、電磁波が、流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収されるのであれば、電磁波は、特に限定されず、例えば、可視光（７８０ｎｍ以上４００ｎｍ未満の波長を備えた電磁波）であってもよい。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、電磁波は、エックス線である。この場合には、流動性媒体に対するエックス線の透過率は、一般的に、流動性媒体に対する紫外線の透過率よりも高い傾向があるため、より多くの種類の流動性媒体の流れを評価することができる。また、エックス線のエネルギーは、ガンマ線のエネルギーよりも低いため、エックス線の使用は、ガンマ線の使用よりも安全である。また、ガンマ線の使用は、ガンマ線を密閉する空間を管理する必要があるため、エックス線の使用は、ガンマ線の使用よりも低コストである。なお、エックス線は、特性エックス線であってもよく、フィルタなどを使用して得られる（特性エックス線を含んでもよい）連続エックス線の一部であってもよい。エックス線が、（特性エックス線を含んでもよい）連続Ｘ線の一部である場合には、高い輝度のエックス線を比較的容易に得ることができる。特性エックス線を使用する場合には、十分な輝度のエックス線を得るために、放射光施設によって発生させた高輝度の特性エックス線を用いることが好ましい。

ここで、好ましくは、エックス線の波長は、少なくとも一種類のトレーサーの少なくとも一つの吸収端の波長を含む。トレーサーの吸収端は、トレーサーの材料に固有である。この場合には、その少なくとも一つの吸収端の波長を備えたエックス線が、その少なくとも一つの吸収端を備えた少なくとも一種類のトレーサーに高い効率で吸収される。少なくとも一種類のトレーサーを含む流動性媒体のエックス線による画像において、吸収端の波長を備えたエックス線を透過させる流動性媒体の像と、その吸収端の波長を備えたエックス線を高い効率で吸収するその吸収端を備えたトレーサーの像とのコントラスト（画像の濃度分布のＳ／Ｎ）を向上させることができる。その結果、流動性媒体の流れと共に移動するトレーサーの対応付けの精度を向上させることができる。また、吸収端を備えたトレーサーが、その吸収端の波長を備えたエックス線を高い効率で吸収するので、流動性媒体に照射するその吸収端の波長を含むエックス線の強度を、低減させることができる。なお、少なくとも一種類のトレーサーの少なくとも一つの吸収端の波長を含むエックス線は、その吸収端を備えたトレーサーの材料と同じ材料を、エックス線の発生源のターゲットとして用いればよい。

加えて、流動性媒体の流れを評価することは、流動性媒体の流れを定性的に観察すること及び流動性媒体の流れを定量的に計測することの両方を含む。

このように、本発明による流動性媒体の流れ評価方法等においては、粉体を含む流動性媒体の流れを追跡する、流動性媒体と異なる少なくとも一種類のトレーサーに対して、流動性媒体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波を適用する。このとき、電磁波は、流動性媒体を透過するが、少なくとも一種類のトレーサーに吸収されるので、流動性媒体及び少なくとも一種類のトレーサーと相互作用した電磁波によって、少なくとも一種類のトレーサーの流れを得ることができる。そして、電磁波によって得られた少なくとも一種類のトレーサーの流れを、流動性媒体の流れとして評価することができる。

また、本発明による流動性媒体の流れ評価方法等によれば、少なくとも一種類のトレーサー及び流動性媒体を透過する電磁波を使用することによって、可視光によって内部を観察することができない粉体を含む流動性媒体の流れを評価することができる。例えば、電子写真装置内における有色の現像剤の流れを計測することができる。さらに、流動性媒体の搬送、輸送、攪拌、分散等が行われる装置内における流動性媒体の流れを評価することができる。例えば、電子写真装置における現像剤の現像挙動、電子写真装置の現像器内における現像剤の攪拌及び分散の挙動等を計測することができる。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、少なくとも一種類のトレーサーを含む流動性媒体に電磁波を照射し、電磁波によって形成される少なくとも一種類のトレーサーを含む流動性媒体の画像を取得し、画像を用いて粉体を含む流動性媒体の流れを評価することを含む。よって、本発明による流動性媒体の流れ評価装置は、好ましくは、少なくとも一種類のトレーサーを含む流動性媒体に電磁波を照射する手段、電磁波によって形成される少なくとも一種類のトレーサーを含む流動性媒体の画像を取得する手段、及び画像を用いて粉体を含む流動性媒体の流れを評価する手段を含む。また、本発明による流動性媒体の流れ評価プログラムは、好ましくは、少なくとも一種類のトレーサーを含む流動性媒体に電磁波を照射する手順、電磁波によって形成される少なくとも一種類のトレーサーを含む流動性媒体の画像を取得する手順、及び画像を用いて粉体を含む流動性媒体の流れを評価する手順をコンピュータに実行させることを含む。

上記の電磁波を照射する手段としては、使用する電磁波の波長に依存して、少なくとも電磁波の発生源を含む公知の電磁波を照射する手段を用いることができる。電磁波の発生源としては、電磁波が、紫外線である場合には、低圧水銀ランプ、キセノンランプなどを用いることができる。また、電磁波が、エックス線である場合には、公知のエックス線管などを用いることができる。さらに、電磁波が、ガンマ線である場合には、そのガンマ線を放出する放射性同位体を用いることができる。

画像を取得する手段としては、使用する電磁波の波長に依存して、公知の撮像手段を用いることができる。電磁波が、紫外線である場合には、銀塩フィルム、光電子増倍管などを用いることができる。また、電磁波が、エックス線である場合には、Ｘ線フィルム、イメージングプレート、イメージインテンシファイアーなどを用いることができる。なお、イメージインテンシファイアーで得られた可視化された画像を、銀塩フィルム及びＣＣＤなどに受光してもよい。さらに、電磁波が、ガンマ線である場合には、放射線検出器を用いることができる。

流動性媒体の流れを評価する手段としては、市販のコンピュータなどを用いることができる。また、流動性媒体の流れを評価する手段としてのコンピュータには、本発明による流動性媒体の流れ評価プログラムを実行させてもよい。

この場合には、少なくとも一種類のトレーサーを含む流動性媒体の画像に基づいて、粉体を含む様々な流動性媒体の流れを評価することができる。よって、得られた流動性媒体の画像を目視によって観察するのみならず、得られた流動性媒体の画像を適切に処理することによって、流動性媒体の流れの定量的な評価をすることができる。また、得られた流動性媒体の画像の誤差（歪み、移動誤差など）を修正することによって、流動性媒体の流れを、より正確に評価することができる。

図１は、本発明による流動性媒体の流れ評価装置及び方法の例を説明する図である。図１に示す流動性媒体の流れ評価装置及び方法においては、電磁波としてエックス線を用いる。エックス線を発生させる手段としてのエックス線発生装置１は、エックス線管であり、エックス線管内のターゲット材料を変更することによって、種々の特性エックス線を発生させることができる。ここでは、ターゲット材料として銅（Ｃｕ）を用いる。エックス線発生装置１から発生する特性エックス線（Ｋ線）のエネルギーは、約８．９８ｋｅＶである。

この特性エックス線を流動性媒体３に照射する。この流動性媒体３は、粉体４及びトレーサー５を含む。粉体４は、カラー樹脂トナーなどのような、可視光を全く透過せず、その内部の状態を観察することができない材料からなる。また、トレーサー５は、銅からなる粒子のトレーサーであり、流動性媒体３内で粉体４と共に流れている。粒子のトレーサー５は、粉体４よりも十分小さい。主成分が樹脂である粉体４は、ほぼすべてのエックス線の波長領域における電磁波を透過させる。一方、粒子のトレーサー５は、ほぼすべてのエックス線の波長領域における電磁波を吸収する。よって、流動性媒体３を透過した特性エックス線の像において、粉体４は、特性エックス線の像の背景として撮像され、粒子のトレーサ５は、特性エックス線の像の背景における明確に黒色の像として撮像される。

また、特性エックス線の照射方向において、流動性媒体３に対してエックス線発生装置１の反対側にイメージインテンシファイヤー６が設けられている。流動性媒体３における粒子状のトレーサー５の位置情報を有する流動性媒体３を透過した特性エックス線の像は、イメージインテンシファイヤー６に入力される。イメージインテンシファイヤー６は、内部における電界の強さを制御することができるガラスの真空管を有する。そして、ガラスの真空管の一端には、入射するエックス線に対して電子を放出するでに電子放出器が設けられ、ガラスの真空管の他端には、管内で発生した電子のエネルギーによって蛍光を放出する蛍光板が設けられている。流動性媒体３を透過した特性エックス線の像は、イメージインテンシファイヤー６の電子放出器に入射し、一旦、電子像に変換される。そして、その電子像は、ガラスの真空管内の電界によって集束及び増幅させられ、蛍光板上に投射される。蛍光板に投射された電子像は、蛍光板に塗布された蛍光材によって、可視光の像に変換される。この可視光の像は、回転鏡７等の光学系を通じて、ＣＣＤカメラ８の受光面に導かれる。なお、図１においては、イメージインテンシファイヤー６、回転鏡７、及びＣＣＤカメラ８が、画像を取得する手段に対応する。

さらに、流動性媒体３を透過した特性エックス線に対応する可視光の像は、時系列のデジタル化された画像として、パーソナルコンピュータなどの演算装置９に送られる。演算装置９は、流動性媒体の流れを評価する手段に対応する。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、トレーサーは、粒子であり、流動性媒体の流れと共に移動する粒子の位置の対応付けを行うことを含む。この場合には、トレーサーが、粒子であるため、複数種類のトレーサーを比較的容易に用いることができる。また、粉体のみからなる流動性媒体の流れを比較的容易に評価することができる。ここで、粒子は、粉体（粉末）を含む。また、流動性媒体の流れと共に移動する粒子の位置の対応付けを行うため、時間と共に変動する流動性媒体の流れの情報を得ることができる。より具体的には、複数の時刻で得られる少なくとも一種類の粒子のトレーサーを含む流動性媒体の複数の画像を取得し、複数の画像の間における粒子のトレーサーの位置を対応付ける。複数の画像の間で対応付けられた粒子のトレーサーの位置における変動から、粒子のトレーサーの流れと実質的に一致する流動性媒体の流れを評価することができる。

複数の画像の間における粒子のトレーサーの位置の対応付けは、例えば、公知のＰＴＶ（Ｐａｒｔｉｃｌｅ Ｔｒａｃｋｉｎｇ Ｖｅｌｏｃｉｍｅｔｒｙ）の方法によって達成される。

図２は、複数の画像の間における粒子のトレーサーの位置の対応付けを説明する図である。まず、複数の時刻における粒子のトレーサーを含む流動性媒体の複数の画像を取得する。図２に示すように、二つの画像Ｉｍａｇｅ Ａ_ｎ及びＩｍａｇｅ Ａ_ｎ＋１は、所定の時間間隔だけ異なる二つの時刻に得られた粒子のトレーサを含む流動性媒体の画像である。ある時刻での画像Ｉｍａｇｅ Ａ_ｎにおける一つの粒子のトレーサーに対応する、別の時刻での画像Ｉｍａｇｅ Ａ_ｎ＋１における一つの粒子のトレーサーを探索し、ある時刻での画像Ｉｍａｇｅ Ａ_ｎにおける一つの粒子のトレーサーを、別の時刻での画像Ｉｍａｇｅ Ａ_ｎ＋１における一つの粒子のトレーサーに対応付ける。その結果、所定の時間間隔における一つの粒子のトレーサーの位置変動の情報を得ることができる。同様に、この操作を、二つの時刻に得られた複数の画像における複数の粒子のトレーサーの各々に対して行うことによって、所定の時間間隔における複数の粒子のトレーサーの流れの位置変動の情報を得ることができ、結果として、流動性媒体の流れを評価することができる。ある時刻での画像における複数の粒子のトレーサーと、別の時刻での画像における複数の粒子のトレーサーとの対応付けには、様々なアルゴリズムが用いられる。

まず、複数の粒子のトレーサーの対応付けに関する前処理として、トレーサーを含む流動性媒体の画像に含まれる複数の粒子のトレーサーの像について、粒子分離処理を行う。図３は、トレーサーを含む流動性媒体の画像に含まれる複数の粒子のトレーサーの像に対する粒子分離処理を説明する図である。図３に示すように、トレーサーを含む流動性媒体の画像には、複数の粒子のトレーサーの像が、二次元的に互いに重なり合っている場合がある。このような互いに重なり合った複数の粒子のトレーサーの像を、粒子のトレーサーの形状に基づいて、複数の個別の粒子のトレーサーの像に分離する。例えば、図３に示すように、複数の粒子のトレーサーが、球形の形状を有するとすれば、互いに重なり合った粒子のトレーサーの像を、別個の球形のトレーサーの像に分離することができる。

図４は、ＰＴＶによる複数の画像の間における粒子のトレーサーの対応付けに関するアルゴリズムの例を説明する図である。図４に示すように、二つの画像Ｉｍａｇｅ Ｂ_ｎ及びＩｍａｇｅ Ｂ_ｎ＋１は、所定の時間間隔だけ異なる二つの時刻に得られた粒子のトレーサを含む流動性媒体の画像である。まず、ある時刻での画像Ｉｍａｇｅ Ｂ_ｎにおける一つの粒子のトレーサーＰｘに注目する。ある時刻での画像Ｉｍａｇｅ Ｂ_ｎにおける一つの粒子のトレーサーＰｘの周囲における、粒子のトレーサーＰｘとの距離が短い四つの粒子のトレーサーＰ１、Ｐ２、Ｐ３、Ｐ４を探索する。ここで、粒子のトレーサーＰｘとの距離は、探索される粒子のトレーサＰ１、Ｐ２、Ｐ３、Ｐ４の重心と粒子のトレーサーＰｘの重心との距離である。次に、別の時刻での画像Ｉｍａｇｅ Ｂ_ｎ＋１において、粒子のトレーサーＰｘに対応する粒子のトレーサの候補として、一つの粒子のトレーサーＰｘ’を選択する。そして、一つの粒子のトレーサーＰｘ’を中心とした所定の探索領域に含まれる、一つの粒子のトレーサーＰｘ’との距離が短い四つの粒子のトレーサーＰ１’、Ｐ２’、Ｐ３’、Ｐ４’を探索する。なお、所定の探索領域は、任意に設定することができる。例えば、流動性媒体の流れが、層流であり、その流れの乱れの最大幅が、Ｒであるときには、一つの粒子のトレーサーＰｘ’の位置を中心とした半径Ｒの円であってもよい。次に、式

に従って、数値ＥＲを計算する。ここで、ＰｘＰｉは、Ｐｘの重心とＰｉ（ｉ＝１，２，３，４）の重心との間の距離であり、Ｐｘ’Ｐｊ’は、Ｐｘ’の重心とＰｊ’（ｊ＝１，２，３，４）の重心との間の距離である。そして、別の時刻での画像Ｉｍａｇｅ Ｂ_ｎ＋１における探索領域内に含まれる全ての候補粒子について数値ＥＲを計算し、最小のＥＲを与える粒子のトレーサーを、粒子のトレーサーＰｘに対応する粒子のトレーサーとする。同様に、この操作を、画像Ｉｍａｇｅ Ｂ_ｎにおける全ての粒子のトレーサーについて行う。

このアルゴリズムは、二つの画像が取得された時間間隔が十分に小さく、一つの画像における複数の粒子のトレーサーにおける相対的な位置が、ほとんど変化しないという仮定のもとで有効である。二つの時刻に得られた二つの画像における粒子のトレーサーが、対応付けられると、二つの画像にの間における粒子のトレーサーの位置に関する変位ベクトルが得られる。そして、この変位ベクトルを、二つの画像が取得された時間間隔で割ることによって、粒子のトレーサーの速度ベクトルを得ることができる。すなわち、その粒子のトレーサーにおける流動性媒体の流れの速度を求めることができる。異なる時刻に取得されった複数の画像における必要な複数の粒子のトレーサーについて速度ベクトルを求めることで、それらの複数の粒子のトレーサーを含む流動性媒体の流れの速度場を得ることができる。

本発明による流動性媒体の流れ評価方法等においては、トレーサーは、液体又は粒子集合体であり、流動性媒体の流れと共に移動する液体又は粒子集合体の対応付けを行うことを含む。この場合には、トレーサーとして、液体又は粒子集合体を用いても、流動性媒体の流れと共に移動する液体又は粒子集合体の対応付けを行うため、時間と共に変動する流動性媒体の流れの情報を得ることができる。

なお、トレーサーとしての液体は、一種類の液体であることが好ましいが、互いに混合しない複数の液体を用いることも可能である。液体のトレーサーとしては、例えば、水銀などが挙げられる。また、粒子集合体は、トレーサーを含む流動性媒体の画像の分解能よりも小さい粒子（粉末）の集合体であり、擬似的に液状に観察される。粒子集合体のトレーサーとしては、例えば、現在のマイクロフォーカスエックス線観察装置の実用分解能が、１μｍ程度であるため、１μｍ以下の（平均）直径を備えた鉄の粉末などが挙げられる。

また、流動性媒体の流れと共に移動する液体又は粒子集合体の位置の対応付けを行うため、時間と共に変動する流動性媒体の流れの情報を得ることができる。より具体的には、複数の時刻で得られる少なくとも一種類の液体又は粒子集合体のトレーサーを含む流動性媒体の複数の画像を取得し、複数の画像の間における液体又は粒子集合体のトレーサーを対応付ける。複数の画像の間で対応付けられた液体又は粒子集合体のトレーサーの位置における変動から、液体又は粒子集合体のトレーサーの流れと実質的に一致する流動性媒体の流れを評価することができる。

複数の画像の間における液体又は粒子集合体のトレーサーの対応付けは、例えば、公知のＰＩＶ（Ｐａｒｔｉｃｌｅ Ｉｍａｇｉｎｇ Ｖｅｌｏｃｉｍｅｔｒｙ）の方法によって達成される。

図５は、本発明による流動性媒体の流れ評価装置及び方法の別の例を説明する図である。図５に示す流動性媒体の流れ評価装置及び方法は、図１に示す流動性媒体の流れ評価装置及び方法と同様であるが、トレーサーが、銅からなる粒子のトレーサー５ではなく、水銀からなる液体のトレーサー１３である点で図１に示す流動性媒体の流れ評価装置及び方法と異なっている。図５において、水銀からなる液体のトレーサー１３もまた、ほぼすべてのエックス線の波長領域における電磁波を吸収する。よって、流動性媒体３を透過した特性エックス線の像において、粉体４は、特性エックス線の像の背景として撮像され、液体のトレーサ１３は、特性エックス線の像の背景における明確に黒色の像として撮像される。

図６は、ＰＩＶによる複数の画像の間における液体のトレーサーの対応付けに関するアルゴリズムの例を説明する図である。図６に示すように、二つの画像Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂは、所定の時間間隔だけ異なる二つの時刻に得られた液体のトレーサを含む流動性媒体の画像である。まず、ある時刻での画像Ｉｍａｇｅ Ａ及び別の時刻での画像Ｉｍａｇｅ Ｂにおける所定の範囲の領域を、それぞれ、横方向にｎ個の領域に分割し、縦方向にｍ個の領域に分割する。すなわち、画像Ｉｍａｇｅ Ａ及び画像Ｉｍａｇｅ Ｂにおける所定の範囲の領域は、ｎ×ｍ個の領域に分割される。画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊ（ｉ＝１，２，…ｎ：ｊ＝１，２，…ｍ）及び画像Ｉｍａｇｅ Ｂにおける分割された領域の一つＱ_ｉｊ（ｉ＝１，２，…ｎ：ｊ＝１，２，…ｍ）の大きさは、それぞれ、流動性媒体の流れの最小単位程度の大きさであることが好ましいが、十分な数の画素を含むことも必要である。そして、ある時刻での画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊに注目する。次に、流体のトレーサーの流れの特性を考慮して、Ｉｍａｇｅ Ｂ中に探索領域を設定する。次に、画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊを、画像Ｉｍａｇｅ Ｂにおける分割された領域の一つＱ_ｉｊに対応付ける。具体的には、画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊの画像の濃度Ｐ_ｉｊ（ｉ，ｊ）及び画像Ｉｍａｇｅ Ｂにおける任意の分割された領域の一つＱ_ｉｊの画像の濃度Ｑ_ｉｊ（ｉ，ｊ）を取得し、式

に従って、数値ＥＲを計算する。そして、別の時刻での画像Ｉｍａｇｅ Ｂにおける探索領域内に含まれる全ての分割された領域Ｑ_ｉｊについて数値ＥＲを計算し、最小のＥＲを与える画像Ｉｍａｇｅ Ｂにおける分割された領域を、画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊに対応する画像Ｉｍａｇｅ Ｂにおける分割された領域とする。このアルゴリズムは、二つの画像が取得された時間間隔が十分に小さく、一つの画像における液体のトレーサーにおける形状及び濃度分布が、ほとんど変化しないという仮定のもとで有効である。二つの時刻に得られた二つの画像における液体のトレーサーが、対応付けられると、画像Ｉｍａｇｅ Ａにおける分割された領域の一つＰ_ｉｊの重心と、Ｐ_ｉｊに対応付けられた画像Ｉｍａｇｅ Ｂにおける分割された領域の一つＱ_ｉｊの重心を結ぶ、液体のトレーサーにおける一つの領域Ｐ_ｉｊの変位ベクトルが得られる。そして、この変位ベクトルを、二つの画像が取得された時間間隔で割ることによって、液体のトレーサーにおける一つの領域Ｐ_ｉｊの重心における速度ベクトルを得ることができる。すなわち、その液体のトレーサーにおける一つの領域Ｐ_ｉｊの重心における流動性媒体の流れの速度を求めることができる。この操作を、画像Ｉｍａｇｅ Ａにおける分割された領域の全てについて行うと、液体のトレーサーが存在する流動性媒体の領域における流動性媒体の流れの速度場を得ることができる。さらに、異なる時刻に取得された複数の画像における液体のトレーサーについて速度場を求めることで、異なる時刻における液体のトレーサーを含む流動性媒体の流れの速度場を得ることができる。なお、異なる時刻に取得された複数の画像を任意の形状に分割することが可能である。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、少なくとも一種類のトレーサーは、複数種類のトレーサーである。この場合には、電磁波が、複数種類のトレーサーによって吸収される。ここで、トレーサーによって吸収される電磁波の吸収率は、一般的に、トレーサーの種類によって異なる。よって、複数種類のトレーサーに電磁波を適用した場合、トレーサーの種類によって異なるトレーサーによる電磁波の吸収率の差によって、流動性媒体に含まれるトレーサーの種類を区別することができる。例えば、複数種類のトレーサーを含む流動性媒体の画像を取得した場合には、トレーサーの種類によって、トレーサーの像の濃度が異なる。流動性媒体に含まれるトレーサーの種類を区別することができるため、異なる時間に取得された複数種類のトレーサーを含む流動性媒体の画像の間において、トレーサーの総数を変更することなく、得られる流動性媒体の流れに関する情報の数を維持し、トレーサーの対応付けをより容易にすると共により正確にすることができる。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、電磁波は、互いに異なる波長を備えた複数の電磁波を含む。この場合には、互いに異なる波長を備えた複数の電磁波が、少なくとも一種類のトレーサーによって吸収される。ここで、トレーサーによって吸収される電磁波の吸収率は、一般的に、電磁波の波長によって異なる。よって、互いに異なる波長を備えた複数の電磁波をトレーサーに適用した場合、電磁波の波長によって、トレーサーによる電磁波の吸収を調整することができる。特に、電磁波の波長が、トレーサーの吸収端を含むときには、その吸収端を備えたトレーサーの吸収を向上させることができる。例えば、少なくとも一種類のトレーサーを含む流動性媒体の画像を取得した場合には、電磁波の波長によって、トレーサーの像の濃度が異なる。

特に、本発明による流動性媒体の流れ評価方法等においては、望ましくは、少なくとも一種類のトレーサーは、複数種類のトレーサーであり、且つ、電磁波は、互いに異なる波長を備えた複数の電磁波を含む。この場合には、互いに異なる波長を備えた複数の電磁波を複数種類のトレーサーに適用することで、複数種類のトレーサーにおける特定のトレーサーによる電磁波の吸収を選択的に調整することができる。

さらに、互いに異なる波長の複数の電磁波の一部又は全てが、複数種類のトレーサーにおける特定のトレーサーの吸収端を含むことが望ましい。例えば、互いに異なる波長の複数の電磁波の一部が、複数種類のトレーサーにおける特定のトレーサーの吸収端を含み、且つ、互いに異なる波長の複数の電磁波の残りが、複数種類のトレーサーにおける特定のトレーサーの吸収端を含まなくてもよい。また、互いに異なる波長の複数の電磁波の全部が、それぞれ、複数種類のトレーサーにおける別個のトレーサーの吸収端を含んでもよい。これらの場合には、複数種類のトレーサーにおける特定のトレーサーによる電磁波の吸収を選択的に向上させることができる。すなわち、流動性媒体に含まれるトレーサーの種類を、より明確に区別することができるため、異なる時間に取得された複数種類のトレーサーを含む流動性媒体の画像の間において、トレーサーの総数を変更することなく、得られる流動性媒体の流れに関する情報の数を維持し、トレーサーの対応付けをさらに容易にすると共にさらに正確にすることができる。

例えば、二種類のトレーサー、それらトレーサーの一方の吸収端を含む電磁波及びそれらトレーサーの両方の吸収端を含まない電磁波を用いてもよい。それらトレーサーの一方の吸収端を含む電磁波を、二種類のトレーサーに照射すれば、その吸収端を備えたトレーサーの吸収のみを選択的に向上させることができ、二種類のトレーサーを明確に区別することができる。その結果、トレーサーの種類毎にトレーサーの対応付けを容易に行うことができ、トレーサーの対応付けの精度を向上させることができる。

同様に、二種類のトレーサー、それらトレーサーの一方の吸収端を含む電磁波及びそれらトレーサーの他方の吸収端を含む電磁波を用いてもよい。それらの電磁波を、二種類のトレーサーに別個に照射すれば、照射した電磁波に対応する吸収端を備えたトレーサーの吸収のみを選択的に向上させることができ、二種類のトレーサーを明確に区別することができる。その結果、トレーサーの種類毎にトレーサーの対応付けを容易に行うことができ、トレーサーの対応付けの精度を向上させることができる。

なお、互いに異なる波長を備えた複数の電磁波を、単一の電磁波の発生源から選択的に発生させてもよい。互いに異なる波長を備えた複数の電磁波を、単一の電磁波の発生源から選択的に発生させるためには、例えば、単一の電磁波の発生源におけるターゲットの種類を変更すればよい。この場合には、互いに異なる波長を備えた複数の電磁波を発生させるために、複数の電磁波の発生源を必要とせず、流動性媒体の流れ評価装置を小型化することができる。

また、互いに異なる波長を備えた複数の電磁波を、複数の電磁波の発生源から発生させてもよい。この場合には、互いに異なる波長を備えた複数の電磁波を発生させるために、単一の電磁波の発生源から発生する電磁波を選択するための操作をする必要が無い。その結果、複数の電磁波の発生源による電磁波の発生を切り替えるだけで、互いに異なる波長を備えた複数の電磁波を、相対的に短い時間間隔で発生させることができる。

図７は、複数種類のトレーサーの材料としての金属によって吸収されるエックス線の質量吸収係数を説明する図である。図７に示すグラフにおいて、横軸は、トレーサーに照射されるエックス線のエネルギー（ｋｅＶ）の対数であり、縦軸は、エックス線に対するトレーサーの質量吸収係数（ｃｍ^２／ｇ）の対数である。なお、質量吸収係数は、エックス線に対するトレーサーの線減弱係数（１／ｃｍ）／トレーサーの密度（ｇ／ｃｍ^３）に等しい。また、トレーサーの材料としての金属は、鉄Ｆｅ、タングステンＷ、鉛Ｐｂ、又は銀Ａｇである。図７に示すように、金属によるエックス線の吸収は、エックス線の波長（又はエネルギー）によって異なり、金属の種類によって、エックス線の吸収端が異なる。このように、複数種類の金属トレーサーを用いると共に特定の金属トレーサーの吸収端を含む波長を備えたエックス線を金属トレーサーに照射することで、複数種類の金属のトレーサーにおいてその特定の金属トレーサーの吸収を選択的に向上させることができる。その結果、複数種類の金属トレーサーにおいてその特定の金属トレーサーを明確に区別することができる。なお、複数種類の金属トレーサーを含む流動性媒体の画像においては、その特定の金属トレーサーの濃度を、他の金属トレーサーの濃度に対して大きく変化させることになる。その結果、異なる時刻に取得された画像における金属トレーサーの総数を増加させても、特定の金属トレーサーのみの対応付けを行うことによって、金属トレーサーの対応付けの誤差は、増加せず、金属トレーサーの対応付けを容易に且つ正確にすることができる。

図８（ａ）は、複数種類の直径１００μｍの球形粒子である複数種類のトレーサーによって吸収されるエックス線の相対吸収率を説明する図である。図８（ａ）に示すグラフにおいて、横軸は、トレーサーに照射されるエックス線のエネルギー（ｋｅＶ）の対数であり、縦軸は、直径１００μｍの球形粒子のトレーサーによって吸収されるエックス線の相対吸収率の対数である。また、トレーサーとしての直径１００μｍの球形粒子の金属材料は、鉄Ｆｅ、タングステンＷ、鉛Ｐｂ、又は銀Ａｇである。図８（ａ）に示すグラフにおいて、縦軸の相対吸収率の基準については、鉛Ｐｂの直径１００μｍの球形粒子によって吸収される１０ｋｅＶのエックス線の吸収量を１とした。

図８（ｂ）は、複数種類のトレーサーとしての複数種類の直径１００μｍの球形粒子の中心を通過するエックス線の相対吸収率を説明する図である。図８（ｂ）に示すグラフにおいて、横軸は、トレーサーに照射されるエックス線のエネルギー（ｋｅＶ）であり、縦軸は、トレーサーとしての直径１００μｍの球形粒子の中心を通過するエックス線の相対吸収率の対数である。また、トレーサーとしての直径１００μｍの球形粒子の金属材料は、鉄Ｆｅ、タングステンＷ、鉛Ｐｂ、銀Ａｇである。図８（ｂ）に示すグラフにおいて、縦軸の相対吸収率の基準については、鉛Ｐｂの直径１００μｍの球形粒子によって吸収される２０ｋｅＶのエックス線の吸収量を１とした。

図８（ａ）及び（ｂ）に示すように、金属によるエックス線の吸収は、エックス線の波長（又はエネルギー）によって異なり、金属の種類によって、エックス線の吸収端が異なる。このように、複数種類の金属トレーサーとしての直径１００μｍの球形粒子を用いると共に特定の金属からなる直径１００μｍの球形粒子の吸収端を含む波長を備えたエックス線を、複数種類の金属トレーサーとしての直径１００μｍの球形粒子に照射することで、複数種類の金属トレーサーとしての直径１００μｍの球形粒子において、その特定の金属からなる直径１００μｍの球形粒子の吸収を選択的に向上させることができる。その結果、複数種類の金属のトレーサーとしての直径１００μｍの球形粒子において、その特定の金属からなる直径１００μｍの球形粒子を明確に区別することができる。なお、複数種類の金属トレーサーとしての直径１００μｍの球形粒子を含む流動性媒体の画像においては、その特定の金属からなる直径１００μｍの球形粒子の濃度を、他の金属からなる直径１００μｍの球形粒子の濃度に対して大きく変化させることになる。その結果、異なる時刻に取得された画像における金属トレーサーとしての直径１００μｍの球形粒子の総数を増加させても、特定の金属からなる直径１００μｍの球形粒子のみの対応付けを行うことによって、直径１００μｍの球形粒子の対応付けの誤差は、増加せず、直径１００μｍの球形粒子の対応付けを容易に且つ正確にすることができる。

このように、エックス線の波長（又はエネルギー）を、複数種類のトレーサーにおける特定の金属材料からなるトレーサーの吸収端付近で変化させることによって、複数種類のトレーサーを含む流動性媒体の画像における特定の金属材料からなるトレーサーの濃度を大きく変化させることができる。そして、流動性媒体の画像において、この特定の金属材料からなるトレーサーの濃度と大きく異なる濃度を備えた他の金属材料のトレーサーによる像を消去して、特定の金属材料からなるトレーサーの像のみを残す画像処理を行うことは、非常に容易である。例えば、鉄Ｆｅ、タングステンＷ、鉛Ｐｂ、又は銀Ａｇからなる四種類のトレーサーを流動性媒体に混入し、流動性媒体に照射するエックス線の波長（又はエネルギー）を適切に変化させて、異なる時刻で四種類のトレーサーを含む流動性媒体の複数の画像を得る。そして、複数の画像の間で流動性媒体におけるトレーサーの対応付けを行うときには、画像におけるトレーサーの濃度差を用いた画像処理によって、特定の金属からなるトレーサー以外のトレーサーの像を消去し、特定のトレーサーの像のみを残す。得られた特定のトレーサーの像についてのみ、対応付けを行う。結果として、流動性媒体の流れの情報の総数に対応する、画像における四種類のトレーサーの総数を増加させても、より少ない数の特定のトレーサーの像に関する対応付けによって、トレーサーの像の対応付けの誤差は、増加せず、トレーサーの像の対応付けを容易に且つ正確にすることができる。

図９（ａ）及び（ｂ）は、互いに異なるエネルギーを備えたエックス線を照射することによって得られた二種類の金属トレーサーを含む流動性媒体の画像における対応付けを説明する図である。図９（ａ）は、二種類の金属トレーサーを含む流動性媒体に６０ｋｅＶのエネルギーを備えたエックス線を照射することによって得られた異なる時刻における二つの画像を示す。図９（ｂ）は、二種類の金属トレーサーを含む流動性媒体に８０ｋｅＶのエネルギーを備えたエックス線を照射すると共に相対的な濃度値０．５以上のトレーサー像を抽出することによって得られた異なる時刻における二つの画像を示す。図９（ａ）及び（ｂ）において、二種類の金属トレーサーの材料は、銀（Ａｇ）及びタングステン（Ｗ）である。『Ｗの１００μｍの吸収率を１』と正規化した場合、銀の１００μｍの吸収率は、６０ｋｅＶにおいて１．１３（Ｗとほぼ同じ）、８０ｋｅＶにおいて０．３１６（Ｗの半分以下）である。よって、図９（ａ）に示すように、６０ｋｅＶのエネルギーを備えたエックス線を照射した場合には、銀のトレーサーによって吸収される６０ｋｅＶのエネルギーを備えたエックス線の吸収率と、タングステンのトレーサーによって吸収される６０ｋｅＶのエネルギーを備えたエックス線の吸収率がほぼ同じであるため、銀及びタングステンのトレーサーの両方は、画像において同程度の濃度を示す。よって、銀及びタングステンのトレーサーを正確に区別することは、容易ではなく、銀及びタングステンのトレーサーの両方を含むトレーサーの全部に対して対応付けをする必要があり、対応付けの誤差は、大きくなる。これに対して、図９（ｂ）に示すように、８０ｋｅＶのエネルギーを備えたエックス線を照射した場合には、銀のトレーサーによって吸収される８０ｋｅＶのエネルギーを備えたエックス線の吸収率と、タングステンのトレーサーによって吸収される８０ｋｅＶのエネルギーを備えたエックス線の吸収率が異なるため、銀及びタングステンのトレーサーの両方は、画像において異なる濃度を示す。すなわち、銀のトレーサーは、０．５を超える相対的な濃度値を示す一方で、タングステンのトレーサーが、０．５未満の相対的な濃度値を示す。よって、画像において銀及びタングステンのトレーサーを容易に且つ正確に区別することができる。そして、画像において０．５以上の相対的な濃度値を示すトレーサーの像のみを、適当な公知の画像処理方法を用いて選択する（画像において０．５未満の相対的な濃度値を示すトレーサーを消去する）。すなわち、画像において０．５を超える相対的な濃度値を示す銀のトレーサーが、抽出される。その結果、銀のトレーサーに対してのみ対応付けをすることができるため、対応付けの誤差を、小さくすることができる。逆に、画像において０．５以下の相対的な濃度値を示すトレーサーの像のみを、適当な公知の画像処理方法を用いて選択すれば、タングステンのトレーサーに対しても対応付けをすることができる。

なお、銀及びタングステンのトレーサーを含む流動性媒体が、他の金属（例えば、Ｐｂ及び／又はＦｅ）のトレーサーをさらに含む場合にも、互いに異なるエネルギーを備えた複数のエックス線を照射することによって得られるトレーサーの像の濃度差を用いて、トレーサーの種類を区別することができる。

本発明による流動性媒体の流れ評価方法等においては、好ましくは、互いに異なる波長を備えた複数の電磁波は、それぞれ、異なる時間に少なくとも一種類のトレーサーを含む流動性媒体に照射される。この場合には、互いに異なる波長を備えた複数の電磁波が、それぞれ、異なる時間に少なくとも一種類のトレーサーによって吸収される。ここで、トレーサーによって吸収される電磁波の吸収率は、一般的に、電磁波の波長によって異なる。よって、互いに異なる波長を備えた複数の電磁波を、それぞれ、異なる時間にトレーサーに適用した場合、異なる時間に適用される電磁波の波長によって、トレーサーによる電磁波の吸収を調整することができる。特に、複数の電磁波における一部の電磁波の波長が、トレーサーの吸収端を含むときには、その吸収端を含む波長を備えた電磁波及びその吸収端を含まない波長を備えた電磁波が、異なる時間に、その吸収端を備えたトレーサーに適用される。その結果、その吸収端を含む波長を備えた電磁波が吸収されるときには、その吸収端を備えたトレーサーによる電磁波の吸収は、向上する。一方、その吸収端を含まない波長を備えた電磁波が吸収されるときには、その吸収端を備えたトレーサーによる電磁波の吸収は、低減される。よって、複数の電磁波が、適用される間に、その吸収端を備えたトレーサーによる電磁波の吸収を最適化することが可能となる。例えば、少なくとも一種類のトレーサーを含む流動性媒体の画像を取得した場合には、異なる時間に照射される電磁波の波長によって、トレーサーの像の濃度を、最適化することが可能となる。

特に、本発明による流動性媒体の流れ評価方法等においては、望ましくは、少なくとも一種類のトレーサーは、複数種類のトレーサーであり、且つ、互いに異なる波長を備えた複数の電磁波は、それぞれ、異なる時間に少なくとも複数種類のトレーサーを含む流動性媒体に照射される。この場合には、互いに異なる波長を備えた複数の電磁波を、異なる時間に複数種類のトレーサーに適用することで、複数種類のトレーサーにおける各々のトレーサーによる電磁波の吸収を、異なる時間に、選択的に最適化することが可能となる。例えば、複数種類のトレーサーを含む流動性媒体の画像を取得した場合には、異なる時間に照射される電磁波の波長によって、複数種類のトレーサーにおける各々のトレーサーの像の濃度を、異なる時間に、選択的に最適化することが可能となる。

さらに、互いに異なる波長の複数の電磁波の一部又は全てが、複数種類のトレーサーにおける特定のトレーサーの吸収端を含むことが望ましい。例えば、互いに異なる波長の複数の電磁波の一部が、複数種類のトレーサーにおける特定のトレーサーの吸収端を含み、且つ、互いに異なる波長の複数の電磁波の残りが、複数種類のトレーサーにおける特定のトレーサーの吸収端を含まなくてもよい。また、互いに異なる波長の複数の電磁波の全部が、それぞれ、複数種類のトレーサーにおける別個のトレーサーの吸収端を含んでもよい。これらの場合には、複数種類のトレーサーにおける特定のトレーサーによる電磁波の吸収を、異なる時間に、十分に最適化することが可能となる。例えば、複数種類のトレーサーを含む流動性媒体の画像を取得した場合には、異なる時間に照射される電磁波の波長によって、複数種類のトレーサーにおける特定のトレーサー像の濃度を、異なる時間に、十分に最適化することが可能となる。

さらに好ましくは、互いに異なる波長を備えた複数の電磁波は、それぞれ異なる時間に且つ少なくとも一回以上繰り返して、少なくとも一種類のトレーサーを含む流動性媒体に照射される。ここで、一回繰り返すことは、互いに異なる波長を備えた複数の電磁波の各々を、一回ずつ、少なくとも一種類のトレーサーを含む流動性媒体に適用することを意味する。この場合には、互いに異なる波長を備えた複数の電磁波の各々が、異なる時間に少なくとも一種類のトレーサーを含む流動性媒体に適用されると共に、互いに異なる波長を備えた複数の電磁波が、少なくとも一回以上繰り返して、少なくとも一種類のトレーサーを含む流動性媒体に適用される。よって、例えば、少なくとも一種類のトレーサーを含む流動性媒体の画像を取得した場合には、互いに異なる波長を備えた複数の電磁波の各々一つを少なくとも一種類のトレーサーに照射して得られる複数の画像を、異なる時間に得ることができる。その結果、互いに異なる波長を備えた複数の電磁波の各々一つを吸収する少なくとも一種類のトレーサーの流れに対応する流動性媒体の流れを評価することができる。

互いに異なる波長を備えた複数の電磁波が、エックス線である場合、少なくとも一回以上繰り返してエックス線を放出する手段としては、公知の様々な手段を用いることができる。例えば、エックス線のターゲットから放出されるエックス線を、鉛からなる物体によって繰り返し機械的に遮断する方法、溝が設けられたエックス線のターゲットを回転させる方法、エックス線を反射する反射膜を振動させてエックス線の照射方向を制御する方法、エックス線を反射する結晶の反射膜に電圧を印加して反射膜の結晶の原子間距離を変化させる方法などが挙げられる。

図１０は、本発明による流動性媒体の流れ評価装置及び方法のまた別の例を説明する図である。図１０に示す流動性媒体の流れ評価装置及び方法は、図１に示す流動性媒体の流れ評価装置及び方法と同様であるが、複数のエックス線発生装置及び複数種類のトレーサーを採用すると共に複数のエックス線発生装置から放出されるエックス線を遮断することが可能なエックス線遮断機構を用いる点で、図１に示す流動性媒体の流れ評価装置及び方法と異なっている。

図１０において、流動性媒体３は、粒子のトレーサー５として、Ｆｅ、Ａｇ、Ｐｂ又はＷからなる四種類の直径１００μｍの粒子を含む。また、三つのエックス線発生装置１Ａ、１Ｂ、１Ｃは、それぞれ、ターゲットとしてＷ、Ａｇ、Ｐｂを用いるエックス線管である。これら三つのエックス線発生装置１Ａ、１Ｂ、１Ｃにおけるターゲットの大きさは、流動性媒体３に対して十分に小さく且つ十分に遠方にあるので、三つのエックス線発生装置１Ａ、１Ｂ、１Ｃから放出されるエックス線は、概略平行光として、流動性媒体３に照射される。ターゲットとしてＷ、Ａｇ、Ｐｂを用いたエックス線管であるエックス線発生装置１Ａ、１Ｂ、１Ｃに十分に大きい管電圧を印加すると、エックス線発生装置１Ａ、１Ｂ、１Ｃにおけるターゲットの各々から、特性エックス線を含む連続エックス線が発生する。そして、三つのエックス線発生装置１Ａ、１Ｂ、１Ｃから発生する各Ｗ、Ｐｂ、ＡｇのＫ吸収端の波長を備えた特性エックス線を、Ｓｉ結晶等を用いたモノクロメータによって、連続エックス線から分光し、それらの特性エックス線を、流動性媒体３に照射する。

また、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する特性エックス線を遮断する又は通過させるエックス線遮断機構１１は、三つの穴を有する十分原子番号の大きい材料（例えば、鉛）からなる遮蔽板を有している。遮蔽板は、ソレノイド機構を用いて一定の時間間隔（例えば、１／６０秒、１／３０秒）で振動しており、遮蔽板及び遮蔽板の穴は、それぞれ、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する特性エックス線を、遮蔽板の振動の時間間隔で、遮蔽する及び通過させる。また、遮蔽板の振動に従って、遮蔽板の三つの穴は、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する特性エックス線の一つを通過させると共に残りの二つを遮断するように、遮蔽板に形成される。なお、エックス線遮断機構１１の振動は、演算装置９によって制御される。

図１１は、図１０に示す流動性媒体の流れ評価方法等における特性エックス線の照射及び特性エックス線による撮像の時間スケジュールを示す図である。図１１の横軸は、時間を表す。図１１に示すように、三つのエックス線発生装置１Ａ、１Ｂ、１Ｃから発生する各Ｗ、Ｐｂ、ＡｇのＫ吸収端（Ｗ：６９．５ｋｅＶ、Ｐｂ：８８．０ｋｅＶ、Ａｇ：２５．５ｋｅＶ）の波長を備えた特性エックス線は、交互に繰り返して照射される。また、各特性エックス線は、繰り返しで、流動性媒体３に１／６０秒間照射され、１／３０秒間遮断される。一方、ＣＣＤカメラ８は、１／１２０秒の一定の時間間隔で、流動性媒体３の画像を撮像する。なお、ＣＣＤカメラ８による流動性媒体３の撮像のタイミングは、流動性媒体３に対する特性エックス線の照射のタイミングと同期している。

また、図１２（ａ）、（ｂ）及び（ｃ）は、図１０に示す流動性媒体の流れ評価方法を用いて、図１１に示す時間スケジュールに従って撮像された粒子のトレーサーの画像を示す図である。図１２における数字は、画像における対応するトレーサーの像の吸収率である。画像として得られる各トレーサーの濃度は、装置の条件によって異なるが、透過率（＝１−正規化された吸収率）と一般的に表すことができる。ここでは、吸収率の差によって、各トレーサーを判定する概念について説明する。図１１に示す時間スケジュールに従って撮像された粒子のトレーサーの画像を、それぞれ、Ｎ＝ｎ０，ｎ１，ｎ２，…と番号付けると、Ｎ＝３ｍ（ｍ：整数）番目の画像におけるＦｅ、Ａｇ、Ｐｂ、及びＷのトレーサーの吸収率を図１２（ａ）に示す。また、Ｎ＝３ｍ＋１番目におけるＦｅ、Ａｇ、Ｐｂ、及びＷのトレーサーの吸収率を図１２（ｂ）に示し、Ｎ＝３ｍ＋２番目におけるＦｅ、Ａｇ、Ｐｂ、及びＷのトレーサーの吸収率を図１２（ｃ）に示す。

まず、図１２（ａ）に示すように、Ｎ＝３ｍ番目の画像においては、Ｗのトレーサーの吸収率が、他の金属材料のトレーサーの吸収率よりも十分に大きい。よって、Ｎ＝３ｍ番目の画像においては、Ｗのトレーサーを、他の金属材料のトレーサーと容易に区別することができる。そして、Ｎ＝３ｍ番目の画像におけるＷのトレーサーの像の流れを、対応付けによって評価することで、Ｗのトレーサーの位置に対応する流動性媒体の流れを評価することができる。

次に、図１２（ｂ）に示すように、Ｎ＝３ｍ＋１番目の画像においては、ＰｂとＷのトレーサーを、他の金属材料のトレーサーと容易に区別することができる。また、Ｎ＝３ｍの流れでＷのトレーサーが特定できるので、その情報をあわせて、Ｗのトレーサーを特定することができる。よって、そして、Ｎ＝３ｍ＋１番目の画像におけるＷのトレーサーの像の流れを、対応付けによって評価することで、Ｗのトレーサーの位置に対応する流動性媒体の流れを評価することができる。

最後に、図１２（ｃ）に示すように、Ｎ＝３ｍ＋２番目の画像においては、特にＦｅのトレーサーは、それぞれ、他より低い吸収率を有する。従って、Ｎ＝３ｍ、Ｎ＝３ｍ＋１番目の画像から得られた情報により、Ｆｅのトレーサー及びＡｇのトレーサーを識別することができる。そして、Ｎ＝３ｍ＋２番目の画像におけるＦｅのトレーサーの像の流れ及びＡｇのトレーサーの像の流れを、対応付けによって評価することで、Ｆｅのトレーサー及びＡｇのトレーサーの位置に対応する流動性媒体の流れを評価することができる。

ここで、ＣＣＤカメラ８による撮像は、１／６０秒間隔で行われており、例えば、Ｎ＝３ｍ＋１の画像の撮像時刻は、Ｎ＝３ｍの画像の撮像時刻と１／６０秒だけ異なる。図１２（ａ）に示すように、Ｎ＝３ｍの画像のみで、Ｐｂのトレーサーを識別することができる。しかしながら、１／６０秒後のＮ＝３ｍ＋１の画像においては、Ｗのトレーサーの吸収率及びＰｂのトレーサーの吸収率は、類似するため、Ｗのトレーサー及びＰｂのトレーサーの識別が、十分ではないことがある。そして、流動性媒体におけるＷ及びＰｂのトレーサーの流れが速く、Ｎ＝３ｍ＋１の画像の撮像時刻とＮ＝３ｍの画像の撮像時刻との差が小さい場合には、Ｎ＝３ｍの画像及びＮ＝３ｍ＋１の画像間におけるＷのトレーサー及びＰｂのトレーサーの変位が大きくなる。その結果、類似の画像の吸収率を備えたＷのトレーサーの対応付け及びＰｂのトレーサーの対応付けの誤差が大きくなる可能性もある。

そこで、互いに異なる複数の波長を備えた電磁波を発生させる複数の電磁波の発生源の各々に対応する、電磁波による画像を取得する手段を設けることが好ましい。この場合には、複数の電磁波の発生源の各々から放射される複数の電磁波による画像を、所定の時間間隔で、一つの電磁波による画像を取得する手段で取得する必要がなく、複数の電磁波による画像の各々を、複数の電磁波の各々に対応する、電磁波による画像を取得する手段によって、同時に取得することができる。その結果、複数のトレーサーの流れが速い場合であっても、複数のトレーサーの各々について、画像における複数のトレーサーの像の変位は、小さく、複数の電磁波の各々を適用して得られた画像ごとに、複数のトレーサーの識別及び対応付けの精度を向上させることができる。

このように、流動性媒体が、高速で移動している場合に、トレーサーの変位が大きいため、画像におけるトレーサーの濃度差を用いても、トレーサーの対応付けが、困難であることもある。図１３は、本発明による流動性媒体の流れ評価装置及び方法のさらに別の例を説明する図である。図１３に示す流動性媒体の流れ評価装置及び方法は、図１０に示す流動性媒体の流れ評価装置及び方法と同様であるが、複数のエックス線発生装置及び複数種類のトレーサーを採用すると共に、エックス線遮断機構を設けずに複数のエックス線発生装置の各々に対応する、エックス線発生装置から発生するエックス線を取得する複数の光学系を有する点で、図１０に示す流動性媒体の流れ評価装置及び方法と異なっている。

図１３においては、エックス線発生装置１Ａ、１Ｂ、１Ｃに対応する光学系１２Ａ、１２Ｂ、１２Ｃが設けられている。各々の光学系は、例えば、イメージインテンシファイヤー及びＣＣＤカメラ等を有してもよい。ここでは、エックス線発生装置１Ａ、１Ｂ、１Ｃの各々に対応して、各々の光学系１２Ａ、１２Ｂ、１２Ｃが設けられているので、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する特性エックス線を、エックス線遮断機構を使用して、遮断する又は通過させる必要はない。光学系１２Ａ、１２Ｂ、１２Ｃの各々は、それぞれ、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する三つの特性エックス線による流動性媒体３の画像を同時に取得する。その結果、エックス線発生装置１Ａ、１Ｂ、１Ｃから発生する三つの特性エックス線による流動性媒体３の画像の間で、流動性媒体３に含まれるトレーサー５の変位は、ほとんどなく、複数種類のトレーサーの識別及び対応付けの精度を向上させることができる。

本発明による流動性媒体の流れ評価方法等においては、少なくとも一種類のトレーサーは、複数種類のトレーサーであり、電磁波は、互いに異なる波長を備えた複数の電磁波を含み、互いに異なる波長を備えた複数の電磁波は、少なくとも一種類のトレーサーの一つの吸収端の波長よりも短い波長を備えた第一の電磁波及びその吸収端の波長よりも長い波長を備えた第二の電磁波を含み、第一の電磁波によって形成される複数種類のトレーサーを含む流動性媒体の第一の画像及び第二の電磁波によって形成される複数種類のトレーサーを含む流動性媒体の第二の画像に少なくとも減算処理を施して、複数種類のトレーサーを識別することを含む。

この場合には、二つの画像に減算処理を施すことによって、複数種類のトレーサーの識別を、より容易にすると共により正確にすることができる。特に、複数種類のトレーサーの像が、互いに重なり合っている場合にも、複数種類のトレーサーを識別することができる。さらに、複数種類のトレーサーの形状によらずに、複数種類のトレーサーの識別をすることができる。従って、複数種類のトレーサーの形状は、球形である必要はない。

図１４は、エネルギーサブトラクション法によるトレーサーの識別を説明する図である。ここでは、エネルギーサブトラクションは、互いに異なる複数の波長を備えた複数のエックス線によって、トレーサーの画像を撮像し、複数のエックス線によって得られた複数の画像に減算処理を施すことによって、トレーサーの像を画像上で強調する技術である。例えば、図１４に示すように、トレーサーを構成する材料が、材料Ｘであり、周囲のノイズ源となる物体（例えば、容器）の材料が、材料Ｙであるとする。材料Ｘ及び材料Ｙは、それぞれ、図１４に示すような、エックス線のエネルギーに対するエックス線の吸収率を有するものとする。ここで、材料Ｘのトレーサー及び材料Ｙの物体に対して、材料Ｘの吸収端よりも小さいエネルギーＡのエックス線及び材料Ｘの吸収端よりも大きいエネルギーＢのエックス線を別個に照射すると共に、材料Ｘのトレーサー及び材料Ｙの物体の画像を取得する。Ｉｍａｇｅ Ａは、材料Ｘの吸収端よりも小さいエネルギーＡのエックス線によって得られた材料Ｘのトレーサー及び材料Ｙの物体の画像であり、Ｉｍａｇｅ Ｂは、材料Ｘの吸収端よりも大きいエネルギーＢのエックス線によって得られた材料Ｘのトレーサー及び材料Ｙの物体の画像である。Ｉｍａｇｅ Ａにおいては、エネルギーＡのエックス線に対する材料Ｘのトレーサー及び材料Ｙの物体の吸収率は、ほぼ同等であるため、材料Ｘのトレーサー及び材料Ｙの物体の像が、ほぼ同等の濃度を有する。一方、Ｉｍａｇｅ Ｂにおいては、エネルギーＢのエックス線に対する材料Ｘのトレーサーの吸収率は、エネルギーＢのエックス線に対する材料Ｙの物体の吸収率よりも顕著に大きいため、材料Ｘのトレーサーの像の濃度も、材料Ｙの物体の像の濃度よりも顕著に高くなる。ここで、Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂに対して、デジタル画像処理演算によって減算処理（サブストラクション）を行うと、ＩｍａｇｅＢにおける濃度分布からＩｍａｇｅＡの濃度分布を差し引いたような画像ＩｍａｇｅＢ−ＩｍａｇｅＡが得られる。すなわち、ＩｍａｇｅＢ−ＩｍａｇｅＡにおいては、材料Ｙの物体の像の濃度は、非常に小さく、材料Ｘのトレーサーの像の濃度のみが高くなる。すなわち、ＩｍａｇｅＢ−ＩｍａｇｅＡでは、ＩｍａｇｅＡ及びＩｍａｇｅｂにおけるよりも材料Ｘのトレーサーの像が強調されている。このようにして、トレーサーの像を強調することで、トレーサーの対応付けの精度を向上させることができる。

本発明による流動性媒体の流れ評価方法等において、好ましくは、電磁波を、少なくとも一種類のトレーサーを含む流動性媒体に対して複数の異なる方向から照射し、複数の方向に対して、電磁波によって形成される少なくとも一種類のトレーサーを含む流動性媒体の画像を取得し、複数の方向に対して取得された画像の対応付けを行うことを含む。この場合には、複数の方向に対して取得された画像の対応付けを行うことによって、流動性媒体の三次元的な流れを評価することが可能となる。

ここで、複数の異なる方向から照射される電磁波を、複数の電磁波の発生源から発生させてもよい。この場合には、電磁波の発生源を移動させる必要が無く、複数の方向に対して取得される複数の画像を、比較的短い時間差で取得することができ、比較的速い流動性媒体の三次元的な流れを評価することができる。

また、複数の異なる方向から照射される電磁波を、移動可能な電磁波の発生源から発生させてもよい。 さらに、少なくとも一種類のトレーサーを含む前記流動性媒体の全体を移動させ、少なくとも一種類のトレーサーを含む流動性媒体の移動の前後に、電磁波によって形成される少なくとも一種類のトレーサーを含む流動性媒体の画像を取得し、移動の前後に取得された画像の対応付けを行うことを含む。これらの場合には、複数の電磁波の発生源を必要とせず、より小型の流動性媒体の流れ評価装置を提供することができる。

図１５は、本発明による流動性媒体の三次元的な流れの評価装置及び方法の例を説明する図である。図１５に示す流動性媒体の三次元的な流れの評価装置及び方法は、図１に示す流動性媒体の流れ評価装置及び方法と同様であるが、二つの同一のエックス線発生装置１を用いる点で、図１に示す流動性媒体の流れ評価装置及び方法と異なっている。また、流動性媒体の画像を取得する手段としては、エックス線による流動性媒体の画像を可視化するために、イメージングプレート１０を使用する。イメージングプレートの表面には、輝尽性蛍光体が塗布されており、この輝尽性蛍光体にエックス線が照射されると、エックス線のエネルギーが蓄積される。その後、イメージングプレートの輝尽性蛍光体にレーザー光等を照射すると、輝尽性蛍光体に蓄積されたエネルギーが、可視光として放射される。その結果、エックス線による流動性媒体の画像を可視化することができる。図１５においては、イメージングプレート１０は、（ｘ，ｙ，ｚ）座標で、ｚ＝０のｘｙ平面に対して平行に設置されている。また、二つのエックス線発生装置１におけるターゲットの位置は、それぞれ、（Ｘ，０，Ｌ）及び（−Ｘ，０，Ｌ）と設定する。ここで、ターゲットは、点状であり、ターゲットから放出される特性エックス線は、十分に平行光であるとする。このような図１５に示す流動性媒体の三次元的な流れの評価装置及び方法を用いて、流動性媒体３の画像を取得する。

図１６は、流動性媒体の三次元的な流れの評価装置及び方法を用いて得られる二つの画像を説明する図である。図１６に示すように、二つのエックス線発生装置１から放射されるエックス線によって、流動性媒体３に含まれる粒子のトレーサー５の二つの像Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂが、イメージングプレート１０に得られる。ここで、Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂにおけるトレーサー５の像は、前述の対応付けにより、同一のトレーサーである。すなわち、Ｉｍａｇｅ Ａは、二つのエックス線発生装置１の一方から放射されるエックス線によって得られたトレーサーの像であり、そのトレーサーの像の二次元座標は、Ｐ１（ｘ１，ｙ１）であるとする。また、Ｉｍａｇｅ Ｂは、二つのエックス線発生装置１の他方から放射されるエックス線によって得られたトレーサーの像であり、そのトレーサーの像の二次元座標は、Ｐ２（ｘ２，ｙ２）であるとする。これらの二つの画像における同一のトレーサー５の座標Ｐ１（ｘ１，ｙ１）及びＰ２（ｘ２，ｙ２）並びに三角測量の原理を用いて、このトレーサー５の三次元的な真の位置（ｘ，ｙ，ｚ）を、以下のような計算式
ｘ＝Ｘ（ｘ１＋ｘ２）／（２Ｘ−（ｘ１−ｘ２））
ｙ＝２Ｘ×ｙ１／（２Ｘ−（ｘ１−ｘ２））
ｚ＝（ｘ１−ｘ２）Ｌ／（（ｘ１−ｘ２）−２Ｘ）
から得ることができる。このようにして、流動性媒体３におけるトレーサー５の三次元的な位置を得ることができる。そして、異なる時間において流動性媒体３の複数の画像を取得し、トレーサー５の三次元的な位置を用いて、トレーサー５の三次元的な流れを評価し、トレーサー５の位置に対応する流動性媒体３の位置における三次元的な流れを評価することができる。

また、二つの同一のエックス線発生装置１を用いる代わりに、単一の移動可能なエックス線装置を用いてもよい。まず、ある時刻に単一の移動可能なエックス線発生装置を、座標（Ｘ，０，Ｌ）に設置して、エックス線発生装置から放出されるエックス線による流動性媒体３の第一の画像を取得する。次に、その単一の移動可能なエックス線発生装置を、流動性媒体３の流れの速度に対して十分速い速度で、座標（−Ｘ，０，Ｌ）に移動させ、エックス線発生装置から放出されるエックス線による流動性媒体３の第二の画像を取得する。このようにして、上述したように、流動性媒体３におけるトレーサー５の三次元的な位置を得ることができる。

図１７は、本発明による流動性媒体の三次元的な流れの評価装置及び方法の別の例を説明する図である。図１７に示す流動性媒体の三次元的な流れの評価装置及び方法は、図１５に示す流動性媒体の三次元的な流れの評価装置及び方法と同様であるが、単一のエックス線発生装置を用いると共にトレーサー５を含む流動性媒体３の全体を移動させる点で、図１５に示す流動性媒体の三次元的な流れの評価装置及び方法と異なっている。より具体的には、トレーサー５を含む流動性媒体３の全体を含有する容器の全体を移動させる。ここで、トレーサー５を含む流動性媒体３の全体を移動させる速度は、流動性媒体３の流れの速度に対して十分速い速度で移動させる。また、トレーサー５を含む流動性媒体３の全体を移動させる際の加速度が、流動性媒体の流れに影響を与えないことが必要である。

図１７においては、単一のエックス線発生装置１は、（ｘ，ｙ，ｚ）座標において、（０，０，Ｌ）の位置に設置されている。そして、単一のエックス線発生装置１から放出されるエックス線を、トレーサー５を含む流動性媒体３の全体を含有する容器の移動の前後に、流動性媒体３に照射する。ここで、トレーサー５を含む流動性媒体３又はトレーサー５を含む流動性媒体３の全体を含有する容器における代表点を設定する。ただし、容器の移動の前後におけるその代表点のｘ座標が、Ｘ及び−Ｘであるように、代表点は、決定される。そして、トレーサー５を含む流動性媒体３の全体を含有する容器の移動の前後における流動性媒体３の二つの画像を、イメージングプレート１０を用いて取得する。

イメージングプレート１０上に得られた、容器の移動の前後における流動性媒体３の二つの画像が、図１６に示すような二つの画像Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂであるとする。なお、二つの画像Ｉｍａｇｅ Ａ及びＩｍａｇｅ Ｂにおけるトレーサー５の像は、対応付けによって、同一のトレーサーであるとする。Ｉｍａｇｅ Ａにおけるトレーサー５の像の二次元座標Ｐ１（ｘ１，ｙ１）及びＩｍａｇｅ Ｂにおけるトレーサー５の像の二次元座標Ｐ２（ｘ２，ｙ２）並びに三角測量の原理を用いて、このトレーサー５の三次元的な真の位置（ｘ，ｙ，ｚ）は、以下のような計算式
ｘ＝ｘ１×Ｘ／（ｘ２−ｘ１）
ｙ＝ｙ１×Ｘ／（ｘ２−ｘ１）
ｚ＝−Ｘ×Ｌ／（ｘ２−ｘ１）
から得ることができる。このようにして、流動性媒体３におけるトレーサー５の三次元的な位置を得ることができる。そして、異なる時間において流動性媒体３の複数の画像を取得し、トレーサー５の三次元的な位置を用いて、トレーサー５の三次元的な流れを評価し、トレーサー５の位置に対応する流動性媒体３の位置における三次元的な流れを評価することができる。

なお、本発明の発明者は、粉体（を含む流動性媒体）における、粉体の密度と大きく異なる密度を備えたトレーサーの流れが、粉体（を含む流動性媒体）の流れと実質的に一致することを発見し、本発明の第一、第二、及び第三の態様である流動性媒体の流れ評価方法、装置、及びプログラムを発明するに至った。

例えば、有色のトナーのような樹脂からなる有色の粉体の流れを計測するために、有色の粉体の内部に到達することが可能なＸ線及びＸ線を吸収することが可能な金属のトレーサー粒子を用いることが考えられるが、樹脂からなる粉体の密度と金属のトレーサー粒子の密度は、大きく（数十倍）異なる。このため、一般的には、樹脂からなる粉体の流れと金属のトレーサー粒子の流れは、大きく異なると予想される。

また、流体の流れの可視化については、流体の流れに対するトレーサーの追随性のガイドラインに関する研究（例えば（社）可視化情報学会編：ＰＩＶハンドブック（２００２），ｐｐ．３０−３５．）は、行われてきたが、粉体の流れに対するトレーサーの追随性に関する研究は、粉体におけるトレーサーの流れが、粉体がトレーサーに及ぼす機械的力場によって、粉体の流れと大きく異なると予想されるため、現在までのところ行われていない。

驚くべきことに、発明の発明者は、樹脂からなる有色の粉体であるトナーの流れを、Ｘ線及び金属のトレーサー粒子を用いて計測し、樹脂からなる粉体であるトナーの流れと、金属のトレーサー粒子の流れとが、実質的に一致することを確認した。すなわち、本発明の発明者は、粉体を含む流動性媒体の流れを、流動性媒体と異なる少なくとも一種類のトレーサー及び粉体を透過すると共に少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価することが可能であることを新たに見出した。

具体的には、本発明による流動性媒体の流れ評価方法等が適用され得る、電子写真装置内の現像器における（トナーの粉体である）現像剤の搬送を、複数種類の金属トレーサーを用いて評価した。ここで、複数種類の金属トレーサーは、それぞれ、タングステンからなる粒径８μｍのトレーサー及び鉄からなる粒径５０μｍのトレーサーである。

図１８（ａ）は、現像剤の搬送の評価に使用したスクリュー搬送式の現像器の概略図である。現像剤の搬送の評価に使用したスクリュー搬送式の現像器による現像剤のスクリュー搬送速度は、４．８４ｍｍ／秒とした。また、現像剤は、平均粒子径８μｍ及び平均密度１２００ｋｇ／ｍ^３のトナーの粉体であった。

図１８（ｂ）は、スクリュー搬送式の現像器によって搬送される、現像剤の搬送方向におけるタングステン及び鉄のトレーサーの移動度に関する測定結果を示す図である。図１８（ｂ）に示すグラフの横軸は、時間（秒）であり、縦軸は、トレーサーの移動距離（ｍｍ）である。本発明による流動性媒体の流れ評価方法等に従って、タングステン及び鉄のトレーサーの移動度を測定した。なお、現像剤（トナー）の移動度は、現像剤のトナーと異なる色のトナー（異色トナー）を現像剤に混入させ、現像剤の表面に露出した異色トナーの位置を測定して、現像剤の移動度を推定した。

図１８（ｂ）に示すように、タングステン及び鉄のトレーサーの移動度は、それぞれ、現像剤（トナー）の移動度とほぼ同一の傾向を示している。よって、少なくともスクリュー搬送方式の現像器における粉体の現像剤の搬送に関しては、トレーサーの流れ（搬送）は、粉体の現像剤の流れ（搬送）と実質的に一致することが確認された。よって、（一定の振動を与えるような場のような）特殊な機械的作用場を除き、トレーサーの流れは、粉体を含む流動性媒体の流れと概略一致し、本発明による流動性媒体の流れ評価方法等は、粉体を含む流動性媒体の流れに広く適用することができると言える。

次に、本発明による流動性媒体の流れ評価プログラムのアルゴリズムの例を説明する。この流動性媒体の流れ評価プログラムのアルゴリズムは、事前計測部、画像取得部、及び画像解析部を含み得るが、事前計測部を除く画像取得部及び画像解析部からなるアルゴリズムであってもよい。

＜事前計測部＞
（ｉ）トレーサーとして使用するｉ種類のトレーサー（ｉは、トレーサーの種類の数であり、自然数である）を用意し、それぞれを質量順にＴ（１）〜Ｔ（ｉ）とインデックスする。例えば、直径１００μｍのＦｅ、Ｗ及びＰｂの３種類の粒子（ｉ＝３）を用意し、それぞれの粒子を、Ｔ（１）、Ｔ（２）、Ｔ（３）とする。

（ｉｉ）流動性媒体の流れの評価に使用する電磁波に対して、画像におけるＴ（１）〜Ｔ（ｉ）のトレーサーの濃度を検出しておき、それらのトレーサーの濃度を、それぞれ、Ｄ（１）〜Ｄ（ｉ）とする。トレーサーの濃度は相対値である。また、複数の波長を備えた複数の電磁波を用いる場合には、複数の電磁波の各々に対して、トレーサーの濃度を検出する。例えば、画像におけるＴ（１）（直径１００μｍのＦｅ粒子）の濃度を計測し、その濃度の値をＤ（１）とする。

＜画像取得部＞
（ｉ）粉体を含む流動性媒体の流れに、Ｔ（１）〜Ｔ（ｉ）のトレーサーを混入させ、電磁波を照射させる。また、複数の波長を備えた複数の電磁波を用いる場合には、複数の電磁波の各々を照射させる。さらに、粉体を含む流動性媒体の流れを三次元的に評価する場合には、同一の電磁波を、流動性媒体に対して異なる方向から照射させる。

（ｉｉ）画像を取得する適当な手段に、Ｔ（１）〜Ｔ（ｉ）のトレーサーを含む流動性媒体に電磁波を照射して得られる静止画像を、適当な時間間隔ｄｔ毎に取得させる。ここで、適当な時間間隔ｄｔは、流動性媒体の流れの性質に依存して決定される。時間間隔ｄｔが長いときには、画像解析部における全体の計算量は少なくなり、計算の高速化が達成され望る。しかしながら、流動性媒体の流れの変動が激しい場合、又は流動性媒体の流れが、複雑な流れ場を示す場合には、時間間隔ｄｔを短くする必要がある。

（ｉｉｉ）時間間隔ｄｔ毎に取得された静止画像を、時系列にＩ（ｎ）とインデックスする。ｎは、最初の静止画像に対する静止画像の番号である。例えばＩ（１００）は、１００番目に取得された静止画像である。その静止画像が取得された時刻は、最初の静止画像が取得された時刻が０であるとすると、１００×ｄｔ経過後の時刻となる。

＜画像解析部＞
（ｉ）流動性媒体の流れの情報を必要とする時間領域ｔを、ｔ０≦ｔ≦ｔ１とする。

（ｉｉ）以下の１）〜４）の操作を、ｔ０≦ｔ≦ｔ１の間に時間間隔ｄｔ毎に得られた静止画像に対して、繰り返す。

１）静止画像Ｉ（ｔ）及びＩ（ｔ＋ｄｔ）に着目する。

２）静止画像Ｉ（ｔ）及びＩ（ｔ＋ｄｔ）内におけるトレーサーの像に対してトレーサー粒子認識を行う。トレーサー粒子認識は、画像内にトレーサーとして認識される画像領域を定義することを意味する。トレーサー粒子認識には、公知の粒子認識アルゴリズムを利用することができる。例えば、８近傍連結成分のラベリング処理が挙げられる。また、トレーサーの形状及び大きさを、予め規定することができれば、粒子認識の精度を向上させることができる。トレーサー粒子認識によって、静止画像Ｉ（ｔ）にｘ個のトレーサーが認識され、静止画像Ｉ（ｔ＋ｄｔ）にｙ個のトレーサーが認識されたとする。静止画像Ｉ（ｔ）にｘ個のトレーサーを、Ｐ（１）〜Ｐ（ｘ）とラベリングし、静止画像Ｉ（ｔ＋ｄｔ）におけるｙ個のトレーサーを、Ｑ（１）〜Ｑ（ｙ）とラベリングする。

３）トレーサーＰ（１）〜Ｐ（ｘ）及びＱ（１）〜Ｑ（ｙ）の粒子画像平均濃度を算出し、粒子画像平均濃度のそれぞれをＤＰ（１）〜ＤＰ（ｘ）、ＤＱ（１）〜ＤＱ（ｙ）とする。ここで、粒子画像平均濃度は、例えば、粒子認識で定義された粒子空間内に含まれる画像ピクセルの濃度の算術平均値である。

４）以下の操作ａ）〜ｄ）を、１≦ｓ≦ｘとなるｓに対して行う。画像Ｉ（ｔ）内の粒子Ｐ（ｓ）に着目し、粒子Ｐ（ｓ）に対応する画像Ｉ（ｔ＋ｄｔ）内に対応する粒子を以下の操作ａ）〜ｄ）で見つける。ここで、粒子Ｐが粒子Ｑに対応するとは、画像Ｉ（ｔ）内で認識された粒子Ｐが、流動性媒体の流れに従って移動し、画像（Ｉ＋ｄｔ）において粒子Ｑとして撮影されることを意味する。

ａ）粒子Ｐ（ｓ）のトレーサー種を判定する。

ＤＰ（ｓ）−ＤＰ（ｉ）の絶対値が最小となるｉを備えたＴ（ｉ）を、Ｐ（ｓ）のトレーサーの種類とする。このｉをａとする。

ｂ）画像Ｉ（ｔ＋１）内に粒子Ｐ（ｓ）に対応する粒子候補を選択する。そのひとつをＱ（ｕ）とする。粒子候補の選択は、例えば、画像Ｉ（ｔ＋ｄｔ）内の座標（Ｐ（ｓ）ｘ、Ｐ（ｓ）ｙ）を中心として、一定の距離しきい値である距離ｄｓ内に含まれる全粒子を候補とする。しきい値ｄｓは、流れの性質によって定める。流動性媒体における解析空間内の流れの予想最大速度が、ｖである場合には、ｄｓは、ｖ×ｄｔ程度とする。次に、粒子Ｑ（ｕ）のトレーサーの種類を判定する。ＤＱ（ｕ）−ＤＱ（ｉ）の絶対値が最小となるｉを備えたＴ（ｉ）を、Ｑ（ｕ）のトレーサーの種類とする。このｉをｂとする。

ｃ）Ｑ（ｕ）のＰ（ｓ）に対する粒子対応指数ＥＲを下記のように計算する。

ＥＲ＝ｋ×位置対応指数＋（１−ｋ）トレーサー種対応指数（ｋは、０以上１以下の重み付け係数とする。）
ここで、位置対応指数には、公知のＰＴＶアルゴリズムを利用することができ、例えば、図４と共に説明した近傍４粒子距離判定アルゴリズムが挙げられる。周囲の粒子の状況が異なるほど、位置対応指数は大きくなる（候補の判定の可能性が下がる）。また、トレーサー種対応指数は、例えばａ−ｂの絶対値と定義する。すなわちトレーサーの種類が一致すれば、トレーサー種対応指数は、０である。トレーサーの種類が異なると判定された候補粒子に対しては、そのトレーサーの質量の差が大きいほど、対応指数は、大きくなる（候補の判定の可能性が下がる）。

そして、最も小さいＥＲを備えたＱ（ｕ）をＰ（ｓ）の対応粒子と定義する。

ｄ）Ｐ（ｓ）及びＱ（ｕ）の重心座標を結ぶ位置ベクトルをもって、流動性媒体の座標（Ｐ（ｓ）ｘ、Ｐ（ｓ）ｙ）における流れのベクトルとする。ここで、Ｐ（ｓ）ｘ及びＰ（ｓ）ｙは、それぞれ、粒子Ｐ（ｓ）の重心のｘ座標及びｙ座標である。

（ｉｉｉ）適当な過誤ベクトル消去処理を行う。過誤ベクトル消去処理には、公知のアルゴリズムを利用することができる。例えば、流れの連続性を仮定した消去アルゴリズムである。

以上のようなアルゴリズムによって、時間ｔ０≦ｔ≦ｔ１における流動性媒体の流れ場のベクトル情報を得ることができる。

なお、複数の波長を備えた複数の電磁波を用いる場合には、複数の電磁波の各々に対して画像解析部の操作を行う。また、流動性媒体の三次元的な流れを評価する場合には、二つの異なる方向から照射された電磁波による二つの画像の間で、トレーサー粒子の対応付けをさらに行う。

なお、本発明による流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムは、粉体及び／又は流体を含む流動性媒体の全体に適用され得る。すなわち、本発明による流動性媒体の流れ評価方法、流動性媒体の流れ評価装置、及び流動性媒体の流れ評価プログラムは、粉体のみからなる流動性媒体並びに粉体及び流体の両方を含む流動性媒体のような、粉体を含む流動性媒体のみならず、粉体を含まない流動性媒体、すなわち、液体及び気体並びに液体及び気体の混合物のような、流体のみからなる流動性媒体などにも適用され得る。

以上、本発明の実施の形態及び実施例を具体的に説明してきたが、本発明は、これらの実施の形態及び実施例に限定されるものではなく、これら本発明の実施の形態及び実施例を、本発明の主旨及び範囲を逸脱することなく、変更又は変形することができる。

本発明による流動性媒体の流れ評価装置及び方法の例を説明する図である。 複数の画像の間における粒子のトレーサーの位置の対応付けを説明する図である。 トレーサーを含む流動性媒体の画像に含まれる複数の粒子のトレーサーの像に対する粒子分離処理を説明する図である。 ＰＴＶによる複数の画像の間における粒子のトレーサーの対応付けに関するアルゴリズムの例を説明する図である。 本発明による流動性媒体の流れ評価装置及び方法の別の例を説明する図である。 ＰＩＶによる複数の画像の間における液体のトレーサーの対応付けに関するアルゴリズムの例を説明する図である。 複数種類のトレーサーの材料としての金属によって吸収されるエックス線の質量吸収係数を説明する図である。 （ａ）は、複数種類の直径１００μｍの球形粒子である複数種類のトレーサーによって吸収されるエックス線の相対吸収率を説明する図であり、（ｂ）は、複数種類のトレーサーとしての複数種類の直径１００μｍの球形粒子の中心を通過するエックス線の相対吸収率を説明する図である。 （ａ）及び（ｂ）は、互いに異なるエネルギーを備えたエックス線を照射することによって得られた二種類の金属トレーサーを含む流動性媒体の画像における対応付けを説明する図である。 本発明による流動性媒体の流れ評価装置及び方法のまた別の例を説明する図である。 図１０に示す流動性媒体の流れ評価方法等における特性エックス線の照射及び特性エックス線による撮像の時間スケジュールを示す図である。 （ａ）、（ｂ）及び（ｃ）は、図１０に示す流動性媒体の流れ評価方法を用いて、図１１に示す時間スケジュールに従って撮像された粒子のトレーサーの画像を示す図である。 本発明による流動性媒体の流れ評価装置及び方法のさらに別の例を説明する図である。 エネルギーサブトラクション法によるトレーサーの識別を説明する図である。 本発明による流動性媒体の三次元的な流れの評価装置及び方法の例を説明する図である。 流動性媒体の三次元的な流れの評価装置及び方法を用いて得られる二つの画像を説明する図である。 本発明による流動性媒体の三次元的な流れの評価装置及び方法の別の例を説明する図である。 （ａ）は、現像剤の搬送の評価に使用したスクリュー搬送式の現像器の概略図であり、（ｂ）は、スクリュー搬送式の現像器によって搬送される、現像剤の搬送方向におけるタングステン及び鉄のトレーサーの移動度に関する測定結果を示す図である。

符号の説明

１ エックス線発生装置
２ エックス線
３ 被写体
４ 粉流体媒体
５ 粒子状トレーサー
６ イメージインテンシファイヤー
７ 回転鏡
８ ＣＣＤカメラ
９ 演算装置
１０ イメージングプレート
１１ 遮蔽版
１２ 光学系
１３ 液体状トレーサー

Claims (20)

1. 粉体を含む流動性媒体の流れを、該流動性媒体と異なる少なくとも一種類のトレーサー及び該流動性媒体を透過すると共に該少なくとも一種類のトレーサーに吸収される電磁波を用いて、評価することを特徴とする流動性媒体の流れ評価方法。
2. 前記少なくとも一種類のトレーサーを含む前記流動性媒体に前記電磁波を照射し、
   該電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、
   該画像を用いて前記粉体を含む流動性媒体の流れを評価することを含むことを特徴とする請求項１に記載の流動性媒体の流れ評価方法。
3. 前記流動性媒体は、粉体であることを特徴とする請求項１又は２に記載の流動性媒体の流れ評価方法。
4. 前記粉体は、トナーを含む現像剤であることを特徴とする請求項１乃至３のいずれか一項に記載の流動性媒体の流れ評価方法。
5. 前記トレーサーは、粒子であり、
   前記流動性媒体の流れと共に移動する該粒子の位置の対応付けを行うことを含むことを特徴とする請求１乃至４のいずれか一項に記載の流動性媒体の流れ評価方法。
6. 前記トレーサーは、液体又は粒子集合体であり、
   前記流動性媒体の流れと共に移動する該液体又は粒子集合体の対応付けを行うことを含むことを特徴とする請求項１乃至４のいずれか一項に記載の流動性媒体の流れ評価方法。
7. 前記少なくとも一種類のトレーサーは、複数種類のトレーサーであることを特徴とする請求項１乃至６のいずれか一項に記載の流動性媒体の流れ評価方法。
8. 前記電磁波は、互いに異なる波長を備えた複数の電磁波を含むことを特徴とする請求項１乃至７のいずれか一項に記載の流動性媒体の流れ評価方法。
9. 前記少なくとも一種類のトレーサーは、複数種類のトレーサーであり、
   前記電磁波は、互いに異なる波長を備えた複数の電磁波を含み、
   該互いに異なる波長を備えた複数の電磁波は、前記少なくとも一種類のトレーサーの一つの吸収端の波長よりも短い波長を備えた第一の電磁波及び該吸収端の波長よりも長い波長を備えた第二の電磁波を含み、
   前記第一の電磁波によって形成される前記複数種類のトレーサーを含む前記流動性媒体の第一の画像及び前記第二の電磁波によって形成される前記複数種類のトレーサーを含む前記流動性媒体の第二の画像に少なくとも減算処理を施して、前記複数種類のトレーサーを識別することを含むことを特徴とする請求項１乃至６のいずれか一項に記載の流動性媒体の流れ評価方法。
10. 前記互いに異なる波長を備えた複数の電磁波は、それぞれ、異なる時間に前記少なくとも一種類のトレーサーを含む前記流動性媒体に照射されることを特徴とする請求項８又は９に記載の流動性媒体の流れ評価方法。
11. 前記互いに異なる波長を備えた複数の電磁波は、単一の電磁波の発生源から選択的に発生させられることを特徴とする請求項８乃至１０のいずれか一項に記載の流動性媒体の流れ評価方法。
12. 前記互いに異なる波長を備えた複数の電磁波は、複数の電磁波の発生源から発生させられることを特徴とする請求項８乃至１０のいずれか一項に記載の流動性媒体の流れ評価方法。
13. 前記電磁波は、エックス線であることを特徴とする請求項１乃至１２のいずれか一項に記載の流動性媒体の流れ評価方法。
14. 前記エックス線の波長は、前記少なくとも一種類のトレーサーの少なくとも一つの吸収端の波長を含むことを特徴とする請求項１３に記載の流動性媒体の流れ評価方法。
15. 前記電磁波を、前記少なくとも一種類のトレーサーを含む前記流動性媒体に対して複数の異なる方向から照射し、
    該複数の方向に対して、該電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、
    前記複数の方向に対して取得された該画像の対応付けを行う
    ことを含むことを特徴とする請求項１乃至１４のいずれか一項に記載の流動性媒体の流れ評価方法。
16. 前記複数の異なる方向から照射される前記電磁波は、複数の電磁波の発生源から発生させられることを特徴とする請求項１５に記載の流動性媒体の流れ評価方法。
17. 前記複数の異なる方向から照射される前記電磁波は、移動可能な電磁波の発生源から発生させられる請求項１５に記載の流動性媒体の流れ評価方法。
18. 前記少なくとも一種類のトレーサーを含む前記流動性媒体の全体を移動させ、
    前記少なくとも一種類のトレーサーを含む前記流動性媒体の移動の前後に、前記電磁波によって形成される前記少なくとも一種類のトレーサーを含む前記流動性媒体の画像を取得し、
    該移動の前後に取得された該画像の対応付けを行う
    ことを含むことを特徴とする請求項１乃至１４のいずれか一項に記載の流動性媒体の流れ評価方法。
19. 請求項１乃至１８のいずれか一項に記載の流動性媒体の流れ評価方法を実行する手段を含むことを特徴とする流動性媒体の流れ評価装置。
20. 請求項１乃至１８のいずれか一項に記載の流動性媒体の流れ評価方法をコンピュータに実行させることを特徴とする流動性媒体の流れ評価プログラム。

Images