JP6824205B2 - コンピュータ実装方法、データ処理システム、及びデータストレージデバイス - Google Patents
コンピュータ実装方法、データ処理システム、及びデータストレージデバイス Download PDFInfo
- Publication number
- JP6824205B2 JP6824205B2 JP2018019074A JP2018019074A JP6824205B2 JP 6824205 B2 JP6824205 B2 JP 6824205B2 JP 2018019074 A JP2018019074 A JP 2018019074A JP 2018019074 A JP2018019074 A JP 2018019074A JP 6824205 B2 JP6824205 B2 JP 6824205B2
- Authority
- JP
- Japan
- Prior art keywords
- reservoir
- matrix
- well
- data
- oil layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 148
- 238000012545 processing Methods 0.000 title claims description 64
- 238000013500 data storage Methods 0.000 title claims description 16
- 239000011159 matrix material Substances 0.000 claims description 188
- 239000012530 fluid Substances 0.000 claims description 94
- 238000004364 calculation method Methods 0.000 claims description 65
- 238000005553 drilling Methods 0.000 claims description 58
- 229930195733 hydrocarbon Natural products 0.000 claims description 12
- 150000002430 hydrocarbons Chemical class 0.000 claims description 12
- 239000004215 Carbon black (E152) Substances 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 150000002632 lipids Chemical class 0.000 claims 3
- 239000013598 vector Substances 0.000 description 74
- 238000004088 simulation Methods 0.000 description 38
- 230000008569 process Effects 0.000 description 27
- 238000004422 calculation algorithm Methods 0.000 description 20
- 238000004891 communication Methods 0.000 description 20
- 238000010168 coupling process Methods 0.000 description 19
- 238000005859 coupling reaction Methods 0.000 description 19
- 230000008878 coupling Effects 0.000 description 18
- 230000008602 contraction Effects 0.000 description 15
- 238000005094 computer simulation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 239000000284 extract Substances 0.000 description 8
- 239000011148 porous material Substances 0.000 description 8
- 239000002131 composite material Substances 0.000 description 6
- 230000006399 behavior Effects 0.000 description 5
- 230000000052 comparative effect Effects 0.000 description 5
- 229910000614 lithium tin phosphorous sulfides (LSPS) Inorganic materials 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000013626 chemical specie Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 125000001475 halogen functional group Chemical group 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000013067 intermediate product Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Operations Research (AREA)
- Geophysics And Detection Of Objects (AREA)
- Complex Calculations (AREA)
Description
導入部として、上記で確認されたタイプの既知のコンピュータ油層シミュレーションのより詳細な説明が提供される。地下油層からの石油およびガスの採取のために、掘削孔は炭化水素流体の採取用にこれらの層に掘削される。採取プロセスの間、水および/またはガスなどの流体は、圧入井に注入され、孔隙内の流体混合物は生産井から生産される。これらの油層の将来の挙動を予測し、代替の開発計画を評価するために、油層シミュレータはシミュレーションモデルを実行するために使用される。これらのモデルは、最初に、既存の生産データを使用して、履歴整合ステップで較正される。較正されたモデルは、次いで、将来の操業シナリオを評価するために使用される。たとえば、履歴整合モデルは、所定の位置に残っている炭化水素のうちの多くを採取するために、さらなる坑井をいつどこで掘削するかを判断するために使用される場合がある。
下記に記載される方程式内のシンボルは以下の意味を有する。
p=圧力
q=生産率
xi=モル分率
Vj=相体積
Sj=相飽和
ci=化学種iの全体濃度
φ=多孔度
ρ=密度
μ=粘度
ω=質量分率
R=均一反応率
D=分散係数
u=速度
Vi=岩孔体積
ni t=モル総数
上付き文字:
ref=基準
p=流体相
t=合計
下付き文字:
i=成分インデックス
j=相インデックス
化学種成分i用の一般的な化学種保存方程式は、
によって与えられる。
ここで、
である。
と簡略化される。
として表すことができる。
ここで、細孔体積Vφは圧力のみの関数であり、
として記述される。
である。
を有する。
WIは、通常、レイヤ坑井インデックスと呼ばれる。WIは、坑井の近くの定常状態の圧力プロファイルを介する坑井の坑底圧力に対する坑井セル有限差分グリッド圧力および坑井セルについての相当半径roの決定に関する。Peaceman数式として知られているものは、しばしば、油層シミュレータにおいてレイヤ坑井インデックスを決定するために使用され、これは、後に、異方性矩形グリッドおよび他のより複雑な状況に拡張されている。
陰解坑井方程式および陰解油層方程式を有する完全連成系は、以下の通りである。
ここで、非線形坑井方程式および非線形油層方程式のヤコビ導関数は、
であり、XWおよびXRは、坑井変数および油層変数であり、RWおよびRRは、初期坑井残差および初期油層残差である。各グリッドセルは複数の方程式を有することができ、坑井グリッド用の方程式の数と油層グリッド用の方程式の数は異なるので、nvarはグリッド当たりの油層方程式の数として、nwvarはグリッド当たりの坑井方程式の数として使用される。nvarは、一般に、nwvarと等しくないことを留意されたい。したがって、AWW内の各非ゼロ値は(nwvar*nwvar)の高密度ブロックであり、ARR内の各非ゼロ値は(nvar*nvar)の高密度ブロックであり、AWR内の各非ゼロ値は(nwvar*nvar)の高密度ブロックであり、ARW内の各非ゼロ値は(nvar*nwvar)の高密度ブロックである。
この方法では、複合行列が直接解かれた。すなわち、行列方程式は
と設定された。
ここで、
である。
この解法は、最初に油層方程式から坑井方程式を切り離して縮約系を形成する。
図5は、式(16)による縮約系行列50を有する、図4の油層モデルM−2の同じ例を示す。図5は、完全陰解坑井結合項によって誘導された追加の非ゼロ高密度ブロック導関数を有する縮約系行列の小さな例である。モデル内の坑井が非常に少ないグリッドセルを貫通するとき、この方法は許容できる可能性がある。しかしながら、油層内の坑井が多くのグリッドセルを貫通する場合、縮約系行列50内の追加の非ゼロ値の数は過剰になり、この方法の魅力を失う。これは、縮約系行列
が多くの坑井誘導非ゼロ項を有していたためであって、坑井項がない油層行列と比較して、解くことをはるかに高価にする。
現在のシミュレータにおいて使用される完全陰解結合系解法用の最新の技術の始まりの主要な方法は、RowsumまたはColsumのいずれかの事前調整方法である。この方法は、簡単で、適当な数の坑井が穴を開けたグリッドセルを有する坑井に適しているが、シミュレーションモデルが、複雑な異成分を有し、非常に高い穿孔多相フローレートを有する数千のグリッドセルに穴を開ける場合がある多くの複合マルチラテラル長距離坑井を含んでいるとき、それは全ソルバーに対して収束が遅くて正常に動かなくなるか、または収束に失敗する場合がある。これらの複雑な状況では、結合ソルバーは、セル間のフロー項に対してわずかに事前調整されているだけの過剰な坑井影響係数の結合に起因して、失速するか、または全く収束できない可能性がある。
は、元のARR行列よりもかなり多くの非ゼロ値を有することができるので、それは明示的に計算されない。事前調整のために、最新技術は、列和(Colsum)または行和(Rowsum)のいずれかを行って坑井影響行列を対角化する。ネスト化因数分解(NF)事前調整アルゴリズムでは、NFは、行列:
に適用される。
に適用される。
行列−ベクトル積は、積が必要なときはいつでも、一連の行列−ベクトル乗数として計算される。
前述の困難を克服するために、本発明は、縮約系行列用の高度に並列化可能なプレコンディショナに基づく、堅牢で、さらに計算するのに効率的である、コンピュータ実装方法論を提供する。本発明は、炭化水素貯留層の挙動を分析し、それらの開発を評価するために使用される、油層シミュレーションの既存の技術的プロセスに対する改善点を提供する。本発明はまた、シミュレータ処理において不十分な収束または収束の失敗に起因して失われる処理時間を削減することにおいて、およびデータ処理システムのネットワーク内の通信の複雑さを低減することにおいても、油層シミュレーションの挙動におけるコンピュータの機能化を改善する。式(16)において、
である。
[命名法]
ARR=油層グリッドヤコビ行列
AWR=坑井から油層への結合ヤコビ行列
ARW=油層からの坑井への結合ヤコビ行列
AWW=マルチセグメント坑井ヤコビ行列
AP=圧力係数ヤコビ行列
C=圧力係数セレクタ
CT=圧力係数セレクタの移項
ER=油層グリッドヤコビ行列の残余成分
EW=坑井影響係数ヤコビ行列
M−1=全系統プレコンディショナ
PR=油層グリッドヤコビ行列の区分化成分
R,r=残差ベクトル
W=圧力減結合演算子行列
X=解ベクトル
添字:
C=補正
P=圧力
R=油層
W=坑井
T=合計
坑井結合および油層グリッド非ゼロ値充填項用の可変次数を用いた並列1レベル線解べき級数(LSPS)事前調整
均一1レベル事前調整方法は、以下のように記載される。
この方法では、以下の置換が行われる。
ARR=[PR+ER](21)
である。
である。
を与える。
である。
である。
である。
の結果は簡単に計算することができる。たとえば、PRの1つの適切な選択は、Z線で順序付けられた油層連立方程式用のブロック三重対角行列である。
を使用する近似逆プレコンディショナは、
と書くことができる。
である。
である。
になる。
である。
ri=b−[A]xi
が行われ得るように推定値xiが選択される。
が方程式(28)に従い、その成分は、完全連成完全陰解坑井影響ヤコビ行列を伴う方程式(21)から(27)によって詳述されている。近似逆プレコンディショナは、記載されるように、クリロフ部分空間反復法内で適用される。
[並列CPR型線解べき級数事前調整]
拘束圧残差CPR事前調整方法として知られているものは、Wallis,J.R.らの「Constrained Residual Acceleration of Conjugate Residual Methods」、SPE 13563、Proceedings of the 8th SPE Reservoir Simulation Symposium、ダラス、米国、1985年2月10〜13日に最初に記載された。CPR法の変形形態は、FungおよびDogruの「Parallel Unstructured Solver Methods for Simulation of Complex Giant Reservoirs」、SPE(2008年12月)、ページ440〜446においてさらに説明された。
と書くことができる。
ここで、
である。
である。
によって与えられる。
W=DIAG−1(A)(47)
として計算される。
がもう少しで正定値になるように選択される必要があることに留意されたい。これは、LSPS、またはILU変形形態のうちのいずれかに対する要件ではない。
は以下のように略述することができる。
1.全系統残差を圧力系に限定する。
rp=CTWr (48)
2.圧力系を繰り返し解く。
Apxp=rp (49)
3.圧力解を全系統に拡張する。
s=Cxp (50)
4.圧力残差から全系統残差を補正する。
rc=r−A・s (51)
5.第2段階の事前調整を使用して全系統を解く。
Mx=rc (52)
6.圧力解から全系統解を補正する。
xc=x+s (53)
7.圧力系事前調整は、並列代数マルチグリッド法(parallel algebraic multi−grid)すなわちPAMG法、または、上記の実施形態1に記述されたILU(k)法もしくはLSPS法の使用を含む、いくつかの方法で遂行することができる。
を伴う。油層解ベクトルが計算された後、行列方程式(17)を使用して、坑井解ベクトルが計算される。図7に示された処理は、図15のデータ処理システムにおいて実施される。
であり、逆数
は、方程式(29)から(33)におけるように、[U]−1[L]−1として計算される。
[並列結合事前調整]
この不均一事前調整方法は、坑井影響係数行列に対して上記の実施形態1による事前調整を使用するが、油層係数行列用の異なる方法を含む。これは、油層項と坑井影響項に対して異なるプレコンディショナを適用するために必要とされる場合の助けになる。結合プレコンディショナは、
と書かれる。
と書くことができる。
は、油層項向けの別の効果的なプレコンディショナであり得る。たとえば、並列アプリケーションのための適切な領域分解法を有するA行列の油層部向けのILU(k)またはILUTが、ここで選ばれる場合がある。
がILU(k)プレコンディショナである場合、
である。
を伴う近似解ベクトルを計算するための近似逆事前調整方法を含む。近似逆プレコンディショナは、クリロフ部分空間反復アルゴリズム内で適用される。
図8および図9は、それぞれ、図2および図3の20および22に示されたような結合陰解掘削孔用のデータを処理するために、図1の領域Dの油層グリッドの非構造化部分領域の分散並列処理のためのグリッドデータの構成を概略的に示す。油層領域向けの分散並列処理用のデータは、構造化グリッドおよび非構造化グリッドのために複数の部分領域に区分化され、異なる計算プロセスに割り当てられた各部分領域は、それぞれ、米国特許第8,433,551号および米国特許第8,386,227号に記載されており、それらの各々において、出願人は発明者と呼ばれる。完全連成完全陰解坑井−油層解法のための分散並列処理は、以下の方式で遂行される。
Claims (15)
- 地下の炭化水素貯留層内のマルチラテラル坑井の掘削孔内の多相流体のフローのシミュレーションのコンピュータ実装方法であって、前記マルチラテラル坑井が、入力された油層データに基づいて油層セルのグリッドに編成された地下の油層のセルを有する前記油層内の長さ方向の範囲に沿った複数の位置で流体の交差を有し、前記油層セルが、その中で起こる多相流体フローを有し、
前記コンピュータにおいて、前記マルチラテラル坑井のうちのいくつかとのフロー交差を有する油層セル用の、油層データを有する陰解油層方程式、圧力方程式、および流動方程式の完全連成非線形セットに基づき、油層からの坑井への結合ヤコビ行列である油層計算行列ARWを編成するステップと、
前記コンピュータにおいて、前記油層セルのうちのいくつかとのフロー交差を有する掘削孔セル用の、坑井データを有する陰解坑井方程式および流動方程式の完全連成非線形セットに基づき、坑井から油層への結合ヤコビ行列である掘削孔計算行列AWRを編成するステップと、
前記コンピュータにおいて、前記掘削孔セルと前記油層セルの前記フロー交差に基づいて坑井影響行列EWを編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRを含む全系統計算行列であって、XWを坑井に関する変数とし、XRを油層に関する変数とし、RWを坑井に関する初期残差とし、RRを油層に関する初期残差とし、AWWを、マルチセグメント坑井ヤコビ行列であるマルチセグメント坑井行列とし、ARRを、油層グリッドヤコビ行列である油層グリッド行列としたときに、以下が成立するものとして定義される全系統計算行列を編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRの圧力係数を抽出するステップと、
全系統残差から圧力残差を抽出するステップと、
前記抽出された圧力残差を最小化することにより、前記全系統計算行列の前記油層セルおよび掘削孔セル内の圧力についての近似圧力解を解くステップと、
前記近似圧力解に基づいて前記全系統計算行列の、前記油層セルについての流体圧力および前記残差を更新するステップと、
前記全系統計算行列、前記坑井影響行列、ならびに前記更新された圧力および残差についての近似全系統更新値を計算するステップと、
前記近似全系統更新値を、前記更新された流体圧力と組み合わせるステップと、
前記全系統残差を更新するステップと、
保存方程式の前記完全連成非線形セットおよび前記更新された全系統残差を使用して前記全系統計算行列を解くことにより、前記多相流体フローを決定するステップと
を備える、コンピュータ実装方法。 - 前記コンピュータ実装方法が、各々が並行して動作する複数のコンピューティングコアを備える複数のコンピュータノードを有するプロセッサ内で実施され、
並列に動作するコンピューティングコアの数に従って、前記油層データおよび前記坑井データをいくつかの並列データ下位領域に区分化するステップと、
前記区分化されたデータ下位領域を割り当てて、並列データ下位領域を形成するステップと
をさらに含む、請求項1に記載のコンピュータ実装方法。 - 前記油層が百万セルよりも大きいグリッドブロックの領域に編成される、請求項1に記載のコンピュータ実装方法。
- 前記決定された多相流体フローの出力表示を形成するステップ
をさらに含む、請求項1に記載のコンピュータ実装方法。 - 前記決定された多相流体フローの記録を形成するステップ
をさらに含む、請求項1に記載のコンピュータ実装方法。 - 地下の炭化水素貯留層内のマルチラテラル坑井の掘削孔内の多相流体のフローをシミュレートするプロセッサを備えたデータ処理システムであって、前記マルチラテラル坑井が、入力された油層データに基づいて油層セルのグリッドに編成された地下の油層のセルを有する前記油層内の長さ方向の範囲に沿った複数の位置で流体の交差を有し、前記油層セルが、その中で起こる多相流体フローを有し、前記プロセッサが、
前記データ処理システムにおいて、前記マルチラテラル坑井のうちのいくつかとのフロー交差を有する油層セル用の、油層データを有する陰解油層方程式、圧力方程式、および流動方程式の完全連成非線形セットに基づき、油層からの坑井への結合ヤコビ行列である油層計算行列ARWを編成するステップと、
前記データ処理システムにおいて、前記油層セルのうちのいくつかとのフロー交差を有する掘削孔セル用の、坑井データを有する陰解坑井方程式および流動方程式の完全連成非線形セットに基づき、坑井から油層への結合ヤコビ行列である掘削孔計算行列AWRを編成するステップと、
前記データ処理システムにおいて、前記掘削孔セルと前記油層セルの前記フロー交差に基づいて坑井影響行列EWを編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRを含む全系統計算行列であって、XWを坑井に関する変数とし、XRを油層に関する変数とし、RWを坑井に関する初期残差とし、RRを油層に関する初期残差とし、AWWを、マルチセグメント坑井ヤコビ行列であるマルチセグメント坑井行列とし、ARRを、油層グリッドヤコビ行列である油層グリッド行列としたときに、以下が成立するものとして定義される全系統計算行列を編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRの圧力係数を抽出するステップと、
全系統残差から圧力残差を抽出するステップと、
前記抽出された圧力残差を最小化することにより、前記全系統計算行列の前記油層セルおよび掘削孔セル内の圧力についての近似圧力解を解くステップと、
前記近似圧力解に基づいて前記全系統計算行列の、前記油層セルについての流体圧力および前記残差を更新するステップと、
前記全系統計算行列、前記坑井影響行列、ならびに前記更新された圧力および残差についての近似全系統更新値を計算するステップと、
前記近似全系統更新値を、前記更新された流体圧力と組み合わせるステップと、
前記全系統残差を更新するステップと、
保存方程式の前記完全連成非線形セットおよび前記更新された全系統残差を使用して前記全系統計算行列を解くことにより、前記多相流体フローを決定するステップと
を実施するものであるデータ処理システム。 - プロセッサが、各々が並行して動作する複数のコンピューティングコアから構成される複数のコンピュータノードを備え、前記コンピューティングコアが、
並列に動作するコンピューティングコアの数に従って、前記油層データおよび前記坑井データをいくつかの並列データ下位領域に区分化するステップと、
前記区分化されたデータ下位領域を割り当てて、並列データ下位領域を形成するステップと
を実施する、請求項6に記載のデータ処理システム。 - 前記油層が百万セルよりも大きいグリッドブロックの領域に編成される、請求項6に記載のデータ処理システム。
- 前記決定された多相流体フローの出力表示を形成するワークステーション
をさらに含む、請求項6に記載のデータ処理システム。 - 前記決定された多相流体フローの記録を記憶するデータメモリ
をさらに含む、請求項6に記載のデータ処理システム。 - 地下の炭化水素貯留層内のマルチラテラル坑井の掘削孔内の多相流体のフローをプロセッサにシミュレートさせるためのコンピュータで実行可能な命令を、非一時的でコンピュータで読み取り可能な媒体に記憶したデータストレージデバイスであって、前記マルチラテラル坑井が、入力された油層データに基づいて油層セルのグリッドに編成された地下の油層のセルを有する前記油層内の長さ方向の範囲に沿った複数の位置で流体の交差を有し、前記油層セルが、その中で起こる多相流体フローを有し、前記データストレージデバイスに記憶された前記命令が、
前記コンピュータにおいて、前記マルチラテラル坑井のうちのいくつかとのフロー交差を有する油層セル用の、油層データを有する陰解油層方程式、圧力方程式、および流動方程式の完全連成非線形セットに基づき、油層からの坑井への結合ヤコビ行列である油層計算行列ARWを編成するステップと、
前記コンピュータにおいて、前記油層セルのうちのいくつかとのフロー交差を有する掘削孔セル用の、坑井データを有する陰解坑井方程式および流動方程式の完全連成非線形セットに基づき、坑井から油層への結合ヤコビ行列である掘削孔計算行列AWRを編成するステップと、
前記コンピュータにおいて、前記掘削孔セルと前記油層セルの前記フロー交差に基づいて坑井影響行列EWを編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRを含む全系統計算行列であって、XWを坑井に関する変数とし、XRを油層に関する変数とし、RWを坑井に関する初期残差とし、RRを油層に関する初期残差とし、AWWを、マルチセグメント坑井ヤコビ行列であるマルチセグメント坑井行列とし、ARRを、油層グリッドヤコビ行列である油層グリッド行列としたときに、以下が成立するものとして定義される全系統計算行列を編成するステップと、
前記油層計算行列ARWおよび前記掘削孔計算行列AWRの圧力係数を抽出するステップと、
全系統残差から圧力残差を抽出するステップと、
前記抽出された圧力残差を最小化することにより、前記全系統計算行列の前記油層セルおよび掘削孔セル内の圧力についての近似圧力解を解くステップと、
前記近似圧力解に基づいて前記全系統計算行列の、前記油層セルについての流体圧力および前記残差を更新するステップと、
前記全系統計算行列、前記坑井影響行列、ならびに前記更新された圧力および残差についての近似全系統更新値を計算するステップと、
前記近似全系統更新値を、前記更新された流体圧力と組み合わせるステップと、
前記全系統残差を更新するステップと、
保存方程式の前記完全連成非線形セットおよび前記更新された全系統残差を使用して前記全系統計算行列を解くことにより、前記多相流体フローを決定するステップと
を前記プロセッサに実施させる、データストレージデバイス。 - 前記命令が、各々が並行して動作する複数のコンピューティングコアを備える複数のコンピュータノードを有するプロセッサ内で実施され、前記命令が、
並列に動作するコンピューティングコアの数に従って、前記油層データおよび前記坑井データをいくつかの並列データ下位領域に区分化するステップと、
前記区分化されたデータ下位領域を割り当てて、並列データ下位領域を形成するステップと
を前記コンピューティングコアに実施させる命令をさらに備える、請求項11に記載のデータストレージデバイス。 - 前記油層が百万セルよりも大きいグリッドブロックの領域に編成される、請求項11に記載のデータストレージデバイス。
- 前記データストレージデバイスに記憶された前記命令が、
前記決定された多相流体フローの出力表示を形成するステップ
を前記コンピューティングコアにさらに実施させる、請求項12に記載のデータストレージデバイス。 - 前記データストレージデバイスに記憶された前記命令が、
前記決定された多相流体フローの記録を形成するステップ
を前記コンピューティングコアにさらに実施させる、請求項12に記載のデータストレージデバイス。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562164083P | 2015-05-20 | 2015-05-20 | |
US62/164,083 | 2015-05-20 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017560232A Division JP2018518613A (ja) | 2015-05-20 | 2016-05-17 | 油層シミュレーションにおける並列解法または完全連成完全陰解掘削孔モデル化 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018119967A JP2018119967A (ja) | 2018-08-02 |
JP6824205B2 true JP6824205B2 (ja) | 2021-02-03 |
Family
ID=56084413
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017560232A Pending JP2018518613A (ja) | 2015-05-20 | 2016-05-17 | 油層シミュレーションにおける並列解法または完全連成完全陰解掘削孔モデル化 |
JP2018019074A Active JP6824205B2 (ja) | 2015-05-20 | 2018-02-06 | コンピュータ実装方法、データ処理システム、及びデータストレージデバイス |
JP2020076641A Active JP6887042B2 (ja) | 2015-05-20 | 2020-04-23 | コンピュータ実装方法、データ処理システム、およびデータストレージデバイス |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017560232A Pending JP2018518613A (ja) | 2015-05-20 | 2016-05-17 | 油層シミュレーションにおける並列解法または完全連成完全陰解掘削孔モデル化 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020076641A Active JP6887042B2 (ja) | 2015-05-20 | 2020-04-23 | コンピュータ実装方法、データ処理システム、およびデータストレージデバイス |
Country Status (6)
Country | Link |
---|---|
US (4) | US10242136B2 (ja) |
EP (2) | EP3399141A1 (ja) |
JP (3) | JP2018518613A (ja) |
CN (2) | CN107849910B (ja) |
CA (1) | CA2986034C (ja) |
WO (1) | WO2016187175A1 (ja) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9164191B2 (en) | 2011-02-09 | 2015-10-20 | Saudi Arabian Oil Company | Sequential fully implicit well model for reservoir simulation |
US10113400B2 (en) * | 2011-02-09 | 2018-10-30 | Saudi Arabian Oil Company | Sequential fully implicit well model with tridiagonal matrix structure for reservoir simulation |
EP2839409B1 (en) * | 2012-06-15 | 2017-02-08 | Landmark Graphics Corporation | Systems and methods for solving a multi-reservoir system with heterogeneous fluids coupled to common gathering network |
CN107368455A (zh) * | 2017-06-22 | 2017-11-21 | 东南大学 | 一种gpu加速的电力潮流上三角方程组回代方法 |
US11499412B2 (en) | 2017-11-24 | 2022-11-15 | Total Se | Method and device for determining hydrocarbon production for a reservoir |
CN109002629B (zh) * | 2018-08-01 | 2023-02-03 | 苏州慧德仿真技术有限公司 | 一种多相流模拟仿真的卷积神经网络及快速可视化方法 |
CN110863818B (zh) * | 2018-08-08 | 2023-08-29 | 中国石油化工股份有限公司 | 一种剩余油/气分布的描述方法及装置 |
US11461514B2 (en) * | 2018-09-24 | 2022-10-04 | Saudi Arabian Oil Company | Reservoir simulation with pressure solver for non-diagonally dominant indefinite coefficient matrices |
CN109108168B (zh) * | 2018-10-30 | 2019-12-27 | 上汽大众汽车有限公司 | 冲压模具合模间隙的计算方法 |
WO2020123698A1 (en) * | 2018-12-11 | 2020-06-18 | Resonant, Inc. | Fast, highly accurate, full-fem surface acoustic wave simulation |
WO2020160490A1 (en) | 2019-02-01 | 2020-08-06 | Lightelligence, Inc. | Processing matrix operations for rate limited systems |
EP3956699B1 (en) * | 2019-04-16 | 2023-03-29 | TotalEnergies OneTech | A method for upscaling of relative permeability of the phase of a fluid |
CN110146903B (zh) * | 2019-05-24 | 2020-11-13 | 国网浙江省电力有限公司信息通信分公司 | 一种基于反馈调整惯性权重的粒子群北斗卫星选择方法 |
CN110366126B (zh) * | 2019-06-17 | 2020-08-04 | 上海交通大学 | 移动式自组网络中脱离节点等待时间最优决策方法及系统 |
CN110348158B (zh) * | 2019-07-18 | 2021-07-27 | 中国水利水电科学研究院 | 一种基于分区异步长解的地震波动分析方法 |
US11499397B2 (en) | 2019-10-31 | 2022-11-15 | Saudi Arabian Oil Company | Dynamic calibration of reservoir simulation models using flux conditioning |
US11501038B2 (en) | 2019-10-31 | 2022-11-15 | Saudi Arabian Oil Company | Dynamic calibration of reservoir simulation models using pattern recognition |
CN111859766B (zh) * | 2020-07-28 | 2024-01-23 | 福建省拳石科技发展有限公司 | 可变计算域的拉格朗日积分点有限元数值仿真系统及方法 |
CN112949238A (zh) * | 2021-03-19 | 2021-06-11 | 梁文毅 | 一种基于迭代法的电气仿真方法 |
CN113324840B (zh) * | 2021-05-31 | 2022-04-22 | 西南石油大学 | 一种非均质地层井壁渐进坍塌过程的流-固-热耦合模拟方法 |
CN117290642A (zh) * | 2022-10-28 | 2023-12-26 | 国家电投集团科学技术研究院有限公司 | 基于Newton-Raphson求解器的热工水力模型的耦合方法、装置及设备 |
Family Cites Families (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5321612A (en) * | 1991-02-26 | 1994-06-14 | Swift Energy Company | Method for exploring for hydrocarbons utilizing three dimensional modeling of thermal anomalies |
JP2956800B2 (ja) * | 1991-09-19 | 1999-10-04 | 株式会社日立製作所 | 連立一次方程式に関する計算装置 |
FR2801710B1 (fr) * | 1999-11-29 | 2002-05-03 | Inst Francais Du Petrole | Methode pour generer un maillage hybride permettant de modeliser une formation heterogene traversee par un ou plusieurs puits |
AU2001270265A1 (en) * | 2000-06-29 | 2002-01-14 | Object Reservoir, Inc. | Method and system for coordinate transformation to model radial flow near a singularity |
CN100489558C (zh) | 2004-06-07 | 2009-05-20 | 埃克森美孚上游研究公司 | 用于求解隐式储层仿真矩阵的方法 |
US7526418B2 (en) * | 2004-08-12 | 2009-04-28 | Saudi Arabian Oil Company | Highly-parallel, implicit compositional reservoir simulator for multi-million-cell models |
US7596480B2 (en) | 2005-04-14 | 2009-09-29 | Saudi Arabian Oil Company | Solution method and apparatus for large-scale simulation of layered formations |
US7516056B2 (en) | 2005-04-26 | 2009-04-07 | Schlumberger Technology Corporation | Apparatus, method and system for improved reservoir simulation using a multiplicative overlapping Schwarz preconditioning for adaptive implicit linear systems |
FR2918179B1 (fr) * | 2007-06-29 | 2009-10-09 | Inst Francais Du Petrole | Methode pour estimer la permeabilite d'un reseau de fractures a partir d'une analyse de connectivite |
GB2455077A (en) | 2007-11-27 | 2009-06-03 | Polyhedron Software Ltd | Estimating the state of a physical system using generalized nested factorisation |
CN101896690B (zh) | 2007-12-13 | 2015-02-18 | 埃克森美孚上游研究公司 | 使用非结构化网格的储层模拟上的并行自适应数据分区 |
US8190414B2 (en) * | 2008-03-26 | 2012-05-29 | Exxonmobil Upstream Research Company | Modeling of hydrocarbon reservoirs containing subsurface features |
BRPI0919456A2 (pt) * | 2008-09-30 | 2015-12-22 | Exxonmobil Upstream Res Co | método para modelar escoamento de fluido em um reservatório de hidrocarboneto |
BRPI0919457A2 (pt) * | 2008-09-30 | 2015-12-01 | Exxonmobil Upstream Res Co | método para simular escoamento de fluido em um reservatório de hidrocarboneto |
GB2477901A (en) * | 2008-12-03 | 2011-08-17 | Chevron Usa Inc | System and method for predicting fluid flow characteristics within fractured subsurface reservoirs |
RU2549147C2 (ru) * | 2009-05-07 | 2015-04-20 | Сауди Арабиан Ойл Компани | Системы, компьютерно-реализуемые способы и компьютерно-считываемые программные продукты для расчета приближенного давления дренирования скважины для имитатора коллектора |
AU2010282773B2 (en) * | 2009-08-14 | 2014-10-23 | Bp Corporation North America Inc. | Reservoir architecture and connectivity analysis |
WO2011028197A1 (en) * | 2009-09-02 | 2011-03-10 | Landmark Graphics Corporation | System and method of hydrocarbon formation modeling |
WO2011056666A2 (en) * | 2009-10-28 | 2011-05-12 | Chevron U.S.A. Inc. | Multiscale finite volume method for reservoir simulation |
CA2785569A1 (en) * | 2010-02-02 | 2011-08-11 | Hector Klie | Multilevel percolation aggregation solver for petroleum reservoir simulations |
US8463586B2 (en) | 2010-06-22 | 2013-06-11 | Saudi Arabian Oil Company | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids |
US9396162B2 (en) | 2010-07-22 | 2016-07-19 | John APPLEYARD | Method and apparatus for estimating the state of a system |
BR112012031776A2 (pt) * | 2010-07-26 | 2016-11-01 | Exxonmobil Upstream Res Co | método e sistema para simulação multinível paralela |
US9187984B2 (en) * | 2010-07-29 | 2015-11-17 | Exxonmobil Upstream Research Company | Methods and systems for machine-learning based simulation of flow |
US8433551B2 (en) | 2010-11-29 | 2013-04-30 | Saudi Arabian Oil Company | Machine, computer program product and method to carry out parallel reservoir simulation |
CN103097657B (zh) * | 2010-09-07 | 2015-05-27 | 沙特阿拉伯石油公司 | 用于模拟由储集层的网格界定的多个油气井的生产特性的机器及计算机实施方法 |
US8386227B2 (en) | 2010-09-07 | 2013-02-26 | Saudi Arabian Oil Company | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation |
KR101908278B1 (ko) * | 2010-09-27 | 2018-10-17 | 엑손모빌 업스트림 리서치 캄파니 | 풀 파동장 반전을 위한 실제적인 해결책으로서의 동시 소스 인코딩 및 소스 분리 |
US8775144B2 (en) | 2010-12-13 | 2014-07-08 | Chevron U.S.A. Inc. | Constrained pressure residual preconditioner for efficient solution of the adjoint equation |
US8583411B2 (en) | 2011-01-10 | 2013-11-12 | Saudi Arabian Oil Company | Scalable simulation of multiphase flow in a fractured subterranean reservoir as multiple interacting continua |
US8731891B2 (en) * | 2011-07-28 | 2014-05-20 | Saudi Arabian Oil Company | Cluster 3D petrophysical uncertainty modeling |
WO2013059279A1 (en) * | 2011-10-18 | 2013-04-25 | Saudi Arabian Oil Company | Reservoir modeling with 4d saturation models and simulation models |
EP2761529B1 (en) * | 2011-12-16 | 2023-01-18 | Landmark Graphics Corporation | System and method for flexible and efficient simulation of varying fracture density in a reservoir simulator |
US9208268B2 (en) | 2012-02-14 | 2015-12-08 | Saudi Arabian Oil Company | Giga-cell linear solver method and apparatus for massive parallel reservoir simulation |
US8977524B2 (en) * | 2012-03-06 | 2015-03-10 | Siemens Aktiengesellschaft | Interior point method for reformulated optimal power flow model |
US9753180B2 (en) | 2012-03-28 | 2017-09-05 | Exxonmobil Upstream Research Company | Method for multiphase flow upscaling |
US9858369B2 (en) * | 2012-10-18 | 2018-01-02 | Helic, Inc. | Large-scale power grid analysis on parallel architectures |
CN104533370B (zh) * | 2014-11-06 | 2017-03-15 | 中国石油大学(北京) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 |
-
2015
- 2015-11-23 US US14/948,605 patent/US10242136B2/en active Active
-
2016
- 2016-05-17 EP EP18178356.4A patent/EP3399141A1/en active Pending
- 2016-05-17 JP JP2017560232A patent/JP2018518613A/ja active Pending
- 2016-05-17 CA CA2986034A patent/CA2986034C/en active Active
- 2016-05-17 EP EP16725707.0A patent/EP3298232B1/en active Active
- 2016-05-17 WO PCT/US2016/032819 patent/WO2016187175A1/en active Application Filing
- 2016-05-17 CN CN201680043093.3A patent/CN107849910B/zh not_active Expired - Fee Related
- 2016-05-17 CN CN201810448806.5A patent/CN108691521B/zh not_active Expired - Fee Related
-
2017
- 2017-01-23 US US15/412,507 patent/US10229237B2/en active Active
-
2018
- 2018-02-06 JP JP2018019074A patent/JP6824205B2/ja active Active
-
2019
- 2019-01-31 US US16/263,689 patent/US10769326B2/en active Active
- 2019-01-31 US US16/263,701 patent/US10762258B2/en active Active
-
2020
- 2020-04-23 JP JP2020076641A patent/JP6887042B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
WO2016187175A1 (en) | 2016-11-24 |
US10229237B2 (en) | 2019-03-12 |
US20190163856A1 (en) | 2019-05-30 |
JP6887042B2 (ja) | 2021-06-16 |
EP3298232A1 (en) | 2018-03-28 |
CN107849910A (zh) | 2018-03-27 |
US20170147728A1 (en) | 2017-05-25 |
US10242136B2 (en) | 2019-03-26 |
US20160341015A1 (en) | 2016-11-24 |
EP3399141A1 (en) | 2018-11-07 |
JP2018119967A (ja) | 2018-08-02 |
CN107849910B (zh) | 2021-05-28 |
CA2986034C (en) | 2021-10-26 |
JP2018518613A (ja) | 2018-07-12 |
US10769326B2 (en) | 2020-09-08 |
EP3298232B1 (en) | 2019-03-06 |
CA2986034A1 (en) | 2016-11-24 |
CN108691521A (zh) | 2018-10-23 |
CN108691521B (zh) | 2020-11-06 |
US20190163855A1 (en) | 2019-05-30 |
JP2020143572A (ja) | 2020-09-10 |
US10762258B2 (en) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6887042B2 (ja) | コンピュータ実装方法、データ処理システム、およびデータストレージデバイス | |
US8433551B2 (en) | Machine, computer program product and method to carry out parallel reservoir simulation | |
US8386227B2 (en) | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation | |
US10839114B2 (en) | Method and system for stable and efficient reservoir simulation using stability proxies | |
US6928399B1 (en) | Method and program for simulating a physical system using object-oriented programming | |
EP2614460B1 (en) | Machine, computer program product and method to generate unstructured grids and carry out parallel reservoir simulation | |
CN109964151B (zh) | 利用加速的含水层计算的并行储层模拟 | |
US10534877B2 (en) | Adaptive multiscale multi-fidelity reservoir simulation | |
NO339000B1 (no) | Framgangsmåte og datasystem for simulering av lagdelte grunnformasjoner | |
NO20121123A1 (no) | Forbehandlingsenhet for reservoarsimulering | |
NO20130073A1 (no) | Fremgangsmate og system for stabilisering av formuleringsmetoder | |
AU2015409650A1 (en) | Simulating fractured reservoirs using multiple meshes | |
CN112668251A (zh) | 一种模拟裂缝性油藏多重嵌套介质内流体流动的方法 | |
CN112752894B (zh) | 具有用于非对角占优的不定系数矩阵的压力求解器的储层模拟 | |
AU2011271175B2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
CN112699617A (zh) | 一种模拟裂缝性油藏多重嵌套介质内流体流动的装置 | |
de Menezes | Development of a multipurpose reservoir simulator based on a plugin architecture |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180206 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180426 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20181114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181120 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20190219 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190405 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20191001 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191205 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200121 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200204 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200225 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200326 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20200825 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200827 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20201222 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210112 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6824205 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |