JP2014202722A - 混相状態分布計測装置 - Google Patents

Abstract

【課題】混相媒体分布の再構成画像の精度を向上させる装置を提供する。【解決手段】管路１０の周辺に配置されたｎ個（ただし、ｎは３以上）のセンシング電極２０を有するセンサ１２と、センサ１２を用いて静電容量又は電気抵抗を計測する計測手段１４と、を備え、センシング電極２０のうち、ｐ個（ただし、１≰ｐ≰ｎ−１）の印加電極に測定用電圧を印加し、ｑ個（ただし、１≰ｑ≰ｎ−ｐ，３≰ｐ＋ｑ）の接地電極を接地した状態で計測手段１４により静電容量又は電気抵抗を計測する混相状態分布計測装置とする。【選択図】図１

Description

本発明は、複数の媒体の混相状態を測定する混相状態分布計測装置に関する。

固体、液体、気体等の混相状態を測定する混相状態分布計測装置は、パイプラインや他の機器における管路の閉塞状態の確認、輸送の効率化、混合の効率化、抵抗低減、腐食防止などのために用いられている。

混相状態分布計測装置では、ｒ個の電極を含むセンサを管路周辺に配置し、ｒ個の電極のうち接地電極と測定電圧を印加する印加電極とを１つずつ設け、接地電極と印加電極との組み合わせを変更しつつ静電容量を測定し、これらを利用して管内の状態を再現する。このときの静電容量の測定可能数は、ｒ個の電極の順列組み合わせの半分であるｒ（ｒ−１）／２である。例えば、８個の電極を用いた場合はｎ＝２８通りの静電容量の測定が可能となる。

混相状態分布を求める方法として、既知のｎ×１行列の静電容量行列Ｃ、既知のｎ×ｍ行列のセンシティビティマップＳｅ、未知のｍ×１行列の静電率分布行列Ｅとが、Ｃ＝ＳｅＥの関係にあり、この式の反復計算により誘電率分布行列Ｅを求める技術が開示されている。このとき、γｉ＝Ｃ^ＴＳｅｉとおき、反復回数をｋとしたときに評価関数をｆ（Ｅ^（ｋ））＝γ・Ｅ^（ｋ）とし、この評価関数が１に近づく誘電率分布行列Ｅを求めることが記載されている（特許文献１）。

特開２００３−１３０８３５号公報

ところで、ｎ個の静電容量の測定結果から解像度ｍの画像を構成する場合、通常はｎ＜ｍであるために画像の再構成時における画像推定問題の解が不定となり、画像が一意に定まらない。ｎがｍに近づくほど解は求め易くなり、ｎ＝ｍであれば画像は一意に定まる。

そこで、静電容量の測定結果の数を解像度に近づけて、画像を安定に再現できる混相状態分布計測装置が望まれている。

本発明の１つの態様は、管路の周辺に配置されたｎ個（ただし、ｎは３以上）のセンシング電極を有するセンサと、前記センサを用いて静電容量又は電気抵抗を計測する計測手段と、を備え、前記センシング電極のうち、ｐ個（ただし、１≦ｐ≦ｎ−１）の電極に測定用電圧を印加し、ｑ個（ただし、１≦ｑ≦ｎ−ｐ，３≦ｐ＋ｑ）の電極を接地した状態で前記計測手段により静電容量又は電気抵抗を計測することを特徴とする混相状態分布計測装置である。

ここで、ｐ又はｑは、少なくとも２以上に変化させることが好適である。

また、前記管路内の媒体の分布を示す画像を生成する画像再構成手段と、前記画像を表示する画像表示手段と、を備えることが好適である。

また、前記計測手段で得られた静電容量又は電気抵抗と一致する前記管路内部の媒体の誘電率分布を反復計算により算出する際に、誘電率分布の変化が静電容量又は電気抵抗に与える影響を誘電率分布の更新毎に算出し、算出して得られた感度を用いて反復計算を実行することが好適である。

また、ｎ＝８，ｐ＝１，ｑ＝７の状態を含むことが好適である。

本発明によれば、混相媒体分布の再構成画像の精度を向上させることができる。

本発明の実施の形態における混相状態分布計測装置の構成を示す図である。 本発明の実施の形態におけるセンサの構成を示す断面図である。 本発明の実施の形態における演算処理フローチャートである。 従来の演算処理フローチャートである。 本発明の実施の形態におけるた従来の演算処理を適用した混相状態分布計測装置で得られた誘電率分布の例を示す図である。 本発明の実施の形態における混相状態分布計測装置で得られた誘電率分布の例を示す図である。

本発明の実施の形態における混相状態分布計測装置１００は、図１に示すように、センサ１２、静電容量計測手段１４、誘電率分布再構成手段１６及び再構成画像出力手段１８を含んで構成される。混相状態分布計測装置１００は、センサ１２が測定対象となる管路１０の外周部に巻かれるように配置される。

センサ１２は、図２の断面図に示すように、複数のセンシング電極２０を組み合わせて環状に構成される。本実施の形態では、センサ１２は、８つのセンシング電極２０−１〜２０−８により構成される例を示している。ただし、センサ１２に含まれるセンシング電極２０の数は、これに限定されるものでなく、ｒ個（ｒは、３以上）であればよい。各センシング電極２０は、静電容量計測手段１４から測定電圧が印加される、又は、接地されて、その組み合わせによって異なるコンデンサを形成する。

静電容量計測手段１４は、リード線２２を介してセンシング電極２０により形成されるコンデンサの静電容量を計測する。本実施の形態では、センサ１２は８つのセンシング電極２０−１〜２０−８を含むので、静電容量計測手段１４も８つのリード線２２−１〜２２−８を含む例を示している。ただし、リード線２２の数は、これに限定されるものでなく、ｒ個（ｒは、３以上）のセンシング電極２０を備える場合には、センシング電極２０の数に合わせてｒ本設ければよい。

静電容量計測手段１４は、ｒ個のセンシング電極２０のうち、ｐ個（ただし、１≦ｐ≦ｒ−１）を印加電極として選択して測定用電圧を印加し、測定用電圧を印加しなかったセンシング電極２０のうちｑ個（ただし、１≦ｑ≦ｒ−ｐ，３≦ｐ＋ｑ）を接地電極として選択して接地する。

例えば、ｒ＝８の場合、ｐは１以上７以下の数から選択される。ここで、ｐを１に設定した場合、ｑは２以上７以下の数から選択される。表１は、これらの条件を満たすセンシング電極２０の組み合わせの例を示す。電極数をｒ、印加電極の数をｐ、接地電極の数をｑとしたとき、電極の組み合わせの数は_ｒＣ_ｐ×_ｒ−ｐＣ_ｑとなる。ただし、印加電極の数ｐ＝１、接地電極の数ｑ＝１のときは、印加電極と接地電極とを交換しても情報量は増えないので計測できる静電容量の数はｎ＝_ｒＣ_１×_ｒ−１Ｃ_１／２＝ｒ（ｒ−１）／２となる。一方、印加電極の数ｐ＝１としても、接地電極の数ｑ＝２であれば、計測できる静電容量の数はｎ＝_ｒＣ_１×_ｒ−１Ｃ_２＝ｒ（ｒ−１）（ｒ−２）となり、接地電極の数ｑ＝１のときより大幅に増加する。表１では、センシング電極２０の数ｒ＝８のときにおいて、センシング電極２０の組み合わせの数及び測定可能な静電容量の数を併せて示している。

このように、本実施の形態では、センシング電極２０のうち印加電極及び接地電極となる電極をそれぞれ１つに限定しないことにより、より多くの静電容量の測定結果に基づいて管路１０内の混相状態の再現を行うことができる。

誘電率分布再構成手段１６は、_ｒＣ_ｐ×_ｒ−ｐＣ_ｑ個の電極の組み合わせからなるコンデンサの静電容量に基づいて管路１０内の計測空間１０ａにおける誘電率分布を演算し、計測空間１０ａの誘電率分布ＣＴの再構成画像信号を生成する。誘電率分布再構成手段１６は、演算装置１６ａを含んで構成することができる。演算装置１６ａは、例えば、プロセッサ、メモリ、入出力インタフェースを備えるコンピュータで構成される。また、静電容量計測手段１４、誘電率分布再構成手段１６、及び演算装置１６ａを全てコンピュータで構成してもよく、再構成画像出力手段１８を当該コンピュータのディスプレイで構成してもよい。演算装置１６ａは、画像再構成プログラムを実行することによって静電容量正規化手段１６ｂ、誘電率明確化手段１６ｃ及び誘電率空間勾配最小化手段１６ｄを実現する。

画像再構成方法では、数値シミュレーションを用いて計測した静電容量と一致する誘電率分布を求める。例えば、数値シミュレーションに有限要素法に基づくものを用いた場合、管路１０内の計測空間１０ａをｍ個の有限要素（メッシュ）に分割し、各メッシュにおける誘電率分布を算出する。メッシュの個数は、再構成画像の解像度に対応する。ここで、計測された静電容量の数をｒ個、求める誘電率の個数（メッシュの個数と等しい）をｍ個とすると、静電容量行列Ｃ（ｒ×１行列）、誘電率分布行列Ｅ（ｍ×１行列）、センシティビティマップ行列Ｓｅ（ｒ×ｍ行列）の関係は、Ｃ＝ＳｅＥで表わされる。

誘電率分布行列Ｅを求めるには、センシティビティマップ行列Ｓｅの逆行列を求めればよいが、センシティビティマップ行列Ｓｅは正方行列でないので逆行列が存在しない。したがって、ＬＢＰ（Ｌｉｎｅｒ Ｂａｃｋ Ｐｒｏｊｅｃｔｉｏｎ）法やＮＲ（Ｎｅｗｔｏｎ−Ｒａｐｈｓｏｎ）法により誘電率分布行列Ｅを求めることができる。

本実施の形態では、例として、以下の方法を示す。誘電率分布再構成手段１６は、静電容量計測手段１４で計測されたｎ個の静電容量Ｃｊ^０に基づいて誘電率分布を再構成する。

図３に、本実施の形態の演算処理をフローチャートで示す。まず、初期値、すなわち管路１０内の初期誘電率分布を設定する（Ｓ１０１）。

次に、初期誘電率分布から順解析を行って静電容量を算出する（Ｓ１０２）。ここで、「順解析」とは、誘電率分布から支配方程式（ポアソン方程式）を用いて電位（静電ポテンシャル）を求めることをいい、電位（静電ポテンシャル）から誘電率分布を求める「計測」に対する概念である。

次に、順解析により得られた静電容量を用いて目的関数を計算する（Ｓ１０３）。目的関数は、具体的には、順解析により得られた静電容量と、計測して得られた静電容量との最小二乗誤差である。但し、本実施の形態では、目的関数を定義する際に、正規化、明確化及び平滑化を行うことで画像の精度を向上させる。

次に、目的関数の感度解析を行う（Ｓ１０４）。感度解析とは、反復計算の終了判定、及び設定された誘電率の更新量を決定するために目的関数の感度を求めることをいう。設定された誘電率は、Ｓ１０１で設定された初期誘電率及び反復計算により順次更新されたその都度の誘電率をいう。

次に、目的関数が収束したか否かを判定する（Ｓ１０５）。この判定は、具体的にはＳ１０４で算出された目的関数の感度が最小値、つまり０ないしその近傍であるか否かにより行われる。目的関数の感度が０ないしその近傍である場合には収束したと判定し（Ｓ１０５でＹＥＳ）、反復計算を終了してそのときの誘電率分布を最終的な結果とする。他方、目的関数の感度が０ないしその近傍でない場合には、未だ収束していないと判定し（Ｓ１０５でＮＯ）、目的関数の感度に基づく設計変更、すなわち新たな誘電率分布を再設定（Ｓ１０６）して再びＳ１０２以降の処理を繰り返す。繰り返しループの中に感度解析の処理が含まれていることから、誘電率分布が更新される毎に、これに応じて感度も再計算されて順次変化していく点に留意されたい。

以下、各処理についてより詳細に説明する。
＜初期値設定＞
まず、適当な初期値として、初期誘電率分布εを設定する。そして、この初期誘電率分布εは、電位φ（静電ポテンシャル）を求めるために用いられる。
＜順解析＞
静電容量Ｃ_ｉは、電位φ、誘電率ε、センサ１２に印加する電位差Ｖ_ｃを用いて、センシング電極２０の境界Γ_ｉでの境界積分として算出される。すなわち、
である。

電位φ_ｉは、連続変数で表した媒体の誘電率εを用いて、コンピュータにより静電場の数値シミュレーションで求めることができる。静電場の支配方程式は、
である。（数２）式をｎ回の計測条件と同じ境界条件の下で数値シミュレーションし、その結果得られる電位φ_ｉから（数１）式を用いて静電容量Ｃ_ｉを求める。数値シミュレーションに必要なプログラムは、予めコンピュータのメモリに記憶される。

なお、数値シミュレーションでは、公知の有限要素法（ＦＥＭ）等で離散化して計算する。離散化した支配方程式は、誘電率分布から決まるマトリクスＥ、電位の分布を表すベクトルΦ、及び境界条件を離散化して生じるベクトルＱを用いて、行列式として、
と表現される。
＜目的関数計算＞
次に、複数のセンシング電極２０を有するセンサ１２で、その組合せの数だけｎ回計測した場合、計測した静電容量Ｃ_ｉ ^０と、数値シミュレーションで計算した静電容量Ｃ_ｉの最小二乗誤差ｆ_１を目的関数として定義する。
（数１）式を考慮して、境界の誘電率で構成されたベクトルＫを用いて離散化すると、
となる。ここで、Φは上記のように電位の分布を表すベクトルである。
＜正規化＞
静電容量Ｃ_ｉ ^０、Ｃ_ｉは、数値シミュレーションで基準となる所定の誘電率ε^１、ε^２を与えて計算した静電容量Ｃ_ｉ ^１、Ｃ_ｉ ^２を用いて、
のように正規化する。ここで、＊は正規化された量であることを示す。

計測又は計算で求めた静電容量Ｃ_ｉは、静電容量の計測に利用する電極の位置に応じてその大きさが異なるため、管路内の誘電率分布の変化が静電容量に明確に表せない。そこで、静電容量Ｃ_ｉから所定の誘電率分布ε^１を与えたときの静電容量Ｃ_ｉ ^１を差し引き、さらに、その値を静電容量Ｃ_ｉ ^１ともう一つの所定の誘電率分布ε^２から計算した静電容量Ｃ_ｉ ^２の差で除することで、すべての静電容量の大きさを同程度にでき、管路内の誘電率分布の変化を静電容量で明確に表わすことができる。

（数６）式のように静電容量Ｃ_ｉ ^０、Ｃ_ｉを正規化した場合、(数４)式は、
となる。

(数１)式を考慮して、境界の誘電率と正規化の効果を含んだベクトルＫ^＊を用いて離散化すると、
となる。
＜感度解析＞
次に、反復計算の収束判定、及び誘電率の更新量を決定するために、ｆ_１ ^＊の感度を算出する。ｆ_１ ^＊の感度は、随伴変数法を用いて算出される。ｆ_１ ^＊に含まれる静電容量が誘電率分布εと電位分布φに依存することから、
となる。この式において、∂Φ／∂εは、計算により直接求めることは困難である。そこで、（数３）式に示した支配方程式にベクトルΛを乗じた項を、（数６）式に加える。
この式において、ＥＱ＋Φは本来その値は０であるので、Ｆとｆ_１ ^＊は同値であり、（数１０）で定義したＦの感度は、ｆ_１ ^＊の感度と同値になる。

