(57)【特許請求の範囲】
【請求項１】
流体が圧入又は産出される坑井内の複数の計算点における前記流体の圧力を、積分差分法を用いて計算する流体圧力計算方法であって、
前記坑井が存在する三次元格子モデルとは別に、一つの前記計算点に一つの格子が対応する一次元格子モデルを設定する一次元格子設定ステップと、
前記計算点が存在する前記三次元格子モデルの１格子の格子点圧力と、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との距離と、前記計算点における流体の相の流量と、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との接続面積と、前記流体の絶対浸透率と、前記流体の相の相対浸透率と、に基づいて、有限体積法によって前記計算点における前記流体の圧力を計算する圧力計算ステップと、
を含むことを特徴とする流体圧力計算方法。
【請求項２】
前記圧力計算ステップにおいて、
前記三次元格子モデルの層α、前記計算点が存在する前記三次元格子モデルの１格子の格子点圧力Ｐｉ、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との接続面積Ａ、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との距離Ｌ、流体の絶対浸透率ｋ、流体の相βの相対浸透率ｋｒβ、及び、前記計算点における流体の相βの流量Ｑβに基づいて、式（１）
【数１】
0006768341.tif 000010
によって、前記計算点における前記流体の圧力Ｐｂｈを計算する
ことを特徴とする請求項１に記載の流体圧力計算方法。
【請求項３】
前記圧力計算ステップにおいて、
前記層αに対する井戸指数ＷＩα及び格子内の流動性λαに基づいて、
Ｌ＝１［ｍ］
ｋ＝１［ｍ２］
Ａ＝ＷＩα
と設定することによって、式（２）
【数２】
0006768341.tif 000011
によって、前記計算点における前記流体の圧力Ｐｂｈを計算する
ことを特徴とする請求項２に記載の流体圧力計算方法。
【請求項４】
前記流体の圧力Ｐｂｈは、前記坑井に圧入される前記流体の圧力であり、
前記井戸指数ＷＩαは、
坑井と直交する方向の有効浸透率ｋ、坑井方向の格子厚さｈ、格子内の総流動性λＴ、坑井の半径ｒｗ、三次元格子モデルにおける１格子の有効半径ｒｅ、スキンファクターＳに基づいて、式（３）
【数３】
0006768341.tif 000012
によって計算されている
ことを特徴とする請求項３に記載の流体圧力計算方法。
【請求項５】
前記流体の圧力Ｐｂｈは、前記坑井に圧入される前記流体の圧力であり、
前記流動性λαは、
前記相対浸透率ｋｒβ及び流体の相βの粘性係数μβに基づいて、式（４）
【数４】
0006768341.tif 000013
によって計算されている
ことを特徴とする請求項３又は請求項４に記載の流体圧力計算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地層中に掘削された坑井に圧入される、又は、坑井から生産される流体の圧力を計算する技術に関する。
【背景技術】
【０００２】
従来、地層中に掘削された坑井に二酸化炭素を圧入して貯留したり、地層中に掘削された坑井から石油や天然ガスを生産したりすることが行われている。
坑井に圧入され、又は、坑井から生産されたりする流体の数値解析には、坑井モデルが用いられる。
坑井モデルは、坑井が掘削された地層に三次元格子モデルを設定し、数値解析対象となる格子の圧力と、坑井内の局所的な圧力とを関係づけるモデルである。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
従来の流体圧力を計算するシミュレータには、前記した坑井モデルが予め組み込まれている。すなわち、坑井モデルを用いて坑井内における流体の圧力を計算するには、坑井モデルが予め組み込まれたプログラムを用いる必要があった。
【０００４】
本発明は、前記した事情に鑑みて創案されたものであり、坑井モデルが組み込まれていない汎用の流体シミュレータにおいて、流体の圧力を好適に計算することが可能な流体圧力計算方法を提供することを課題とする。
【課題を解決するための手段】
【０００５】
前記課題を解決するため、本発明は、流体が圧入又は産出される坑井内の複数の計算点における前記流体の圧力を、積分差分法を用いて計算する流体圧力計算方法であって、前記坑井が存在する三次元格子モデルとは別に、一つの前記計算点に一つの格子が対応する一次元格子モデルを設定する一次元格子設定ステップと、前記計算点が存在する前記三次元格子モデルの１格子の格子点圧力と、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との距離と、前記計算点における流体の相の流量と、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との接続面積と、前記流体の絶対浸透率と、前記流体の相の相対浸透率と、に基づいて、有限体積法によって前記計算点における前記流体の圧力を計算する圧力計算ステップと、を含むことを特徴とする。
【０００６】
かかる構成によると、坑井モデルが組み込まれていない汎用の流体シミュレータに対して一次元格子モデルを設定することで模擬的な坑井モデルを適用し、流体の圧力を好適に計算することができる。
【０００７】
前記流体圧力計算方法は、前記圧力計算ステップにおいて、前記三次元格子モデルの層α、前記計算点が存在する前記三次元格子モデルの１格子の格子点圧力Ｐｉ、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との接続面積Ａ、前記計算点が存在する前記三次元格子モデルの１格子と当該１格子に対応する前記一次元格子モデルの１格子との距離Ｌ、流体の絶対浸透率ｋ、流体の相ｊの相対浸透率ｋｒβ、及び、前記計算点における流体の相βの流量Ｑβに基づいて、式（１）
【数１】
0006768341.tif 000002
によって、前記計算点における前記流体の圧力Ｐｂｈを計算する構成であってもよい。
【０００８】
前記流体圧力計算方法は、前記圧力計算ステップにおいて、前記層αに対する井戸指数ＷＩα及び格子内の流動性λαに基づいて、
Ｌ＝１［ｍ］
ｋ＝１［ｍ２］
Ａ＝ＷＩα
と設定することによって、式（２）
【数２】
0006768341.tif 000003
によって、前記計算点における前記流体の圧力Ｐｂｈを計算する構成であってもよい。
【０００９】
前記流体の圧力Ｐｂｈは、前記坑井に圧入される前記流体の圧力であり、前記井戸指数ＷＩαは、坑井と直交する方向の有効浸透率ｋ、坑井方向の格子厚さｈ、格子内の総流動性λＴ、坑井の半径ｒｗ、三次元格子モデルにおける１格子の有効半径ｒｅ、スキンファクターＳに基づいて、式（３）
【数３】
0006768341.tif 000004
によって計算されている構成であってもよい。
【００１０】
前記流体の圧力Ｐｂｈは、前記坑井に圧入される前記流体の圧力であり、前記流動性λαは、前記相対浸透率ｋｒβ及び流体の相βの粘性係数μβに基づいて、式（４）
【数４】
0006768341.tif 000005
によって計算されている構成であってもよい。
【発明の効果】
【００１１】
本発明によると、坑井モデルが組み込まれていない汎用の流体シミュレータに対して模擬的な坑井モデルを適用し、流体の圧力を好適に計算することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態に係る流体圧力計算装置を示すブロック図である。
【図２】三次元格子モデル及び一次元格子モデルを示す模式図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施形態について、坑井内に圧入される流体の圧力を計算する場合を例にとり、適宜図面を参照しながら説明する。
【００１４】
図１に示すように、本発明の実施形態に係る流体圧力計算装置1は、操作部１０と、表示部２０と、これらを制御する制御部３０と、を備える。
第一の実施形態において、流体圧力計算装置１は、流体（例えば、二酸化炭素）が圧入される坑井内の複数の計算点における流体の圧力を、積分差分法（有限体積法ともいう）を用いて計算する。
【００１５】
＜操作部＞
操作部１０は、キーボード、マウス等によって構成されている。操作部１０は、ユーザによる当該操作部１０の操作内容に基づく入力信号を制御部３０へ出力する。
【００１６】
＜表示部＞
表示部２０は、ディスプレイ等によって構成されている。表示部２０は、制御部３０から出力された表示用信号を取得し、取得された表示用信号に基づいて画面表示を行う。
【００１７】
＜制御部＞
制御部３０は、図示しないＣＰＵ（Central Processing Unit）、ＲＯＭ（Read-Only Memory）、ＲＡＭ（Random Access Memory）、入出力回路等によって構成されている。
制御部３０は、流体が圧入される坑井内の複数の計算点における流体の圧力を、積分差分法を用いて計算する
制御部３０は、機能部として、記憶部３１と、数値解析部３２と、一次元格子モデル設定設定部３３と、圧力計算部３４と、を備える。
【００１８】
＜記憶部＞
記憶部３１には、数値解析部３２及び圧力計算部３４によって用いられる数値として、以下に示す定数が記憶されている。かかる定数は、ユーザによる操作部１０の操作内容に基づいて、記憶部３１に記憶される。
（１）ＷＩα（Well injectivity Index）［−］：三次元格子モデル２（図２参照）の層αに対する流体の相ｊの井戸指数
（２）ｋ［ｍ２］：坑井３と直交する方向における、流体の浸透率（絶対浸透率）
（３）ｋｒβ［−］：坑井３と直交する方向における、流体の相βの相対浸透率
（４）ｈ［ｍ］：坑井３が延びる方向における三次元格子モデル２の格子２ａの厚さ
（５）ｒｗ［ｍ］：坑井３の半径
（６）ｒｅ［ｍ］：三次元格子モデル２における１格子２ａの有効半径（立方体又は直方体である１格子２ａを円柱であるとみなした場合の半径であって、ｒｅ＝ＣＣ（Ａｉ／ｗｆｒａｃ）０．５によって計算される）
（７）Ｓ［−］：スキンファクター（坑井３の壁面に生じる泥の膜の影響）
（８）λα［Ｐａ−１・ｓ−１］：格子２ａ内の流動性
（９）μβ［Ｐａ・ｓ］：流体の相βの粘性係数
層α（α＝１，２，３，…）は、三次元格子モデル２の格子に対して、鉛直方向に順に付与された通し番号である。同一水平面内の（同じ深さに位置する）格子は、同じ層αに属する。
ｉ（ｉ＝１，２，３，…）は、三次元格子モデル２の全ての格子に対して、個別に付与された通し番号である。
本実施形態において、相βは、ｗ（水、海水）、ｇ（ガス、二酸化炭素）、ｏ（油）等といった物質としての違いと、同一物質における固相、液相、気相といった状態の違いと、の両方を含む。すなわち、複数の流体が界面を形成して区分されつつ共存する場合に、これらは別の相として区別される。
また、各浸透率ｋ，ｋｒβは、深さ方向によって異なる値となる。
【００１９】
＜数値解析部＞
数値解析部３２は、坑井モデルが組み込まれていない汎用の流体シミュレータである。
数値解析部３２は、坑井３が掘削された地層に三次元格子モデル２を模擬した三次元格子モデル２について積分差分法（有限体積法ともいう）を用いた数値解析を行うことによって、三次元格子モデル２内の複数の格子２ａの圧力Ｐｉ、圧入流量Ｑβ等を計算する（図２参照）。
【００２０】
≪三次元格子モデル≫
図２に示すように、三次元格子モデル２は、立方体又は直方体形状を呈する複数の格子２ａを三次元方向に配列したモデルである。図２は、三次元格子モデル２を平面的に表現したものであり、実際には、紙面と直交する方向にも格子２ａが配列されている。
かかる三次元格子モデル２には、流体が圧入される坑井３が設定されている。
坑井３は、海底の地層に掘削されており、流体としての二酸化炭素は、坑井３内の水（海水）に圧入される。
坑井３には、深さ方向に配列された複数の計算点３ａが設定されている。すなわち、複数の計算点３ａは、それぞれ異なる深さに設定されている。
【００２１】
図１に戻り、本実施形態において、数値解析部３２は、三次元格子モデル２において、積分差分法（有限体積法ともいう）を用いた公知の数値解析手法によって、各格子２ａの中央部の圧力Ｐｉ、坑井３の計算点３ａにおける流体の相βの圧入流量Ｑβを計算し、計算結果を記憶部３１に記憶させる。
【００２２】
＜一次元格子モデル設定部＞
一次元格子モデル設定部３３は、坑井３が存在する三次元格子モデル２とは別に、一つの計算点３ａに一つの格子４ａが対応する一次元格子モデル４を設定する（一次元格子設定ステップ）。
一次元格子モデル設定部３３は、自動的、又は、ユーザによる操作部１０の操作に基づいて、一次元格子モデル４を設定する（図２参照）。
設定された一次元格子モデル４は、流体の圧力Ｐｂｈを計算するための模擬的な坑井モデルとして機能する。
【００２３】
≪一次元格子モデル≫
図２に示すように、一次元格子モデル４は、立方体又は直方体形状を呈する複数の格子４ａを一次元方向（深さ方向）に配列したモデルである。
一次元格子モデル４は、複数の格子４ａが鉛直方向に配列されるように設定されている。
一つの格子４ａは、坑井３の一つの計算点３ａに対応している。また、一つの格子４ａと、対応する計算点３ａが含まれる一つの格子２ａとは、同じ深さに設定されている。すなわち、互いに対応する一つの格子４ａと一つの格子２ａとを結ぶ線分（接続線）は、水平方向に延びている。ここで、一次元格子モデル４の複数の格子４ａは、平面視で同じ形状を呈するとともに、深さ方向の寸法は、適宜異なる値に設定されている。
【００２４】
ここで、一次元格子モデル設定部３３は、一次元格子モデル４の格子４ａの空間座標を設定する必要はなく、以下に示すＡ，Ｌ、及び格子２ａ，４ａの接続線の重力方向からの角度を設定すればよい。
Ａ：計算点３ａが存在する三次元格子モデル２の１格子２ａと当該１格子２ａに対応する一次元格子モデル４の１格子４ａとの接続面積（ここでは、Ａ２）
Ｌ：計算点３ａが存在する三次元格子モデル２の１格子２ａと当該１格子２ａに対応する一次元格子モデル４の１格子４ａとの距離
本実施形態において、互いに接続された格子２ａ，４ａの接続線の重力方向からの角度は、９０度である。
【００２５】
＜圧力計算部＞
圧力計算部３４は、計算点３ａが存在する三次元格子モデル２の１格子２ａの格子点圧力Ｐｉと、計算点３ａが存在する三次元格子モデル２の１格子２ａと当該１格子２ａに対応する一次元格子モデル４の１格子４ａとの距離Ｌと、計算点２ａにおける流体の相βの流量Ｑβと、に基づいて、積分差分法（有限体積法ともいう）を用いて、計算点における流体の圧力を計算する（圧力計算ステップ）。
圧力計算部３４は、記憶部３１に記憶された数値を用いて、下記式（１）によって計算点３ａにおける流体の圧力Ｐｂｈを計算する。
【数５】
0006768341.tif 000006
すなわち、圧力計算部３４は、一次元格子モデル４の格子４ａの中心における流体の圧力を求めることによって、坑井３内の計算点３ａにおける流体の圧力Ｐｂｈを計算することができる。
【００２６】
本実施形態において、圧力計算部３４は、
Ｌ＝１［ｍ］
ｋ＝１［ｍ２］
【数６】
0006768341.tif 000007
ただし、β＝ｇ（坑井３内の液体に圧入される気体、例えば二酸化炭素），ｗ（気体が圧入される液体、例えば、水、海水）
Ａ＝ＷＩα
【数７】
0006768341.tif 000008
と設定することによって、下記式（２）によって、計算点３ａにおける流体の圧力Ｐｂｈを計算する。
【数８】
0006768341.tif 000009
ここで、λαは、格子２ａ内の流動性であって、流体の動きやすさを表し、格子２ａ内における流体の飽和度に基づいて計算され、予め記憶されている。λαは、層αごとに異なる値となる。
【００２７】
また、Ａは、三次元格子モデル２内の格子２ａにおける一次元格子モデル４内の格子４ａと対向する面の面積Ａ１（記憶部３１内に記憶）であってもよく、一次元格子モデル４内の格子４ａにおける三次元格子モデル２内の格子２ａと対向する面の面積Ａ２（一次元格子設定部３３によって設定）であってもよい。
【００２８】
本発明の実施形態に係る流体圧力計算装置１及び流体圧力計算装置１による流体圧力計算方法は、坑井３が存在する三次元格子モデル２とは別に一次元格子モデル４を設定し、これらの格子モデル２，４による積分差分法を用いて坑井３内の計算点３ａにおける流体の圧力Ｐｂｈを計算するので、坑井モデルが組み込まれていない汎用の流体シミュレータに対して模擬的な坑井モデルを適用し、流体の圧力Ｐｂｈを好適に計算することができる。
【００２９】
以上、本発明の実施形態について説明したが、本発明は前記実施形態に限定されず、本発明の要旨を逸脱しない範囲で適宜変更可能である。例えば、圧力の計算対象となる流体は、圧入される二酸化炭素、生産される石油、天然ガス等に限定されない。また、本発明は、坑井３から産出される流体の圧力を、前記した式（１）（２）を用いて計算する流体圧力計算装置１及び流体圧力計算方法としても具現化可能である。
【符号の説明】
【００３０】
１      流体圧力計算装置
２      三次元格子モデル
２ａ    格子
３      坑井
３ａ    計算点
４      一次元格子モデル
４ａ    格子
３３    一次元格子モデル設定部
３４    圧力計算部
【図１】
0006768341.tif 000014
【図２】
0006768341.tif 000015