(57)【要約】有限要素モデルを解くための系および方法である。ここで、解を支配するマトリクスはマトリクスの重み付け係数を調節することにより変更され、マトリクスの対角線上にある要素は非ネガティブであり、対角線上にない要素は非ポジティブである。一つの実施形態では、系は、各要素をまたぐ流体フローの方向に基づいて重み付けされている区切りに対する各ノードの寄与とともに、有限要素メッシュ上に区分される。各要素の他方のノードから上流のノードは、得られたマトリクスが実質的に対角線とされるようにより重く重み付けされる。このマトリクスは伝統的な方法を使用して解かれる。【特許請求の範囲】
【請求項１】
多相流体フローが存在する系に対応する有限要素モデルを解く方法であって、
前記モデルに対応する有限要素マトリクスを生成するステップであって、ここで、前記マトリクスは複数の係数を含む、ステップと、
前記係数を調節して、対角線上の要素は非ネガティブであり対角線上にない要素は非ポジティブであるマトリクスを獲得するステップと、
有限要素法を使用して前記マトリクスを用いる前記モデルについての解を生成するステップと
を包含する、方法。
【請求項２】
前記係数を調節するステップは、前記要素をまたぐ流体フローの方向に従って各要素のノードに重み付けをするステップを包含する、請求項１に記載の方法。
【請求項３】
前記要素をまたぐ流体フローの方向に従って各要素のノードに重み付けをするステップは、前記要素をまたぐ流体フローの方向を決定するステップと、前記ノードが前記要素の他方のノードから上流にある場合には各ノードをより重く重み付けし、前記ノードが前記要素の他方のノードから下流にある場合には各ノードをあまり重くないように重み付けするステップとを包含する、請求項２に記載の方法。
【請求項４】
前記要素の大きい方の部分が前記要素の他方のノードからより前記ノードから下流にある場合には各ノードをより重く重み付けし、前記要素の小さい方の部分が前記要素の他方のノードからより前記ノードから下流にある場合には各ノードあまり重くないように重み付けする、請求項３に記載の方法。
【請求項５】
前記有限要素マトリクスは少なくとも２つの流体相が存在する系に対応する、請求項１に記載の方法。
【請求項６】
前記有限要素マトリクスは３つ以上の流体相が存在する系に対応する、請求項１に記載の方法。
【請求項７】
前記有限要素マトリクスは４次元有限要素モデルに対応する、請求項１に記載の方法。
【請求項８】
前記系は油層に対応する、請求項１に記載の方法。
【請求項９】
前記マトリクスは任意の時間に物理的に非現実的でなくはない解を生成するように構成される、請求項１に記載の方法。
【請求項１０】
前記マトリクスは振動しない解を生成するように構成される、請求項１に記載の方法。
【請求項１１】
前記系のモデルを分離して有限要素メッシュを生成するステップと、前記有限要素メッシュに基づいて前記マトリクスを生成するステップとをさらに包含する、請求項１に記載の方法。
【請求項１２】
有限要素モデルを解く場合に向上した精度を得るための方法であって、
多相流体フローが存在する系のモデルを分離するステップと、
前記モデルに対応する有限要素マトリクスを生成するステップであって、ここで、前記マトリクスは、単調的でありかつ線形性を保持する解を生成するように構成される、ステップと、
前記マトリクスを使用して前記モデルの解を生成するステップと
を包含する、方法。
【請求項１３】
複数の指示を含むコンピュータ読み取り可能な媒体であって、ここで、前記指示は、多相流体フローが存在する系に対応する有限要素モデルを解く方法をコンピュータに実行させるように構成され、前記方法は、
前記モデルに対応する有限要素マトリクスを生成するステップであって、ここで、前記マトリクスは複数の係数を含む、ステップと、
前記係数を調節して、対角線上の要素は非ネガティブであり対角線上にない要素は非ポジティブであるマトリクスを獲得するステップと、
有限要素法を使用して前記マトリクスを用いる前記モデルについての解を生成するステップと
を包含する、コンピュータ読み取り可能な媒体。
【請求項１４】
前記係数を調節するステップは、前記要素をまたぐ流体フローの方向に従って各要素のノードに重み付けをするステップを包含する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項１５】
前記要素をまたぐ流体フローの方向に従って各要素のノードに重み付けをするステップは、前記要素をまたぐ流体フローの方向を決定するステップと、前記ノードが前記要素の他方のノードから上流にある場合には各ノードをより重く重み付けし、前記ノードが前記要素の他方のノードから下流にある場合には各ノードをあまり重くないように重み付けするステップとを包含する、請求項１４に記載のコンピュータ読み取り可能な媒体。
【請求項１６】
前記要素の大きい方の部分が前記要素の他方のノードからより前記ノードから下流にある場合には各ノードをより重く重み付けし、前記要素の小さい方の部分が前記要素の他方のノードからより前記ノードから下流にある場合には各ノードあまり重くないように重み付けする、請求項１５に記載のコンピュータ読み取り可能な媒体。
【請求項１７】
前記有限要素マトリクスは少なくとも２つの流体相が存在する系に対応する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項１８】
前記有限要素マトリクスは３つ以上の流体相が存在する系に対応する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項１９】
前記有限要素マトリクスは４次元有限要素モデルに対応する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項２０】
前記系は油層に対応する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項２１】
前記マトリクスは任意の時間に物理的に非現実的でなくはない解を生成するように構成される、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項２２】
前記マトリクスは振動しない解を生成するように構成される、請求項１３に記載のコンピュータ読み取り可能な媒体。
【請求項２３】
前記方法は、前記系のモデルを分離して有限要素メッシュを生成するステップと、前記有限要素メッシュに基づいて前記マトリクスを生成するステップとをさらに包含する、請求項１３に記載のコンピュータ読み取り可能な媒体。
【発明の詳細な説明】
【０００１】
（関連出願）
本出願は、米国特許法１１９（ｅ）のもとに、米国特許出願第６０／２１５，６９７号の優先権を主張する。この出願は、Ｓｔｅｐｈｅｎ　Ｒ．ｋｅｎｎｏｎ、Ｋｏｋ　Ｔｈｙｅ　Ｌｉｍ、Ｓｃｏｔｔ　Ａ．Ｃａｎａａｎ、Ｓｔｅｖｅｎ　Ｂ．Ｗａｒｄ、Ｓｔｕａｒｔ　Ｗ．Ｐｏｎｄ，Ｊｒ．、およびＥｄｗａｒｄ　Ｊ．Ｂａｒｒａｇｙらによる「Ｍｅｔｈｏｄ　ａｎｄ　Ｓｙｓｔｅｍ　ｆｏｒ　Ｏｉｌ　Ｒｅｓｅｒｖｏｉｒ　Ｓｉｍｕｌａｔｉｏｎ　ａｎｄ　Ｍｏｄｅｌｉｎｇ」という名称の、２０００年６月２９日に出願された出願であり、あたかもその全文が説明されたかのように、本明細書中で参考として援用される。
【０００２】
（発明の分野）
本発明は、概して、有限要素解析を用いて、物理的システムをモデリングする方法、より具体的には、モデルの解を制御するマトリクスが、単調性および線形性の保護を維持するように構成される、有限要素モデルを解く方法に関する。
【０００３】
（発明の背景）
物理的なシステムは、ある一定の条件下でのふるまいをシミュレートするように、機械的にモデリングされ得る。これらのシステムをモデリングする手段には様々なものがあり、非常に単純なものから、極度に複雑なものにわたる。物理的システムをモデリングする、より高度な手段のうちの１つは、有限解析の使用によるものである。名前が示すように、有限要素解析には、数学的モデルにおける物理的システムの個々の有限の要素の表現、および所定の組の境界条件がある場合におけるこのモデルの解が含まれる。物理的システムをモデリングする他の類似する手段は、差分解析の使用によるものである。差分解析には、モデリングされた空間内の個々の点のモデリング、およびこれらの点の間の差の計算が含まれる。差分解析は、しばしば、流体の動的なふるまいをシミュレートするために用いられる。
【０００４】
合理化された、上流（ｕｐｗｉｎｄｉｎｇ）の方法を用いる、従来の差分技術および有限要素技術は、典型的には、油層での油の生成をシミュレートするために用いられる。これらの技術の各々は、利点を有し、問題点をも有する。概して、差分技術は、物理的に現実的な結果を生成するが、あまり正確ではない。一方、有限要素技術は、より正確であるが、物理的に現実的でない結果を生成する。
【０００５】
それぞれの問題点の結果として、従来用いられていた差分技術、および有限要素技術は、通常、かなりの量のコンピュータリソースを必要としていた。差分技術の場合、妥当な精度が達成され得るが、有限要素モデルにおいて必要な要素より、かなり多いノードを必要とする。これによって、正確な解を計算するために必要なメモリおよびＣＰＵの時間の量が増大する。有限要素技術の場合、より少ないノードで、同等の精度が達成され得るが、解が現実的でないことがある。例えば、濃度が、１００パーセントより多くなったり、透過率が負の値になったりする。従って、解を収束させないか、または非常に遅く収束させる問題点があり得る。従って、これらの技術は、より信頼性が低く、許容できる精度が達成される前に、多数のインタラクションを必要とし得る。
【０００６】
これらの標準的な技術の各々が欠点を有し、欠点が、許容できる解を生成するために必要なコンピュータリソースの量を増大させるので、油層のようなシステムをモデリングして、このようなシステムについて、正確、かつ現実的な解を確実に生成する方法を提供することが望まれる。
【０００７】
（発明の要旨）
上記の問題点のうちの１つ以上は、本発明の様々な実施形態によって解決され得る。概して、本発明は、解を支配するマトリクスが改良されて、単調性および線形性の保持の目的を満たす、有限要素モデルを解くシステムおよび方法を含む。これは、マトリクスの重み付け係数を調節して、マトリクスの対角線上にある要素が負にならないようにし、対角線から離れている要素を正にならないようにすることによって、達成される。これによって、解が、従来の差分技術より正確になり、従来の有限要素技術より現実的になる。
【０００８】
本発明の技術は、従来技術と比較して、様々な利点を提供し得る。例えば、本発明の技術がより正確なので、許容できる結果を生成するために必要なインタラクションの数は低減され得、計算するリソースが節約される。また、本発明の技術は、解が現実的な範囲に制約され、従来の有限要素技術の解と同じ様態で振動しないので、解がより確実に収束されるようにする。従って、本発明の技術は、概して、従来技術と比較すると、より高速であり、より正確であり、リソースを使う量がより少なく、より信頼できる。
【０００９】
本発明の一実施形態は、有限要素技術を用いてモデリングされたシステムについて、解を生成する方法を含む。このシステムは、有限要素メッシュに対して分離され、各ノードの分離への寄与が、各要素にわたる流体の流れの方向に基づいて重み付けされる。寄与は、それぞれの要素の他のノードから上流にあるノードを有利にするように重み付けされる。このようにして寄与の重み付けを行うことによって、マトリクスが実質的に対面マトリクスとなる。すなわち、マトリクスの対角線にある要素が強調され、対角線から離れている要素に重きがおかれなくなる。マトリクスがこのように改良された後、従来の技術を用いて解かれる。
【００１０】
本発明の代替的な実施形態は、油層における複数相の流体の流れをモデリングする方法を含む。この実施形態において、油層に対応する四次元の有限要素メッシュが構成される。メッシュのノードの間の有限要素オペレータの残留物を分配し、その後、従来の有限要素技術を用いて組み立て、解くことによって、解がメッシュについて生成される。この実施形態において、有限要素オペレータの残留物を分配する工程は、各要素について、有限要素オペレータの値を各ノードについて計算する工程と、これらの値を、正の値の組および負の値の組に分割する工程と、対角線の要素を強調し、対角線から離れた要素に重きをおかないように、残留物を分配する工程と、マトリクスについてタンジェントオペレータを計算する工程とを含む。
【００１１】
本発明の別の代替的な実施形態は、コンピュータ可読媒体に組み込まれたソフトウェア製品を含む。このソフトウェア製品は、コンピュータまたは任意の他のデータプロセッサに、本開示による方法の工程を実施させるように構成された複数の命令を含む。このコンピュータ可読媒体は、ＣＤ−ＲＯＭ、ＤＶＤ、フロッピー（Ｒ）ディスク、コンピュータメモリ、磁気テープ等（但し、これらに限定されない）を含む多数の異なる媒体のうちのいくつかを含み得る。本発明のさらに別の代替的な実施形態は、本開示による方法の工程を実施するように構成されたコンピュータシステムを含む。このコンピュータシステムは、ハードウェア、ファームウェア、またはソフトウェアに収納された命令により、本開示による方法を実施するようにプログラムされ得る。
【００１２】
多数のさらなる代替的な実施形態も存在可能である。
【００１３】
（好適な実施形態の詳細な説明）
本発明の他の目的および利点は、下記の詳細な説明を読み、かつ添付の図面を参照することにより明白になり得る。
【００１４】
本発明は、種々の修正が行われ、かつ代替的な形態を有することが仮定される一方、本発明の具体的な実施形態を図面における例示、およびそれに伴う詳細な説明により示す。しかしながら、図面および詳細な説明は、本発明を説明される特定の実施形態に限定することを意図するものでないことが理解される。本開示は、反対に、特許請求の範囲により規定される本発明の範囲内に相当する全ての修正、等価物、および代替物を含むことを意図するものである。
【００１５】
本発明の好適な実施形態を以下に示す。下記に示される好適な実施形態および任意の他の実施形態は例示であり、本発明を制限するのではなく、説明することを意図するものであることに留意されたい。
【００１６】
概して、本発明は有限要素モデルを求めるための方法を含み、その解は、単調であること、および直線性の保持を満足させることの両方を保証する。これにより、その解が正確であり、かつ現実的であることの両方が確実となる。解の単調性および直線性は、対角線上にあるマトリクスの要素が負ではなく、対角線上にない要素が正ではないように、有限要素モデルの解についてのマトリクスを生成することにより確実となる。そして、このマトリクスは、所望の正確で現実的な解を生成するために、従来の有限要素技術を用いて求められ得る。
【００１７】
本発明の方法は、伝統的な有限要素技術、および伝統的な有限差分技術に対して多数の利点を提供する。上述したとおり、本発明の方法は、非現実的で振動する解という不利な点がない、伝統的な有限要素技術の精度を提供する。同様に、本発明の方法は、伝統的な有限差分技術の安定性を提供する一方で、はるかに高い精度を提供する。
【００１８】
本発明の方法により提供される精度および現実性は、コンピュータリソースの利用におけるさらなる利点へと繋がる。例えば、本発明の方法は、物理的システムのための解をより速く生成することを許容し得る。これは、いくつかの事項の結果である。まず、解は物理的に現実的であるため振動せず、その結果、より迅速に収束する傾向がある。さらに、本発明の方法を用いる解の生成は、通常、伝統的な有限要素技術ほど反復を必要とせず、かつ速いため、解を生成するために必要とされるＣＰＵ時間が短い。本発明の方法はまた、伝統的な技術ほどメモリを使用し得ない。本発明の方法は、一般に、伝統的な有限差分技術よりも正確であるため、物理的システムをモデリングするほどの多くのノードを用いる必要がない。その結果、そのモデルを格納するために必要なメモリは少ない（そして、これらのノードに基づいて解を算出するためのＣＰＵ時間が短い）。本発明の方法はまた、伝統的な技術よりもさらに信頼でき得る。伝統的な有限要素技術を用いて生成された解の振動は、いくつかの場合において、解が収束することを妨げ得る。上記のとおり、本発明の方法は、現実的な（例えば、現実的な解と存在し得ない解との間を振動しない）解を提供するため、解は確実に収束する。
【００１９】
有限要素モデリングでは、解析される領域が要素と呼ばれる小領域に分けられる。この領域を小領域に分割するプロセスは、離散化（ｄｉｓｃｒｅｔｉｚａｔｉｏｎ）またはメッシュ生成と呼ばれ得る。この領域は、各要素にわたって定義された関数により表わされる。これにより、領域全体を表わすために必要とされる局所関数よりもはるかに単純である多数の局所関数が生成される。次の工程は、各要素に対する応答を解析することである。これは、その領域内の種々の要素の特性を定義するマトリクス、およびそのドメイン内の各要素に従う残留物（ｒｅｓｉｄｕａｌｓ）を定義するベクトルを構築することにより達成される。全ての要素マトリクスおよびベクトルが作成されると、それらは組み合わせられ、マトリクス式となる。境界条件を適用した後、マトリクス式は、未知のノードの応答を得るために求められ得る。要素内応答が、各要素にわたって定義された関数を用いて、ノードの値から補間され得る。
【００２０】
本発明の一実施形態についての以下の説明は、油層のモデリングに関する。この実施形態は、例示のために示されるものであり、本発明の方法は、各種の物理的システム、例えば、航空機の表面上の流れ、天気予報、内燃機関を流れる流れ、他の産業的な流体の流れに適用され得ることに留意する必要がある。潜在的起こり得る適用例のリストは、限定するものではなく、例示するものである。
【００２１】
好適な実施形態の細部を、以下に示す。しかし、まずいくつかの用語を定義することが有用であり得る。
【００２２】
ノードは、空間内の点である。有限要素モデリングにおいて、ノードは、モデリングされる要素の頂点を形成する。また、ノードは、モデリングされる空間内の要素間の境界を規定する、ノードおよび辺のメッシュの一部を形成する。
【００２３】
辺は、要素の頂点を形成する２つのノード間の線である。辺は、モデリングされる空間内の要素間の境界を規定するメッシュの一部を形成する。
【００２４】
単体は、ｎ次以上の空間における、ｎ＋１個の点によって決定される、ｎ次元の空間形状である。すなわち、単体は、所与の次元の空間を囲む必要がある、最小限の数の境界点を有する、幾何学的、空間的要素である。例えば、二次元において、単体は、その三角形によって境界が定められる内部領域を有する三角形を含む（図１参照）。２点では、その点（および、その点を相互に接続する線）によって境界が定められる領域が存在しないので、二次元空間における単体を形成するためには不十分である。４点は、二次元領域の境界を定めるためには充分であり得るが、二次元領域の境界を定め得る境界点の最小限の数ではない。三次元において、単体は、４つの頂点によって境界が定められる、四面体を含む（図１を参照）。四次元において、単体は、超四面体（ハイパーテットとも呼ばれる）を含む。
【００２５】
メッシュは、空間を満たす要素の集合である。これらの要素は、その空間内に存在するシステムを表す。各要素が、複数のノード、および／またはこれらのノード間の辺によって規定され得るので、メッシュは、ノード、および／またはノード間の辺の集合としても考慮され得る。この開示の様々な箇所で、「メッシュ」は、この用語が用いられる文脈によって、要素の集合またはノード／辺の集合のどちらをも意味するように用いられる。本明細書中、メッシュは、有限要素モデリング、または単にモデリングを意味するようにも用いられ得る。
【００２６】
単一の流体相があるモデリングシステムと複数の相があるモデリングシステムとの間には、根本的な違いがある。複数の流体相がある場合、対流があり、相が混合する可能性がある。（流体を駆動する圧力があり、流体間の面が移動すると仮定する。）従って、有限要素モデルは、対流オペレータまたは要素を含む必要がある。これらのオペレータは、１つの流体相しかモデリングされない場合、または、複数の流体相が完全に混合されるかまたは、乳化される場合は必要ない。
【００２７】
複数の流体相を有するシステムにおいて、システムの特徴は、非常に唐突に変化し得る。この状況を図２に示す。図２は、物理的システムにおける実際の値と、従来の有限要素技術および本発明による技術を用いて求められた解との間の差を示す図である。実際の（「現実世界の」）値１１が、第１の値ｖ１から、距離ｄ１に達した直後に、第２の値ｖ０まで変化する、本質的には階段関数である関数を形成する。この実施例において、値ｖ０およびｖ１は、物理的に実現可能な値の範囲の限界を規定する。
【００２８】
従来の有限要素技術が用いられてこのようなシステムがモデリングされる場合、解は、典型的には、変動する傾向がある。図２に示すように、これらの技術を用いて生成される解１２は、モデリングを試みている、実際の物理的な値の上と下とに変動する。この解の最も大きな問題点のうち２つは、明らかな不正確さ（解は、実際の値を厳密にたどらない）、および、いくつかの点において、解が物理的に実現可能な最大値より大きいか、または、物理的に実現可能な最小値より小さいことである。従って、従来の有限要素技術を用いて生成される解において起こる変動は、非現実的な結果を生成するか、または、全く、収束しなくなることもある。
【００２９】
しかし、この開示内容による方法は、解が現実的になるように制約されるので、これらの不利益を被らない。この制約は、解ベクトルを管理するマトリクス（すなわち、ヤコビのマトリクス）が単調性の性質、および線形性の保護を満たすことを確実にすることによって、達成される。より実際的な用語で記載すると、これは、マトリクスの対角線上の全ての係数が負でないことを確実にし、マトリクスの対角線から離れている全ての係数が正でないことを確実にすることを含む。マトリクスがこれらの特徴を有する場合、解は、全ての点において現実的であり続け、図２における曲線１３に示すように、物理的システムのふるまいをより正確にたどる。従って、モデリングは、正確な結果を生成し、確実に収束する可能性がより高い。
【００３０】
ヤコビのマトリクスは、ダーシーの法則に基づいて、多孔質媒質を流れる複数相の流れを記述し、システムにおける質量および運動量の保存を表す偏微分方程式から得られる。これらの式は、有限要素メッシュに対して分離され、式の係数はマトリクスを形成する。各係数は、式の勾配が正であるか、または負であるか（すなわち、流体の流れが正であるか、または負であるか）に依存して、正または負のいずれかになる。有限要素技術が微分方程式に適用される場合、代数方程式の大きな集合が生成される。この代数方程式の大きな集合は、非線形である。式の集合は、ニュートン法（式の第１の導関数を取ることを含む）を用いることによって、解かれ得る。第１の導関数は、ヤコビのマトリクスである。このマトリクスの係数は、対角線上に負でない項を有し、対角線から離れた場所に正でない項を有するマトリクスを得るように、重み付けされる。項の適切な重み付けは、マトリクスの単調性の性質、および線形性を確実にする。
【００３１】
ヤコビのマトリクスの係数は、実質的には、有限要素メッシュにおける要素の各ノードについての重み付け係数である。従来の有限要素技術において、要素のノードは、均一に重み付けされる。例えば、二次元要素において、３つのノードの各々は、１／３で重み付けされる。三次元要素が用いられる場合、４つのノードの各々は、１／４で重み付けされる。しかし、本発明の方法において、ノードは、モデリングされたシステムにおいて、流体の流れを反映させるように重み付けされる。これは、図３〜６を参照しながら示される。
【００３２】
図３は、１つの二次元要素（三角形）を示す。要素に関する流体の流れを矢印で示す。この矢印は、流体が図の上方右手部分から図の下方左手部分へと流れていることを示す。従来の有限要素法を用いた場合、ノード２１〜２３の各々は、等しく重み付けされている（すなわち、１／３、１／３、１／３）。本発明の技術を用いた場合には、これらのノードは等しく重み付けされない。流体が概してノード２２および２３からノード２１へと流れるため、ノード２２および２３のそれぞれは１／２で重み付けされ、一方ノード２１は０で重み付けされている。図４を参照して、流体の流れる方向が変化した場合、ノードの重み付けが異なることが分かる。図４において、流体は、概して、図の下方左手部分から図の上方右手部分へと流れている。従って、ノード２１は１で重み付けされ、一方ノード２２および２３は０で重み付けされている。
【００３３】
図５および６を参照して、ノードの重み付けは、対応する要素の形状によっても影響を受けることが分かる。図５において、流体の流れは、要素のエッジに対してほぼ垂直である。この状況は、図３の状況と類似している。この場合、要素の前方のエッジを終端させるノードの各々は１／２で重み付けされた。この場合、要素の残りのノードは、要素の前方のエッジにある中点からのダウンストリームであった。図５において、代わりに、残りのノードは第１のノードから別のノードへの直接的なダウンストリームである。この状況では、最後のノード３３からのアップストリームであるノード３２は２／３で重み付けされ、一方他のアップストリームであるノード３１は１／３で重み付けされる。図６には１つのアップストリームノード３３がある。ノード３１および３２はともに要素の後縁上にあるので、それらは共に０で重み付けされている。一方、アップストリームノード３３は１で重み付けされている。
【００３４】
これらの場合のそれぞれで、流体が流れる方向についてのノードの重み付けは、対角要素から外れた要素の影響を減らすかまたはなくしつつ、マトリクスの対角要素を効果的に強調している。従って、このことは、マトリクスが単調性および線形性といった特性を有することを確実にしている。
【００３５】
前述の記載は、有限要素法の分野の当業者が本発明を実行し得るに十分であるが、本技術の簡単な技術的説明を以下に述べる。
【００３６】
時間にわたってモデリングされる物理システムを考えると、有限要素メッシュがまず構成される（すなわち、システムは時空ドメインにおいて区切られる）。図７を参照すると、時空ドメインの二次元投影図が示される。
【００３７】
Ω＝Ωｓｐａｃｅ×［Ｏ，Ｔ］は、モデル化される時空ドメインである。ここで
【００３８】
【数１】
2004502999.tif 000002
であり、Ωを壊して（セグメント化して）、複数の重なることのないサブドメインΩｅにする。これらのサブドメインΩｅの各々および
【００３９】
【数２】
2004502999.tif 000003
は、四次元単体（ハイパー四面体）を含む。線形有限要素の基本は、各要素（Ωｅ）となり、試行関数およびテスト関数が各要素に規定される。
【００４０】
多孔質媒体内での多相流体の流れは以下の式
【００４１】
【数３】
2004502999.tif 000004
に基づく。
【００４２】
この多相の流れの式である基本的なＧａｌｅｒｋｉｎ有限要素区切りは、各要素に対して
【００４３】
【数４】
2004502999.tif 000005
のように与えられる。
【００４４】
本発明の技術を採用する基本的な方法は、
【００４５】
【数５】
2004502999.tif 000006
が規定される様態で規定される。従って、本方法は以下のように規定される：
【００４６】
【数６】
2004502999.tif 000007
図８を参照して、本発明の１実施形態による方法を図示するフローチャートを示す。上述したように、微分方程式が有限要素に対して区切りを与える場合には、これら微分方程式は、基本的には、微分方程式からノードの値間の代数的関係へと変換される。このことは、ノード間に対応する特性（例えば、飽和、圧力、浸透率、前立腺等）の線形変化があると仮定することで達成される。これは、すべての可能な解の空間を線形である解のみに制限する。この仮定を用いる場合、１つの関連性のセットは、各要素にわたって線形である係数の値についての微分方程式から導かれる。これが行われる場合、各要素は、特定量の区切りをノードの各々に付加する。各ノードにおいて合計された付加は、残差となる。非線形代数系を解くためには、残差＝０とすることが望ましい。
【００４７】
有限要素演算子（γｊｉ）の残差を分散させるために、まず、４四面体（ハイパー四面体）の各ノードについて、Ｋｊを計算する。得られるＫｊは正の値（Ｋｊ＋）と負の値（Ｋｊ−）とに分けられる。次に、残差γｊがノードに分散され、正弦演算子（Ｔｊｋ）が計算される。次いで、システムが構築されて、通常の有限要素法を用いた場合と同様に解を求める（有限要素システムを構築し解を求めるプロセスが周知であるので、本明細書中ではこれ以上詳述しない）。
【００４８】
本発明によって提供され得る利益および利点は、特定の実施形態について上述してきた。これらの利益および利点、ならびにこれら利益および利点を生じさせるかまたはより明確にし得る任意の要素または制限は、特許請求の範囲に記載の任意または全ての請求項にとって、重要、必要、または必須の特徴であると解釈されるべきではない。本明細書中で用いられる用語「含む（ｃｏｍｐｒｉｓｅｓ）」「含む（ｃｏｍｐｒｉｓｉｎｇ）」またはこれらの任意の他の変形は、文言に従う要素または制限を含むが限定しないものとして解釈されるように意図される。従って、一連の要素を含むプロセス、方法、部品または装置は、これらの要素だけでなく、請求の範囲に係るプロセス、方法、物品または装置に明示的ではないものの一連のまたは固有の他の要素を含み得る。
【００４９】
本発明は、特定の実施形態を参照して説明してきたが、これら実施形態は例示であって、本発明の範囲をこれらの実施形態に制限しないことを理解されたい。上述の実施形態に対する多くの変更、改変、追加および改良が可能である。特に、これらの変更は、コンピュータまたは他のデータ処理デバイス、コンピュータ読取り可能媒体（例えば、フロッピー（Ｒ）ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ等）、ソフトウェアを格納する格納デバイス、コンピュータメモリ等、ファームウェアまたは上述の方法を採用した他のプログラミングを含み得る。これらの変更、改変、追加および改良は、上掲の特許請求の範囲で規定される本発明の範囲内において行われることを企図する。
【図面の簡単な説明】
【図１】
図１は、二次元および三次元の単体（ｓｉｍｐｌｉｃｅｓ）を示す図である。
【図２】
図２は、物理的システムにおける実際の値、伝統的な有限要素技術を用いて生成された解、および本発明による技術を用いて生成された解の間の差分を示す図である。
【図３】
図３は、第１の二次元単体要素をわたる流体の流れを示す図である。
【図４】
図４は、図３に示す流れと反対方向に第１の二次元単体要素をわたる流体の流れを示す図である。
【図５】
図５は、第２の二次元単体要素をわたる流体の流れを示す図である。
【図６】
図６は、図５に示す流れと反対方向に第２の二次元単体要素をわたる流体の流れを示す図である。
【図７】
図７は、時間−空間ドメインの二次元投影を示す図である。
【図８】
図８は、本発明の１つの実施形態による方法を示すフローチャートである。
【図２】
2004502999.tif 000008
【図８】
2004502999.tif 000009
【国際公開パンフレット】
2004502999.tif 000010
2004502999.tif 000011
2004502999.tif 000012
2004502999.tif 000013
2004502999.tif 000014
2004502999.tif 000015
2004502999.tif 000016
2004502999.tif 000017
2004502999.tif 000018
2004502999.tif 000019
2004502999.tif 000020
2004502999.tif 000021
2004502999.tif 000022
2004502999.tif 000023
2004502999.tif 000024
2004502999.tif 000025
2004502999.tif 000026
2004502999.tif 000027
2004502999.tif 000028
2004502999.tif 000029
2004502999.tif 000030
2004502999.tif 000031
2004502999.tif 000032
2004502999.tif 000033
2004502999.tif 000034
2004502999.tif 000035
2004502999.tif 000036
2004502999.tif 000037
2004502999.tif 000038
2004502999.tif 000039
2004502999.tif 000040
2004502999.tif 000041
【国際公開パンフレット（コレクトバージョン）】
2004502999.tif 000042
【国際公開パンフレット（コレクトバージョン）】
2004502999.tif 000043
【国際調査報告】
2004502999.tif 000044
2004502999.tif 000045
2004502999.tif 000046