Ｆの感度は、
となる。ベクトルΛの任意性から、右辺第２項が０となるようにベクトルΛを決定すると、
である。このときのベクトルΛを用いると感度を計算できる。
(数１２)式と（数１３）式に（数８）式を用いれば、
となる。上式で用いる値は、既知量あるいは数値シミュレーションの結果として求められる値である。したがって、１回のマトリクス演算で感度を求めることができ、感度計算の時間が短縮化できる。離散化した要素１つ１つの値を変更し、その都度数値シミュレーションを実施して感度を計算する場合と比べると、計算時間を大幅に短縮できる。
＜収束判定＞
（数１１）式で求めた感度の値が０ないしその近傍となったときに、そのときの誘電率分布がｆ_１ ^＊を最小とする分布であるとし、収束したものとみなして反復計算を終了する。
目的関数ｆ_１ ^＊の感度が０ないしその近傍であれば目的関数ｆ_１ ^＊の極値となる誘電率分布ε^ｏｐｔを求めたことになる。ε^ｏｐｔは、目的関数ｆ_１ ^＊が凸関数であれば全ての誘電率分布の中で目的関数ｆ_１ ^＊を最小とする値であり、非凸関数であれば誘電率分布ε^ｏｐｔ近傍において最小となる値である。従って、本実施形態では、目的関数ｆ_１ ^＊の感度と数学的に等価なＦの感度が０ないしその近傍となることを反復計算の終了条件とすることで、反復計算で得られた誘電率分布εがその近傍において目的関数ｆ_１ ^＊を最小とする値であることを保証する。
＜感度に基づく設計変更＞
他方、反復計算の終了条件を満足しないとき、反復計算のｋ回目の誘電率ε^ｋは、
のように更新する。ここで、αは更新量を決めるパラメータである。αを十分小さくとれば確実に目的関数ｆ_１ ^＊を小さくするように誘電率を更新できる。反復計算を効率的に行うためには、数理計画法に従ってαを決定する。また、感度解析を反復計算中に毎回行うことで、正確な感度に基づいて目的関数値を減少させるように誘電率分布が更新できる。
＜明確化＞
誘電率εを連続変数で扱う場合、計測対象の媒体の誘電率とは異なる中間的な誘電率が許容されることになる。これをできる限り回避するために、管路１０内の媒体が２種類の媒体１、媒体２（例えば気体と液体）からなるときに、これらの媒体の誘電率ε^１、ε^２と計算で求められる誘電率εの差を評価する関数ｆ_２
を（数７）式に加えて、
で定義される関数ｆが最小となる誘電率εを求める。但し、ε^１＜ε^２とする。

目的関数ｆ_１ ^＊に加えて、連続変数で表された誘電率分布εと管路１０内部に分布する２つの媒体１，２の誘電率の差を小さくするような関数ｆ_２を付加することで、誘電率分布を２つの媒体１，２の誘電率のいずれかに近づけた値を求めるようになり、誘電率分布を明確化することができる。関数ｆ_２が誘電率分布を明確化する関数である。
＜平滑化＞
一般的に、計測するｎ個の静電容量が、誘電率分布を表すベクトルΦの要素の数よりも少ないため、誘電率分布が一意に定まらない。そこで、誘電率の空間勾配の二乗の領域積分量ｆ_３ ^＊
を（数７）式に加えて、
で定義される関数ｆが最小となる誘電率分布εを求める。

ｎ回計測された静電容量から、数値シミュレーションで計算する領域に設定したｍ個の誘電率を決定する場合、ｎ＝ｍであれば誘電率は一意に定まるが、通常ｎ＜ｍであるため解が不定となり、上記のように誘電率が一意に定まらない。従って、不足する情報を補う必要がある（正則化）。目的関数ｆ_１ ^＊を正則化するために、特異な誘電率分布が得られないような関数ｆ_３を目的関数ｆ_１ ^＊に付加することで、近似的な解として誘電率分布を求めることができるようになり平滑化できる。関数ｆ_３ ^＊が誘電率分布を平滑化する関数である。

