(57)【要約】【課題】多重材料物体を電気ネットワークの形で表現し解析する方法を提供すること。【解決手段】多重材料物体の特性を決定するための方法（１２）を提供する。本方法は、多重材料物体を取り囲む電極に個々の電気信号パターンをもつ印加電気信号組を提供することによって回転電界を発生させるステップ（１２２）を含む。本方法はさらに、印加した各電気信号パターンに対応する電極から電気信号の計測電気信号を取得するステップ（１２４）を含む。電気ネットワークは、印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて決定される。本方法はさらに、電気ネットワークを解析することによって多重材料物体の特性を決定するステップ（１２８）を含む。【選択図】図２
2011141275.tif 000002
【特許請求の範囲】
【請求項１】
多重材料物体の特性を決定する方法（１２０）であって、
回転電界を発生させるために物体を取り囲む電極に個々の電気信号パターンをもつ印加電気信号組を提供するステップ（１２２）と、
各印加電気信号パターンに対応する電極において電気信号の計測電気信号組を取得するステップ（１２４）と、
印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて物体の電気ネットワーク表現を決定するステップ（１２６）と、
電気ネットワークを解析することによって多重材料物体の特性を決定するステップ（１２８）と、
を含む方法。
【請求項２】
前記個々の電気信号パターンの最大数が多くとも電極の半数に等しく、かつ前記電気ネットワークが純抵抗性要素からなる、請求項１に記載の方法。
【請求項３】
電気ネットワークを決定する前記ステップは、印加電気信号組をフル印加電気信号組に変換するステップと、計測電気信号組をフル計測電気信号組に変換するステップと、を含む、請求項２に記載の方法。
【請求項４】
印加電気信号組を変換する前記ステップは既知の電気信号パターンを利用するステップを含む、請求項３に記載の方法。
【請求項５】
計測電気信号組を変換する前記ステップは電気ネットワークの純抵抗的性質を利用するステップを含む、請求項３に記載の方法。
【請求項６】
特性を決定する前記ステップは個々の電気信号パターンの最大数が多くとも電極の半数に等しいときに新たな電気ネットワークを決定するステップを含む、請求項１に記載の方法。
【請求項７】
新たな電気ネットワークを決定する前記ステップは電気ネットワークの抵抗値に基づいて電気ネットワークを修正するステップを含む、請求項６に記載の方法。
【請求項８】
特性を決定する前記ステップは、印加電気信号組を用いてシミュレーションし実際の計測電気信号組とシミュレーションの計測電気信号組の間の誤差を最小化することによって新たな電気ネットワークを完成させるステップを含む、請求項６に記載の方法。
【請求項９】
多重材料物体の電気ネットワーク表現を決定する方法であって、
物体内部に回転電界を発生させるために物体を取り囲む電極に対して個々の電気信号パターンをもつ印加電気信号組を提供するステップと、
各印加電気信号パターンに対応する電極において電気信号の計測電気信号組を取得するステップと、
印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて多重材料物体の電気ネットワーク表現を決定するステップと、
を含む方法。
【請求項１０】
回転電界を発生させるように多重材料物体を取り囲む電極に個々の電気信号パターンをもつ印加電気信号組を提供するための電源（６２）と、
各印加電気信号パターンに対応する電極において電気信号の計測電気信号組を取得するための計測ユニット（６２）と、
印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて物体の電気ネットワーク表現を計算し多重材料物体の特性を評価するための処理回路（５０）と、
を備える多重材料検知システム（４０）。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、全般的には多重材料物体内の個々の材料の組成及び特性を決定するための方法に関し、また具体的には、多重材料物体を電気ネットワークの形で表現し解析することに関する。
【背景技術】
【０００２】
電気インピーダンス断層撮像法（ＥＩＴ）は、多重材料物体内部の材料分布に関する定量的マッピングに使用可能な侵害性を最小とした計測技法の１つである。ＥＩＴでは、導電率と誘電率のマップを使用し当該物体内部の様々な材料からなる分布を推定する。物体を取り囲む電極を介して物体に対して様々な電流パターンまたは電圧パターンが印加され、これに対応する電圧や電流が計測される。この電流−電圧の関係に基づいて、内部インピーダンスまたは内部アドミッタンス分布が決定される。
【０００３】
多重材料物体の一例は、パイプやコンジットの内部を少なくとも２つの材料または相が一緒に流れているような混相流（ｍｕｌｔｉｐｈａｓｅ  ｆｌｏｗ）である。混相流処理は、例えば石油業界、薬品業界、食品業界、化学業界を含む多種多様な産業分野にとって重要である。これらの種類の混相流処理では既存の処理装置や新規の処理装置に対する設計の改善及び動作効率の上昇を可能にするためにその内部特性に関する直接的な知見が必要である。混相処理のパフォーマンスの予測に使用される特性には例えば、空間的な相分布（空間体積相比率）、フローレジーム（ｆｌｏｗ  ｒｅｇｉｍｅ）、界面面積、並びに絶対速度及び相間や材料間での相対速度を含むことがある。材料の分布が不均一であると化学変化や変換に利用可能な材料間の界面面積を低減させる傾向となると共に空間的に不均一な反応ゾーンや濃度を生成させるような再循環フローを生じさせることがあるため、材料の空間分布を知ることは特に有用である。さらに体積相比率及び速度が、混相流の適正かつ適時の制御を可能にする重要なパラメータである。
【０００４】
ＥＩＴ技法の１つでは、一方の電極に電流が流入しかつもう一方の電極では電流が流出しているような対となった境界電極対に１度に１対ずつ電流が印加されると共に、これらの電極のすべてにおいて電圧が計測される。類似の一技法では、境界電極対に対して１度に１対ずつ電圧が印加されると共に、すべての電極において電流が計測される。これらの技法に関連する問題点の１つは、信号対雑音比が低いことである。さらに、比誘電率や導電率の変動が小さい材料を対象とする場合、分解能が低下する。
【０００５】
別のＥＩＴ技法では、計測完了に必要なデータを生成するように電極のすべてに同時に電流または電圧が印加される。電極に印加される電流または電圧はすべて、互いに電気的に同相でありかつ異なる振幅を有している。しかしこの技法では、電極に印加する電流または電圧パターンの数が多くなる（典型的には、電極数より１つだけ少ない）ためそれだけ時間がかかる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許第７１６７００９号
【発明の概要】
【発明が解決しようとする課題】
【０００７】
したがって、上述の問題に対処するような方法及びシステムを提供することが望ましい。
【課題を解決するための手段】
【０００８】
本発明の一実施形態では、多重材料物体の特性を決定する方法を提供する。本方法は、多重材料物体内部に回転電界を発生させるために多重材料物体を取り囲む電極に個々の電気信号パターンをもつ印加電気信号組を提供するステップと、各印加電気信号パターンに対応する電極から電気信号の計測電気信号組を取得するステップと、を含む。本方法はさらに、印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて電気ネットワークを決定するステップを含む。本方法はさらに、電気ネットワークを解析することによって多重材料物体の特性を決定するステップを含む。
【０００９】
別の本発明の実施形態では、多重材料物体の電気ネットワークを決定するための方法を提供する。本方法は、回転電界を発生させるために個々の電気信号パターンをもつ印加電気信号組を電極に提供するステップを含む。本方法はさらに、各印加電気信号パターンに対応する電極から電気信号の計測電気信号組を取得するステップと、印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて電気ネットワークを決定するステップと、を含む。
【００１０】
さらに別の本発明の実施形態では、多重材料検知システムを提供する。本システムは、物体内部に回転電界を発生させるように多重材料物体を取り囲む電極に個々の電気信号パターンをもつ印加電気信号組を提供するための電源と、各印加電気信号パターンに対応する電極から電気信号の計測電気信号組を取得するための計測ユニットと、を含む。本システムはさらに、印加電気信号組、計測電気信号組及び印加電気信号組の逆行列に基づいて多重材料物体の特性を評価するための処理回路を含む。
【００１１】
本発明に関するこれらの特徴、態様及び利点、並びにその他の特徴、態様及び利点については、同じ参照符号が図面全体を通じて同じ部分を表している添付の図面を参照しながら以下の詳細な説明を読むことによってより理解が深まるであろう。
【図面の簡単な説明】
【００１２】
【図１】例示の一実施形態による石油産生施設のブロック図である。
【図２】本発明の一実施形態で使用される電気インピーダンス断層撮像（ＥＩＴ）式の混相流量計の概要図である。
【図３】本発明の一実施形態による電圧印加型ＥＩＴシステムのブロック図である。
【図４】本発明の一実施形態による回転磁場パターンを表した図である。
【図５】本発明の一実施形態によるコンジット内の多重材料物体を取り囲む電極の概要並びに基準正弦波形を表した図である。
【図６】本発明の一実施形態による電気抵抗断層撮像（ＥＲＴ）技法を用いて多重材料物体の特性を決定する方法を表した流れ図である。
【図７】本発明の一実施形態によるＥＩＴ技法を用いて多重材料物体の特性を決定する方法を表した流れ図である。
【図８】本発明の一実施形態による図６のＥＲＴ技法を用いて多重材料物体の特性を決定する別の方法を表した流れ図である。
【発明を実施するための形態】
【００１３】
以下で詳細に検討するように、本発明の実施形態は多重材料物体の特性を計測するためのシステム及び方法を提供するように機能する。混相流の一例における物体の特性には、コンジット内を流れる石油、水及び気体（このコンテキストでは、気体状の炭化水素について言及）に関する体積率及び流量を含むことがある。本発明について石油／気体／水の計測値の使用に関連して記載しているが、本発明はこうした用途に限定されるものではなく、本発明の態様は例えば、がん診断や水処理などの広く多種多様な産業上の処理、ヘルスケア処理及び化学処理に用途を見出すことができる。こうした多重材料物体の特性は混相流に関して得られる例とかなり異なることがある。
【００１４】
図１は、例示の一実施形態に従った石油産生施設１０を表している。この石油産生施設は典型的には、その各々が配管系１４に相互接続された複数の油井１２を含む。配管系１４は混相流量計（ＭＰＦＭ）１８に結合させた産生物分岐管１６を含む。混相流量計は、油井の非常に近くにおける未処理の油井ストリームの計測を可能にし、これによりリザーバ管理の向上に使用し得る油井パフォーマンスの連続監視の提供が可能となる。油井１２から汲み出された流体は、産生物分岐管１６を通して不純物分離装置２０に送られる。施設１０内において追加的な試験分離装置（図示せず）がＭＰＦＭに加えて用いられること、あるいはこれに代えて用いられることがあることに留意すべきである。試験分離装置に対するＭＰＦＭの利点の１つは、計測実行に要する時間が短縮されることである。試験分離装置では油井が変化したときに充填及び安定化をさせなければならないが、ＭＰＦＭでは油井流体の変化により迅速に応答できると共に安定化に要する時間がより短い。
【００１５】
不純物分離装置２０は、油井から汲み出された石油、気体及び水を分離させる。不純物分離装置２０は、１つまたは複数の計測機器を含むことがある。この計測機器には例えば、油井から抽出される水の量または流量を計測するための水メーターや、油井から抽出される石油の量を計測するためのエマルジョン（ｅｍｕｌｓｉｏｎ）メーターを含むことがある。さらに計測機器は、油井ヘッド圧力センサ、温度計、塩分計、ｐＨ計など油井パフォーマンスの監視に典型的に利用されるような別の機器を含むことがある。
【００１６】
図２は、電気インピーダンス断層撮像（ＥＩＴ）式のＭＰＦＭシステム４０を表している。ＥＩＴでは、物体を取り囲む電極を介して実施した電気計測から多重材料物体内部の導電率または誘電率分布が推定される。多重材料物体の周りに導電性の電極が取り付けられ、この電極の幾つかまたはすべてに交番電流または電圧が印加される。得られる電気ポテンシャルまたは電流が計測されると共に、印加する電流及び／または電圧の構成またはパターンを異ならせて処理が反復される。
【００１７】
図２のＭＰＦＭシステム４０は、コンジット４６の周りに分布させた多数の電極４４を備えた電極アレイ４２を含む。このコンジットは複数の材料やその内部の相を保持するパイプやタンクなどの容器を成すことがあり、またこのコンジットは人体の一部や人の全身など別の受け容れ体を成すことがある。より具体的な実施形態ではその電極数は、コンジットのサイズ及び必要とする確度に応じて８個、１２個または１６個とすることがある。これらの電極は、良好な電気的接触を保証するのに適したコーティングの必要に応じた使用を含むような一実施形態によるコンジットの内部壁に直接取り付けられることがある。電極とコンジット壁の間には、適当な電気絶縁材料を設けることがある。これらの電極は、電流または電圧源、Ｄ−Ａ変換器、Ａ−Ｄ変換器、差動増幅器、濾過器、ディジタルマルチプレクサ、アナログマルチプレクサ、クロック、及び／またはコンピュータ５０に結合させたディジタルＩ／Ｏユニットなどの構成要素を含み得る電子コンディショニング回路４８に接続されている。一実施形態ではコンピュータ５０は、画像再構成処理で使用されるディジタル信号プロセッサカード及び画像の表示に適したディスプレイ５２を装備したパーソナルコンピュータを含む。別の実施形態では、現場プログラム可能ゲートアレイ（ｆｉｅｌｄ−ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙ：ＦＰＧＡ）やプログラム可能複合論理デバイス（ｃｏｍｐｌｅｘ  ｐｒｏｇｒａｍｍａｂｌｅ  ｌｏｇｉｃ  ｄｅｖｉｃｅ：ＣＰＬＤ）などの別の処理回路を使用することもある。本発明の一実施形態では電極４４は、コンジット内の多重材料物体内部に回転電界を生成するような電圧パターンからなる電圧の組を電極の両端に印加することによって付勢を受ける。ある任意の時点において、１つまたは複数の電圧源を用いて電極にある電圧パターンを印加しており、電極の各々において対応する電流信号の組が計測される。一実施形態では、電圧の印加に代えて１つまたは複数の電流源を用い電流を電極内に注入することによりこれを付勢させており、その電極の両端において対応する電圧が計測される。
【００１８】
図３は、Ｌ個の電極を有する電圧印加型ＥＩＴシステムのブロック図６０である。各電極６１は、印加電圧を発生させるための電圧源６２、電流を計測する電流計、並びに印加電圧を直接計測する電圧計を含む回路に接続されている。システム全体に対する単一の基準に合わせた較正を可能とさせるために、切替ネットワーク６３によって、電圧源／電流計／電圧計回路のいずれに対しても単一の較正回路６４を接続することを可能としている。電圧源（電流計及び電圧計を伴う）並びにスイッチに対してディジタル制御器（図示せず）がインタフェースされることがあり、また較正回路によってシステム構成が設定されると共に電圧及び電流のディジタル計測値が収集される。電極に対して電圧に代えて電流が印加されるような別の実施形態では、例えば直流源または電圧対電流変換器系を備え得る電流源を用いて電流を発生させている。
【００１９】
電流源と電圧源のいずれの実施形態においても、得られた計測値はコンピュータ５０（図２）により処理されると共に、多重材料物体内部の電気インピーダンスまたはアドミッタンス分布に関する表現が決定される。この電気インピーダンスまたはアドミッタンス分布は次いで、多重材料物体の特性を提供するようにコンピュータ５０によってさらに解析を受ける。多重材料物体の特性は例えば、システムの材料の組成や分布を含むことがある。またこの電気インピーダンスまたはアドミッタンス分布は、多重材料物体の個々の材料のフローレジーム、相比率及び速度を決定するためにさらに解析を受けることもある。フローレジームは、例えばバブリー（ｂｕｂｂｌｙ）フロー、チャーン（ｃｈｕｒｎ）フロー、スラグ（ｓｌｕｇ）フロー、アニュラー（ａｎｎｕｌａｒ）フロー（ただし、これらに限らない）を含むことがある。
【００２０】
電圧源の実施形態では、ある電極から別の電極まで流れる電流は、すべての電極両端に印加された相対電圧並びにすべての電極同士の間に存在する材料の導電率及び誘電率の関数である。この材料は例えば、石油だけのことや、石油と気体の混合体のこともある。材料とその分布に応じて、すべての電極同士の間のインピーダンスまたはアドミッタンスは異なると共に、電極間を流れる電流も異なる。したがって、印加した電圧組と計測した電流組から、各電極対間のインピーダンス組または該インピーダンスからなるインピーダンス行列を計算することが可能である。同様に別法としてあるいは追加として、印加した電圧組と計測した電流組から各電極対間のアドミッタンス組またはインピーダンスからなるアドミッタンス行列を計算することが可能である。インピーダンスやアドミッタンスは電極間の材料の導電率や誘電率の関数であり、インピーダンスまたはアドミッタンスネットワークの解析によって、材料分布及びその特性を決定することができる。
【００２１】
電圧パターンからなる電圧組によってコンジット内に回転電界が生成される。回転電界は、印加信号周波数またはその整数倍に対応する一定の角速度でその方向が変化する電界である。このことは、その回転磁界が一定の角速度で方向を変化させる多相交流モータを動作させる基本原理である回転磁界と同様である。
【００２２】
図４は、様々な空間高調波パターンに関してパイプ内部のシミュレーションを行った幾つかの例示的な電界を表している。プロット５６、５７、５８及び５９はそれぞれ、第１、第２、第４及び第８高調波に関する電界を示している。第１高調波パターンの電界は、平行の直線からなり、第２高調波パターンの電界は直角の双曲線の群からなり、第４及び第８高調波パターンの電界曲線は、双曲線ではないがこれと幾らか類似したものとなる。
【００２３】
電圧パターンからなる電圧組を電極に印加することによって回転電界を生成させることがある。パイプ全体にわたりＬ個の電極があると仮定すると、これら電極間に印加される電圧パターンは次式で表すことが可能である。
【００２４】
【数１】
2011141275.tif 000003
【００２５】
【数２】
2011141275.tif 000004
【００２６】
【数３】
2011141275.tif 000005
上式において、
【００２７】
【数４】
2011141275.tif 000006
は
【００２８】
【数５】
2011141275.tif 000007
の励起パターンで第ｋ番目の電極に印加される電圧であり、
【００２９】
【数６】
2011141275.tif 000008
はピーク印加電圧であり、ｔは時間であり、またωは正弦波の角周波数である。式（２）及び（３）は式（１）についてプラスとマイナスの関係を分離して展開させたものである。式（２）は物体内部に時計回りの回転電界を発生させる順電圧パターンであり、一方式（３）は物体内部に逆時計回りの回転電界を発生させる逆電圧パターンである。各励起パターンは、式（１）及び（２）で与えられるすべてのパターン組が完了するまで順次電極に印加される。その励起パターンに対応する変数
【００３０】
【数７】
2011141275.tif 000009
は１からＮｐまで変化する（ここで、Ｌが偶数の場合はＮｐ＝Ｌ／２であり、またＬが奇数の場合はＮｐ＝（Ｌ−１）／２である）。一実施形態では電気抵抗断層撮像法（ＥＲＴ）の場合、全部でＮｐ個の順電圧パターンまたはＮｐ個の逆電圧パターンが印加される。別の実施形態では、多重材料物体に関するより完全な情報を取得するために、順パターンと逆パターンの両方が電極に印加される。順電圧パターンと逆電圧パターンが一体となって、多重材料物体のネットワーク表現内の電気的要素の実数成分と虚数成分に関する情報が提供される。式（１）、（２）及び（３）と同様の方式で、順電流パターンと逆電流パターンを記述し計算することも可能である。
【００３１】
図５は、コンジットの周囲に沿って配置させた電極の概要図８０並びに振幅
【００３２】
【数８】
2011141275.tif 000010
の波形組８２を表している。波形組８２の内部には互いからの位相シフト範囲を等しくした８つの正弦波が存在していることに留意されたい。これらの波形にはＶ１〜Ｖ８の記号を付けており、その対応電極に印加される電圧信号に対応している。例えばＶ１は電極１に印加される電圧信号に対応し、またＶ７は電極７に印加される電圧信号を意味している。波形組８２は１つの電圧パターンに対応する。こうしたパターンの幾つかを、式（１）に従って計算し電極に印加することがある。８個の電極の使用は単に一例であり、別に所望の任意の数の電極を使用し得る。一実施形態ではそのパイプ８４は、石油８６と水８８などの材料または相を含む。
【００３３】
本発明の一実施形態では、多重材料物体を成した材料のうちの幾つかまたはすべての誘電率または透磁率の効果は導電率と比べてかなり小さいことがアプリオリに知られることがあり、またこうしたケースにおいて多重材料物体を表現する電気ネットワークは抵抗性要素のみから構成されるものと云うことが可能である。これは、電気抵抗断層撮像法（ＥＲＴ）と呼ぶＥＩＴの１つの特異的ケースである。ＥＲＴのみが実行されるような実施形態では、相分布に関する完全な情報を取得するために印加される電圧または電流パターンの総数はＮｐである。このＮｐ個のパターンは時計回りの回転電界を発生させる順パターンと逆時計回りの回転電界を発生させる逆パターンのいずれにもなり得る。
【００３４】
図６は、ＥＲＴにおけるインピーダンスまたはアドミッタンス行列の決定のための方法１００を表している。ステップ１０２では、Ｎｐ個のパターンからなる第１の電気信号組が電極に印加される。印加した各第１の電気信号パターンに関するそれぞれの第２の電気信号パターンを有する対応する第２の電気信号組がステップ１０４で計測される。一実施形態では、第１の電気信号パターンが電圧組を成すことがあり、かつ第２の電気信号パターンが電流組を成すことがある。別の実施形態では、第１の電気信号パターンが電流組を成すことがあり、かつ第２の電気信号パターンが電圧組を成すことがある。例えば第１の電気信号組が電圧組でありかつ第２の電気信号組が電流組であるような実施形態では、Ｎｐ個の電圧パターンだけが電極に印加されるため、印加された電圧パターン組を表現する電圧行列と計測された対応する電流パターン組を表現する電流行列とを構築することができる。この電圧行列及び電流行列は、Ｌ個の横列とＮｐ個の縦列を有することになる。第ｎ番目のパターンに関する第ｍ番目の電極における電圧または電流は、電圧または電流行列の要素（ｍ，ｎ）によって与えられることになる。さらに同じ実施形態に関する一例として、印加されたＮｐ個の電圧パターンからなる組が式（２）で与えられる順回転パターンであると仮定してみる。Ｌが奇数であるケースでは、電圧行列の残りのエントリーを式（３）から決定することが可能であり、また得られた新たな電圧行列はＬ個の横列とＬ−１個の縦列を有することになる。Ｌが偶数であるケースでは、式（３）を使うことによって１つの追加の電圧パターンが得られると共に、新たな電圧行列がＬ個の横列とＬ個の縦列を有することになる。ここで、式（３）を用いて算定した第１のパターンは電圧パターンの印加電圧組における最終パターンである式（２）によって算定した最終パターンと事実上同じである。したがって、式（３）を用いて算定した第１のパターンに対応する縦列は除去され得る、またしたがって新たな電圧行列はＬ個の横列とＬ−１個の縦列を有する。同様に新たな電流行列の残りのエントリーは、ＥＲＴのケースではインピーダンスまたはアドミッタンスネットワーク内の要素がすべて純抵抗性であるため第ｉ番目の印加電圧信号パターンに由来する計測電流パターン内の個々の電流が第（Ｌ＋１−ｉ）番目の電流パターンに対応する同じ電極における電流の複素共役になるということによって決定することができる。したがってステップ１０４では、第１と第２の電気信号行列が修正される。
【００３５】
ステップ１０８ではフル印加電気信号組、得られたフル電気信号組及びフル印加電気信号組の逆行列に基づいて、多重材料物体内部の材料分布を表現するインピーダンス行列またはアドミッタンス行列が決定される。一実施形態では、印加信号組の逆行列を決定するために擬似逆行列（Ｐｓｅｕｄｏ  ｉｎｖｅｒｓｅ）計算が利用されることがある。こうして計算された行列は、物体内部の材料分布に関する電気ネットワーク表現となる。これは、複素エントリーからなる対称行列である。換言すると、インデックスｉ及びｊのすべてについて、行列ａの第ｉ横列で第ｊ縦列の要素は該行列ａの第ｊ横列で第ｉ縦列の要素に等しい。
【００３６】
【数９】
2011141275.tif 000011
上式において、ａｉ，ｊは電極番号ｉ及びｊの間の等価インピーダンスまたはアドミッタンスである。次いでステップ１１０において、多重材料物体の特性を決定するために計算したインピーダンスまたはアドミッタンス行列が解析される。アドミッタンス行列では、電流を通過させる漏洩経路が存在ないとき、各横列ごとの横列総和はゼロでありかつ各縦列ごとの縦列総和はゼロである。このことを利用して、電極並びにそれ自体の付属回路から離れているシステム内に存在する任意の電流経路を検出することができる。
【００３７】
図７は、電圧源の実施形態のＥＩＴにおいてインピーダンスまたはアドミッタンス行列を決定するための方法１２０を表している。一実施形態ではＥＩＴの場合に、多重材料物体内部の材料に関する誘電率が既知でないとき、ステップ１２２において電極に対して順電圧パターンと逆電圧パターンの両方が印加される。ステップ１２４では、その電極から得られる電流信号が計測される。電圧パターンの印加ではなく、電流パターンを印加することも可能であり、電圧パターンが計測されることもあり得ることに留意すべきである。印加する電圧パターンの総数は電極数より１少ない数である（残りのパターンは事実上冗長性であるため）。したがって、電極がＬ個の場合、印加される総電圧パターンはＬ−１となる。次いでステップ１２６において、第１の電気信号組、第２の電気信号組並びに第１の電気信号組の逆行列に基づいて電気ネットワークが決定される。この電気ネットワークは、多重材料物体の特性を決定するためにステップ１２８においてさらに解析を受ける。
【００３８】
図８は、ＥＩＴにおいてインピーダンスまたはアドミッタンス行列を決定するための別の方法１４０を表している。この実施形態では、図６で説明したようなＥＲＴ手順が先ず実行される。ステップ１４２においてＮｐ個の電圧パターンからなる電圧組が電極に印加されると共に、ステップ１４４において得られる電流組が計測される。ステップ１４６では、サイズが（Ｌ，Ｌ−１）のフルの電圧及び電流行列にさせるように図６に関連して上述したようにして行列のうちの残り半分を満たすことによって電圧組と電流組が修正される。ステップ１４８では、印加電圧組、計測電流組並びに印加電圧組の逆行列に基づいて、アドミッタンス行列またはインピーダンス行列が決定される。さらに、ＥＲＴにより決定した電気ネットワーク内の純抵抗性要素の値に基づいて、物体内部の異なる材料の導電率及び分布は、数値モデルからの事前算定結果並びに純抵抗性要素のみを有するインピーダンスまたはアドミッタンス行列を用い、このステップ（行列内の実数）が容量性または誘導性要素（行列内の複素数）を有する新たなインピーダンスまたはアドミッタンス行列になるように新たに修正され終わるまで評価される。例えばＥＲＴが物体内部の３つの異なる箇所において３つの異なる導電率値を提供していれば、導電率値が１０−６ジーメンス／メートルの範囲にあれば、その材料は石油と同様でありしたがってその比誘電率は１〜５の範囲となり得ることを評価し得る。同様に導電率値が１０−３ジーメンス／メートル以上の範囲にあれば、その材料は水でありしたがってその比誘電率は７５と８０の間の値とし得ると決定し得る。比誘電率を評価し終えると、電気ネットワーク要素のうちの容量性部分を評価することが可能である。次いでこの新たなインピーダンスまたはアドミッタンス行列を既知の印加電圧パターンを用いてシミュレーションし、シミュレーション電流が決定される。この新たな電気ネットワークは、ステップ１４９において計測電流組とシミュレーション電流組の間の誤差に対する反復最小化によって完成させる。上で説明したように、完成させたこの電気ネットワークを次いで、ステップ１５０において解析し多重材料物体の特性が決定される。導電率値に基づいて、並びに別法または追加として例えば密度計からの情報の利用など別のセンサ情報の利用によって材料の比誘電率を評価する別の方法も存在し得る。さらに、例えば人の頭部内など多重材料物体内部の異なる材料の空間分布に関して事前の知見が存在することがある。ステップ１４９はこうした情報を用いてインピーダンスまたはアドミッタンス行列を、その要素に容量性または誘導性部分が追加されるように適当に修正している。
【００３９】
コンピュータ５０（図２）は、印加電圧パターン、計測電流信号並びに計算したインピーダンスまたはアドミッタンス行列に関する情報の保存、処理並びに通信を行う。電圧パターンは全体で１つの電圧組または印加信号組を形成しており、また計測電流は１つの電流組または計測信号組を形成している。上述したように別の実施形態では、電極内に電流が注入されることがあり、かつ電圧が計測されることがある。印加信号組が電圧組でありかつ計測信号組が電流組である場合、行列Ｐによって印加信号組を表現し得、また行列Ｑによって計測信号組を表現し得る。次いでそのインピーダンスネットワークが次式のような行列Ｚによって表現される、
Ｚ＝Ｐ＊Ｑ−１        （４）
ここで、Ｑ−１は行列Ｑの逆行列を意味している。しかし行列Ｑは計測信号組であり、また数学的に非適切な行列であるのが一般的である。したがって、この行列の逆行列を見出すのは困難である。したがってこのケースでは、次式で与えられるようにしてアドミッタンス行列を計算可能であるため、アドミッタンス行列の計算がより容易である。
【００４０】
Ｙ＝Ｑ＊Ｐ−１        （５）
印加信号組は計測信号組と比べて忠実度がより良好である。さらに、印加信号組は事前によく分かっているため、これをコンピュータ内に格納しておくことが可能であり、またその逆行列もコンピュータ内に格納しておくことが可能である。したがって、アドミッタンス行列を決定するためには、この利用可能な印加電圧信号の逆行列が計測信号組に対して直接乗算される。一実施形態ではさらに、印加信号組が電流組Ｍでありかつ計測信号組が電圧組Ｓであれば物体はインピーダンス行列Ｚによって表現し計算され得る。このケースではインピーダンス行列Ｚは次式で与えられる。
【００４１】
Ｚ＝Ｓ＊Ｍ−１        （６）
ここでも、Ｍ−１は行列Ｍの逆行列を意味している。ここでも同じく、印加信号組Ｍの逆行列の計算は計測信号組Ｓの逆行列計算と比べてより容易である。
【００４２】
式（５）と（６）のいずれにおいても、一実施形態では行列Ｐ及びＭの逆行列の計算のために擬似逆行列が利用される。擬似逆行列技法は例えば、片側（ｏｎｅ−ｓｉｄｅｄ）逆行列、Ｄｒａｚｉｎ逆行列、グループ（Ｇｒｏｕｐ）逆行列、Ｂｏｔｔ−Ｄｕｆｆｉｎ逆行列、またはＭｏｏｒｅ−Ｐｅｎｒｏｓｅ擬似逆行列を含むことがある。したがって数値方式や反復法を使用することなくインピーダンス行列またはアドミッタンス行列の直接計算を得るように逆行列計算技法が用いられる。次いで、式（５）及び（６）のアドミッタンス行列Ｙとインピーダンス行列Ｚを次式のように計算することができる。
【００４３】
Ｙ＝Ｑ＊（（ＰＨＰ）−１＊ＰＨ）        （６）
Ｚ＝Ｓ＊（（ＭＨＭ）−１＊ＭＨ）        （７）
上式において、ＰＨ及びＭＨは行列Ｐ及びＭの共役転置行列である。
【００４４】
幾つかの処理システムでは、多重材料系の電気ネットワーク表現におけるアドミッタンスまたはインピーダンス行列要素は時間と共に変動することがある。しかし、印加信号組は一定でありかつ計測信号組がリアルタイムで計測されるため、時変動するアドミッタンス行列またはインピーダンス行列も計算することができる。一実施形態では、異なるランダムフローレジームモデルに関する様々な事前算定のテーブルまたは行列をコンピュータ内に保存しておき、システムが電気ネットワーク表現を決定するときに単にこれを保存済みの行列と比較することによって多重材料物体内部のフローレジーム及び材料分布が特定される。別の実施形態では、多重材料物体の材料分布その他の特性を決定するためにヒューリスティックアルゴリズムによってアドミッタンスまたはインピーダンス行列が解析される。ヒューリスティックアルゴリズムは例えば、ファジィ論理アルゴリズムや遺伝的なアルゴリズムあるいはニューラルネットワークアルゴリズムを含むことがある。ファジィ論理アルゴリズムは、漠然とした情報、曖昧な情報、不明瞭な情報、雑音性の情報あるいは散逸した入力情報に基づいて結論に至るための簡便な方法を提供する。ファジィ論理アルゴリズムの結果は厳密ではなく近似である。さらにニューラルネットワークは、物体に関するもっとも確からしい材料分布その他の特性を特定するように多重材料物体の挙動履歴に基づいてトレーニングさせることがある。
【００４５】
本発明のある種の特徴についてのみ本明細書において図示し説明してきたが、当業者によって多くの修正や変更がなされるであろう。したがって添付の特許請求の範囲が、本発明の真の精神の範囲に属するこうした修正や変更のすべてを包含させるように意図したものであることを理解されたい。
【符号の説明】
【００４６】
１０  石油産生施設
１２  油井
１４  配管系
１６  産生物分岐管
１８  混相流量計（ＭＰＦＭ）
２０  不純物分離装置
４０  電気インピーダンス断層撮像方式のＭＰＦＭシステム
４２  電極アレイ
４４  電極
４６  コンジット
４８  電子コンディショニング回路
５０  コンピュータ
５２  ディスプレイ
６０  電圧印加型ＥＩＴシステムのブロック図
６１  電極
６２  電圧源
６３  切替ネットワーク
６４  較正回路
５６  電界プロット
５７  電界プロット
５８  電界プロット
５９  電界プロット
８０  コンジットの周囲に沿って配置させた電極の概要図
８２  電極に印加される波形組
８４  パイプ
８６  石油
８８  水
１００  ＥＲＴにおけるインピーダンスまたはアドミッタンス行列の決定方法
１０２  方法ステップ
１０４  方法ステップ
１０６  方法ステップ
１０８  方法ステップ
１１０  方法ステップ
１２０  ＥＩＴにおけるインピーダンスまたはアドミッタンス行列の決定方法
１２２  方法ステップ
１２４  方法ステップ
１２６  方法ステップ
１２８  方法ステップ
１４０  ＥＩＴにおけるインピーダンスまたはアドミッタンス行列の決定方法
１４２  方法ステップ
１４４  方法ステップ
１４６  方法ステップ
１４８  方法ステップ
１５０  方法ステップ
【図１】
2011141275.tif 000012
【図２】
2011141275.tif 000013
【図３】
2011141275.tif 000014
【図４】
2011141275.tif 000015
【図５】
2011141275.tif 000016
【図６】
2011141275.tif 000017
【図７】
2011141275.tif 000018
【図８】
2011141275.tif 000019