"当該構造46を通る電気信号の伝達に影響を与える様々な導電率又は誘電率を有する多数の領域を含む構造46の内部の画像を形成する方法であって、前記構造46を通る平面に沿って、前記構造46の外部表面に沿って配置される一連の電極44を設けるステップと、前記電極44への入力のための一連の電気信号集合であって、各々が前記電極44の各々への入力のための電気信号パターンを含むような一連の電気信号集合を与えるステップと、前記電極44に前記電気信号集合Viの各々を相次いで印加するステップであって、各々の信号集合の印加が前記構造46の内部に回転電場を発生する、印加するステップと、各々の信号集合の印加に応答して、前記電極44から測定電気信号集合を得るステップと、前記測定電気信号集合により、前記電極44が周囲に配置されている前記構造46の領域での前記構造46の内部を表わす導電率分布又は誘電率分布を算出するステップと、様々な導電率又は誘電率を有する2以上の領域を含む前記構造46の前記内部を表わす画像データを決定するステップとを備えており、前記導電率分布又は誘電率分布を算出する前記工程は、前記導電率分布に基づくエネルギ関数、及び前記測定電気信号集合を導電率に関する電圧値の変化の関数として表現する順モデルを定義することと、前記エネルギ関数と前記順モデルに基づく関数との間の差の偏導関数がゼロに等しくなることを要求することとを含んでいる、方法。前記一連の印加電気信号集合Viは各々電圧入力パターンを含み、前記測定電気信号集合は各々測定される電流レベルを含む、請求項1に記載の方法。前記一連の印加電気信号集合Viは電圧入力パターンを含み、前記測定電気信号集合は測定電流レベルを含む、請求項1に記載の方法。前記構造46は、様々な導電率を有する物質が流れるコンジットである、請求項1に記載の方法。前記構造46は、様々な相を有する物質が流れるコンジットである、請求項1に記載の方法。前記物質は、石油86及び炭化水素ガス88を含んでいる、請求項4に記載の方法。前記コンジット46を流れる液体物質及び気体物質の相対比を決定するために画像情報が用いられる、請求項4に記載の方法。前記印加電気信号集合は、λ個の電圧励起パターン集合を含んでおり、各々のパターンがLL=1~k個の電極44の一つに印加され、前記電圧パターンは、に従って発生され、式中、Vkλは第λの励起パターン集合として第kの電極44に印加される時間変化型電圧正弦波82であり、Vハット^記号は参照正弦波のピーク電圧であり、tは時間であり、ωは前記正弦波の周波数である、請求項1に記載の方法。前記測定電気信号集合は、前記導電率分布を決定するのに用いられる少なくともλL個の測定データ要素を含んでいる、請求項8に記載の方法。前記導電率分布又は誘電率分布を算出する前記工程は、Nが前記導電率分布又は誘電率分布を定義するのに用いられる値σjの総数に等しい場合にj=1~Nについての導電率値又は誘電率値σjの集合について前記エネルギ関数をとして定義するステップであって、前記測定電気信号集合を導電率又は誘電率に関する電圧値の変化の関数として表現する前記順モデルは、として表現され、式中、Uiは、i=1,2,3,...,Mである場合に前記電極44において測定される前記電気信号集合において総数Mのデータ要素を含むベクトルである、定義するステップと、である場合に∂F/∂σj=0となるような関数Fσjを、前記エネルギ関数と前記順モデルに基づく関数との間の差の偏導関数がゼロに等しくなるように評価し、これにより前記構造の前記内部の部分を表わす導電率分布又は誘電率分布を決定するステップとを含んでいる、請求項1に記載の方法。各電極が被検物体を切断する平面に沿って位置するときに、これらのアルゴリズムは同じ平面に沿った導電率分布を与えることができる。一般的には、広範な数学的方法及び数値的手法を適用して、被検物体の導電率分布に似た導電率分布を決定する。有用な画像分解能を生成するためには膨大な計算が必要とされる。EIT方式の画像再構成を商用の相対的に大きいアレイに導入するために、相対的に短時間で相対的に少ない計算を行なって満足な画像形成を提供し得る手法を見出すことが望ましい。本発明の一実施形態によれば、方法が、様々な導電率を有する多数の領域を含む構造の内部の画像を形成する。一実施形態では、各々が構造の内部に回転電場を発生するような一連の電気信号集合が一連の電極に相次いで印加される。測定電気信号集合が電極から得られ、これらの集合を用いて、電極が周囲に配置されている領域の構造の内部を表わす導電率分布を算出する。導電率分布を算出する工程は、導電率分布に基づく費用関数例えばエネルギ関数、及び測定電気信号集合を導電率に関する電圧値の変化の関数として表現する順モデルを定義することと、このエネルギ関数と順モデルに基づく関数との間の差の偏導関数がゼロに等しくなることを要求することとを含んでいる。後述する式6~式12を参照されたい。以下、1よりも多い物質がパイプ又はコンジットの内部を流れているような多相流において成分比を決定するために撮像が用いられ得るような応用について本発明を記載する。本発明の各実施形態は、パイプを流れる石油、水、及び炭化水素ガスの分率及び流量を推定するために、製油産業に用いられるもののような多物質系の特性を測定するシステム及び方法を提供するように作用する。本発明を石油/ガス/水測定での利用に関して記載するが、本発明はかかる応用に限定されない。寧ろ、本発明は、例えば医用撮像法、癌診断法及び水処理法のような広範な工業、医療及び化学工程に適用可能である。パイプを流れる物質は、固体、液体及び気体といった物理的状態の多数の状態にあってよい。多相流工程の正確な特性決定は、石油工業、医薬品工業、食品工業、及び化学工業に用いられるシステム及び処理設備の設計を改善し、動作効率を高めることを可能にする。多相工程の性能を予測するのに用いられる関連する流れ特性としては、例えば空間的相分布空間的容積相分率、流れ領域、界面面積、並びに相同士の間又は物質同士の間の絶対速度及び相対速度等が挙げられる。多相流における各物質の非一様な空間的分布を決定すると、さらに望ましい結果を得るために工程を管理し制御することが可能になる。例えば、石油製品を運ぶパイプラインでは、液体成分に対する気体成分の比を最小にすることができる。石炭スラリー輸送では、石炭に対する水の容積比を、単位容積当たり達成可能な最大のエネルギ量を保証するように最適化することができる。製造工程を最適化するのに適用可能なもう一つの例として、化学反応又は転化を起こしている物質の非一様な分布を監視して減少させることが重要である。かかる工程は、例えば再循環流が空間的に非一様な反応帯又は濃度を生成することにより、物質同士の間に減少した界面面積を有する場合がある。このように、容積相分率及び関連する成分の実時間知見が、多相流の時宜を得た実効的な制御を可能にする。しかしながら、かかる特性決定を生成する計算要件は、不適当な複雑さを加え、監視されている工程の時宜を得た制御を具現化するのに必要とされる応答時間を遅延させる場合がある。この理由から、先行システムの多くは、近似に頼って物質比を決定して、多相系での各物質の分離を特性決定していた。従来、EITを用いた画像再構成は、電極の各対に入力信号を相次いで与えて、1又は複数の他の電極において信号を感知することに基づいている。例えば、電流を一度に一対ずつ電極対の間に流して、残りの電極において電圧を測定することができる。代替的には、電圧を一度に一対ずつ電極対の間に印加して、他の電極の1又は複数において電流を測定することができる。実施形態の各例は、電流又は電圧を電極の全てに同時に印加することにより、相対的に多量の電気データを迅速に生成するように構成される。印加される電流又は電圧は、予め画定されている位相シフト又は周波数シフトを有していてもよいし、互いに対して電気的に同相にあって振幅変化を有していてもよい。従来、電極に印加される多数の電流又は電圧パターンに基づいて再構成を実行するためには相対的に長い処理時間が必要とされていた。一方、増大した量の電気データの生成によって、さらに高い信号対雑音比及びさらに高い分解能を実現することが可能である。このことは、被検物体が相対的な誘電率又は導電率のばらつきが小さい成分を含んでいるときには、必要な分解能を与えるために相対的に多数の測定が必要とされるため、特に適当である。図1に示す単純化した概略図を参照して、本発明の実施形態の一例を配管系14に接続された多数の油井12iを有する典型的な製油施設10について説明する。配管系14は、各々の油井12iから流量Fiを受け取って総流量Fを出力するように結合されている製油分岐管16を含んでいる。各々の油井からの流量Fiは多相流量計MPFM18を通過した後に分岐管16に入り、例えば液相及び気相の全体比を調整するように弁系統によって制御され得る。各々の多相流量計は、油井に極く近い処理前の油井流の測定を可能にし、このようにして油井の性能の連続監視を提供し得るシステムの一部である。各々の多相流量計システムによって処理される情報を用いて、貯油槽のさらに十分な管理及び流れ制御を行なうことができる。油井12iから汲み上げられる流体は製油分岐管16を介して製油分離装置20へ送られる。試験分離装置図示されていないを多相流量計システムに組み入れてもよい。多相流量計システムの試験分離装置に勝る一つの利点は、流れの組成を特性決定するのに必要とされる時間が短縮されることである。製油分離装置20は、油井から汲み上げられた石油、ガス、及び水を分離する。製油分離装置20はさらに、1又は複数の測定装置を含み得る。これらの測定装置は、例えば油井から抽出された水の量又は速度を測定する水量計、及び油井から抽出された石油の量を測定するエマルション計を含み得る。油井性能を監視するために坑口圧力センサ及び温度計を含むさらに他の測定装置を含めてもよい。次に、図2の単純化した概略図を参照すると、電気インピーダンス断層写真法EITに基づくMPFM18を含み、また本発明の実施形態の一例によれば導電率再構成に基づく画像再構成システム40が示されている。但し、原則として、EIT画像再構成は導電率及び/又は誘電率permittivity。admittivityとも呼ぶに基づくものであってよいことを理解されたい。システム40によれば、観察されている物質の導電率及び/又は誘電率は、電流信号集合又は電圧信号集合が印加されているときの被検容積の周辺の近くの電極から採取される境界電気測定に基づいて推測される。実施形態の各例によれば、各々の信号集合について全ての電極から測定を相次いで得ることにより、多量のデータが取得される。すなわち、印加電流及び/又は印加電圧の多数の異なる構成又はパターンについて測定を繰り返す。図2のMPFMシステム40のMPFM18は、任意の数L個の電極44と、この応用では全体的に円筒形を有するコンジット46の周りに分配された配線とを含む電極アレイ42を含んでいる。明示的に図示していないが、コンジットは、多相成分がコンジットを通過するように油井12iの一つと分岐管16との間に一列に配置されている。本発明を説明する目的のために図に示す電極の数Lは12であるが、12個よりも少ない又は多い電極をコンジットの周りに配置して測定を行なってよく、この数は部分的には、パイプの寸法及び所望の画像分解能に依存する。一実施形態では、電極の各々がコンジットの内壁に沿って装着されてもよく、外部信号発生源からの電気的接続が外壁を通過して電極と接続される。各々の電極には、電極と、コンジット46を流れる媒体との間の接触抵抗を低減するための皮膜が施工され得る。電極44のアレイは、円筒形コンジット46の中心対称軸に直交する平面Pにおいて構成される。アレイ42の電極は各々、電子回路48への配線を通じて接続され、電子回路48は、例えば各々の電極毎に、電流又は電圧発生源、D/A変換器、A/D変換器、差動増幅器、及び1又は複数のフィルタを含んでいる。回路48はさらに、入力信号集合を各電極に与えるためにコンピュータ50に結合されているアナログマルチプレクサ、クロック回路、及び/又はディジタル入出力ユニットを含んでいる。回路48の特定的な構成は当技術分野で公知であり、本発明の説明は回路設計の詳細な説明を必要とする訳ではないのでここではこれ以上立ち入らない。コンピュータ50は、画像再構成に用いられる形式の従来のディジタル信号プロセッサを含むハイエンドの高速プロセッサと、本発明に従って形成される画像を視認するのに適当な表示器52とを含み得る。フィールドプログラマブルゲートアレイFPGA及びコンプレクスプログラマブルロジックデバイスCPLDのような他の処理回路をコンピュータ50に含めてもよいし、システム40の他の部分に含めてもよい。各々の電気信号集合が電極44に印加されたときに、コンジットの内部に回転電場が発生される。以下の議論は、電圧信号が電極に印加されて、画像再構成の必要データを生成するためにこれらの電圧信号に応答して電流信号が取得されることを仮定している。電圧発生源を用いて各々のパターン集合を電極44に印加し、電極の各々において対応する電流信号集合が測定される。図3は、L個の電極を有する印加電圧型EITシステム60の単純化した概略図である。各々の電極44が、印加電圧を発生する電圧発生源と、電流を測定する電流計と、印加電圧を直接測定する電圧計とを含む回路62に接続されている。スイッチ開閉網63が、単一の較正回路64を電圧発生源/電流計/電圧計回路の任意のものに接続することを可能にして、システム全体を単一の参照発生源に対して較正することを可能にしている。ディジタル制御器図示されていないが、上述の各構成要素にインタフェイスを介して接続されており、較正、信号発生及び測定データの取得を含めた切り換え作用を動作させる。測定データが本発明に従ってコンピュータ50によって処理されて、電極が配置されている平面Pでの導電率分布を決定する。導電率及び誘電率分布はさらに処理されて、液体相、例えば石油及び水における様々な相の分布又は様々な物質の分布を含めた多相流の物性を決定する。導電率分布の解析によって、多物質系の個々の成分の流れ領域、相分率及び速度を決定することができる。流れ領域としては、例えば気泡流、チャーン流、スラグ流、又は環状流等がある。電圧発生源の実施形態では、各電極の間を流れる電流は、全ての電極を横断して印加される相対的な電圧と、各電極が周囲に構成されている平面Pに沿った導電率及び誘電率分布との関数である。例えば、物質は石油のみである場合もあれば、油水液体混合物及び分離した気相である場合もある。物質成分の分布に依存して、各電極の間の導電率及び誘電率は変化し、結果として、各電極の間の電流のレベルも変化する。すると、1又は複数の電圧集合の印加及び対応する測定電流データによって、平面Pに沿って変化する導電率に対応する導電率の値の行列を算出することができる。これらの導電率の値は平面Pを通る物質の導電率及び誘電率に対応しているので、導電率を解析することにより物質組成又は相流分布を決定することができる。各々の電圧パターン集合は時間変化型関数であって、この関数はコンジット44の内部に回転電場を発生し、すなわち電場は印加信号の周波数又はの整数倍の関数である角速度で方向を変化させる。定数の周波数の場合には、電場が回転する角速度も定数になる。L個の電極に電圧パターン集合を相次いで印加することにより回転電場が発生されるときに、各電極を横断して印加される電圧パターンは次式として表わされ得る。式中、Vkλは第λの励起パターン集合として第kの電極に印加される電圧であり、Vハット^記号は参照正弦波のピーク電圧であり、tは時間であり、ωは正弦波の周波数である。式2及び式3は、式1の正符号及び負符号についての二つの代替形態の一方に各々対応する式1の代替形態であり、正符号を選択したものが順電圧パターンを与え、負符号を選択したものが反転電圧パターンを与える。式2による信号集合が、コンジット46の内部での電場の時計回り回転の順電圧パターンを発生し、式3による信号集合が、コンジットの内部での電場の反時計回り回転の反転電圧パターンを発生する。励起パターン集合は、式2による順電圧パターン集合又は式3による反転電圧パターン集合を含み得る。かかる実施形態では、λは1~L/2、又はこれ以上にわたり、結果として合計でλ個の順電圧パターン集合又は反転電圧パターン集合の印加が各電極44に適用される。もう一つの実施形態では、各電極44に印加される信号集合が順パターン及び反転パターンの両方を含むことができ、結果として多物質系についてのさらに完全な情報が得られる。順電圧パターン及び反転電圧パターンは合わせて、システムの電気的要素の実成分及び虚成分についての情報を与える。図4は、L=8すなわち8個の電極44iをi=1~8について示す電極44の概略図80である。L=8個の電極の図は例示的なものである。各電極44iは、円筒形コンジット46の平面Pによって切断される円の周囲の各位置に一様間隔で隔設されている。図4はさらに、正弦波形に配置される波形値Viの集合をi=1~8について示している。所与のiの値について、各々の図示の値Viが対応する電極44iに入力される信号に対応する。すなわち、各々の波形値Viは、Vハット=10Vである場合にλ=1を有する印加パターン集合について式2又は式3に従って発生される8個の時間依存型電圧入力信号の異なる一つ一つに同時に生ずる時間的な値である。8個の電極44iの系列において隣り合った電極に印加される電圧振幅Vkλの間の差は、k=i及びλ=1である場合に2πλk-1/Lずつ与えられる位相シフトに基づく。このように、波形値V1~V8は、所与の時刻に各電極に同時に印加される位相シフトした電圧信号である。この例では、コンジット46は、油水エマルション86、及びエマルションから分離したガス88を含んでおり、ガスはコンジット46の中心領域を流れ、エマルション86がガス流を包囲している。例えばλがそれぞれ2及び/又は3に等しい場合の式2又は式3による第二高調波及び/又は第三高調波の電圧パターンを含めた他の電圧パターン集合を各電極に印加してもよい。図3のシステム60によれば、各々の電極の専用となっている回路62は全て共通のクロック信号の下で動作し、各々が隣接する電極について発生される波形の位相に対して別個の位相を有する波形を発生する。この観点で、波形値Viの集合がプロットされる横軸は、円筒形コンジット46の周りでの平面Pにおける角度位置に対応し、所与の時刻例えば時刻t=0においては、正弦波のV1~V8とラベルされた8個の点の各々が、各々の電極における印加信号の振幅に対応する。このように、2πLλ=1の場合の位相シフトを隣り合った電極同士の間に生じさせながら、多数の位相シフトした正弦波形が、等間隔で隔設された電極の異なるものに同時に印加され、波形に示される値V1~V8は時刻tにおいて8個の電極に印加される瞬間電圧を指す。集合的に、8個の入力信号は、平面Pにおいて回転する時間変化型電場パターンを発生する。一実施形態では、導電率分布のみを決定すればよいときには、相分布についての完全な情報を取得するために各電極に印加される電圧パターン集合の総数はλ=L/2まで少なくすることができる。L/2個のパターンは、時計回り回転電場を発生する順パターン又は反時計回り回転電場を発生する反転パターンの何れかとなり得る。本発明の一実施形態によるデータ取得工程は、各電極44にλ個の電圧パターン集合を相次いで印加する。他の実施形態では、印加される電気信号集合は、1又は複数の電流値集合及び1又は複数の電圧値集合を含み得る。理解し易くする目的のために、第一の電気信号集合は電圧値集合であるものと仮定する。各々の印加電圧パターン集合毎に、電気信号パターン集合が例えば図3のシステム60に従って測定される。各々の印加電圧パターン集合毎に取得されるデータ要素の例示的な数は、電極の数と、所与のλの値について印加される信号のサイクル当たり採取される測定の数との積である。A/D変換器を含むシステム40の回路48によって、各々の電圧励起パターンについて、電流データの時系列がL個の電極の各々について記録される。この時系列データは、全L個の電極における電流の振幅及び位相を抽出するように処理される。各々の電極における時系列データは、当該電極における電流について一つの振幅及び一つの位相値を与える。故に、L個の電極についてはL個の測定が存在する。この工程を全λ個の電圧パターン集合について繰り返すと、λL個の時系列測定データ要素が得られる。λL個の時系列測定データ要素に基づいて、コンジット46における物質の分布を表わす導電率行列が決定される。画像再構成は、導電率及び/又は誘電率行列に基づいて行なわれる。物質分布を示す行列値の決定例についてのさらに詳細な議論については、参照により本出願に援用される米国特許出願第12/652,116号を参照されたい。本発明によれば、導電率行列又は誘電率行列の決定は、順モデルの不整合を最小にすることにより分布を求める従来のアプローチとは異なる。代わりに、導電率分布のエネルギが順モデル不整合の束縛条件を受ける目的関数となるような最大値問題を解く。各々の電極において電圧がUelであり、順モデルがUel=Jσσ式4と定義される場合の例を以下に掲げる。ヤコビアンJは導電率σに関するUelの偏導関数であるが、σは代替的には誘電率を表わしていてもよい。行列形態では、この順モデルは次式として表わされる。式中、Uiは総数M個のデータ要素、例えば測定電流データを含むベクトルであり、i=1,2,3,...,Mである。導電率行列σjは、コンジット46の内部の平面Pに沿った領域について決定されるべき導電率分布を集合的に定義するj個の要素を含んでいる。目的関数は、電場に関係する多様な可能なエネルギ関数から選択され得る。ここでは、一例としてを選択する。すなわち、順モデル不整合を束縛条件として導電率分布のエネルギが最大化される。明確に述べると、最大化はラグランジュの未定乗数法を用いて達成されることができ、λiがラグランジュの未定乗数項を示す。従って、Fσが∂F/∂σj=0式8の条件を受ける。すると、ここから及びが得られ、従って順モデル表現は、となり式中、Uiは実測された電圧データである、この式を解くと、各々のλk毎にλp1-λp=0についてガウス=ザイデル繰り返し法を介してラグランジュの未定乗数λkの値を求めることができる。λp1=λpDTDαRTR1DTgexp-gcomp式13式中、Rは正則化パラメータであり、Dの各々の要素はと定義され、αは正則化パラメータであり、gexp=giであり、gcomp=2Ucompであり、Ucompは順モデルUel=Jσσを用いて算出される電圧データであり、項DT及びRTはそれぞれD及びRの転置行列である。構造の内部の画像を表示する方法について記載した。この画像は、人体の器官又は他の部分の医用画像であってよいが、さらに一般的には、コンジットを流れる多相流体例えば物質は石油及び炭化水素ガスを含むのような異なる電気的特性によって特徴付けられる容積、又は容積を流れる電気信号の伝達に影響を与える様々な導電率を有する多数の領域を含む容積の内部の特徴の画像であってよい。決定された画像情報を用いて、コンジットを流れる液体物質及び気体物質の相対比を推定することができる。本発明の各実施形態の特徴は、平面に沿って配置される各電極に入力される一連の電気信号集合であって、各々の信号集合が電極の各々に入力される電気信号パターンを含んでいるような一連の電気信号集合の提供である。様々な導電率を有する2以上の領域を含む構造によって、本発明は、電気インピーダンス断層写真法による画像形成に対するさらに計算効率のよいアプローチによる態様で構造の内部を表わす画像を決定する。比較のために、ニュートンワンステップエラーリコンストラクタNewtonOneStepErrorReconstructor手法を用いると、繰り返し式逆ソルバは導電率分布のN個の要素の全てについての集約的計算による解決を必要とするが、上述の方法はM個のラグランジュの未定乗数λkについて解く。すなわち、導電率行列における要素の数Nは取得される測定の数Mよりも遥かに大きく、順モデルと測定データとの間の差に基づいて解を繰り返し求める従来の再構成アプローチは、繰り返し式手法によってN個の解を求めることを要求する。本発明の開示された実施形態では、一連の印加電気信号集合が各々電圧入力パターンを含んでおり、測定電気信号集合が各々測定電流レベルを含んでいる。他の実施形態では、一連の印加電気信号集合が各々電流入力パターンを含み、測定電気信号集合が各々測定電圧レベルを含む。さらに他の実施形態では、一連の印加電気信号集合の幾つかが電圧入力パターンを含み、一連の印加電気信号集合の残部が電流入力パターンを含み、測定電気信号集合の一部が測定電流レベルを含み、測定電気信号集合の残部が測定電圧レベルを含む。本発明の方法を全体的に円筒形である容積を画定するコンジットについて説明したが、本発明は、様々な寸法及び形状の本体空洞及び本体実質を含む様々な形状の構造に適用可能である。また、複数の電極から成る単一のアレイについて説明したが、信号集合は、多数の隔設された電極アレイに同時に印加されても相次いで印加されてもよく例えば被検構造を切断する平面に沿って、各電極及び構造が被検容積の様々な部分に関連する画像情報を生成するために互いに対して変位され得るような構成も思量される。本書では発明の幾つかの特徴のみを図示して説明したが、当業者には多くの改変及び変形が想到されよう。従って、特許請求の範囲は、発明の要旨に含まれるような全ての改変及び変形を網羅するものと理解されたい。"