関数ｆ_２及び関数ｆ_３ ^＊は、独立に評価される関数である。そこで、（数１８）式に（数１９）式を加えて、
で定義される関数ｆが最小となる誘電率εを求めるようにすることで、明確化と平滑化をともに実現できる。明確化する関数ｆ_２を付加すると特異な誘電率分布を誘発する可能性があり、他方で、平滑化する関数ｆ_３ ^＊を付加すると特異な誘電率分布を抑えているため媒体の界面が明確にならない可能性がある。これらを同時に目的関数ｆ_１に付加することで、特異な誘電率分布を抑えつつ、明確化した誘電率分布を再構成できる。

（数１８）式、（数２０）式、及び（数２１）式のパラメータβ、κは問題に応じて適宜設定する。例えば、初期誘電率分布をε^０としたときのｆ_１ ^＊に対して、ｆ_２、ｆ_３ ^＊がｆ_１ ^＊の１０％程度の大きさになるようにβ、κを決定する。

このようにして、誘電率分布Ｅを再構成することができる。再構成画像出力手段１８は、誘電率分布再構成手段１６で再構成された誘電率分布Ｅを画像として表示する。

なお、目的関数に対して明確化の関数ｆ_２及び平滑化の関数ｆ_３ ^＊を付加しているが、必要に応じて、明確化の関数ｆ_２のみを目的関数に付加する、あるいは平滑化の関数ｆ_３ ^＊のみを目的関数に付加することも可能である。明確化の関数ｆ_２のみを付加して再構成した画像と、平滑化の関数ｆ_３ ^＊のみを付加して再構成した画像と、関数ｆ_２及び関数ｆ_３ ^＊をともに付加して明確化及び平滑化して再構成した画像をいずれもディスプレイ上に表示し、相互に比較可能に並列表示することも可能である。再構成画像を表示する際に、ユーザが適宜、明確化、平滑化を画面上で選択できるようにし、ユーザにより選択された演算のみを実行してその結果を表示してもよい。

比較のため、図４に従来の誘電率分布Ｅの再構成の解析のフローチャートを示す。従来の解析手法においては、まず感度行列、つまりセンシティビティマップを作成する（Ｓ２０１）。そして、管路１０内の初期誘電率分布を設定し（Ｓ２０２）、順解析を行って静電容量を算出する（Ｓ２０３）。順解析により得られた静電容量と、計測して得られた静電容量との差（最小二乗誤差）を用いて収束判定し（Ｓ２０４）、収束していなければＳ２０１で作成した感度行列に基づいて誘電率を再設定し（Ｓ２０５）、Ｓ２０３以降の処理を繰り返す。

従来の解析手法では、感度行列は最初に作成され、誘電率分布によらずに反復計算において常に一定である。しかし、本来管路内部の誘電率εを反復計算により求める際には誘電率分布が常に異なり、誘電率εが変化したときのｆ_１への影響度を示すｆ_１の感度は、管路内部の誘電率分布に応じて変化する。したがって、従来の解析手法では必ずしも適切な感度に基づいて誘電率分布を推定していない。その結果、最終的に得られる誘電率分布の画像は本実施の形態で説明した手法に比べて不鮮明になりやすい。ただし、本実施の形態では、ｒ個のセンシング電極２０のうち、ｐ個（ただし、１≦ｐ≦ｒ−１）を印加電極として選択して測定用電圧を印加し、測定用電圧を印加しなかったセンシング電極２０のうちｑ個（ただし、１≦ｑ≦ｒ−ｐ，３≦ｐ＋ｑ）を選択して接地することによって、静電容量の計測ためのセンシング電極２０の組み合わせの数を増やしており、これによって従来手法を用いた場合においても誘電率分布の画像をより鮮明にすることができる。

図５は、従来の解析手法でセンシング電極２０の数ｒ＝８個のときに印加電極の数ｐ＝１、接地電極の数ｑ＝７として、組み合わせ５６通りの計測を行って得られたｎ＝５８個の静電容量Ｃ_ｉ ^０を用いて誘電率分布を再構成した結果を示す。計測対象として（ａ）、（ｂ）２つのパターンをそれぞれ示す。なお、図５では、比較のために、センシング電極２０の数ｒ＝８個のときに印加電極の数ｐ＝１、接地電極の数ｑ＝１として、組み合わせ２８通りの計測を行って得られたｎ＝２８個の静電容量Ｃ_ｉ ^０を用いて誘電率分布を再構成した結果も示している。

