(57)【特許請求の範囲】
【請求項１】
地下媒体内の特定規模の空間における炭化水素（ＨＣ）の地質学的時間スケールでの移動をシミュレートするための、コンピュータによって実施される並列処理方法であって、
複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩から、ＨＣが排出されたことを特定するステップであって、前記複数のサブドメインは、分散メモリ高性能計算（ＨＰＣ）クラスタ内のそれぞれの計算コアに割り当てられる、ステップと；
潜在的なトラップピークを前記複数のグリッドセル内で識別するステップであって、前記潜在的なトラップピークは、前記ＨＣが到達したときにそれ以上移動しなくなる位置である、ステップと；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行するステップであって、前記侵入プロセスは、前記複数のサブドメインの一のサブドメインと隣接するサブドメインとの間の通信を含み、前記ＨＣは、前記計算コア間を並行に移動する、ステップと；
ＨＣ類を含むグリッドセルが過剰量のＨＣ類を含んでいるかどうかを特定するステップと；
前記グリッドセルが過剰量のＨＣ類を含んでいると判断した場合に、前記識別された潜在的なトラップピークに関連付けられたトラップにおける、グリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行するステップであって、前記集積プロセスは、非隣接サブドメイン間の通信を含み、前記侵入プロセスと前記集積プロセスは、定常状態に達するまで交互に行われる、ステップと；
前記グリッドセルが過剰量のＨＣ類を含んでいないと判断した場合に、前記トラップをＨＣのポテンシャル値で更新するステップと；を備える、
コンピュータによって実施される並列処理方法。
【請求項２】
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項３】
スピルポイントに達したことを特定するステップであって、前記スピルポイントは、前記トラップに充填された前記ＨＣが前記トラップ外に溢流するポイントである、ステップと；
前記侵入プロセスを実行するステップに戻るステップと；を更に備える、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項４】
前記侵入プロセスは、
前記ＨＣがサブドメイン境界に達したことを特定するステップと；
前記ＨＣを、隣接するサブドメインに移動するステップと；
前記ＨＣが同じサブドメインで充填トラップに達した場合、前記ＨＣを前記充填トラップと併合するステップと；を備える、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項５】
前記サブドメイン境界においてＨＣ量を通信するステップと；
サブドメイン境界を越えた前記ＨＣを、別のサブドメインに関連付けられたトラップと併合するステップと；を更に備える、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項６】
前記ＨＣがサブドメイン境界に達していることを特定するステップと；
前記サブドメイン境界に関連付けられた隣接するサブドメインに、ＨＣ量、最小の油相ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信するステップと；を更に備える、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項７】
前記ＨＣがトラップ境界に達したことを特定するステップと；
複数のトラップが集積境界を共有するかどうかを特定するステップと；を更に備える、
請求項１に記載のコンピュータによって実施される並列処理方法。
【請求項８】
地下媒体内の特定規模の空間における炭化水素（ＨＣ）の地質学的時間スケールでの移動をシミュレートするための、コンピュータによって実施される並列処理方法の１又は複数の命令を格納する非一時的なコンピュータ読取り可能媒体であって、前記１又は複数の命令は所定の操作を実行するためにコンピュータシステムによって実行可能であり、前記操作は、
複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作であって、前記複数のサブドメインは、分散メモリ高性能計算（ＨＰＣ）クラスタ内のそれぞれの計算コアに割り当てられる、操作と；
潜在的なトラップピークを前記複数のグリッドセル内で識別する操作であって、前記潜在的なトラップピークは、前記ＨＣが到達したときにそれ以上移動しなくなる位置である、操作と；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作であって、前記侵入プロセスは、前記複数のサブドメインの一のサブドメインと隣接するサブドメインとの間の通信を含み、前記ＨＣは、前記計算コア間を並行に移動する、操作と；
ＨＣ類を含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；
前記グリッドセルが過剰量のＨＣ類を含んでいると判断した場合に、前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作であって、前記集積プロセスは、非隣接サブドメイン間の通信を含み、前記侵入プロセスと前記集積プロセスは、定常状態に達するまで交互に行われる、操作と；
前記グリッドセルが過剰量のＨＣ類を含んでいないと判断した場合に、前記トラップをＨＣのポテンシャル値で更新する操作と；を備える、
非一時的なコンピュータ読取り可能媒体。
【請求項９】
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１０】
スピルポイントに達していることを特定する操作であって、前記スピルポイントは、前記トラップに充填された前記ＨＣが前記トラップ外に溢流するポイントである、操作と；
前記侵入プロセスを実行する操作に戻る操作と；を実行する１又は複数の命令を更に備える、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１１】
前記侵入プロセスは、
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記ＨＣを隣接するサブドメインに移動する操作と；
前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する操作と；を実行する１又は複数の命令を備える、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１２】
前記サブドメイン境界でＨＣ量を通信する操作と；
サブドメインの境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する操作と；を実行する１又は複数の命令を更に備える、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１３】
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油相ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信する操作と；を実行する１又は複数の命令を更に備える、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１４】
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
複数のトラップが集積境界を共有するかどうかを特定する操作と；を実行する１又は複数の命令を更に備える、
請求項８に記載の非一時的なコンピュータ読取り可能媒体。
【請求項１５】
地下媒体内の特定規模の空間における炭化水素（ＨＣ）の地質学的時間スケールでの移動をシミュレートするための、コンピュータによって実施される並列処理システムであって、
コンピュータメモリと；
前記コンピュータメモリと相互作動可能に結合され、所定の操作を実行するように構成されたハードウェアプロセッサであって、前記操作は、
複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作であって、前記複数のサブドメインは、分散メモリ高性能計算（ＨＰＣ）クラスタ内のそれぞれの計算コアに割り当てられる、操作と；
潜在的なトラップピークを前記複数のグリッドセル内で識別する操作であって、前記潜在的なトラップピークは、前記ＨＣが到達したときにそれ以上移動しなくなる位置である、操作と；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作であって、前記侵入プロセスは、前記複数のサブドメインの一のサブドメインと隣接するサブドメインとの間の通信を含み、前記ＨＣは、前記計算コア間を並行に移動する、操作と；
ＨＣを含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；
前記グリッドセルが過剰量のＨＣ類を含んでいると判断した場合に、前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作であって、前記集積プロセスは、非隣接サブドメイン間の通信を含み、前記侵入プロセスと前記集積プロセスは、定常状態に達するまで交互に行われる、操作と；
前記グリッドセルが過剰量のＨＣ類を含んでいないと判断した場合に、前記トラップをＨＣのポテンシャル値で更新する操作と；を備える、前記ハードウェアプロセッサと；を備える、
コンピュータによって実施される並列処理システム。
【請求項１６】
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
請求項１５に記載のコンピュータによって実施される並列処理システム。
【請求項１７】
スピルポイントに達したことを特定する操作であって、前記スピルポイントは、前記トラップに充填された前記ＨＣが前記トラップ外に溢流するポイントである、操作と；
前記侵入プロセスを実行する操作に戻る操作と；を実行するように更に構成される、
請求項１５に記載のコンピュータによって実施される並列処理システム。
【請求項１８】
前記侵入プロセスは、
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記ＨＣを、隣接するサブドメインに移動する操作と；
前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する操作と；を実行するように更に構成される、
請求項１５に記載のコンピュータによって実施される並列処理システム。
【請求項１９】
前記サブドメイン境界でＨＣ量を通信する操作と；
サブドメイン境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する操作と；を実行するように更に構成される、
請求項１５に記載のコンピュータによって実施される並列処理システム。
【請求項２０】
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油相ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信する操作と；を実行するように更に構成される、
請求項１５に記載のコンピュータによって実施される並列処理システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、２０１７年６月２３日に出願された米国仮出願第６２／５２４，２２３号及び２０１８年６月１３日に出願された米国特許出願第１６／００７，１７５号の優先権を主張し、その内容全体は参照により本明細書に組み込まれる。
【背景技術】
【０００２】
堆積盆（ｂａｓｉｎ）モデリング／シミュレーションは、石油システムモデリング／シミュレーションとしても知られ、堆積盆モデルのデジタル化と、堆積盆に関連付けられた相互関係のあるプロセスのシミュレーションとにより、数億年という地質学的時間スケールで堆積盆の進化とその流体含有量を追跡する。近年、この堆積盆モデリングは、地球探査科学者が炭化水素（ＨＣ）流体の存在とその種類を予測し、探査坑井を掘削してＨＣ流体を回収する前に行う地質リスク評価にとって重要なツールとなっている。堆積盆シミュレータには、１）逆抽出（ｂａｃｋ−ｓｔｒｉｐｐｉｎｇ）と圧縮（ｃｏｍｐａｃｔｉｏｎ）、２）圧力計算、３）熱流解析と動力学、４）石油の生成、吸着、及び排出、５）ＨＣ流体相の挙動、６）ＨＣのマイグレーション（ｍｉｇｒａｔｉｏｎ、以下「移動」ともいう）とエントラップメント（ｅｎｔｒａｐｍｅｎｔｓ、溜まり、閉じ込め）、を計算する数値モジュールが含まれていることが多い。最近、インベージョンパーコレーション（Ｉｎｖａｓｉｏｎ  Ｐｅｒｃｏｌａｔｉｏｎ、ＩＰ）は、二次ＨＣマイグレーションをモデル化する一般的な方法になっている。堆積盆モデリングソフトウェアの現行の商用リリース版には、シリアル計算専用のＩＰ手法が実装されている。これは、数百万セルのしきい値までの堆積盆モデルに対しては十分に高速である。しかし、このしきい値は、ＨＣマイグレーションの詳細なグリッドモデリングが必要な大規模堆積盆シミュレーションに対しては望ましくない制約となる。
【発明の概要】
【０００３】
本開示は、分散メモリ、高性能計算（ＨＰＣ）クラスタの対象となる２相の油水（ｏｉｌ−ｗａｔｅｒ）システムのための並列処理のインベージョンパーコレーション（ＩＰ）マイグレーション方法について述べる。
【０００４】
実施において、炭化水素（ＨＣ）は、１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩（ｓｏｕｒｃｅ  ｒｏｃｋ）から排出されたものであると特定される。潜在的なトラップピークが、複数のグリッドセル内で識別される。侵入（インベージョン）プロセスは、ＨＣが複数のグリッドセル内で移動を停止するまで実行される。ＨＣを含むグリッドセルに過剰（体積）量のＨＣが含まれているかどうかが特定される。集積プロセスを実行して、識別された潜在的なトラップピークに関連付けられたトラップにおいて、グリッドセルへのＨＣ充填をモデル化する。トラップはＨＣポテンシャル値でもって更新される。
【０００５】
前述の主題は、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータにより実施されるシステムを用いて実施可能である。非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータ読取可能な命令を格納（保存）する。コンピュータにより実施されるシステムは、１又は複数のコンピュータと相互運用可能に結合された１又は複数のコンピュータメモリデバイスと、命令を格納する有形の非一時的な機械読取可能媒体と、を備え、１又は複数のコンピュータによって実行されるとき、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法／コンピュータ読取可能な命令を実行する。
【０００６】
本明細書で述べる主題は、次の利点の１又は複数を実現できるように、特定の実施において実施される。すなわち、第１に、従来の市販の堆積盆シミュレータは、ＩＰ手法を用いてシリアル処理を実行するように構成されている。従来の堆積盆シミュレータの構成は、堆積盆シミュレーションの実行に必要な計算資源のために、データモデルのサイズと結果の出力解像度とを制限する。第２に、ここで述べる方法論は、分散メモリのＨＰＣクラスタ用に構成される。その結果、データモデルのサイズと出力解像度は、従来の堆積盆シミュレータのように制限されることはない。第３に、ここで述べる方法論と、ＨＰＣクラスタの使用とにより、計算パフォーマンスとスケーラビリティが向上する。超大規模堆積盆シミュレーションを、従来の堆積盆シミュレータ他の方法で必要とされる期間よりもはるかに短い期間で実行でき、堆積盆モデリングのための高解像度のマイグレーションシミュレーションを実行できる性能を、地球科学者やエンジニアに提供する。
【０００７】
本願に記載された主題における１又は複数の実施の詳細は、発明の詳細な説明、特許請求の範囲及び添付図面に記載されている。本主題の他の特徴、態様及び利点は、発明の詳細な説明、特許請求の範囲及び添付の図面から明らかになるであろう。
【図面の簡単な説明】
【０００８】
本特許又は本願ファイルには、カラーで作成した少なくとも１枚のカラー図面が含まれている。カラー図面を伴う本願の公開公報のコピーは、米国特許商標庁へ申請し、必要な料金を納付することにより提供される。
【０００９】
【図１Ａ】図１Ａは、一実施に係る、分散メモリ高性能計算（ＨＰＣ）クラスタを対象とする２相油水システムのための並列処理インベージョンパーコレーション（ＩＰ）マイグレーション方法を示すフローチャートである。
【００１０】
【図１Ｂ】図１Ｂは、一実施に係る、図１Ａの枠１２１で描いた方法ステップを改めて説明するフローチャートである。
【００１１】
【図１Ｃ】図１Ｃは、一実施に係る、図１Ａの枠１５１で描いた方法ステップを改めて説明するフローチャートである。
【００１２】
【図２Ａ】図２Ａは、一実施に係る、内側ハロー（ｈａｌｏ）と外側ハローとの間の境界におけるトラップピークの識別を示すブロック図である。
【００１３】
【図２Ｂ】図２Ｂは、一実施に係る、外側ハローの隣接値に基づく内側ハローのトラップピークの識別を示すブロック図である。
【００１４】
【図３】図３は、一実施に係る、流れルーチン中の炭化水素（ＨＣ）交換を示すブロック図である。
【００１５】
【図４Ａ】図４Ａは、一実施に係る、グリッドセルが過剰なＨＣ量を含んでいるとしてマーク付けされた最初のステップを示すブロック図である。
【００１６】
【図４Ｂ】図４Ｂは、一実施に係る、過剰なＨＣ量を含んでいるとしてマーク付けされたトラップピークを示すブロック図であって、最小ポテンシャル値を有する侵入された隣接セルは、集積されたグリッドセルとなり、すべて、より大きいポテンシャル値を有する。
【００１７】
【図５】図５は、一実施に係る、０から３の番号を付された４つの異なるサブドメインにまたがる集積を示すブロック図である。
【００１８】
【図６】図６は、一実施に係る、複数のサブドメインにわたって集積が行われるときに実行される通信（ｃｏｍｍｕｎｉｃａｔｉｏｎ）を示すブロック図である。
【００１９】
【図７】図７は、一実施に係る、１又は複数のサブドメインから情報を受信する１つのサブドメインを示すブロック図である。
【００２０】
【図８】図８は、一実施に係る、単一の計算コアでテストケース1モデルを実行した後の結果を示すプロットである。
【００２１】
【図９】図９は、一実施に係る、複数のトラップピークを含む集積体を示すプロットである。
【００２２】
【図１０】図１０は、一実施に係る、５０個の計算コアを用いて実行された、図８及び図９のテストケースの結果のプロットである。
【００２３】
【図１１】図１１は、一実施に係る、異なるサブドメインに存在する併合のトラップピークを示すプロットである。
【００２４】
【図１２】図１２は、一実施に係る、テストケース２ａの連続実行からの結果を示すプロットである。
【００２５】
【図１３】図１３は、一実施に係る、複数のサブドメインにわたって発生し、互いに直接隣接していないサブドメインにまたがるテストケース２ａの併合を示すプロットである。
【００２６】
【図１４Ａ】図１４Ａは、一実施に係る、単一の計算コアでテストケース２ａモデルを実行した結果を示すプロットである。
【００２７】
【図１４Ｂ】図１４Ｂは、一実施に係る、１，０００個の計算コアを用いてテストケース２ａモデルを実行した結果を示すプロットである。
【００２８】
【図１５】図１５は、一実施に係る、１，０００個の計算コアを用いてテストケース２ｂモデルを実行した結果を示すプロットである。
【００２９】
【図１６】図１６は、一実施に係る、テストケース１についてプロットされた計算コアの数に対する正規化された実行時間のグラフである。
【００３０】
【図１７】図１７は、一実施に係る、テストケース２ａ／２ｂについてプロットされた計算コアの数に対する正規化された実行時間のグラフである。
【００３１】
【図１８】図１８は、一実施に係る、本願で説明される、記述されたアルゴリズム、記述された方法、記述された機能、記述されたプロセス、記述されたフロー、及び記述された手順に関連付けられた計算機能を提供するために使用された例示的なコンピュータシステムを示すブロック図である。
【００３２】
種々の図面における類似の参照番号及び記号は、同様の要素を指す。
【発明を実施するための形態】
【００３３】
以下の詳細な説明は、分散メモリ、高性能計算（ＨＰＣ）クラスタを対象とする２相油水システムの並列処理インベージョンパーコレーション（ＩＰ）マイグレーション方法について述べ、当業者が１又は複数の特定の実施の文脈で開示された主題を作り上げて、これを使用できるようにするために提示される。開示された実施において様々な改変、変更、及び並べ替えを行えることは当業者には容易に明らかであり、また定義された一般的な原理は、本開示の範囲から逸脱することなく、他の実施及び用途に適用するようにしてもよい。本開示は、記載又は図示の実施に限定されることを意図してはおらず、記載する原理及び特徴に整合する最も広い範囲が付与されることを意図している。
【００３４】
石油システム／堆積盆は、炭化水素（ＨＣ）集積の存在に必要な地質要素とプロセスとからなる地質システムである。地質要素として含まれるのは、通常は、根源岩、貯留岩、シール岩、及び表土岩である。ＨＣ集積に必要なプロセスとして、トラップの形成、ＨＣ類の連続生成、排出、移動、及び集積が含まれる。
【００３５】
堆積盆モデリング／シミュレーションは、石油システムモデリング／シミュレーションとしても知られ、堆積盆モデルのデジタル化と堆積盆に関連付けられた相互関係のあるプロセスのシミュレーションとにより、数億年の地質学的時間スケールで堆積盆の進化とその流体含有量を追跡する。近年、これは、地球探査科学者がＨＣ流体の種類と存在を予測し、ＨＣ流体を回収するために探査坑井を掘削する前に地質リスクを評価するための重要なツールとなっている。堆積盆シミュレータには、１）逆抽出と圧縮、２）圧力計算、３）熱流解析及び動力学、４）石油の生成、吸着、及び排出、５）ＨＣ流体相の挙動、並びに、６）ＨＣマイグレーションとエントラップメント、を計算するための数値モジュールが含まれていることが多い。最近、ＩＰは、二次ＨＣマイグレーションをモデル化する一般的な方法になっている。
【００３６】
従来の堆積盆シミュレーションでは、Ｄａｒｃｙの法則の多相流拡張が適用されている。Ｄａｒｃｙの法則を用いる３次元（３Ｄ）の超大規模堆積盆モデルのＨＣマイグレーションを十分詳細にシミュレートするための計算要件は、通常、ＨＣ探査会社及び生産会社がアクセス可能な典型的なコンピュータハードウェア資源を超えており、コンピュータ処理時間のかなりの部分を消費することになる。ルールに基づくＩＰ手法を用いたＨＣマイグレーションのシミュレーションは、計算量の少ない挑戦的な方法であり、従来の慣例よりも高い解像度で二次ＨＣマイグレーションシミュレーションを実行できる。堆積盆モデリングソフトウェアの現行の商用リリースには、ルールベースのＩＰ法をシリアルに実行するために用いられるシリアル計算構成を備えている。通常は、最高数百万セルのしきい値までの堆積盆モデルにシリアルな方法が用いられるが、このしきい値は、ＨＣマイグレーションの細かいグリッドモデリングにより得られる高解像度のために望まれる大規模な堆積盆シミュレーションにとって望ましくない制約である。
【００３７】
パーコレーション理論は、一定の印加圧力における毛管力の影響下での多孔質媒体内の流動現象を説明するものである。パーコレーション理論と、毛管力の作用下での、多孔質媒体中の或る流体の別の流体による置き換え（ｄｉｓｐｌａｃｅｍｅｎｔ）との関係が研究されてきた。用語ＩＰは、一定圧力ではなく一定流量での、或る流体の別の流体による置き換えに関する用語である。加えて、選択されたしきい値抵抗へのインタフェースの前進とは対照的に、最小の抵抗ポイントへのインタフェースの前進が導入された。その違いが、ここで述べるＩＰ手法へと導いたのである。
【００３８】
この新しいＩＰ手法では、媒体が、サイト（ｓｉｔｅ）と結合（ｂｏｎｄ）とのネットワークとして特徴付けられる。シミュレーションでは、この特徴はグリッドセルと隣接するグリッドセルとのグリッドセル接続に変換される。新しいＩＰ手法は、ネットワークのすべてのサイトが常駐（防御（ｄｅｆｅｎｄｉｎｇ））流体で充填されていることを前提としている。システムに導入される侵入流体が防御流体と置き換わる。毛管力は、ＨＣがシステム内でどのように移動するかという点で最大の役割を果たす。この移動は、浮力により促進され、侵入流体と常駐流体との間の毛管圧により阻害されることになり、ＨＣ流体（例えば、油）は毛管圧に打ち勝って岩の細孔から隣接する細孔スロートを通らねばならなくなる。
【００３９】
本開示のために、記号とその意味を以下に列挙する。
ｇ  重力定数
Ｐｃｅ  毛管入口圧力
Ｐｃｏｗ  油水毛管圧
Ｐｏ  油圧
Ｐｗ  水圧
ＰＢ  浮力圧力
Ｓｃｏ  臨界油飽和度
Ｓｗｃ  遺留水飽和度
Ｖ  （体積）量
ｚ  基準からの深さ
Φ０  油ポテンシャル
Φｗ  水ポテンシャル
ｐ０  油密度
ｐｗ  水密度
φ  気孔率
上に列挙した記号と通常の意味との間に矛盾がある限りにおいて、列挙した意味が優先するものとする。記号の意味が依然として明確でない限りにおいて、当業者が理解する記号の意味を前提とする。
【００４０】
高レベルでは（概括すれば）、移動する前線を備える、非標準で非典型的な、並列処理による新しいＩＰ手法によって、堆積盆シミュレーションでのＨＣマイグレーションが加速され、もっと大規模で解像度の高いモデルを用いることが可能になる。この新しいＩＰ手法は、ＨＣマイグレーションと集積（ａｃｃｕｍｕｌａｔｉｏｎ）の両方に対して分散メモリの並列処理を実行する。堆積盆ドメインは、この新しいＩＰ手法によっていくつかの計算サブドメインに分割されるが、各サブドメインに必要な作業量は、他のサブドメインとおおよそバランスが取れている。各サブドメインは、ＨＰＣクラスタ内（１又は複数の計算ノード内）の中央処理装置（ＣＰＵ）の割り当てられたワークグループの計算コアに割り当てられる。ＨＰＣクラスタの各計算ノードには、１又は複数のＣＰＵ／プロセッサが含まれていることに留意されたい。各プロセッサには、１又は複数の計算コアが含まれている。
【００４１】
分散したＨＣは、計算プロセス全体で根源（ｓｏｕｒｃｅ）からトラップまで並列に移動する。移動後に集積プロセスが行われ、そこではトラップが動的に成長し、複数のサブドメインを横断して併合する可能性がある。スピル（ｓｐｉｌｌ、溢流）や漏出に達すると、又は、過剰なＨＣ（体積）量が排出されると、集積は停止する。移動と集積処理は定常状態に達するまで交互に行われる。さらに、この新しいＩＰ手法はプロセス間通信を最小限に抑え、優れたスケーラビリティ（ｓｃａｌａｂｉｌｉｔｙ、拡張性）を実現する。シリアルＩＰ手法は、計算ノードで使用可能なメモリストレージと、移動及び集積をシリアルに処理するために必要な実行時間とによって制約を受ける。新しいＩＰ手法は、より多くの計算ノードを用いてモデルのストレージ要件を分散でき、個別の計算プロセスに存在する複数の情報（ｔｉｐ）の同時処理を用いて、マイグレーションステップと集積ステップとを交互に実行（ｉｎｔｅｒｌｅａｖｉｎｇ）することで処理のスピードアップが図れる。スピルと漏出も管理され、その後に三次マイグレーションが行われる場合がある。ＩＰ手法を用いるマイグレーションのすべてのコンポーネントは並列化される。記載する新しいこのＩＰ手法は汎用的であり、あらゆるタイプの潜在的なフィールドに適用できる。
【００４２】
堆積盆モデルには、さまざまなタイプの複雑性と岩石タイプとが含まれてもよい。この方法は、石油システムを複数のサブドメインに分散することと、集積段階で必要な非隣接サブドメイン間の通信（ｃｏｍｍｕｎｉｃａｔｉｏｎ）とから生じる複雑さを扱う。さらに、この方法は、２つのトラップ間で併合が発生した場合に、又は移動する細脈（ｓｔｒｉｎｇｅｒ）が集積された油（オイル）を含むトラップと併合した場合に、必要な通信を扱う。この通信では、移動細脈が１つのサブドメインでトラップに入る可能性があるが、同じトラップの漏出又はスピルポイントは異なるサブドメインに存在する可能性があるので、非隣接サブドメイン間の通信が必要になる場合もある。
【００４３】
移動（マイグレーション）は、排出されたＨＣ流体（例えば、油）の、根源岩からトラップ又は集積（トラップ境界で発生）への動きをモデル化する。移動中に、排出されたＨＣ流体はキャリア岩を通って動くことができる。堆積盆モデルには複数の根源岩が含まれる場合があり、各根源岩には任意の数の根源ポイントがあり、複数のサブドメインに広がる場合がある。根源が１つのサブドメインに存在し、トラップが別のサブドメインに存在する可能性がある。
【００４４】
移動経路は複数のサブドメインを通過して、離れた集積トラップに接続する場合がある。埋め戻し（ｂａｃｋｆｉｌｌｉｎｇ）プロセスでは、トラップが複数のサブドメインにまたがって成長する場合がある。複数の集積トラップが合体してより大きなトラップを形成する一方で、スピルと漏出は正確にモデル化され、トラップの位置と集積量を正しく予測する。集積プロセスは、スピル又は漏出が発生したとき、又は過剰なＨＣ流体が排出されたときに停止する。最終結果に至るために、すべてのプロセスは複数のサブドメインにまたがる場合があり、そして、さまざまな形式のプロセス間通信（例えば、送信又は受信）又は集合プロセス（例えば、ルートプロセスではなく複数のプロセスから結果を集める）を伴う場合がある。
【００４５】
典型的な実施において、この方法は、トラップ集積、併合、及びスピルと漏出を管理するための、一意の並列分散データシステム及び関連する方法を含む。データシステムには、各トラップが一意に識別されるトラップ配列が含まれている。石油システムのトラップが推定され、トラップピークのグローバルセル番号が識別子として用いられる。トラップは、ＨＣがトラップへ移動するまでは非アクティブである。このように、格納データは、もしトラップがサブドメインに存在すれば、そのトラップに関連するそのサブドメインに存在する指標のみが境界リストに格納されるデータである。トラップピークのグローバルセル指標を用いて、どのリストがどのトラップに属するかを識別する。集積トラップが成長し、その境界がサブドメインの外側ハロー（ｈａｌｏ）に入ると、リストは分割され、外側ハローの指標が、対応する内側ハロー値にマッピングされ、それらのセル指標が実際に存在するサブドメインへ送信される。さらに、侵入に対する抵抗が最も少ないセルの指標を特定するために、対応する指標のポテンシャル値も配列に保存される。すべてのアクティブトラップをグローバルに追跡し、集積プロセスの実行中に正しいリストを参照するために、必要な情報を用いてアクティブなトラップ配列が構成される。ＨＣには、貯留層内の異なる流体（水攻法の水及びＨＣなど）の間を移動する前線（又は境界）がある。例えば、水とＨＣに関しては、より多くの水が押し込まれ、ＨＣが抽出されるにつれて、前線が地層を通過する。本開示では、この前線を、侵入及び集積の先端（ｉｎｖａｓｉｏｎ  ａｎｄ  ａｃｃｕｍｕｌａｔｉｏｎ  ｔｉｐ）と呼ぶこともある。特定のトラップの侵入先端は、常にアクティブなトラップ配列に維持される。２つのトラップ間の併合が発生すると、トラップピークのグローバル値もアクティブなトラップ配列に併合される。
【００４６】
さらに、典型的な実施では、侵入プロセスと集積プロセスは別々であると見なされ、この方法は、これら２つのメインプロセス間で交互に実行される。侵入プロセスには、ＨＣがサブドメインの外側ハローに達するまで局所的に最小抵抗の経路をたどって高ポテンシャル部位から低ポテンシャル部位へ移動する、侵入メインルーチンを含むドライバが含まれる。この侵入ドライバルーチンは、このサブドメインの外側ハローから隣接サブドメインの内側ハローまでで起こる通信を処理する。通信が完了すると、ＨＣが新しい侵入先端から、より低いポテンシャルの隣接セルへ移動できる限り、そのサブドメインで侵入ルーチンが続行される。侵入プロセスは、ＨＣの流れが発生しなくなるまで、通信及び流れのルーチンを交互に繰り返す。この段階で、根源ポイントからのすべてのＨＣがトラップピークへの移動を終えている。侵入プロセスが完了すると、方法は集積段階に入る。ここでも、集積ドライバが必要な通信を処理し、集積ルーチンがローカルな集積を処理する。集積プロセスは、ＨＣ又は境界リストが外側ハローに達したときに、ＨＣのローカル（局所的）な埋め戻しと通信とを交互に繰り返す。スピル又は漏出が発生すると、対応するトラップ集積プロセスは停止し、プロセスは集積がそれ以上できなくなるまで他のトラップのＨＣを埋め戻し続ける。その後、集積プロセスは終了し、スピル段階又は漏出段階にＨＣがある場合、本手法は再び侵入プロセスに入り、ＨＣの移動が始まる。この方法は、ＨＣが石油システムへ均等に分配されるまで、これら２つのメインプロセスを交互に実行する。
【００４７】
メッセージパッシングインタフェース（Ｍｅｓｓａｇｅ  Ｐａｓｓｉｎｇ  Ｉｎｔｅｒｆａｃｅ、ＭＰＩ）標準は、分散メモリＨＰＣクラスタを対象とする並列処理方法を実施するために用いられる。ここで述べる方法論とＨＰＣクラスタとを用いることにより、計算パフォーマンスとスケーラビリティが向上する。超大規模な堆積盆シミュレーションは、従来の又は他の方法で必要とされるよりもはるかに短い期間で実行でき、堆積盆モデリングのための高解像度のマイグレーションシミュレーションを地球科学者及びエンジニアが実行できる性能をもたらす。
【００４８】
この開示では、新しいＩＰ手法について説明した後、この新しいＩＰ手法における新しい並列処理の実施の数値成分について順を追って説明する。先に述べたように、実施は２相油水システムであり、１層又は複数層の複数の根源から排出されるＨＣ流体の移動をシミュレートすることができる。実施について、１４億セルの堆積盆モデルを含む数値例を用いて更に説明する。この数値例では、並列処理の新しいＩＰ手法の実施の正確性とスケーラビリティとが検証される。
【００４９】
新しいＩＰ手法
【００５０】
ＩＰを二次ＨＣマイグレーションに適用する場合、グリッドセルの防御流体は湿潤流体（例えば水）と見なされ、侵入する非湿潤流体（例えば油）によって置き換えられる。毛管入口圧力Ｐｃｅは湿潤相（水）に対する非湿潤相（油）のための毛管しきい圧力としても知られ、Ｓｗ＝１−Ｓｃｏでの毛管圧力値である。ここでＳｃｏは臨界油飽和度である（つまり、ＨＣは流れない）。毛管圧は、非湿潤相が湿潤相で占められている細孔に侵入するための入口圧力である。２相油水システムでは、毛管圧は次のように定義される。
【数１】
0006966577.tif 000002
【００５１】
新しいＩＰ手法において、水相は連続相と見なされる。水相の圧力とポテンシャルは、堆積盆シミュレーションに割り当てられた時間ステップで解決される。水相の圧力解から、油相のＨＣ相ポテンシャルを計算できる。
【００５２】
水相ポテンシャルは、次の方法で堆積盆の時間ステップで水圧解から並列に計算できる。
【数２】
0006966577.tif 000003
【００５３】
すべてのグリッドセルの油相ポテンシャルは以下のように定義できる。
【数３】
0006966577.tif 000004
ＣＯ２隔離問題（ｓｅｑｕｅｓｔｒａｔｉｏｎ  ｐｒｏｂｌｅｍ）のためにインベージョンパーコレーションの別の派生が示されており、ここで水相ポテンシャルは一定であり、油の動きは浮力と毛管現象との相互作用の結果である。式（３）で、zはデータ基準（ｄａｔｕｍ）の深さである。隣接するすべてのグリッドセルと比較したグリッドセルの油相ポテンシャルの違いにより、新しいＩＰ手法で油がとる移動経路が特定される。油は、ポテンシャルの大きいセルからポテンシャルの低いセルへ移動する。セル内の余分な油は、ポテンシャル値が最小の隣接セルへ移動する。細脈としても知られる移動経路の残りの油飽和は、臨界油飽和であり、油の除去量は次のように定義できる。
【数４】
0006966577.tif 000005
【００５４】
ポテンシャル値がすべての隣接セルよりも小さいグリッドセルに油が達すると、集積が始まる。トラップは、最初に集積されたセルから、油のポテンシャルが最も低い隣接セルまで満たされる。充填油量は以下のように計算できる。
【数５】
0006966577.tif 000006
ここで、ＳＷＣは、遺留（ｃｏｎｎａｔｅ）水飽和度である。
【００５５】
トラップ内に集積された油は、垂直平衡条件を満たさねばならない。したがって、トラップ内の充填されたグリッドセルのすべては、同じポテンシャル値に達する。この充填プロセス中におけるトラップされた油のポテンシャルは上昇していく必要があり、ポテンシャルは少なくとも最後に充填されたセルのポテンシャルでなければならない。油の集積は離散的に発生し、各グリッドセルの厚さと深さは異なる可能性があるので、トラップのポテンシャル値は、最後に充填されたセルのポテンシャルで更新され、増分浮力ポテンシャルはセルの厚さの半分に等しくなる。
【数６】
0006966577.tif 000007
式（６）は、水と油の接触深さは、それが充填前のトラップポテンシャルよりも大きい値である場合において、最後に充填されたグリッドセルの底に近いことを意味する。
【００５６】
新しいＩＰ手法を用いるマイグレーションと集積に関連するプロセスには、多くの並列化の課題がある。分散メモリ並列処理では、堆積盆ドメインはいくつかのサブドメインに分割され、各サブドメインは、ＨＰＣクラスタのＣＰＵの計算コアによって処理される計算プロセスに割り当てられる。したがって、先に述べたように、トラップは複数のサブドメインにまたがることができ、根源ポイントは複数のサブドメインにまたがることもできる。根源岩からトラップへのマイグレーションプロセスをシミュレートする際、ＨＣはある計算プロセスから別の計算プロセスへ移動する場合がある。これには、１つのサブドメインから隣接サブドメインへの移動情報を交換するためのデータ通信が必要となる。
【００５７】
並行して、隣接するグリッドセルのリストを作成する必要があるため、埋め戻しは計算上の難しい課題となる。そのグリッドセルの最小抵抗は、対応する油ポテンシャルから特定される。これらのリストは、侵入する次善（２番目の）のグリッドセルを特定するときに用いられる。トラップは複数のサブドメインにまたがることができるため、侵入する次善のグリッドセルは、非隣接サブドメインに存在し得る。埋め戻しの並列化により、ＨＰＣクラスタの異なる計算ノードに割り当てられる可能性のある非隣接サブドメイン間で必要な通信が行われる。
【００５８】
典型的な実施では、プロセス間通信に、先に述べたＭＰＩ標準が用いられる。外側ハローは、隣接サブドメインに属するグリッドセルの外側の境界層であるが、情報は現在のサブドメインで維持され、並列プロセスの同時計算を容易にする。内側ハローは、現在のサブドメインに属するグリッドセルの内側の境界層である。
【００５９】
図１Ａは、一実施に係る、分散メモリ及びＨＰＣクラスタを対象とする２相油水システムのための並列処理ＩＰ方法１００ａを示すフローチャートである。提示を明確にするために、以下の説明は、本明細書の他の図の文脈で方法１００ａを一般的に説明する。しかしながら、方法１００ａは、例えば、任意の適切なシステム、環境、ソフトウェア、及びハードウェア、又は、システム、環境、ソフトウェア、及びハードウェアの適切な組み合わせによって実行できることは言うまでもない。実施によっては、方法１００ａのさまざまなステップは、並行して、組み合わせて、ループで、又は任意の順序で実行することもできる。
【００６０】
図１Ａでは、方法１００ａは、並列侵入プロセス（例えば、図１Ａの上半分）と並列集積プロセス（例えば、図１Ａの下半分）とを交互に実行する（インターリーブする）。方法１００ａが対象とするのは並列での実施であって、そこでは、課題を解決するためにいくつかの同時並行プロセスが機能している。例えば、説明のように、侵入（入り込む（ｉｎｖａｄｅ））（１２０）の間、各プロセスは、すべての細脈（ここで細脈はそのプロセスのために存在する）がトラップに併合（マージ）される（１４５）まで、又は境界で失われるまで、又は現プロセスに属するサブドメインから移動して出る（１４０）まで、細脈の移動を処理する。集積（１５５）の間、すべての過剰なＨＣを用いて、複数のサブドメイン（プロセス）にまたがるトラップを構築する。ステップ１６５は、一プロセスから別のプロセスへの集積先端のシフトを処理する。トラップが大きくなると、別のトラップと結合する場合がある。ステップ１７５はこの状況でのトラップの併合を処理する。
【００６１】
ステップ１１０では、ＨＣが排出されたかどうか（例えば、根源岩のケロジェン（ｋｅｒｏｇｅｎ、油母）がＨＣに反応して排出するのに十分な温度変化を受けた場合）を特定する。ＨＣが排出されたことが特定された場合（例えば、１つのサブドメインに、排出されたＨＣ／根源ポイントが含まれている場合）、集合的な通信サブルーチンを用いて、すべてのサブドメインが連帯してステップ１１５へ進む必要があるすべてのプロセスと通信する。そうではなく、ＨＣが排出されていないことが特定された場合、方法１００ａはステップ１９０へ進む。
【００６２】
ステップ１１５では、すべての潜在的なトラップピークが識別／マーク付けされる。ここで、方法１００ａは、絶対油相ポテンシャル値をチェックし、特定のセルを見つけるが、ここでは、すべての隣接セルのポテンシャル値は、特定のセルのポテンシャル値より大きい。
【００６３】
図２Ａに移ると、図２Ａは、一実施に係る、内側ハローと外側ハローとの間の境界におけるトラップピーク識別（凡例２０５を参照）のブロック図２００ａである。ここでは、グリッドセル２１０に極小値（５．８）が示されている。グリッドセル２１０がトラップピークとして識別されると、グリッドセル２１０は（図２Ｂに示すように）トラップピークとしてマーク付けされるが、ＨＣがそこへ移動するまでは非アクティブのままである。識別には、すべてのサブドメインの外側ハローのグリッドセルに油相ポテンシャル値を入力するために並列化が必要である。トラップピークは内側ハローに存在する可能性があるため、その特定を行うには外側ハロー油相のポテンシャル値が必要である。
【００６４】
図２Ｂに移ると、図２Ｂは、一実施に係る外側ハローの隣接値に基づく内側ハローのトラップピークの識別を示すブロック図２００ｂである。図示されているように、外側ハローにあるグリッドセル２１５の値（６．２）は、５．８が実際にその隣接セルの最小油相ポテンシャル値であることを特定するために必要である。説明したアプローチでは、現在のグリッドセルとセルの面を共有する隣接グリッドセルを考慮することに留意されたい。したがって、グリッドセル２１０（５．８の値）に対して斜めのグリッドセルは考慮されない。当業者によって理解されるように、他の実施において、説明されたアプローチを、対角グリッドセル又は３次元（又はそれ以上）のグリッドセルの比較を可能にするために修正できる。したがって、グリッドセル２１０はトラップピークとして識別される。
【００６５】
識別されたトラップピーク（グリッドセル２１０）のグローバル指標は、ローカル指標からマッピングされ、保存される。これらのトラップピーク指標は、先に述べたトラップロジックに用いられる。グローバル指標は、トラップとそれに関連する埋め戻されたセル、境界リスト、及び油相ポテンシャル値を区別及び識別するために用いられる。必要なトラップ情報にアクセスするために、すべてのアクティブなトラップのグローバル配列のローカルコピーと関連情報とが、すべてのプロセッサに保存されるように構成される。
【００６６】
典型的な実施では、グローバルアクティブトラップ配列には表１に示す情報が含まれる。
【表１】
0006966577.tif 000008
他の実施では、より多くのまたはより少ない情報を表１に含めることができる。図１Ａへ戻り、方法１００ａはステップ１１５からステップ１２０へ進む。
【００６７】
ステップ１２０で、侵入プロセスが始まる。侵入プロセスには、方法１００ａの複数のステップが含まれる（破線の枠１２１で示す）。ここで、排出サイトからのＨＣ類は、隣接するグリッドセルのより低い油相ポテンシャル値を探すというＩＰルールに従って移動する。より小さい油相ポテンシャル値が見つかった場合、油相ポテンシャル値が、最小のグリッドセルが次に侵入するセルとして選択される。ＨＣ量はそのサイトを離れて新しい侵入部位に入るが、式（４）のように計算されるクリティカルパスの飽和が残っているため、その量だけ減少する。移動経路にあるグリッドセルが新たに侵入され、既に残留できない量のＨＣが含まれている場合、その量を、移動される過剰なＨＣに追加する必要がある。根源ポイントから次の侵入サイトへのフローは、すべてのプロセッサでローカルに始まる。セル内にＨＣ量を有するすべてのプロセッサは侵入ルーチンを通過し、ＨＣはポテンシャル勾配に沿って、ＨＣの移動が停止するまで、高から低に移動する。方法１００ａは、ステップ１２０からステップ１２５へ進む。
【００６８】
ステップ１２５では、ＨＣが移動中であるかどうか判定される。方法１００ａは、ＨＣの移動を継続的にチェックする。ＨＣが移動中であると特定された場合、方法１００ａはステップ１２０へ戻り、侵入が継続する。そうではなく、ＨＣが移動中ではないと特定された場合、方法１００ａはステップ１３０へ進む。
【００６９】
ステップ１３０では、ＨＣの動きが停止しているかどうかの特定がなされる（例えば、ＨＣ／移動細脈が、充填されたトラップに達しているか？）。ＨＣの移動が停止していると特定された場合、方法１００ａはステップ１４５へ進む。そうではなく、ＨＣの移動は停止していないと特定された場合、方法１００ａはステップ１３５へ進む。
【００７０】
ステップ１４５では、移動細脈がトラップに併合される。方法１００ａは、ステップ１４５からステップ１２０へ戻る。
【００７１】
ステップ１３５では、ＨＣが特定のサブドメインの境界に達したことによりＨＣの移動が停止しているどうかの特定がなされる。参照される境界は、堆積盆システムの内部境界であるが、サブドメインの外部境界である。これは、サブドメインの外側ハローとも呼ばれる。ＨＣが特定のサブドメインの境界に達したことにより移動が停止していると特定された場合、方法１００ａはステップ１４０へ進む。そうではなく、ＨＣが特定のサブドメインの境界に達したことにより移動が停止しているのではないと特定された場合、方法１００ａはステップ１５０へ進み、集積プロセスが始まる（破線の枠１５１で示す）。
【００７２】
ステップ１４０では、一旦、ＨＣの動きがすべてのプロセッサで停止し、ＨＣ量が外側ハローに存在する場合、方法１００ａは外側ハローのすべてのＨＣを隣接サブドメインの内側ハローに通信する（送る）。
【００７３】
図３に移ると、図３は、一実施に係る、フロールーチン中のＨＣ交換を示すブロック図３００である。図３の内側及び外側ハローグリッドセルの識別については、凡例２０５を参照されたい。ここで、ＨＣがサブドメイン１  ３１０の外側ハローグリッドセル３０５に達していることがわかる。ステップ１４０で説明したように、方法１００ａは、外側ハローグリッドセル（例えば、サブドメイン１のグリッドセル３０５）上のすべてのＨＣを隣接サブドメインの内側ハローグリッドブロック（例えば、サブドメイン２  ３２０の内側ハローグリッドセル３１５）へ送る。説明の実施では、ＨＣは斜めには横断しないと見なされる。当業者にとっては言うまでもなく、他の実施では、説明するアプローチは、ＨＣが斜めに横断できるように改変することができる。
【００７４】
方法１００ａでは、ＨＣ類がどのプロセス上でも移動しなくなるポイントにシステムが達するまで、ＨＣ類をローカルに移動させること及びサブドメイン境界でＨＣ量を通信することを交互に行う。これ以上の移動がないということは、すべてのＨＣ類が、それ以上の低い油ポテンシャルが存在しない位置まで達し、これ以上移動できないことを示している（これらは以前に特定されたトラップピークである）。
【００７５】
実施によっては、ＨＣ類によっては外部境界を介して堆積盆システムから移動できるが、ＨＣ損失と見なされる。侵入プロセスが終了すると、方法１００ａは集積プロセスへ進む。並行して、方法１００ａは、すべてのプロセッサがマイグレーションプロセスを完了し、すべての排出されたＨＣが関連するトラップピークに達し、それ以上流れなくなるのを待つ。図１Ａへ戻り、方法１００ａは、ステップ１４０からステップ１２０に戻る。
【００７６】
ステップ１５０では、ＨＣ量を含むセルを介する検索が、ＨＣ量が過剰であることを示しているかどうか特定される。埋め戻しは、セル内のＨＣ量が移動する場所がもはや無く、集積飽和値より多い量が含まれているときに開始される。最初は、過剰状態であるセルが見つかるとするなら、それは以前に特定されたトラップピークにあるであろう。ＨＣがトラップピークに入ると、トラップピークがアクティブになり、埋め戻しが開始される。埋め戻しは、式（５）のように、充填されたセル内のすべての可動水を移動させ、ＨＣに置き換える。過剰条件の存在が特定された場合、方法１００ａはステップ１５５へ進む。そうではなく、過剰条件が存在しないと特定された場合、方法１００ａはステップ１８０へ進む。
【００７７】
スピルは、新たに侵入されたセルが、そのセルより油相ポテンシャル値が低い直接隣接するセルを有する場合に、埋め戻し中に発生する。ＨＣは、その後に移動可能で、埋め戻しプロセスの中断が発生することがある。集積プロセスは、すべてのトラップが充填し終えるまで、又はスピルポイントに達するまで継続し、方法１００ａは侵入プロセスに切り替わる。侵入プロセスに切り替わるまで、移動すべき過剰なＨＣ量を含むスピルポイント指標にフラグを付すことにより、ルーチンが残りのすべてのアクティブトラップの処理を完了するまで、埋め戻しプロセスによってセルをスキップできるようにする。
【００７８】
漏出（ｂｒｅａｋｔｈｒｏｕｇｈ）は、集積体の油ポテンシャルが、侵入するべき次善のグリッドセルがシールポイントの１つに属するポイントまで増加した場合に発生する。漏出ポイントは、シール上の最も弱いポテンシャルポイントになるであろう。シールの漏出ポイントに過剰なＨＣ量が侵入した後、埋め戻しルーチンは、スピルが処理されると、漏出を処理する。シールグリッドセルは直接隣接するセルを検索し、ポテンシャル値のより低い隣接セルが見つかった場合に、油が移動する。コードは埋め戻しルーチンを終了させ、マイグレーションルーチンを再び始める。
【００７９】
スピルの後、ＨＣは移動を始める。ＨＣは、移動損失としてモデルから出てアクティブでないトラップに移動することができる、又は、以前に集積されたトラップへ移動することもできるが、この場合には併合を生じるであろう。
【００８０】
ＨＣが既に集積されたセルを除き、より低いポテンシャルが存在しないポイントまで移動した場合、ステップ１３５で実行された併合がなされることに留意されたい。集積されたセルは集積段階でタグ付けされていたため、アクティブなトラップ配列を検索することによって、侵入先端とそのプロセッサ番号を知ることができる。その情報が分かれば、過剰なＨＣ量がトラップ侵入先端へ送られる。侵入先端は、別のサブドメインに存在する可能性があるが、アクティブトラップ配列の第２のエントリから特定できる。ＭＰＩは、過剰なＨＣ量をそのプロセッサへ送るために用いられるであろう。
【００８１】
ステップ１８０では、埋め戻し中にスピルが発生したかどうかが特定される。スピル条件が存在すると特定された場合、方法１００ａはステップ１２０へ戻る。そうではなく、スピル条件が存在しないと特定された場合、方法１００ａはステップ１８５へ進む。
【００８２】
ステップ１８５では、更新されたＨＣポテンシャル値が対応するトラップに割り当てられ、その値がグローバルアクティブトラップ配列に保存される。ポテンシャル配列のＨＣポテンシャル値は、効率化の目的で新しいＩＰ手法の最後に限って更新される。新しいＩＰ手法を終了しようとする（ステップ１９０）と、グローバルトラップピークでタグ付けされた集積指標の配列が、アクティブトラップ配列とともに検索され、対応するＨＣの更新されたトラップポテンシャルは、集積されたセルの絶対ＨＣポテンシャル値を置き換える。更新プロセスは、集積体が静的平衡状態にあり、集積内に流れがないというルールに従って行われる。方法１００ａは、ステップ１８５からステップ１９０へ進む。
【００８３】
ステップ１９０で方法１００ａが終了する。実施によっては、方法１００ａが終了する前に、本開示と矛盾しない任意のデータに対して後処理又は他の機能を実行することができる。例えば、終了通知を生成及び送信することができ、終了通知には分析用のレポートデータを含めることができる。
【００８４】
ステップ１５５では、小さいポテンシャル値から大きいポテンシャル値へ順番に、グリッドセルへのＨＣの集積が生ずる。隣接するグリッドセル間でＨＣの移動が発生する。隣接するグリッドセルから、ポテンシャル値が最小のグリッドセルを検索する。このグリッドセルが、侵入すべきグリッドセルとして選択される。
【００８５】
図４Ａに移ると、図４Ａは、一実施に係る、（図１Ａの符号１１５で先に特定された）潜在的なトラップピークが過剰なＨＣ量を含むとしてマーク付けされた最初のステップを示すブロック図４００ａである。図４Ａに関して、潜在的なトラップピークであるグリッドセル４０５ａは、過剰なＨＣ量を含むとしてマーク付けされている。図４Ｂに移ると、図４Ｂは、一実施に係る、過剰なＨＣ量を含むとしてマーク付けされたトラップピークを示すブロック図４００ｂであり、図において、最小のポテンシャル値を有する侵入された隣接セルは集積グリッドセルとなり、その隣接セルはすべて、より大きいポテンシャル値を有する。図４Ｂに関して、トラップピークのグリッドセル４０５ｂは、過剰なＨＣ量を含むとしてマーク付けされ、隣接するグリッドセル４１０ｂは、集積されたグリッドセルとしてマーク付けされる。
【００８６】
隣接するグリッドセル指標は、ポテンシャルリストに加えられているであろう対応するポテンシャルとともに、境界リストに追加される。これらのリストは、集積が行われている際に侵入すべき次善のグリッドセルを検索するときに一緒に用いられる。対応する最小のポテンシャル値を有するグリッドセル指標が、次に侵入すべきグリッドセルとして選択される。
【００８７】
例えば、図４Ａのグリッドセル３が、その特定の近隣で最小のポテンシャルを含むならば、それが次の侵入ターゲットになるように選択され、グリッドセル４０５ａからの過剰なＨＣで充填される。トラップグリッドセルの油量は式（５）に従って設定され、トラップの過剰ＨＣ量はその量だけ減少した。充填されたセルはトラップの一部になり、式（６）を用いてトラップポテンシャルが更新される。
【００８８】
新しいＩＰ手法では、トラップごとに１つのＨＣ更新値が保持される。次に、新たに侵入されたグリッドセルの新しい隣接セルがチェックされ、その新しい隣接セルのポテンシャルが新しく侵入されたセルよりも大きい場合は、埋め戻しが続行され、その新しい隣接セルが境界リストに追加される。
【００８９】
図４Ｂは、グリッドセル３（４０５ｂ）（例えば、図４Ａのグリッドセル３）がグリッドセル４０５ａから（例えば、図４Ａから）侵入されたことを示している。新しいグリッドセルが境界リストに追加される（ここでは、グリッドセル５、６、７）。トラップピークに達した後に埋め戻しされるすべてのグリッドセルは、トラップピークのグローバル指標値でマーク付けされる。このマーキングプロセスは、どの埋め戻しセルがどのトラップに属するかを識別するために用いられる。シリアル処理方式では、プロセスはＨＣ量が使い果たされるまで、又は、スピルポイント若しくは漏出に達するまで続行される。
【００９０】
図１Ａへ戻ると、方法１００ａは、ステップ１５５からステップ１６０へ進む。
【００９１】
ステップ１６０では、ＨＣの集積が別の計算プロセッサで発生しているかどうかについて特定がなされる（ＨＣがサブドメインの外側ハローに侵入したかどうか、又は異なるサブドメインに存在する指標が最小のポテンシャル指標として選択された場合）。並列処理では、トラップの埋め戻しが複数のサブドメインにまたがることがある。
【００９２】
図５に移ると、図５は、一実施に係る、０から３の番号が付された４つの異なるサブドメインにまたがる集積を示すブロック図５００である。斜線で示されたグリッドセルは埋め戻されたグリッドセルを表し、グリッドセル１０１は過剰なＨＣ量を含む侵入先端を表し、灰色の影付きグリッドセルはトラップの境界セルを表す。
【００９３】
図１Ａへ戻ると、ＨＣ集積が別の計算プロセッサで発生していると特定された場合、方法１００ａはステップ１６５へ進む。そうではなく、ＨＣ集積が別の計算プロセッサで発生していないと特定された場合、方法１００ａはステップ１７０へ進む。
【００９４】
ステップ１６５では、並列通信が行われる。外側ハローの指標が選択された場合、その指標は最初に、そのセルが実際に存在するサブドメインの対応する内側ハロー指標値にマッピングされる。通信負荷（ｃｏｍｍｕｎｉｃａｔｉｏｎ  ｏｖｅｒｈｅａｄ）を減らすために、境界リストはすべてのサブドメインでローカルに作成され、そこに存在する。すべてのサブドメインは、ローカルに最も抵抗の少ない境界指標を見つける。グローバルな最小位置が見つかる。もし集積中のサブドメインの外側ハローにあるリストに境界グリッドセルが存在するならば、この方法は外側ハロー指標をフィルタリングし、ローカルな境界リストからそれらのセルを削除する。外側ハローのセルは、そのセルが存在する正しいプロセッサへ送られ、そのサブドメインの対応する内側ハローにトラップ識別子でタグ付けされた境界リストが作成される。新しいＩＰ手法は、境界リストにローカルな指標のみが含まれることを保証する。侵入先端を有するプロセッサは、他のサブドメインに存在するポテンシャル値とともに、ポテンシャルが最も低い単一の境界指標を受け取る。
【００９５】
図６に移ると、図６は、一実施に係る、複数のサブドメインにわたって集積が行われるときに実行される通信を示すブロック図６００である。トラップに関連する境界リストを含む侵入先端を含むサブドメインとは別のサブドメインは、最小抵抗指標を見出し、その最小抵抗指標とそれらの間のポテンシャル値とを、侵入先端を有するサブドメインへ送信する。ＨＣがサブドメインを離れると、過剰なＨＣ量、最小抵抗境界指標、及びその対応するポテンシャル値が、ＨＣが侵入した新しいサブドメインへ送られる。網掛けされたグリッドセル６０５はトラップ識別子であり、これは以前に見出されたトラップピーク指標の対応するグローバル番号である。方法１００ａは、ステップ１６５からステップ１７０へ進む。
【００９６】
ステップ１７０では、トラップが集積境界を共有しているかどうか特定される。埋め戻しの段階では、もし、あるトラップの侵入先端が別のトラップによって埋め戻されるために選択されたならば、併合が発生するであろう。ステップ１７０はこの状態をチェックし、もしそうしたことが起きた場合には、最初のトラップの侵入先端にある過剰なＨＣと第２のトラップの過剰なＨＣ量とが合計され、２つのトラップは併合されることになる。トラップが集積境界を共有していると特定された場合、方法１００ａはステップ１７５へ進む。そうではなく、トラップが集積境界を共有していないと特定された場合、方法１００ａはステップ１５０へ戻る。
【００９７】
ステップ１７５ではトラップが併合される。この併合プロセスは、先に述べた細脈併合とは異なる。埋め戻しはまだ進行中なので、侵入先端は２つのトラップで共有されたポイントになり、２つのトラップが１つのエンティティ（実体）になったときに境界とポテンシャルリストとは併合されなければならない。さらに、各トラップのアクティブトラップピーク配列は併合される必要があり、ここで、その併合されたトラップのＨＣポテンシャルは同じポテンシャル値を持ち、その２つのトラップのグローバル指標が同じ配列に追加される。
【００９８】
表２は、併合が発生したときにアクティブなトラップピーク配列がどのように拡張するかを示す例である。
【表２】
0006966577.tif 000009
【００９９】
集積プロセスは、１つのプロセッサ上で一回に１つのトラップで生じる。ただし、１つのサブドメインは複数のサブドメインから情報を受信できる。これは、１つのサブドメインに２つ以上の隣接するサブドメインがあり、そのサブドメインが拡大する集積をともなう場合に発生する。
【０１００】
図７は、一実施に係る、１又は複数のサブドメインから情報を受信する１つのサブドメインを示すブロック図７００である。図示のように、２つのトラップ（例えば、それぞれトラップ７１０ａ及び７１０ｂ）からの２つの侵入先端（例えば、７０５ａ及び７０５ｂ）は、１つのサブドメイン（ここでは、サブドメイン１）に同時に集積している。侵入先端と埋め戻しされたグリッドセルに関連付けられた色については、凡例７１５を留意されたい。複数のトラップ情報を１つのプロセッサに通信する場合、その通信される情報は特定のトラップに関連付けられる。先に述べたように、境界リストは、第１の指標でトラップ識別子によってタグ付けされるが、その指標はグローバルトラップピーク指標である。過剰なＨＣはアクティブなトラップ配列から特定できる侵入先端指標へ送られる。集積プロセスは、スピル又は漏出が発生するまで、又は、過剰なＨＣがなくなるまで続く。新しいＩＰ手法は、すべてのプロセッサが埋め戻し完了ポイントに達するまで、すべてのプロセッサで埋め戻しを続ける。
【０１０１】
図１Ａへ戻り、方法１００ａは、ステップ１７５からステップ１５０へ戻る。
【０１０２】
図１Ｂに移ると、図１Ｂは、一実施に係る、図１Ａの破線の枠１２１で描いた方法ステップを改めて述べるフローチャート１００ｂを示す。
【０１０３】
ステップ１０２ｂでは、ＨＣの移動が停止するまで侵入が発生する。方法１００ｂは、ステップ１０２ｂからステップ１０４ｂへ進む。
【０１０４】
ステップ１０４ｂでは、ＨＣがトラップ又はサブドメインの境界に達したかどうかがチェックされる。ＨＣがサブドメインの境界に達した場合、方法１００ｂはステップ１０６ｂへ進む。ＨＣがトラップに達した場合、方法１００ｂはステップ１０８ｂへ進む。
【０１０５】
ステップ１０６ｂでは、ＨＣサブドメインの境界に達していれば、ＨＣは近隣のサブドメインへ送られ、フローはステップ１１０ｂへ進む。そうでなければ、フローはステップ１０２ｂへ戻る。
【０１０６】
ステップ１０８ｂでは、ＨＣが別のサブドメインにピークがあるトラップに達すると、通信と併合が発生し、フローはステップ１０２ｂへ戻る。そうでない場合、それ以外のフローはステップ１０２ｂへ戻る。
【０１０７】
ステップ１１０ｂでは、ＨＣが同じサブドメインにピークを有するトラップに達すると、併合が発生し、フローは１０２ｂへ戻る。
【０１０８】
侵入プロセスは、システム内でＨＣの移動がなくなるまで続く。もし、いずれかのグリッドセル内のＨＣの体積がグリッドセルの容積を上回る場合、システムは非平衡状態にあり、集積プロセスが行われる。
【０１０９】
図１Ｃに移ると、図１Ｃは、一実施に係る、図１Ａの枠１５１によって描いた方法のステップを再度説明するフローチャート１００ｃを示す。
【０１１０】
ステップ１０２ｃでは、ＨＣはシステムが平衡状態になるまで埋め戻しを受ける。方法１００ｃは、ステップ１０２Ｃからステップ１０４Ｃへ進む。
【０１１１】
ステップ１０４ｃでは、ＨＣがサブドメインの境界、トラップの境界、又はスピルポイントに達したかどうかのチェックがなされる。ＨＣがサブドメインの境界に達した場合、方法１００ｃはステップ１０６ｃへ進む。ＨＣがトラップの境界に達した場合、方法１００ｃはステップ１０８ｃへ進む。スピルポイントに達した場合、方法１００ｃはステップ１１０ｃへ進む。
【０１１２】
ステップ１０６ｃでは、ＨＣサブドメインの境界に達すると、ＨＣ量、最小ポテンシャル値と指標、及びリストが隣接サブドメインへ送られる。方法１００ｃはステップ１０６ｃからステップ１０２ｃへ戻る。
【０１１３】
ステップ１０８ｃでは、トラップが集積境界を共有している場合、複数のトラップを併合する必要がある。トラップが異なるサブドメインにある場合、トラップ情報の通信が行われる。方法１００ｃは、ステップ１０８ｃからステップ１０２ｃへ戻る。
【０１１４】
ステップ１１０ｃでは、図１Ａ及び図１Ｂで説明したように、終了し、侵入プロセスが始まる。
【０１１５】
テストケース
【０１１６】
ここで述べる新しいＩＰ手法は、数百の岩相と複数の根源岩とを含むモデルを含んだ、さまざまなモデルサイズの複数の堆積盆モデルで（直列及び並列に）試された。試された最大のモデルサイズは、面積グリッドサイズが２５０メートル（ｍ）の１４億個のセルの堆積盆モデルであった。本書には２つのテストケースの結果が含まれている。これらのテスト例を用いて、並列実施の検証、及びこの手法の計算パフォーマンスとスケーラビリティを示す。この方法は、超大規模堆積盆シミュレーションを数時間で実行できるため、地球科学者やエンジニアに、堆積盆モデリングのための高解像度のマイグレーションシミュレーションを実行するための性能の強化を提供する。
【０１１７】
テストケース１
【０１１８】
使用モデルは、ｘ方向に１３５グリッドセル、ｙ方向に１３５グリッドセル、層序が８２層で構成され、モデルサイズは１，４９４，４５０グリッドセルである。面積グリッドサイズは、ｘ方向とｙ方向が共に２キロメートル（ｋｍ）であり、７３，０００平方キロメートル（ｋｍ２）の表面積をカバーする。シミュレーションの実行は、５億５，０００万年前の最初の層序堆積から始まり、現在まで、既知の反応速度モデルを用いる。このモデルの岩相は２つのタイプのみを用いて単純化される。一方はシール層として、他方は残りの層序層に対して用いられる。
【０１１９】
図８に移ると、図８は、一実施に係る、単一の計算コアでテストケース１のモデルを実行した後の結果を示すプロット８００である。図示のように、変換率（ｔｒａｎｓｆｏｒｍａｔｉｏｎ  ｒａｔｉｏ）は、赤（例えば符号８０５）から水色（例えば符号８１０）の色の範囲であり、それぞれ高い所から低い所への変換を示す。変換率は、ケロジェンのどれだけが反応し成熟してＨＣを生成したかを表す。濃い緑色と薄い緑色のグリッドセルは、集積飽和を表す。図示のように、成熟した根源岩は集積体の真下ではなく、モデルの北東にあり、ＨＣは根源岩からシールの下の最小ポテンシャル点に移動したことを示している。移動後に集積プロセスを実行した。集積体を図８に示す。符号８１５が集積体例である。
【０１２０】
図９に移ると、図９は、一実施に係る、複数のトラップピークを含む集積体を示すプロット９００である。集積体９０５内では、トラップピークは黄色のグリッドセルとして示されている（例えば、トラップピーク９１０）。シリアル処理では、この方法は一回に１つのトラップピークを埋め戻し、３つのトラップピークのそれぞれに関連付けられた集積が１つのエンティティに併合されるが、それは併合方法が正しく機能していることを示している。集積体９０５に隣接する２つのグリッドセル９１５ａ／９１５ｂは、ＨＣの集積のない孤立したトラップピークであることに留意されたい。
【０１２１】
図１０に移ると、図１０は、一実施に係る、５０個の計算コアを用いて実行された図８及び図９のテストケースの結果のプロット１０００である。図示のように、さまざまな色の長方形（例えば、長方形１００５）は異なるサブドメインを表す。結果は同一であることが示され（例えば、図９の集積体９０５及び図８の集積体８１５と比較した集積体１０１０を参照）、新しいＩＰ方法の並列化が検証されている。言い換えれば、図１０の集積体１０１０は（図８の集積体８１５の拡大切取図である）図９の集積体９０５と同一である。
【０１２２】
並列処理の併合に関しては、結果もシリアル処理の結果と一致することが示される。例えば、図１１に移ると、図１１は、一実施に係る、異なるサブドメインに存在する併合におけるトラップピークを示すプロット１１００である。図示するように、サブドメインには０から３のラベルが付されている。このケースの場合、並列方式はトラップ１（１１０５）とトラップ２（１１１０）を同時に集積する。これらは別々のサブドメイン（０から３）にあるからである。トラップ２（１１１０）が集積するためにサブドメイン３を離れてサブドメイン１に入ると、トラップ３（１１１５）はサブドメイン３で集積を開始する。異なるサブドメインの異なるトラップからの集積が同時に処理され、シリアル実行と同じ集積の幾何形状（ａｃｃｕｍｕｌａｔｉｏｎ  ｇｅｏｍｅｔｒｉｅｓ）が得られる。これは、並列処理の併合方法も正しく機能していることを示す。
【０１２３】
テストケース２ａ：粗いグリッドモデル
【０１２４】
第２のモデルは、ｘ方向に５３５グリッドセル、ｙ方向に５０５グリッドセルで構成され、層序が８２層で構成され、モデルサイズは２２,１５４,３５０グリッドセルのものである。面積グリッドサイズは、ｘ方向とｙ方向の両方ともに２キロメートル（ｋｍ）で、１，０８０，７００平方キロメートル（ｋｍ２）の表面積をカバーする。テストケース１と同様に、シミュレーションの実行は、５億５，０００万年前での最初の層序堆積から始まり現在まで、同じ反応速度モデルを用いる。繰り返すが、このモデルの岩相は２つのタイプ、すなわち１）一方をシール層として、及び２）他方を残りの層序層に対してのみ用いるように単純化されている。
【０１２５】
図１２に移ると、図１２は、一実施に係る、テストケース２ａのシリアル実行による結果を示すプロットである。濃い緑色のグリッドセル（例えば、符号１２０５及び符号１２１０のグリッドセル）は、集積飽和（ボディ）を示す。赤（例えば、符号１２１５）から青（例えば、符号１２２０）までの範囲の色は、先に述べた変換率を示す。
【０１２６】
ＨＣがモデル境界から流出せずにより長い距離を流れることができる、横方向でより大規模な例を用いて、侵入ルーチンの性能を試すことができる。
【０１２７】
いくつかの集積体は変換率の範囲から遠く離れた領域で見られるが、ここでも、排出されたＨＣがモデル全体におけるポテンシャルの低い領域へ移動していることが示されている。並行して、このテストケースにおいて、ＨＣは更に多くのサブドメインをわたって移動する。
【０１２８】
図１３に移ると、図１３は、一実施に係る、複数のサブドメインにわたって発生し、互いに直接隣接していないサブドメインにまたがるテストケース２ａの併合を示すプロットである。図示のように、表示された黄色のグリッドセル（例えば、グリッドセル１３０５）は、トラップ１３１０のトラップピークを示す。トラップ１３１０は、６つの異なるサブドメイン０、１、４、５、７、及び８と通信するポテンシャルを有する。トラップ１３１０では、サブドメイン１のグリッドセルで埋め戻しが発生した後、侵入すべき次のセルがサブドメイン８に存在する可能性があって、その結果、隣接しないドメイン通信が発生する。このような、より複雑なケースを実行することで、新しいＩＰ手法のロバスト（堅牢）性が更に試される。
【０１２９】
図１４Ａに移ると、図１４Ａは、一実施に係る、単一の計算コアでテストケース２ａモデルを実行した結果を示すプロット１４００ａである。マイグレーションプロセスから生成された集積体（１４０５ａ、１４１０ａ、１４１５ａ等）を表す濃い緑色のパッチに留意されたい。
【０１３０】
図１４Ｂに移ると、図１４Ｂは、一実施に係る、１，０００個の計算コアでテストケース２ａモデルを実行した結果を示すプロット１４００ｂである。図１４Ｂの図示は、新しいＩＰ手法を更に検証するものである。例えば、濃い緑色のパッチ例１４０５ｂ、１４１０ｂ、及び１４１５ｂは、マイグレーション処理で生成された集積体であり、シミュレーションの実行に用いられる計算コアの数に関して不変であるべきである。図示のように、（１，０００個の計算コアを用いて生成される）集積体の例１４０５ｂ、１４１０ｂ、及び１４１５ｂは、（単一の計算コアを用いて生成される）図１４Ａの集積体の例１４０５ａ、１４１０ａ、及び１４１５ａに対応する。
【０１３１】
テストケース２ｂ：細かいグリッドモデル
【０１３２】
本方法を用いて超大規模な高解像度の堆積盆モデルを実行できることを示すために、テストケース２ｂは、ｘ方向に４，２７３グリッドセル、ｙ方向に４，０３３グリッドセルのグリッドサイズで、８２層の層序に細分され、サイズ１，４１３，１０６，７３８グリッドセルのモデルとなった。モデルの他のすべての特性を同じに保ち、計算ノードごとに１４個のプロセッサを用いて１，０００個の計算コアでモデルを実行した。
【０１３３】
図１５に移ると、図１５は、一実施に係る、１，０００個の計算コアを用いてテストケース２ｂモデルを実行した結果を示すプロット１５００である。図に見られるように、示されるトラップ位置は、粗い（ｃｏａｒｓｅ）グリッドのテストケース２ａモデル（図１４Ｂ）と同質である。しかし、テストケース２ａにおいて低解像度の粗いグリッドゆえにアーチファクトであると疑われた、小さく示されたトラップの多くは、テストケース２ｂでは消失していた。実施では、テストケース２ｂの並列処理のシミュレーション処理時間は、テストケース２ｂのモデル（例えば、１０億個のセルを含む）では１３時間という実時間（つまり、壁時間（ｗａｌｌ  ｔｉｍｅ））である。これと比べて、通常は１，０００万から５，０００万セルのモデルサイズを用いる従来のシリアルＩＰ方式にとって、テストケース２ｂのモデルは大きすぎる。テストケース２ｂのモデルで従来のシリアルＩＰ手法を用いようとすると、並列処理のシミュレーション処理を用いたテストケース２ｂ処理の何倍もの（例えば、５から１０倍）の壁時間がかかってしまう。
【０１３４】
テストケース１及びテストケース２ａ／２ｂのスケーラビリティ
【０１３５】
報告されたスケーラビリティの結果は、ＱＬＯＧＩＣ  ＩＮＦＩＮＩＢＡＮＤ  ｉｎｔｅｒｃｏｎｅｃｔ（ＱＬｏｇｉｃ  Ｃｏｒｐｏｒａｔｉｏｎ，アリソビエホ，カリフォルニア州，米国）に接続されたＬＩＮＵＸ  ＰＣクラスタを用いたシミュレーション実行のためのものである。クラスタ内の各計算ノードは、２つのＸＥＯＮ  Ｅ５−２６８０ＣＰＵ（Ｉｎｔｅｌ  Ｃｏｒｐｏｒａｔｉｏｎ，サンタクララ，カリフォルニア州，米国）で構成され、２．８０ＧＨｚのクロック速度で実行され、各ＣＰＵには１０個の計算コアがある。両テストケースとも、計算ノードあたり１６個のプロセッサを用いて実行された。
【０１３６】
図１６に移ると、図１６は、一実施に係る、テストケース１についてプロットされた計算コアの数に対する正規化された実行時間のグラフ１６００である。同様に、図１７に移ると、図１７は、一実施に係る、テストケース２ａ／２ｂについてプロットされた計算コアの数に対する正規化された実行時間のグラフ１７００である。図１６及び図１７に示す結果は、並列化された新しいＩＰ手法が実際に非常にスケーラブルであることを示している。グラフ１６００及びグラフ１７００の理想ライン１６０２及び１７０２は、それぞれ完全なスケーラビリティを表していることに留意されたい。つまり、理想的には、１コアに対して１０コアのコンピュータを用いる場合、並列化によるオーバーヘッドがなく、コンピュータに処理のボトルネックがないとしたならば、コードの実行速度は１０倍に及ぶと予想される。実際には、実際の高速化因子は理想よりも小さくなるのが普通である。グラフ１６００及びグラフ１７００は、記載されたアプローチが約５７／６４（グラフ１６００）及び図１７では約５１／６４の高速化因子に達していることを示している。これは、実際の場合のこの方法の並列処理の効率を示しているが、達成されるスケーラビリティは、実際のモデルにおける過剰なＨＣの分布、及び結果として得られる集積体にも依存する。
【０１３７】
図１８は、一実施に係る、本願で説明されるような記述されたアルゴリズム、方法、機能、プロセス、フロー、及び手順に関連付けられた計算機能を提供するために使用された例示的なコンピュータシステム１８００を示すブロック図である。図示されたコンピュータ１８０２は、サーバ、デスクトップコンピュータ、ラップトップ／ノートブックコンピュータ、無線データポート、スマートフォン、パーソナルデータアシスタント（ＰＤＡ）、タブレットコンピューティングデバイス、これらのデバイスのうちの１つ若しくは複数のプロセッサ、又は任意の他の適切な処理デバイスといった任意のコンピューティングデバイスを包含することを意図しており、該コンピューティングデバイスの物理インスタンス又は仮想インスタンス（又はその両方）を含む。加えて、コンピュータ１８０２は、キーパッド、キーボード、タッチスクリーン、又はユーザ情報を受け入れ可能な他のデバイスといった入力デバイスと、コンピュータ１８０２の動作に関連付けられた情報を伝達する出力デバイスとを含むコンピュータを包含していてもよく、デジタルデータ、視覚情報、または音声情報（若しくは情報の組み合わせ）、又はグラフィカルユーザインタフェース（ＧＵＩ）を含む。
【０１３８】
コンピュータ１８０２は、クライアント、ネットワークコンポーネント、サーバ、データベースあるいは他の持続性のもの、又は本願に説明された主題を実行するためのコンピュータシステムの任意の他のコンポーネント（又は役割の組み合わせ）としての役割を果たすことができる。例示のコンピュータ１８０２は、ネットワーク１８３０と通信可能に結合される。いくつかの実施では、コンピュータ１８０２の１又は複数のコンポーネントは、クラウドコンピューティングベースの、ローカルな、グローバルな、又は他の環境（又は環境の組み合わせ）を含む環境内で動作するように構成できる。
【０１３９】
高レベルでは、コンピュータ１８０２は、説明された主題に関連付けられたデータ及び情報を受け、送出し、処理し、保存し、又は管理するように動作可能な電子的計算デバイスである。いくつかの実施によれば、コンピュータ１８０２は、アプリケーションサーバ、電子メールサーバ、ウェブサーバ、キャッシングサーバ、ストリーミングデータサーバ、又は他のサーバ（又はサーバの組み合わせ）を含んでいても、或いは通信可能に結合されていてもよい。
【０１４０】
コンピュータ１８０２は、（例えば、別のコンピュータ６０２上で動作する）クライアントアプリケーションからネットワーク１８３０を介してリクエストを受けることができ、受信したリクエストを適切なソフトウェアアプリケーションを用いて処理することによって受信したリクエストに応答することができる。加えて、リクエストは、コンピュータ１８０２へ、（例えば、コマンドコンソールから又は他の適切なアクセス方法によって）内部ユーザ、外部ユーザ又はサードパーティ、他の自動化されたアプリケーション、更には、その他の適切なエンティティ、個人、システム、又はコンピュータから送信されるようにしてもよい。
【０１４１】
コンピュータ１８０２の各コンポーネントは、システムバス１８０３を用いて通信できる。いくつかの実施では、コンピュータ１８０２のコンポーネントのいずれか又は全ては、ハードウェア又はソフトウェアの両方（又はハードウェア及びソフトウェアの組み合わせ）が、アプリケーションプログラミングインタフェース（ＡＰＩ）１８１２又はサービスレイヤ１８１３（又は、ＡＰＩ１８１２及びサービスレイヤ１８１３の組み合わせ）を用いて、システムバス１８０３を介して、互いに又はインタフェース１８０４（又は両方の組み合わせ）とインタフェース接続することができる。ＡＰＩ１８１２は、ルーチン、データ構造、及びオブジェクトクラスの仕様を含むことができる。ＡＰＩ１８１２は、コンピュータ言語に依存しないか又は依存するかのいずれかであり、また完全なインタフェース、単一の関数、又は一組のＡＰＩを指していてもよい。サービスレイヤ１８１３は、コンピュータ１８０２へソフトウェアサービスを、又はコンピュータ１８０２に通信可能に結合された他のコンポーネントを（図示されているかどうかに関わらず）提供する。コンピュータ１８０２の機能は、このサービスレイヤを用いて全てのサービスコンシューマにアクセス可能であってもよい。サービスレイヤ１８１３によって提供されたサービスといったソフトウェアサービスは、定義済みのインタフェースを介して再利用可能な定義済み機能を提供する。例えば、インタフェースは、ＪＡＶＡ（登録商標）、Ｃ＋＋、又は、拡張マークアップ言語（ＸＭＬ）フォーマット若しくは他の適切なフォーマットでデータを提供する他の適切な言語で書かれたソフトウェアであり得る。コンピュータ１８０２の統合コンポーネントとして図示される一方で、代替の実施は、ＡＰＩ１８１２又はサービスレイヤ１８１３を、コンピュータ１８０２の他のコンポーネント又はコンピュータ１８０２に通信可能に結合された他のコンポーネントに（図示されているかどうかに関わらず）関連してスタンドアロンコンポーネントとして例示することができる。さらには、ＡＰＩ１８１２又はサービスレイヤ１８１３のいずれか又は全ての部分は、本願の範囲から逸脱することなく、別のソフトウェアモジュール、エンタープライズアプリケーション、又はハードウェアモジュールの子モジュール又はサブモジュールとして実施されるようにしてもよい。
【０１４２】
コンピュータ１８０２は、インタフェース１８０４を含む。図１８には単一のインタフェース１８０４として図示されているが、２つ又はそれより多くのインタフェース１８０４を、コンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従って使用することができる。インタフェース１８０４は、分散環境においてネットワーク１８３０に（図示されているかどうかに関わらず）接続される他のシステムと通信するために、コンピュータ１８０２によって使用される。一般にインタフェース１８０４は、ソフトウェア又はハードウェア（又はソフトウェア及びハードウェアの組み合わせ）でエンコードされたロジックを含み、またネットワーク１８３０と通信するように動作可能である。より具体的には、インタフェース１８０４は、ネットワーク１８３０又はインタフェースのハードウェアが図示のコンピュータ１８０２の内側及び外側において物理信号を通信するように動作可能であるように、通信に関連付けられた１又は複数の通信プロトコルをサポートするソフトウェアを含むことができる。
【０１４３】
コンピュータ１８０２は、プロセッサ１８０５を含む。図１８においては単一のプロセッサ１８０５として示されるが、２又はそれより多くのプロセッサを、コンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従って使用することができる。一般に、プロセッサ１８０５は、命令を実行し、また、本願に記載されるように、コンピュータ１８０２の動作、及び任意のアルゴリズム、方法、機能、プロセス、フロー、若しくは手順を実行するためにデータを処理する。
【０１４４】
また、コンピュータ１８０２は、コンピュータ１８０２又はネットワーク１８３０に（図示されているかどうかに関わらず）接続可能な他のコンポーネント（又はこの両方の組み合わせ）のためのデータを保持することができるデータベース１８０６を含む。例えば、データベース１８０６は、インメモリの、従来型の、又は他のタイプのデータベースであることができ、該データベースは、本願に整合するデータを格納する。いくつかの実施では、データベース１８０６は、説明された機能及びコンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従う、２以上の異なるデータベースタイプの組み合わせ（例えば、インメモリと従来型とのハイブリッド型データベース）であり得る。図１８において単一のデータベース１８０６として示されるが、（同じタイプ又はタイプの組み合わせの）２以上のデータベースは、説明された機能及びコンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従って、使用できる。データベース１８０６は、コンピュータ１８０２の不可欠な構成要素として示されるが、代替の実施では、データベース１８０６は、コンピュータ１８０２の外部にあってもよい。
【０１４５】
また、コンピュータ１８０２は、ネットワーク１８３０に（図示されているかどうかに関わらず）接続できるコンピュータ１８０２又は他のコンポーネント（又は両方の組み合わせ）のためのデータを保持可能なメモリ１８０７を含む。例えば、メモリ１８０７は、ランダムアクセスメモリ（ＲＡＭ）、読み取り専用メモリ（ＲＯＭ）、光学的なもの、磁気的なもの等であることができ、本願に整合するデータを格納する。いくつかの実施では、メモリ１８０７は、説明された機能及びコンピュータ１８０２の特定のニーズ、要望又は特定の実施に従う、２又はより多くの異なるタイプのメモリの組み合わせ（例えば、ＲＡＭと磁気記憶装置との組み合わせ）であることができる。図１８には単一のメモリ１８０７として示されるが、（同じタイプ又はタイプの組み合わせの）２以上のメモリ１８０７が、説明された機能及びコンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従って使用できる。メモリ１８０７は、コンピュータ１８０２の不可欠な構成要素として示されるが、代替の実施では、メモリ１８０７はコンピュータ１８０２の外部にあってもよい。
【０１４６】
アプリケーション１８０８は、コンピュータ１８０２の特定のニーズ、要望、又は特定の実施に従う機能を、特に本願において説明された機能に関して提供するアルゴリズムソフトウェアエンジンである。例えば、アプリケーション１８０８は、１又は複数のコンポーネント、モジュール、又はアプリケーションとして働くことができる。さらに、単一のアプリケーション１８０８として示されているが、アプリケーション１８０８は、コンピュータ１８０２上において複数のアプリケーション１８０８として実施されていてもよい。加えて、コンピュータ１８０２と一体として図示されているが、代替の実施では、アプリケーション１８０８は、コンピュータ１８０２の外部にあってもよい。
【０１４７】
また、コンピュータ１８０２は、電源１８１４を含むことができる。電源１８１４は、ユーザによる交換可能又はユーザによる交換不可能に構成できる充電式又は非充電式バッテリを含むことができる。いくつかの実施では、電源１８１４は、（再充電、スタンバイ、又は別の電力管理の機能を含む）電力変換又は管理回路を含むことができる。いくつかの実施では、電源１８１４は、コンピュータ１８０２を壁のソケットや別の電源に接続でき電源プラグを含み、これにより例えばコンピュータ１８０２に電力を供給したり充電式バッテリを再充電したりすることができる。
【０１４８】
任意の数のコンピュータ１８０２が、このコンピュータ１８０２を含むコンピュータシステムに関連付けられ、又は該コンピュータシステムの外部にあり、各コンピュータ１８０２はネットワーク１８３０を介して通信する。さらに、用語「クライアント」、「ユーザ」、又は他の適切な用語は、本願の範囲から逸脱することなく、必要に応じて、交換可能に使用される。さらに、本願では、多くのユーザが１台のコンピュータ１８０２を使用できることや、１人のユーザが複数のコンピュータ１８０２を使用できることを想定している。
【０１４９】
いくつかの実施では、既述された方法論は、コンピュータにより実施されるコントローラ、データベース、又はコンピュータにより実施される他のシステムに、メッセージ、命令、又は他の通信を送って、コンピュータにより実施される操作又は他の機能の操作をコンピュータにより実施されるシステムの制御を動的に開始し、コンピュータにより実施される操作又は他の機能の操作をコンピュータにより実施されるシステムを制御し、又はコンピュータにより実施される操作又は他の機能／操作を別のコンピュータにより実施されるシステムに実行させるように構成することができる。例えば、データ、操作、出力、又はＧＵＩとのインタラクション（対話、相互作用）に基づく操作が送信されて、コンピュータ、データベース、ネットワーク、又はコンピュータに基づく他のシステムに関連付けられた操作に、記憶効率の操作、データ検索の操作、又はこの開示に整合する他の操作を実行させることができる。別の例では、図示されたＧＵＩ（例えば、図８−１３、１４Ａ−１４Ｂ、１５−１７）とインタラクションすると、自動的に、一又は複数の命令がＧＵＩから送信されて、データのためのリクエスト、データの保存、データの分析、又は本願に整合する他の操作がトリガされるようになる。
【０１５０】
場合によっては、送信された命令により、計算機器又は他の機器という現実世界の有形品に関する制御、動作、修正、強化、又は他の動作が生じ得る。例えば、説示されたＧＵＩは、コンピュータデータベースの磁気／光学ディスクドライブの速度を遅く若しくは速くしたり、コンピューティングシステムを起動若しくは停止したり、ネットワーク接続にわたって許可されるデータ帯域幅をネットワークインタフェースデバイスに無効化させ、調整させ、若しくは増加させたり、あるいは、記述された方法論に関連付けられたコンピューティングシステムに関する若しくは記述された方法論に関連付けられた該コンピューティングシステムとのインタラクションの結果、動作、特定、若しくは分析の通知として（機械的アラーム／発光デバイスといった）可聴／可視アラームを発したりする。
【０１５１】
実施によっては、先の方法論の出力を使って、炭化水素の生産、分析、及び回収に関連する、又は本開示に沿った他の目的のための、具体的な実体面での機器に、動的に影響を及ぼす、指示する、制御する又は、それらを管理することができる。例えば、進行中の掘削作業から受信したリアルタイムデータを、先の方法論を用いてなされる分析に入れ込むことができる。前述した並列処理のインバージョンパーコレーション（ＩＰ）マイグレーション方法の出力は、さまざまな目的に使える。例えば、先の方法論で生成された結果に応じて、掘削軌道を変更したり、ドリル速度を加減したり、ドリルを停止したり、（視覚、聴覚、又は音声アラームで）アラームを有効／無効にしたりして、（例えば、停止、再起動、加速、又は減少させるように）製油所やポンプの操作に影響を及ぼすことができる。他の例としては、ＨＣトラップの検出に基づいて、ジオステアリング及び掘進方向スタッフへの警告（視覚、聴覚、音声アラームなど）を含めることができる。実施によっては、先の方法論は、本開示と整合する炭化水素関連の、又はその他の有形である実体面での機器を制御し、影響を及ぼし、又は使用するための動的なコンピュータにより実施される制御システムの一部として統合できる。
【０１５２】
本主題で説明する実施は、１つ以上の特徴を単独で又は組み合わせて含むことができる。
【０１５３】
例えば、第１の実施では、地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理方法であって：１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩から、ＨＣが排出されたことを特定するステップと；潜在的なトラップピークを前記複数のグリッドセル内で識別するステップと；前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行するステップと；ＨＣ類を含むグリッドセルが過剰量のＨＣ類を含んでいるかどうかを特定するステップと；前記識別された潜在的なトラップピークに関連付けられたトラップにおける、グリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行するステップと；前記トラップをＨＣのポテンシャル値で更新するステップと；を含む。
【０１５４】
前述した及び他に説明された実施は、オプションとして、以下に示す機能の１又は複数を含むことができる。
【０１５５】
以下の機能のいずれかと組み合わせ可能な第１の機能として、前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする。
【０１５６】
前述の又は以下の機能のいずれかと組み合わせ可能な第２の機能として、スピルポイントに達したことを特定するステップと；戻って前記侵入プロセスを実行するステップと；を更に含む。
【０１５７】
前述の又は以下の機能のいずれかと組み合わせ可能な第３の機能として、前記侵入プロセスは：前記ＨＣがサブドメイン境界に達したことを特定するステップと；前記ＨＣを、隣接するサブドメインに移動するステップと；前記ＨＣが同じサブドメインで充填トラップに達した場合、前記ＨＣを前記充填トラップと併合するステップと；を含む。
【０１５８】
前述の又は以下の機能のいずれかと組み合わせ可能な第４の機能として、前記サブドメイン境界においてＨＣ量を通信するステップと；サブドメイン境界を越えた前記ＨＣを、別のサブドメインに関連付けられたトラップと併合するステップと；を更に含む。
【０１５９】
前述の又は以下の機能のいずれかと組み合わせ可能な第５の機能として、前記ＨＣがサブドメイン境界に達していることを特定するステップと；前記サブドメイン境界に関連付けられた隣接するサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信するステップと；を更に含む。
【０１６０】
前述の又は以下の機能のいずれかと組み合わせ可能な第６の機能として、前記ＨＣがトラップ境界に達したことを特定するステップと；複数のトラップが集積境界を共有するかどうかを特定するステップと；を更に含む。
【０１６１】
第２の実施では、地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理の方法の１又は複数の命令を格納する非一時的なコンピュータ読取り可能媒体であって、前記１又は複数の命令は所定の操作を実行するためにコンピュータシステムによって実行可能であり、前記操作は：１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作と；潜在的なトラップピークを前記複数のグリッドセル内で識別する操作と；前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作と；ＨＣ類を含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作と；前記トラップをＨＣのポテンシャル値で更新する操作と；を含む。
【０１６２】
前述した及び他に説明された実施は、オプションとして、以下に示す機能の１又は複数を含むことができる。
【０１６３】
以下の機能のいずれかと組み合わせ可能な第１の機能として、前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする。
【０１６４】
前述の又は以下の機能のいずれかと組み合わせ可能な第２の機能として、スピルポイントに達していることを特定する操作と；戻って前記侵入プロセスを実行する操作と；を実行する１又は複数の命令を更に含む。
【０１６５】
前述の又は以下の機能のいずれかと組み合わせ可能な第３の機能として、前記侵入プロセスは：前記ＨＣがサブドメイン境界に達したことを特定する操作と；前記ＨＣを隣接するサブドメインに移動する操作と；前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する操作と；を実行する１又は複数の命令を含む。
【０１６６】
前述の又は以下の機能のいずれかと組み合わせ可能な第４の機能として、前記サブドメイン境界でＨＣ量を通信する操作と；サブドメインの境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する操作と；を実行する１又は複数の命令を更に含む。
【０１６７】
前述の又は以下の機能のいずれかと組み合わせ可能な第５の機能として、前記ＨＣがサブドメイン境界に達したことを特定する操作と；前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信する操作と；を実行する１又は複数の命令を更に含む。
【０１６８】
前述の又は以下の機能のいずれかと組み合わせ可能な第６の機能として、前記ＨＣがサブドメイン境界に達したことを特定する操作と；複数のトラップが集積境界を共有するかどうかを特定する操作と；を実行する１又は複数の命令を更に含む。
【０１６９】
第３の実施では、地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理システムであって：コンピュータメモリと；前記コンピュータメモリと相互作動可能に結合され、所定の操作を実行するように構成されたハードウェアプロセッサとであって、前記操作は：１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作と；潜在的なトラップピークを前記複数のグリッドセル内で識別する操作と；前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作と；ＨＣを含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作と；前記トラップをＨＣポテンシャル値で更新する操作と；を備える、前記ハードウェアプロセッサと；を含む。
【０１７０】
前述した及び他に説明された実施は、オプションとして、以下に示す機能の１又は複数を含むことができる。
【０１７１】
以下の機能のいずれかと組み合わせ可能な第１の機能として、前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする。
【０１７２】
前述の又は以下の機能のいずれかと組み合わせ可能な第２の機能として、スピルポイントに達したことを特定し；戻って前記侵入プロセスを実行する；ように更に構成される。
【０１７３】
前述の又は以下の機能のいずれかと組み合わせ可能な第３の機能として、前記侵入プロセスは：前記ＨＣがサブドメイン境界に達したことを特定し；前記ＨＣを、隣接するサブドメインに移動し；前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する；ように更に構成される。
【０１７４】
前述の又は以下の機能のいずれかと組み合わせ可能な第４の機能として、前記サブドメイン境界でＨＣ量を通信し；サブドメイン境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する；ように更に構成される。
【０１７５】
前述の又は以下の機能のいずれかと組み合わせ可能な第５の機能として、前記ＨＣがサブドメイン境界に達したことを特定し；前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信するように更に構成される。
【０１７６】
前述の又は以下の特徴のいずれかと組み合わせ可能な第６の特徴として、ＨＣがトラップ境界に達していることを特定し；複数のトラップが集積境界を共有しているかどうかを特定する；ようにさらに構成される。
【０１７７】
コンピュータにより実施される並列処理方法では、地下媒体内の炭化水素（ＨＣ）のマイグレーション（移動）シミュレーションは、堆積盆シミュレータで実行される。堆積盆ドメインは複数のサブドメインに分割され、各サブドメインは堆積盆体積のサブセットを表し、ＨＣ流体に関するデータを有するグリッドセルの集合体に編成される。堆積盆のシミュレーションは、分割された堆積盆のサブドメインが割り当てられている少なくとも１つのマスタノードで構成されるクラスタコンピュータで実行される。グリッドセルの対象となる流体特性パラメータがシミュレートされる。このコンピュータにより実施される方法は：炭化水素（ＨＣ）類が、１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩から排出されたことを特定するステップと；複数のグリッドセル内の潜在的なトラップピークを特定するステップと；ＨＣが複数のグリッドセル内での移動を停止するまで、侵入プロセスを実行するステップと；ＨＣ類を含むグリッドセルに過剰量のＨＣ類が含まれているかどうかを特定するステップと；集積プロセスを実行して潜在的な識別されたトラップピークに関連付けられたトラップにおいてグリッドセルにＨＣを充填することをモデル化するステップと；トラップをＨＣポテンシャル値で更新するステップと；を含む。
【０１７８】
本明細書で記述された主題及び機能的な操作の実施は、デジタル電子回路、有形に具体化されたコンピュータのソフトウェア又はファームウェア、コンピュータハードウェアに、又はそれらの一又は複数の組み合わせで実施でき、これらデジタル電子回路、ソフトウェア、ファームウェア、及びコンピュータハードウェアは、本明細書に開示される構造及びそれらの構造的な同等物を含む。記述された主題におけるソフトウェアの実施は、一又は複数のコンピュータプログラム、つまり、コンピュータプログラム命令の一又は複数のモジュールとして実施されることができ、コンピュータプログラム命令は、コンピュータ若しくはコンピュータにより実施されるシステムによる実行のために、有形の非一時的なコンピュータ読取可能媒体上にエンコードされ、又はコンピュータ若しくはコンピュータにより実施されるシステムの動作を制御する。代替的に又は追加的に、プログラム命令は、人工的に生成された伝播信号、例えば、マシン生成の電気的、光学的、又は電磁気的な信号にエンコードされ、この信号は、コンピュータ又はコンピュータにより実施されるシステムによる実行のために受振器装置への送信用の情報をエンコードするために生成される。コンピュータ記憶媒体は、機械により読取可能な記憶装置、機械により読取可能な記憶基板、ランダム若しくはシリアルアクセスのメモリデバイス、又はコンピュータ記憶媒体の組み合わせであり得る。一又は複数のコンピュータを構成することは、一又は複数のコンピュータがハードウェア、ファームウェア、又はソフトウェア（又はハードウェア、ファームウェア、及びソフトウェアの組み合わせ）をインストールし、一又は複数のコンピュータがソフトウェアを実行するときに特定の計算操作が実行されることを意味する。
【０１７９】
「リアル−タイム」、「リアルタイム」、「リアル（ファースト）タイム（ＲＦＴ）」、「ほぼリアルタイム（ＮＲＴ）」、「準リアルタイム」、又は（当業者が理解する）類似の用語は、個々が、実質的に同時に動作及び応答が生じることを知覚するように、該動作及び応答が時間的に近接していることを意味する。例えば、データにアクセスするための個々の動作に続くデータの表示（又は表示の開始）に対する応答の時間差は、１ミリ秒（ｍｓ）未満、１秒（ｓ）未満、又は５ｓ未満であり得る。要求されたデータが、即時に表示（又は表示開始）される必要はないが、記述されたコンピューティングシステムの処理制限と、データを例えば収集し、正確に評価し、分析し、処理し、格納し、又は送信するために必要な時間とを考慮して、意図的な遅延なしに表示（又は表示開始）される。
【０１８０】
用語「データ処理装置」、「コンピュータ」、又は「電子コンピュータデバイス」（又は当業者によって理解される同等の用語）は、データ処理ハードウェアを指し、あらゆる種類の装置、デバイス、及びマシンを包含し、これらの種類は、データを処理するためのものであり、例示として、プログラマブルプロセッサ、コンピュータ、又は複数のプロセッサ若しくはコンピュータを含む。コンピュータは、例えば、中央処理装置（ＣＰＵ）、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、又は特定用途向け集積回路（ＡＳＩＣ）を包含する専用ロジック回路であるか、又はこれらを更に含むことができる。いくつかの実施では、コンピュータ若しくはコンピュータにより実施されるシステム、又は専用ロジック回路（又はコンピュータ若しくはコンピュータにより実施されるシステムと専用ロジック回路との組み合わせ）は、ハードウェア系又はソフトウェア系（又はソフトウェア系及びハードウェア系の両方の基づく組み合わせ）であることができる。コンピュータは、オプションとして、コンピュータプログラムの実行環境を生成するコード、例えば、プロセッサファームウェア、プロトコルスタック、データベース管理システム、オペレーティングシステム、又は実行環境の組み合わせを構成するコードを含むことができる。本願は、あるタイプのオペレーティングシステムを持つコンピュータ又はコンピュータにより実施されるシステムの使用を想定しており、オペレーティングシステムは、例えば、ＬＩＮＵＸ、ＵＮＩＸ（登録商標）、ＷＩＮＤＯＷＳ（登録商標）、ＭＡＣＯＳ、ＡＮＤＲＯＩＤ（登録商標）、又はＩＯＳである。
【０１８１】
コンピュータプログラムは、プログラム、ソフトウェア、ソフトウェアアプリケーション、ユニット、モジュール、ソフトウェアモジュール、スクリプト、コード、又は他のコンポーネントとして参照され又は記述されるものであって、コンパイル言語若しくはインタープリタ言語を含むプログラミング言語、又は宣言型若しくは手続き型の言語の形式で記載でき、またコンピュータプログラムは、コンピューティング環境で使用するための任意の形式で展開でき、例えばスタンドアロンプログラム、モジュール、コンポーネント、又はサブルーチンとして含む。コンピュータプログラムは、ファイルシステム内のファイルに対応できるが、必ずしも対応する必要はない。プログラムは、他のプログラム又はデータを保持するファイルの一部、例えばマークアップ言語ドキュメントに格納された一又は複数のスクリプトに保存でき、該他のプログラムは、問題のプログラム専用の単一ファイル内に、或いは複数の連携したファイル、例えば一又は複数のモジュール、サブプログラム、又はコードの一部に格納できる。コンピュータプログラムは、一又は複数のコンピュータ上に展開でき、該コンピュータは、一のサイトに位置し、又は複数のサイトに分散されて配置され、これらは通信ネットワークによって相互接続される。
【０１８２】
様々な図面に例示されたプログラムの部分は、ユニット又はモジュールといった個々のコンポーネントとして例示でき、個々のコンポーネントは、様々なオブジェクト、方法、又は他のプロセスを用いて、説明された特徴及び機能を実施しているが、プログラムは、それらの代わりに、いくつかのサブユニット、サブモジュール、サードパーティのサービス、コンポーネント、ライブラリ、及び他のコンポーネントを、必要に応じて含むことができる。逆に、必要に応じて、様々なコンポーネントの機能及び特徴は、単一のコンポーネントに組み合わせることができる。計算上の特定を為すために使用されたしきい値は、静的に、動的に、又は静的と動的との両方で特定できる。
【０１８３】
説明された方法、プロセス、又は論理フローは、本願と整合する機能の一又は複数の例示を表したものであって、説示された又は図示された実施に本願を限定することを意図したものではなく、説示された原理及び特徴に整合する最も広い範囲が与えられるべきである。説示された方法、プロセス、又は論理フローは、一又は複数のプログラム可能なコンピュータによって実行でき、該コンピュータは、入力データを操作して出力データを生成することによって機能を実行する一又は複数のコンピュータプログラムを実行する。方法、プロセス、又は論理フローは、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実行でき、またコンピュータも、特定用途のロジック回路、例えばＣＰＵ、ＦＰＧＡ、又はＡＳＩＣとして実施できる。
【０１８４】
コンピュータプログラムの実行のためのコンピュータは、汎用又は特定用途のマイクロプロセッサ、これらの両方、又は別のタイプのＣＰＵに基づくことができる。一般的には、ＣＰＵは、メモリから命令及びデータを受け取り、メモリに書き込む。コンピュータの重要な要素は、命令を行い又は実行するためのＣＰＵ、及び命令及びデータを保存するための一又は複数のメモリデバイスである。一般に、コンピュータは、データを保存するための一又は複数の大容量記憶装置、例えば磁気、光磁気ディスク、又は光ディスクを含み、或いはこれらの大容量記憶装置に動作可能に結合され、データを受信し、転送し、又はこの両方を行う。しかし、コンピュータは、そのようなデバイスを持つ必要はない。さらには、コンピュータは、別のデバイス、例えば携帯電話、携帯情報端末（ＰＤＡ）、モバイルオーディオ若しくはビデオプレーヤー、ゲームコンソール、全地球測位システム（ＧＰＳ）受信機、又はポータブルメモリストレージデバイスに組み込まれることができる。
【０１８５】
コンピュータプログラム命令及びデータを格納するための非一時的なコンピュータ読取可能媒体は、例示として半導体メモリデバイス、磁気デバイス、及び光学メモリデバイスを含む、あらゆる形態の永続的／非永続的又は揮発性／不揮発性のメモリ、媒体、及びメモリデバイスを含むことができる。あらゆる形態のメモリデバイスは：半導体メモリデバイス、例えばランダムアクセスメモリ（ＲＡＭ）、読み取り専用メモリ（ＲＯＭ）、相変化メモリ（ＰＲＡＭ）、スタティックランダムアクセスメモリ（ＳＲＡＭ）、ダイナミックランダムアクセスメモリ（ＤＲＡＭ）、消去可能プログラマブルリードオンリーメモリ（ＥＰＲＯＭ）、電気的消去可能プログラマブル読み取り専用メモリ（ＥＥＰＲＯＭ）、及びフラッシュメモリデバイスと；磁気デバイス、例えば、テープ、カートリッジ、カセット、内部／リムーバブルディスクと；光磁気ディスクと；光学メモリデバイス、例えば、デジタルバーサタイル／ビデオディスク（ＤＶＤ）、コンパクトディスク（ＣＤ）ＲＯＭ、ＤＶＤ＋／−Ｒ、ＤＶＤ−ＲＡＭ、ＤＶＤ−ＲＯＭ、高精細／密度（ＨＤ）−ＤＶＤ、及びＢＬＵ−ＲＡＹ／ＢＬＵ−ＲＡＹＤＩＳＣ（ＢＤ）、及び他の光学メモリテクノロジーと；を含む。メモリは、様々なオブジェクト又はデータを格納でき、オブジェクト又はデータは、キャッシュ、クラス、フレームワーク、アプリケーション、モジュール、バックアップデータ、ジョブ、Ｗｅｂページ、Ｗｅｂページテンプレート、データ構造、データベーステーブル、動的情報を格納するリポジトリ、又は他の適切な情報を包含しており、他の適切な情報は、任意のパラメータ、任意の変数、任意のアルゴリズム、任意の命令、任意のルール、任意の制約、又は任意の参照を含む。さらに、メモリは、ログ、ポリシー、セキュリティ又はアクセスのデータ、レポートファイルといった他の適切なデータを含むことができる。プロセッサ及びメモリは、特定用途ロジック回路によって補完され又は組み込むことができる。
【０１８６】
ユーザとのインタラクションを提供するために、本明細書で記載された主題の実施は、ユーザに情報を表示するためのディスプレイデバイス、例えば陰極線管（ＣＲＴ）、液晶ディスプレイ（ＬＣＤ）、発光ダイオード（ＬＥＤ）又はプラズマモニターと、ユーザがコンピュータに入力を提供するキーボード及びポインティングデバイス、例えばマウス、トラックボール、トラックパッドとを有するコンピュータ上で実施できる。入力が、タブレットコンピュータの感圧性表面、静電容量式若しくは電気式のセンシングを使用したマルチタッチスクリーンといった、タッチスクリーンを用いてコンピュータに提供することもできる。他のタイプのデバイスが、ユーザとインタラクションするために使用できる。例えば、ユーザに提供されるフィードバックは、（視覚、聴覚、触覚、又は複数のフィードバックタイプの組み合わせといった）あらゆる形態の感覚のフィードバックであることができる。ユーザからの入力は、音響入力、音声入力、触覚入力を包含する、あらゆる形式で受信できる。さらに、コンピュータは、ユーザによって使用されるクライアントコンピューティングデバイスからのドキュメントを受信し及び該ドキュメントをクライアントコンピューティングデバイスに送信することによって（例えば、Ｗｅｂブラウザから受けたリクエストに応答してユーザのモバイルコンピューティングデバイス上においてＷｅｂブラウザにＷｅｂページを送ることにより）ユーザとインタラクションできる。
【０１８７】
用語「グラフィカルユーザインタフェース」又は「ＧＵＩ」は、単数又は複数で使用して、一又は複数のグラフィカルユーザインタフェース及び特定のグラフィカルユーザインタフェースのディスプレイの各々を説明できる。したがって、ＧＵＩは、Ｗｅｂブラウザ、タッチスクリーン、又はコマンドラインインタフェース（ＣＬＩ）を含むが、これらに限定されることなく、任意のグラフィカルユーザインタフェースを提示することができ、コマンドラインインタフェースは、情報を処理すると共に情報の結果をユーザに効率的に提示する。一般的に、ＧＵＩは、いくつかユーザインタフェース（ＵＩ）要素、例えばインタラクティブフィールド、プルダウンリスト、ボタンを含むことができ、これらのうちのいくつか又は全ては、Ｗｅｂブラウザに関連付けられる。これら及び他のＵＩ要素は、Ｗｅｂブラウザの機能に関連付けられまたＷｅｂブラウザの機能を表すことができる。
【０１８８】
本明細書で記述された主題の実施は、バックエンドコンポーネントを例えばデータサーバとして含むコンピューティングシステム、ミドルウェアコンポーネントを、例えばアプリケーションサーバとして含むコンピューティングシステム、フロントエンドコンポーネントを、例えばクライアントコンピュータとして含むコンピューティングシステム、及び、一又は複数のこのようなバックエンド、ミドルウェア、又はフロントエンドコンポーネントの任意の組み合わせを含むコンピューティングシステムにおいて実施でき、クライアントコンピュータは、ユーザが本明細書においで記述される主題の実施とインタラクションできるグラフィカルユーザインタフェース又はＷｅｂブラウザを有する。システムのコンポーネントは、有線又は無線のデジタルデータ通信（又はデータ通信の組み合わせ）、例えば通信ネットワークの任意の形式又は媒体によって相互接続されることができる。通信ネットワークの例示は、ローカルエリアネットワーク（ＬＡＮ）、無線アクセスネットワーク（ＲＡＮ）、メトロポリタンエリアネットワーク（ＭＡＮ）、広域ネットワーク（ＷＡＮ）、マイクロ波アクセスの世界的な相互運用性（ＷＩＭＡＸ）、例えば８０２．１１ａ／ｂ／ｇ／ｎ又は８０２．２０（又は８０２．１１ｘと８０２．２０の組み合わせ、又は本願に整合する他のプロトコル）を用いる無線ローカルエリア（ＷＬＡＮ）、インターネットの全て若しくは一部、別の通信ネットワーク、又は通信ネットワークの組み合わせを含む。通信ネットワークは、例えば、インターネットプロトコル（ＩＰ）パケット、フレームリレーフレーム、非同期転送モード（ＡＴＭ）セル、音声、ビデオ、データ、又はネットワークノード間の他の情報と通信できる。
【０１８９】
コンピューティングシステムは、クライアント及びサーバを含むことができる。クライアント及びサーバは、一般的には、互いに離れており、また典型的には、通信ネットワークを介してインタラクションする。クライアント及びサーバの関係は、それぞれのコンピュータ上で動作すると共に互にクライアント・サーバ関係にあるコンピュータプログラムのおかげで生じる。
【０１９０】
本明細書は多くの特定の実施の詳細を含む一方で、これらは、任意の発明概念の範囲又は請求される得る範囲の制限として解釈されるべきではなく、むしろ、特定の発明概念の特定の実施固有の特徴の説明として解釈されるべきである。個別の実施の観点で本明細書において説明された機能は、組み合わせて又は単一の実施で、実現できる。逆に、単一の実施の観点で既述された様々な特徴は、複数の実施で、個別に、又は任意のサブコンビネーションで実施されることもできる。さらには、既述の特徴は特定の組み合わせで動作するものとして説示され、最初はそのようなものとして請求されているが、一又は複数の特徴は、請求された組み合わせから場合によっては削除され、またサブコンビネーション又はサブコンビネーションの変形であり得る。
【０１９１】
主題の特定の実施が説明された。記載された実施の置換、変更、及び他の実施は、当業者には明らかである以下の請求の範囲内である。操作は特許請求の範囲に又は特定の順序で図面に描かれている一方で、これは、望ましい結果を達成するために、示された特定の順序で又は順番でそのような操作が実行されること又は全ての図示された操作が実行される（いくつかの操作は随意的と見なされる）ことを要求するものとして理解されるべきではない。ある状況では、マルチタスク又は並列処理（又はマルチタスク及び並列処理の組み合わせ）を行うことが有利であり、適切と思われる場合に実行される。
【０１９２】
さらに、既述の実施における様々なシステムモジュール及びコンポーネントの分離又は統合が、全ての実施においてそのような分離又は統合を必要とするものとして理解されるべきではない。また、記載されたプログラムコンポーネント及びシステムは、一般的には、単一のソフトウェア製品に統合でき、又は複数のソフトウェア製品にパッケージ化できることが理解されるべきである。
【０１９３】
これに従って、既述の例示的な実施は、本願を定義し又は制約しない。他の変更、置換、及び変更も、本開示の範囲及び精神から逸脱することなく可能である。
【０１９４】
さらには、請求された任意の実施は、少なくとも、コンピュータにより実施される方法、非一時的なコンピュータ読取可能媒体、及びコンピュータシステムに適用可能であると考えられ、非一時的なコンピュータ読取可能媒体は、コンピュータにより実施される方法を実行するコンピュータにより読取可能な命令を格納し、またコンピュータシステムは、ハードウェアプロセッサに相互に動作可能に結合されたコンピュータメモリを含み、ハードウェアプロセッサは、非一時的なコンピュータ読取可能媒体に格納されたコンピュータにより実施される方法又は命令を実行するように構成される。
以下、本発明の実施の態様の例を列挙する。
［第１の局面］
地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理方法であって：
１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩から、ＨＣが排出されたことを特定するステップと；
潜在的なトラップピークを前記複数のグリッドセル内で識別するステップと；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行するステップと；
ＨＣ類を含むグリッドセルが過剰量のＨＣ類を含んでいるかどうかを特定するステップと；
前記識別された潜在的なトラップピークに関連付けられたトラップにおける、グリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行するステップと；
前記トラップをＨＣのポテンシャル値で更新するステップと；を備える、
コンピュータによって実施される並列処理方法。
［第２の局面］
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第３の局面］
スピルポイントに達したことを特定するステップと；
戻って前記侵入プロセスを実行するステップと；を更に備える、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第４の局面］
前記侵入プロセスは：
前記ＨＣがサブドメイン境界に達したことを特定するステップと；
前記ＨＣを、隣接するサブドメインに移動するステップと；
前記ＨＣが同じサブドメインで充填トラップに達した場合、前記ＨＣを前記充填トラップと併合するステップと；を備える、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第５の局面］
前記サブドメイン境界においてＨＣ量を通信するステップと；
サブドメイン境界を越えた前記ＨＣを、別のサブドメインに関連付けられたトラップと併合するステップと；を更に備える、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第６の局面］
前記ＨＣがサブドメイン境界に達していることを特定するステップと；
前記サブドメイン境界に関連付けられた隣接するサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信するステップと；を更に備える、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第７の局面］
前記ＨＣがトラップ境界に達したことを特定するステップと；
複数のトラップが集積境界を共有するかどうかを特定するステップと；を更に備える、
第１の局面に記載のコンピュータによって実施される並列処理方法。
［第８の局面］
地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理の方法の１又は複数の命令を格納する非一時的なコンピュータ読取り可能媒体であって、前記１又は複数の命令は所定の操作を実行するためにコンピュータシステムによって実行可能であり、前記操作は：
１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作と；
潜在的なトラップピークを前記複数のグリッドセル内で識別する操作と；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作と；
ＨＣ類を含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；
前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作と； 
前記トラップをＨＣのポテンシャル値で更新する操作と；を備える、
非一時的なコンピュータ読取り可能媒体。
［第９の局面］
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１０の局面］
スピルポイントに達していることを特定する操作と；
戻って前記侵入プロセスを実行する操作と；を実行する１又は複数の命令を更に備える、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１１の局面］
前記侵入プロセスは：
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記ＨＣを隣接するサブドメインに移動する操作と；
前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する操作と；を実行する１又は複数の命令を備える、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１２の局面］
前記サブドメイン境界でＨＣ量を通信する操作と；
サブドメインの境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する操作と；を実行する１又は複数の命令を更に備える、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１３の局面］
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信する操作と；を実行する１又は複数の命令を更に備える、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１４の局面］
前記ＨＣがサブドメイン境界に達したことを特定する操作と；
複数のトラップが集積境界を共有するかどうかを特定する操作と；を実行する１又は複数の命令を更に備える、
第８の局面に記載の非一時的なコンピュータ読取り可能媒体。
［第１５の局面］
地下媒体内の炭化水素（ＨＣ）の移動をシミュレートするための、コンピュータによって実施される並列処理システムであって：
コンピュータメモリと；
前記コンピュータメモリと相互作動可能に結合され、所定の操作を実行するように構成されたハードウェアプロセッサであって、前記操作は：
１又は複数のサブドメインに分割された複数のグリッドセルに関連付けられた根源岩からＨＣが排出されたことを特定する操作と；
潜在的なトラップピークを前記複数のグリッドセル内で識別する操作と；
前記ＨＣが前記複数のグリッドセル内で移動を停止するまで侵入プロセスを実行する操作と；
ＨＣを含むグリッドセルに過剰量のＨＣが含まれているかどうかを特定する操作と；
前記識別された潜在的なトラップピークに関連付けられたトラップにおけるグリッドセルへの前記ＨＣの充填をモデル化するために、集積プロセスを実行する操作と；
前記トラップをＨＣポテンシャル値で更新する操作と；を備える、前記ハードウェアプロセッサと；を備える、
コンピュータによって実施される並列処理システム。
［第１６の局面］
前記潜在的なトラップピークを識別するために、隣接するすべてのグリッドセルが特定のグリッドセルよりも高いポテンシャル値を有する、前記特定のグリッドセルを特定するように各セルの絶対油相ポテンシャル値をチェックする、
第１５の局面に記載のコンピュータによって実施される並列処理システム。
［第１７の局面］
スピルポイントに達したことを特定し；
戻って前記侵入プロセスを実行する；ように更に構成される、
第１５の局面に記載のコンピュータによって実施される並列処理システム。
［第１８の局面］
前記侵入プロセスは：
前記ＨＣがサブドメイン境界に達したことを特定し；
前記ＨＣを、隣接するサブドメインに移動し；
前記ＨＣが同じサブドメインで充填されたトラップに達した場合、前記ＨＣを前記充填されたトラップと併合する；ように更に構成される、
第１５の局面に記載のコンピュータによって実施される並列処理システム。
［第１９の局面］
前記サブドメイン境界でＨＣ量を通信し；
サブドメイン境界を越えて前記ＨＣを別のサブドメインに関連付けられたトラップと併合する；ように更に構成される、
第１５の局面に記載のコンピュータによって実施される並列処理システム。
［第２０の局面］
前記ＨＣがサブドメイン境界に達したことを特定し；
前記サブドメイン境界に関連付けられた隣接のサブドメインに、ＨＣ量、最小の油ポテンシャル値を有するグリッドセルのグリッドセル指標、及び境界グリッドセルリストを送信するように更に構成される、
第１５の局面に記載のコンピュータによって実施される並列処理システム。
【符号の説明】
【０１９５】
１８００  コンピュータシステム
１８０２  コンピュータ
１８０５  プロセッサ
１８０７  メモリ
【図１Ａ】
0006966577.tif 000010
【図１Ｂ】
0006966577.tif 000011
【図１Ｃ】
0006966577.tif 000012
【図２Ａ】
0006966577.tif 000013
【図２Ｂ】
0006966577.tif 000014
【図３】
0006966577.tif 000015
【図４Ａ】
0006966577.tif 000016
【図４Ｂ】
0006966577.tif 000017
【図５】
0006966577.tif 000018
【図６】
0006966577.tif 000019
【図７】
0006966577.tif 000020
【図８】
0006966577.tif 000021
【図９】
0006966577.tif 000022
【図１０】
0006966577.tif 000023
【図１１】
0006966577.tif 000024
【図１２】
0006966577.tif 000025
【図１３】
0006966577.tif 000026
【図１４Ａ】
0006966577.tif 000027
【図１４Ｂ】
0006966577.tif 000028
【図１５】
0006966577.tif 000029
【図１６】
0006966577.tif 000030
【図１７】
0006966577.tif 000031
【図１８】
0006966577.tif 000032