図６は、図５（ａ）及び（ｂ）に示した測定対象より複雑な媒体の分布について計測及び誘電率分布の再構成の解析を行った例を示す。図６に示されるように、従来の解析手法と本実施の形態における解析手法のいずれにおいても組み合わせ２８通りの計測よりも５６通りの計測を行った方がより正確かつ鮮明な誘電率分布を再構成することができた。また、同じ組み合わせ数を比較すると、従来の解析手法より本実施の形態における解析手法を適用した方が正確かつ鮮明な誘電率分布を再構成することができた。

以上のように、本実施の形態によれば、従来の混相状態分布計測装置よりもセンシング電極の組み合わせ数を増加させることができ、これに伴って計測できる静電容量の数を増加させることができる。したがって、測定対象の誘電率分布をより正確かつ鮮明に再構成することができる。

また、本実施の形態によれば、感度を一定とするのではなく、感度が誘電率分布に応じて変化し得ることに鑑みて、反復計算毎に得られる誘電率分布に応じてその都度感度を再計算しているため、再構成画像の精度が向上する。また、目的関数の正規化、明確化及び平滑化を行うことで再構成画像の精度をさらに向上させ、鮮明な画像を得ることができる。

なお、本実施の形態において、静電容量に代えて電気抵抗を用いてもよい。また、印加電極の数ｐ及び接地電極の数ｑの組み合わせは一通りに限定されるものでなく、複数通りの組み合わせについて静電容量の計測を行い、その結果に基づいて誘電率分布の解析を行ってもよい。例えば、ｐ＝１及びｑ＝７としてｎ＝５６通りの静電容量を測定すると共に、ｐ＝２及びｑ＝１としてｎ＝１６８通りの静電容量を測定し、すべての測定結果を利用して誘電率分布の解析を行ってもよい。

１，２ 媒体、１０ 管路、１０ａ 計測空間、１２ センサ、１４ 静電容量計測手段、１６ 誘電率分布再構成手段、１６ａ 演算装置、１６ｂ 静電容量正規化手段、１６ｃ 誘電率明確化手段、１６ｄ 誘電率空間勾配最小化手段、１８ 再構成画像出力手段、２０ センシング電極、２２ リード線、１００ 混相状態分布計測装置。

Claims (5)

1. 管路の周辺に配置されたｎ個（ただし、ｎは３以上）のセンシング電極を有するセンサと、
   前記センサを用いて静電容量又は電気抵抗を計測する計測手段と、
   を備え、
   前記センシング電極のうち、ｐ個（ただし、１≦ｐ≦ｎ−１）の電極に測定用電圧を印加し、ｑ個（ただし、１≦ｑ≦ｎ−ｐ，３≦ｐ＋ｑ）の電極を接地した状態で前記計測手段により静電容量又は電気抵抗を計測することを特徴とする混相状態分布計測装置。
2. 請求項１に記載の混相状態分布計測装置であって、
   ｐ又はｑは、少なくとも２以上に変化させることを特徴とする混相状態分布計測装置。
3. 請求項１又は２に記載の混相状態分布計測装置であって、
   前記管路内の媒体の分布を示す画像を生成する画像再構成手段と、
   前記画像を表示する画像表示手段と、
   を備えることを特徴とする混相状態分布計測装置。
4. 請求項１〜３のいずれか１項に記載の混相状態分布計測装置であって、
   前記計測手段で得られた静電容量又は電気抵抗と一致する前記管路内部の媒体の誘電率分布を反復計算により算出する際に、誘電率分布の変化が静電容量又は電気抵抗に与える影響（感度）を誘電率分布の更新毎に算出し、算出して得られた感度を用いて反復計算を実行することを特徴とする混相状態分布計測装置。
5. 請求項１〜４のいずれか１項に記載の混相状態分布計測装置であって、
   ｎ＝８，ｐ＝１，ｑ＝７の状態を含むことを特徴とする混相状態分布計測装置。