JP3837513B2 - 核磁気共鳴法を用いた地層の浸透率の推定方法 - Google Patents
核磁気共鳴法を用いた地層の浸透率の推定方法 Download PDFInfo
- Publication number
- JP3837513B2 JP3837513B2 JP2002220536A JP2002220536A JP3837513B2 JP 3837513 B2 JP3837513 B2 JP 3837513B2 JP 2002220536 A JP2002220536 A JP 2002220536A JP 2002220536 A JP2002220536 A JP 2002220536A JP 3837513 B2 JP3837513 B2 JP 3837513B2
- Authority
- JP
- Japan
- Prior art keywords
- formation
- permeability
- self
- pore
- equation
- 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Measuring Magnetic Variables (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、土木工学、地盤工学、資源探査および物理探査の分野において利用可能な核磁気共鳴法を用いた地層の浸透率の推定方法に関する。
【０００２】
【従来の技術】
従来、地層の浸透率の核磁気共鳴（以下「ＮＭＲ」という。）による推定は、地層中の間隙流体のプロトン緩和データを用いて行っている（Stambaugh,2000,Oil&GasJournal,vol.98,pp.45-52）。
【０００３】
【発明が解決しようとする課題】
しかしながら、従来の方法では、岩石の化学組成がわからないと、正確な浸透率推定はできないという問題があった。これはプロトンの緩和時間が、岩石の化学組成に強く影響されるからである。そのため、化学組成がわからない地層の浸透率を推定するときは、大きな誤差が生じるものであった。
また、Balzarini et al,1996.PhysicaB,vol.226,pp.10-14で,拡散データを用いた浸透率の推定方法が述べられているが、その推定方法は誤ったものであった。
【０００４】
本発明は、核磁気共鳴法を用いて多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測することにより、地層の化学組成がわからくても正確な浸透率を推定可能とすることを目的とする。
【０００５】
【課題を解決するための手段】
上記目的を達成するため、本発明の核磁気共鳴法を用いた地層の浸透率の推定方法は、静磁場磁石、傾斜磁場コイルおよびラジオ波コイルから構成されるセンサーユニットを地層に押し当てて地層中の間隙流体分子の自己拡散係数を計測することにより、地層の浸透率を推定する方法において、浸透率ｋを次の式から求めることを特徴とする核磁気共鳴法を用いた地層の浸透率の推定方法。
【数５】
ただし、φは地層の空隙率、D０はバルク流体分子の自己拡散係数、D∞は空隙サイズよりも十分大きな距離を移動したときの間隙流体分子の自己拡散係数である。また、（S／V）ｐｏｒｅは空隙の比表面積であって、次の式から求める。
【数２】
ただし、Ｄは間隙流体分子の自己拡散係数、ｔは空隙サイズよりも小さな距離を移動したときの間隙流体分子の拡散時間、ｃ１は定数である。
【０００６】
【発明の実施の形態】
以下、本発明による実施の形態を図面に基づき説明する。
図１は、本発明の実施の形態の１例を示したもので、ＮＭＲ計測機器のセンサー部分のみ示してある。
１は、間隙水などの含水素流体の自己拡散係数を計測する対象である流体で満たされた多孔質地層、例えば、岩石、堆積物等である。
２は、センサーユニットであり、静磁場磁石３、傾斜磁場コイル４およびラジオ波コイル５から構成されている。
静磁場磁石３は、原子核が共鳴するように一定の強度をもった磁場を形成するものである。 傾斜磁場コイル４は、流体分子の自己拡散による空間移動を検出するために場所により強さの異なる磁場、すなわち、傾斜磁場を作る。ラジオ波コイル５は、高周波磁場を多孔質地層１に照射し共鳴を起こさせると共に、共鳴の結果、多孔質地層１から発せられる信号を受信するものである。また、図示しないコンピュータにより、傾斜磁場の切り換えや高周波磁場の発生のタイミングを制御したり、受信した信号をもとに計算を行うようになっている。
【０００７】
図１のセンサーユニット２を岩盤の表面やボーリング孔の抗壁の地層１などに押し当てて、地層１の表面のＮＭＲ物性を計測する。その際、岩盤表面から一定の距離、例えば、数cm〜数ｍ奥の位置にあるプロトンに感度を持つよう磁場を調整する。計測対象は、多孔質な地層１中の間隙を満たしている軽水素（１Ｈ）の原子核（プロトン）を含む流体である。たとえば、地下水の水分子中のプロトン、あるいは石油の炭化水素分子中のプロトンが対象核である。
【０００８】
間隙流体は、多孔質岩石中を熱運動によって自己拡散して移動している。その自己拡散係数データをパルス傾斜磁場を用いたＮＭＲ（ＰＦＧ ＮＭＲ）で計測するのが、まず浸透率計測の第一歩である。この自己拡散係数を計測するためには、たとえば、図２で示すパルス系列（ＰＦＧスピンエコー系列）を使用する必要がある。プロトンの歳差運動の回転速度、位相などをラジオ波と傾斜磁場のパルスでコントロールし、エコー時間ＴＥのタイミングで現れるスピンエコー波形を取得する。傾斜磁場強度（Ｇ）の上昇に起因する、エコー波形の高さ（Ｉ）の減少を計測し、それから間隙流体分子の自己拡散係数（Ｄ）を求める。すなわち、
【数１】
Ｉ０は、傾斜磁場強度がゼロのときのスピンエコー強度である。ただし、ｂは次式で定義される量である。
ｂ＝（γＧδ）２（Δ−δ／３） （２）
γは核種に固有の磁気回転比であり、プロトンの場合は2.675×108 rad/Tsである。Δとδは、図２に示すように、それぞれ、傾斜磁場パルスの間隔と持続時間である。ｂは既知量なので、ｂをいろいろ変えてＩ０ で規格化されたスピン・エコーシグナル（Ｉ／Ｉ０）を計測し、最小二乗法で未知量Ｄを決定する。
【０００９】
ＰＦＧ ＮＭＲでは、間隙流体がΔのタイムスケールで拡散できる距離を計測する。Δを変えて拡散係数の変化を調べると、多孔質媒体の比表面積と屈曲度が計測できることを、Mitra,P.P.,P.N.Sen,およびL.M.Schwartzが１９９３年発見した（Physical Review B,47巻,8565-8574頁）。概略は次のとおりである。
多孔質媒体中の間隙流体の拡散係数（Ｄ）という量は、拡散時間・拡散距離に依存する。ある空隙サイズをもつ多孔質媒体内部で間隙流体がランダムウォークしている系を考えると、空隙の中央付近からランダムウォークを開始した流体分子は最初のうちは固体表面にめったに衝突しないのでバルクの流体と同じ自己拡散係数（Ｄ０）をもつ。流体分子の移動距離は、時間が経つにつれて増加する。ランダムウォーカーの根平均二乗変位が空隙サイズとほぼ同じかあるいはより大きくなると固体との衝突が頻繁におこり、流体分子は固体を回り込みながら空隙中を移動していく。この固体部分の幾何学的な障害物効果によって、ＤはＤ０よりすこしだけ低下する。重要なポイントは、その低下量が空隙の比表面積に関係する点である。
【００１０】
【数２】
（Ｓ／Ｖ）poreは空隙の比表面積、ｔは拡散時間、ｃ１は定数、πは円周率である。ＰＦＧ ＮＭＲではΔが拡散時間なので、Δをいろいろ変えてＤの低下量を計測し、（3）式によって比表面積をもとめることができる。ＮＭＲによる世界初の実測例が、Latour, L.L., P.P. Mitra, R.L. Kleinberg, C.H. Sotak著の論文で報告されている（Journal of Magnetic Resonance, A101巻, 342-346頁, 1993年）。
【００１１】
さらに、拡散時間Δを大きくとると、流体分子は空隙サイズよりも十分大きな距離を移動できるようになり、その結果、自己拡散係数はある値, Ｄ ∞, に漸近する。
Ｄ＝Ｄ∞ ａｓ ｔ → ∞ （４）
Ｄ０／Ｄ∞ という比は,多孔質媒体の空隙構造の屈曲度(tortuosity)と呼ばれる量である（たとえばLatour, L.L., R.L. Kleinberg, P.P. Mitra and C.H. Sotak, Journal of Magnetic Resonance, A112巻, 83-91頁, 1995年）。また、地層の空隙率は、従来のＮＭＲ検層技術等で確立しているように、ＣＰＭＧ時系列データの解析から推定できる（たとえば、Straley, C., Rossini, D., Vinegar, H., Tutunjian, P., and Morriss, C.著The Log Analyst, ３−４月号、84−９４頁1997年）。
このように、ＮＭＲを用いれば、地層の空隙率、比表面積、屈曲度を計測できる。
【００１２】
多孔質地層の浸透率は、空隙率、比表面積、屈曲度がわかればカルマン・コゼニー式を用いて、推定できる。ここが、本発明のオリジナル部分である。以下に、その推定方法を述べる。一般に、多孔質媒質の浸透率（ｋ）は、カルマン・コゼニー式で近似できる場合が多い（たとえばGueguen, Y., とV. Palciauskas, Introduction to the Physics of Rocks, Princeton University Press, Princeton, 1994年）。すなわち、
【数３】
φは空隙率、（Ｓ／Ｖ）solid＝ （Ｓ／Ｖ）poreφ／(１ - φ)は固体部分の比表面積、αは１のオーダーの定数、ｌは実在の曲がりくねったパイプ状の空隙の長さ、ｌ０は直線的にのびた仮想的な空隙の長さである。(５)式は、多孔質媒体の複雑な空隙形状を、曲がりくねった円筒状のパイプとして近似し、ポアズイユ流モデルで浸透率を計算したものである。
【００１３】
さて、(５)式において、（ｌ／ｌ０） という比はパイプの伸び率を表す尺度である。本来は長さがｌ０というパイプが多孔質岩石の複雑な空隙構造のためにｌに伸びたことを意味する。このパイプを満たしている流体分子は、自己拡散によって、曲がりくねったパイプに沿って常時移動している。したがって、自己拡散係数の大小は、曲がりくねったパイプ構造の複雑性を反映している。ところで、拡散現象はランダムウォーク過程であり、数学的には放物型の偏微分方程式に分類される。放物型においては、ある距離Ｌを拡散移動するのに必要な時間はＬの２乗に比例する。したがって、ｌの距離を拡散するのに必要な時間はｌ０の距離を移動するのに必要な時間の（ｌ／ｌ０）２倍かかる。この比は、拡散係数の比（屈曲度）そのものである。
したがって、(５)式における（ｌ／ｌ０）２ という比は、（Ｄ０／Ｄ∞）という比で置き換えることができる。
【００１４】
【数４】
（６）式の右辺において、定数αは、１のオーダーの補正係数なので、α＝１と近似する。その結果、(Ｓ／Ｖ)poreに関して（７）式を得る。この（７）式から、拡散データを用いて浸透率を推定する。
【数５】
なお、Balzarini ほか (1996年) Physica B, 226巻, 10-14頁で、拡散データを用いた浸透率の推定方法が述べられているが、そこで導出された式は間違っている。(Ｓ／Ｖ)solidと(Ｓ／Ｖ)poreの混同が誤りの原因である。
【００１５】
（７）式の右辺は、すべてＮＭＲで計測可能な点が特徴である。空隙率φは、従来のＮＭＲ検層技術等で確立しているように、ＣＰＭＧ時系列データの解析から推定できる（たとえば、StraleyほかThe Log Analyst, ３−４月号、84−９４頁1997年）。空隙の比表面積(S/V)poreは、（３）式で示したように、比較的短いΔの条件のもとでＰＦＧ ＮＭＲ法を用いてＤ／Ｄ０ の変化を計測すれば、得られる。Ｄ0は間隙流体の種類（水か石油かなど）と温度がわかれば既存の文献等から知ることができる。Ｄ∞は、（４）式で示したように、比較的長いΔの条件の下でＰＦＧ ＮＭＲ法を用いてDを計測すれば、得られる。
【００１６】
（７）式は、浸透率（ｋ）を拡散データを用いて推定するというオリジナルな方法である。ＮＭＲで地層の浸透率を推定する場合、これまで２つの方法（Timur-Coates モデルとKenyon モデル）が使われてきた（たとえばStambaugh の解説論文(2000年) Oil&Gas Journal, 98巻, 45-52頁を参照）。その２つのモデルは、地層中の間隙流体のプロトン緩和データを用いてｋの推定を行っている。すなわち、ＣＰＭＧ法を用いてプロトンの横緩和現象を計測し、プロトン横緩和時間のヒストグラムから空隙サイズ分布を推定して浸透率を計算している。本発明のような傾斜磁場コイルを用いたＰＦＧ ＮＭＲを実施しているわけではない。プロトン緩和時間を用いた従来の浸透率推定方法には重大な欠点がある。すなわち、プロトン緩和時間が岩石の化学組成（鉄の含有量など）に依存する点である。鉄イオンなどの常磁性物質はプロトンの緩和を促進するので、たとえ同じ空隙サイズでも鉄含有量の多い多孔質岩石は緩和時間が短くなり、少ない岩石では緩和時間が長くなる。したがって、同じ浸透率・同じ空隙構造の地層でも、鉄分含有量が異なるとプロトン緩和時間も異なる。したがって、従来の方法で浸透率を正しく推定するには、岩石の化学組成の情報が必須である。しかし、たとえばＮＭＲ検層でボーリング孔のスキャンを行う場合、地下深部の岩石の種類に関する情報は、獲得困難な場合がしばしばあり、これが浸透率の信頼性を落としている一因になっている。一方、（７）式では、緩和時間の情報は使っていないので、岩石の化学組成の情報は必要ではない。したがって、岩石の化学組成の情報がなくても、浸透率を精度よく推定できる。
【００１７】
【実施例】
(７)式の有効性を実証するために、以下のランダムウォークのコンピューターシミュレーションを行った。サンプルは、直径２．１１ｍｍのガラスビーズ堆積物である（粒度のそろった砂質堆積層のモデル物質として選んだ）。このサンプルの浸透率を（７）式で推定できることを実証する。ビーズ堆積物の空隙率は34％である。ガラスビーズをプラスチック容器に充填して、マイクロフォーカスＸ線ＣＴ（コンピューター断層）画像を取得した。その３次元デジタル画像データを図３に示す。データサイズは、２５６×２５６×２５６画素で、各画素は一辺が０．０５３ｍｍの立方体である。このデジタル画像中の空隙部分を水が満たしていると仮定し、その水分子の自己拡散を離散格子上のランダムウォークでシミュレーションした。ランダムウォークプログラムは、インターネットで無料公開されているマセマティカプログラムＲＷ３Ｄ．ｍ（Watanabe and Nakashiama、 Computers and Geosciences, 28巻, 583-586頁,2002年）を使用した。一般に、ランダムウォークする物質の拡散係数は（８）式で与えられる。
【数６】
【００１８】
<ｒ2>は、ランダムウォーカーの平均自乗変位であり、今回の離散格子上のコンピューターシミュレーションでは（９）式で計算される。
【数７】
ｎ はランダムウォークの総試行回数, ｒｉ はi番目の試行における時刻τでのウォーカーの位置とτ＝０の時のウォーカーの位置との距離である。τは、0, 1, 2, 3, …というように増加する時間ステップ数（無次元）である。
【００１９】
Ｄ0はビーズのない自由空間（バルクの水）中の拡散係数であるが、それは(１０)式で与えられる（たとえば、P. T. Callaghan著Principles of Nuclear Magnetic Resonance Microscopy, オックスフォード大学出版、オックスフォード、１９９１年）。
<ｒ2>＝６Ｄ0ｔ＝ａ２τ （１０）
aは、３次元デジタル単純立方格子の格子定数であり、この場合はＣＴ画像の画素サイズの値（０．０５３ｍｍ）である。（４）（８）（１０）式によれば、ビーズ堆積物の屈曲度（Ｄ０／Ｄ∞）はランダムウォークシミュレーションでは、（１１）式で計算できる。
【数８】
ビーズ堆積物の比表面積 (Ｓ／Ｖ)poreは （３）（８）（１０）式に基づいて（１２）式で計算できる。
【数９】
ただしｃ２は定数である。
【００２０】
ランダムウォークシミュレーションの軌跡の一例を図４に示す。比較的短時間の拡散挙動を図５ (n = 15000)に、比較的長時間の拡散挙動を図６(n = 300)に示す。図５のデータを、（１２）式を用いて最小二乗解析した結果を図７に示す。図７によって(Ｓ／Ｖ)poreは4.5×103 1/mという数値が得られた。一方、直径2.11mm、空隙率34%の均一球の堆積物は(Ｓ／Ｖ)pore = 4.8×103 1/mになるはずである。したがって、今回のランダムウォークシミュレーションでは、誤差６％で比表面積を正しく推定することができた。また、図６を（１１）式を用いて解析すると、Ｄ０／Ｄ∞は平均二乗変位の傾きの比として計算できるので 0.002809/0.001752 = 1.60という値を得た。一方、Latour, L.L., P.P. Mitra, R.L. Kleinberg, C.H. Sotak著の論文によると、均一球の堆積物についてはD０／D∞ = 1.62である（Journal of Magnetic Resonance, A101巻, 342-346頁, 1993年）。したがって、今回のランダムウォークシミュレーションでは、誤差１％で屈曲度（Ｄ０／Ｄ∞）を推定することができた。
【００２１】
以上の結果（φ = 0.34, (Ｓ／Ｖ)pore=4.5×103 1/m, Ｄ０／Ｄ∞ = 1.60）を、(7)式に代入すると、 ｋ = 1.0×10-8 m2という値を得る。一方、このビーズ堆積物の浸透率を通常の室内透水試験方法で実測するとｋ = 0.74×10-8 m2という値を得た。したがって、ランダムウォークシミュレーションによって、誤差35％で浸透率の推定に成功した。このことは、拡散データを用いた浸透率推定の妥当性を例証するものである。さらに、カルマン・コゼニー式は天然の多孔質岩石でも成立する場合が多いので、（７）式の天然の地層への適用の有効性をも示唆するものである。
【００２２】
【発明の効果】
本発明は、以下の効果を奏する。
従来、地層中のプロトンの横緩和現象を計測し、プロトン横緩和時間のヒストグラムから空隙サイズ分布を推定して浸透率を計算しているが、プロトン緩和時間が地層の化学組成（鉄の含有量など）に依存するため、浸透率を正しく推定するには、地層の化学組成の情報が必須である。
これに対し、本発明は、核磁気共鳴法を用いて多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測することにより、地層の化学組成がわからくても正確な浸透率を推定することができる。
【図面の簡単な説明】
【図１】本発明の実施の形態のＮＭＲ計測機器のセンサーユニット部分を説明した図である。
【図２】本発明の実施の形態のＮＭＲ計測機器を使用して多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測する場合に使用するパルス系列の１例である。
【図３】本発明の有効性を実証するためのガラスビーズ堆積物の３次元デジタル画像データを示したものである。
【図４】図３に示すガラスビーズ堆積物中のランダムウォークシミュレーションの軌跡の一例を示したものである。
【図５】ガラスビーズ堆積物中のランダムウォークシミュレーションにおける比較的短時間の拡散挙動を示したものである。
【図６】ガラスビーズ堆積物中のランダムウォークシミュレーションにおける比較的長時間の拡散挙動を示したものである。
【図７】図５のデータを、（１２）式を用いて最小二乗解析した結果を示したものである。
【符号の説明】
１ 多孔質地層
２ センサーユニット
３ 静磁場磁石
４ 傾斜磁場コイル
５ ラジオ波コイル
【発明の属する技術分野】
本発明は、土木工学、地盤工学、資源探査および物理探査の分野において利用可能な核磁気共鳴法を用いた地層の浸透率の推定方法に関する。
【０００２】
【従来の技術】
従来、地層の浸透率の核磁気共鳴（以下「ＮＭＲ」という。）による推定は、地層中の間隙流体のプロトン緩和データを用いて行っている（Stambaugh,2000,Oil&GasJournal,vol.98,pp.45-52）。
【０００３】
【発明が解決しようとする課題】
しかしながら、従来の方法では、岩石の化学組成がわからないと、正確な浸透率推定はできないという問題があった。これはプロトンの緩和時間が、岩石の化学組成に強く影響されるからである。そのため、化学組成がわからない地層の浸透率を推定するときは、大きな誤差が生じるものであった。
また、Balzarini et al,1996.PhysicaB,vol.226,pp.10-14で,拡散データを用いた浸透率の推定方法が述べられているが、その推定方法は誤ったものであった。
【０００４】
本発明は、核磁気共鳴法を用いて多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測することにより、地層の化学組成がわからくても正確な浸透率を推定可能とすることを目的とする。
【０００５】
【課題を解決するための手段】
上記目的を達成するため、本発明の核磁気共鳴法を用いた地層の浸透率の推定方法は、静磁場磁石、傾斜磁場コイルおよびラジオ波コイルから構成されるセンサーユニットを地層に押し当てて地層中の間隙流体分子の自己拡散係数を計測することにより、地層の浸透率を推定する方法において、浸透率ｋを次の式から求めることを特徴とする核磁気共鳴法を用いた地層の浸透率の推定方法。
【数５】
【数２】
【０００６】
【発明の実施の形態】
以下、本発明による実施の形態を図面に基づき説明する。
図１は、本発明の実施の形態の１例を示したもので、ＮＭＲ計測機器のセンサー部分のみ示してある。
１は、間隙水などの含水素流体の自己拡散係数を計測する対象である流体で満たされた多孔質地層、例えば、岩石、堆積物等である。
２は、センサーユニットであり、静磁場磁石３、傾斜磁場コイル４およびラジオ波コイル５から構成されている。
静磁場磁石３は、原子核が共鳴するように一定の強度をもった磁場を形成するものである。 傾斜磁場コイル４は、流体分子の自己拡散による空間移動を検出するために場所により強さの異なる磁場、すなわち、傾斜磁場を作る。ラジオ波コイル５は、高周波磁場を多孔質地層１に照射し共鳴を起こさせると共に、共鳴の結果、多孔質地層１から発せられる信号を受信するものである。また、図示しないコンピュータにより、傾斜磁場の切り換えや高周波磁場の発生のタイミングを制御したり、受信した信号をもとに計算を行うようになっている。
【０００７】
図１のセンサーユニット２を岩盤の表面やボーリング孔の抗壁の地層１などに押し当てて、地層１の表面のＮＭＲ物性を計測する。その際、岩盤表面から一定の距離、例えば、数cm〜数ｍ奥の位置にあるプロトンに感度を持つよう磁場を調整する。計測対象は、多孔質な地層１中の間隙を満たしている軽水素（１Ｈ）の原子核（プロトン）を含む流体である。たとえば、地下水の水分子中のプロトン、あるいは石油の炭化水素分子中のプロトンが対象核である。
【０００８】
間隙流体は、多孔質岩石中を熱運動によって自己拡散して移動している。その自己拡散係数データをパルス傾斜磁場を用いたＮＭＲ（ＰＦＧ ＮＭＲ）で計測するのが、まず浸透率計測の第一歩である。この自己拡散係数を計測するためには、たとえば、図２で示すパルス系列（ＰＦＧスピンエコー系列）を使用する必要がある。プロトンの歳差運動の回転速度、位相などをラジオ波と傾斜磁場のパルスでコントロールし、エコー時間ＴＥのタイミングで現れるスピンエコー波形を取得する。傾斜磁場強度（Ｇ）の上昇に起因する、エコー波形の高さ（Ｉ）の減少を計測し、それから間隙流体分子の自己拡散係数（Ｄ）を求める。すなわち、
【数１】
ｂ＝（γＧδ）２（Δ−δ／３） （２）
γは核種に固有の磁気回転比であり、プロトンの場合は2.675×108 rad/Tsである。Δとδは、図２に示すように、それぞれ、傾斜磁場パルスの間隔と持続時間である。ｂは既知量なので、ｂをいろいろ変えてＩ０ で規格化されたスピン・エコーシグナル（Ｉ／Ｉ０）を計測し、最小二乗法で未知量Ｄを決定する。
【０００９】
ＰＦＧ ＮＭＲでは、間隙流体がΔのタイムスケールで拡散できる距離を計測する。Δを変えて拡散係数の変化を調べると、多孔質媒体の比表面積と屈曲度が計測できることを、Mitra,P.P.,P.N.Sen,およびL.M.Schwartzが１９９３年発見した（Physical Review B,47巻,8565-8574頁）。概略は次のとおりである。
多孔質媒体中の間隙流体の拡散係数（Ｄ）という量は、拡散時間・拡散距離に依存する。ある空隙サイズをもつ多孔質媒体内部で間隙流体がランダムウォークしている系を考えると、空隙の中央付近からランダムウォークを開始した流体分子は最初のうちは固体表面にめったに衝突しないのでバルクの流体と同じ自己拡散係数（Ｄ０）をもつ。流体分子の移動距離は、時間が経つにつれて増加する。ランダムウォーカーの根平均二乗変位が空隙サイズとほぼ同じかあるいはより大きくなると固体との衝突が頻繁におこり、流体分子は固体を回り込みながら空隙中を移動していく。この固体部分の幾何学的な障害物効果によって、ＤはＤ０よりすこしだけ低下する。重要なポイントは、その低下量が空隙の比表面積に関係する点である。
【００１０】
【数２】
【００１１】
さらに、拡散時間Δを大きくとると、流体分子は空隙サイズよりも十分大きな距離を移動できるようになり、その結果、自己拡散係数はある値, Ｄ ∞, に漸近する。
Ｄ＝Ｄ∞ ａｓ ｔ → ∞ （４）
Ｄ０／Ｄ∞ という比は,多孔質媒体の空隙構造の屈曲度(tortuosity)と呼ばれる量である（たとえばLatour, L.L., R.L. Kleinberg, P.P. Mitra and C.H. Sotak, Journal of Magnetic Resonance, A112巻, 83-91頁, 1995年）。また、地層の空隙率は、従来のＮＭＲ検層技術等で確立しているように、ＣＰＭＧ時系列データの解析から推定できる（たとえば、Straley, C., Rossini, D., Vinegar, H., Tutunjian, P., and Morriss, C.著The Log Analyst, ３−４月号、84−９４頁1997年）。
このように、ＮＭＲを用いれば、地層の空隙率、比表面積、屈曲度を計測できる。
【００１２】
多孔質地層の浸透率は、空隙率、比表面積、屈曲度がわかればカルマン・コゼニー式を用いて、推定できる。ここが、本発明のオリジナル部分である。以下に、その推定方法を述べる。一般に、多孔質媒質の浸透率（ｋ）は、カルマン・コゼニー式で近似できる場合が多い（たとえばGueguen, Y., とV. Palciauskas, Introduction to the Physics of Rocks, Princeton University Press, Princeton, 1994年）。すなわち、
【数３】
【００１３】
さて、(５)式において、（ｌ／ｌ０） という比はパイプの伸び率を表す尺度である。本来は長さがｌ０というパイプが多孔質岩石の複雑な空隙構造のためにｌに伸びたことを意味する。このパイプを満たしている流体分子は、自己拡散によって、曲がりくねったパイプに沿って常時移動している。したがって、自己拡散係数の大小は、曲がりくねったパイプ構造の複雑性を反映している。ところで、拡散現象はランダムウォーク過程であり、数学的には放物型の偏微分方程式に分類される。放物型においては、ある距離Ｌを拡散移動するのに必要な時間はＬの２乗に比例する。したがって、ｌの距離を拡散するのに必要な時間はｌ０の距離を移動するのに必要な時間の（ｌ／ｌ０）２倍かかる。この比は、拡散係数の比（屈曲度）そのものである。
したがって、(５)式における（ｌ／ｌ０）２ という比は、（Ｄ０／Ｄ∞）という比で置き換えることができる。
【００１４】
【数４】
【数５】
【００１５】
（７）式の右辺は、すべてＮＭＲで計測可能な点が特徴である。空隙率φは、従来のＮＭＲ検層技術等で確立しているように、ＣＰＭＧ時系列データの解析から推定できる（たとえば、StraleyほかThe Log Analyst, ３−４月号、84−９４頁1997年）。空隙の比表面積(S/V)poreは、（３）式で示したように、比較的短いΔの条件のもとでＰＦＧ ＮＭＲ法を用いてＤ／Ｄ０ の変化を計測すれば、得られる。Ｄ0は間隙流体の種類（水か石油かなど）と温度がわかれば既存の文献等から知ることができる。Ｄ∞は、（４）式で示したように、比較的長いΔの条件の下でＰＦＧ ＮＭＲ法を用いてDを計測すれば、得られる。
【００１６】
（７）式は、浸透率（ｋ）を拡散データを用いて推定するというオリジナルな方法である。ＮＭＲで地層の浸透率を推定する場合、これまで２つの方法（Timur-Coates モデルとKenyon モデル）が使われてきた（たとえばStambaugh の解説論文(2000年) Oil&Gas Journal, 98巻, 45-52頁を参照）。その２つのモデルは、地層中の間隙流体のプロトン緩和データを用いてｋの推定を行っている。すなわち、ＣＰＭＧ法を用いてプロトンの横緩和現象を計測し、プロトン横緩和時間のヒストグラムから空隙サイズ分布を推定して浸透率を計算している。本発明のような傾斜磁場コイルを用いたＰＦＧ ＮＭＲを実施しているわけではない。プロトン緩和時間を用いた従来の浸透率推定方法には重大な欠点がある。すなわち、プロトン緩和時間が岩石の化学組成（鉄の含有量など）に依存する点である。鉄イオンなどの常磁性物質はプロトンの緩和を促進するので、たとえ同じ空隙サイズでも鉄含有量の多い多孔質岩石は緩和時間が短くなり、少ない岩石では緩和時間が長くなる。したがって、同じ浸透率・同じ空隙構造の地層でも、鉄分含有量が異なるとプロトン緩和時間も異なる。したがって、従来の方法で浸透率を正しく推定するには、岩石の化学組成の情報が必須である。しかし、たとえばＮＭＲ検層でボーリング孔のスキャンを行う場合、地下深部の岩石の種類に関する情報は、獲得困難な場合がしばしばあり、これが浸透率の信頼性を落としている一因になっている。一方、（７）式では、緩和時間の情報は使っていないので、岩石の化学組成の情報は必要ではない。したがって、岩石の化学組成の情報がなくても、浸透率を精度よく推定できる。
【００１７】
【実施例】
(７)式の有効性を実証するために、以下のランダムウォークのコンピューターシミュレーションを行った。サンプルは、直径２．１１ｍｍのガラスビーズ堆積物である（粒度のそろった砂質堆積層のモデル物質として選んだ）。このサンプルの浸透率を（７）式で推定できることを実証する。ビーズ堆積物の空隙率は34％である。ガラスビーズをプラスチック容器に充填して、マイクロフォーカスＸ線ＣＴ（コンピューター断層）画像を取得した。その３次元デジタル画像データを図３に示す。データサイズは、２５６×２５６×２５６画素で、各画素は一辺が０．０５３ｍｍの立方体である。このデジタル画像中の空隙部分を水が満たしていると仮定し、その水分子の自己拡散を離散格子上のランダムウォークでシミュレーションした。ランダムウォークプログラムは、インターネットで無料公開されているマセマティカプログラムＲＷ３Ｄ．ｍ（Watanabe and Nakashiama、 Computers and Geosciences, 28巻, 583-586頁,2002年）を使用した。一般に、ランダムウォークする物質の拡散係数は（８）式で与えられる。
【数６】
<ｒ2>は、ランダムウォーカーの平均自乗変位であり、今回の離散格子上のコンピューターシミュレーションでは（９）式で計算される。
【数７】
【００１９】
Ｄ0はビーズのない自由空間（バルクの水）中の拡散係数であるが、それは(１０)式で与えられる（たとえば、P. T. Callaghan著Principles of Nuclear Magnetic Resonance Microscopy, オックスフォード大学出版、オックスフォード、１９９１年）。
<ｒ2>＝６Ｄ0ｔ＝ａ２τ （１０）
aは、３次元デジタル単純立方格子の格子定数であり、この場合はＣＴ画像の画素サイズの値（０．０５３ｍｍ）である。（４）（８）（１０）式によれば、ビーズ堆積物の屈曲度（Ｄ０／Ｄ∞）はランダムウォークシミュレーションでは、（１１）式で計算できる。
【数８】
【数９】
【００２０】
ランダムウォークシミュレーションの軌跡の一例を図４に示す。比較的短時間の拡散挙動を図５ (n = 15000)に、比較的長時間の拡散挙動を図６(n = 300)に示す。図５のデータを、（１２）式を用いて最小二乗解析した結果を図７に示す。図７によって(Ｓ／Ｖ)poreは4.5×103 1/mという数値が得られた。一方、直径2.11mm、空隙率34%の均一球の堆積物は(Ｓ／Ｖ)pore = 4.8×103 1/mになるはずである。したがって、今回のランダムウォークシミュレーションでは、誤差６％で比表面積を正しく推定することができた。また、図６を（１１）式を用いて解析すると、Ｄ０／Ｄ∞は平均二乗変位の傾きの比として計算できるので 0.002809/0.001752 = 1.60という値を得た。一方、Latour, L.L., P.P. Mitra, R.L. Kleinberg, C.H. Sotak著の論文によると、均一球の堆積物についてはD０／D∞ = 1.62である（Journal of Magnetic Resonance, A101巻, 342-346頁, 1993年）。したがって、今回のランダムウォークシミュレーションでは、誤差１％で屈曲度（Ｄ０／Ｄ∞）を推定することができた。
【００２１】
以上の結果（φ = 0.34, (Ｓ／Ｖ)pore=4.5×103 1/m, Ｄ０／Ｄ∞ = 1.60）を、(7)式に代入すると、 ｋ = 1.0×10-8 m2という値を得る。一方、このビーズ堆積物の浸透率を通常の室内透水試験方法で実測するとｋ = 0.74×10-8 m2という値を得た。したがって、ランダムウォークシミュレーションによって、誤差35％で浸透率の推定に成功した。このことは、拡散データを用いた浸透率推定の妥当性を例証するものである。さらに、カルマン・コゼニー式は天然の多孔質岩石でも成立する場合が多いので、（７）式の天然の地層への適用の有効性をも示唆するものである。
【００２２】
【発明の効果】
本発明は、以下の効果を奏する。
従来、地層中のプロトンの横緩和現象を計測し、プロトン横緩和時間のヒストグラムから空隙サイズ分布を推定して浸透率を計算しているが、プロトン緩和時間が地層の化学組成（鉄の含有量など）に依存するため、浸透率を正しく推定するには、地層の化学組成の情報が必須である。
これに対し、本発明は、核磁気共鳴法を用いて多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測することにより、地層の化学組成がわからくても正確な浸透率を推定することができる。
【図面の簡単な説明】
【図１】本発明の実施の形態のＮＭＲ計測機器のセンサーユニット部分を説明した図である。
【図２】本発明の実施の形態のＮＭＲ計測機器を使用して多孔質な地層中の間隙水などの含水素流体の自己拡散係数を計測する場合に使用するパルス系列の１例である。
【図３】本発明の有効性を実証するためのガラスビーズ堆積物の３次元デジタル画像データを示したものである。
【図４】図３に示すガラスビーズ堆積物中のランダムウォークシミュレーションの軌跡の一例を示したものである。
【図５】ガラスビーズ堆積物中のランダムウォークシミュレーションにおける比較的短時間の拡散挙動を示したものである。
【図６】ガラスビーズ堆積物中のランダムウォークシミュレーションにおける比較的長時間の拡散挙動を示したものである。
【図７】図５のデータを、（１２）式を用いて最小二乗解析した結果を示したものである。
【符号の説明】
１ 多孔質地層
２ センサーユニット
３ 静磁場磁石
４ 傾斜磁場コイル
５ ラジオ波コイル
Claims (1)
- 静磁場磁石、傾斜磁場コイルおよびラジオ波コイルから構成されるセンサーユニットを地層に押し当てて地層中の間隙流体分子の自己拡散係数を計測することにより、地層の浸透率を推定する方法において、浸透率ｋを次の式から求めることを特徴とする核磁気共鳴法を用いた地層の浸透率の推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002220536A JP3837513B2 (ja) | 2002-07-30 | 2002-07-30 | 核磁気共鳴法を用いた地層の浸透率の推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002220536A JP3837513B2 (ja) | 2002-07-30 | 2002-07-30 | 核磁気共鳴法を用いた地層の浸透率の推定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004061327A JP2004061327A (ja) | 2004-02-26 |
| JP3837513B2 true JP3837513B2 (ja) | 2006-10-25 |
Family
ID=31941095
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002220536A Expired - Lifetime JP3837513B2 (ja) | 2002-07-30 | 2002-07-30 | 核磁気共鳴法を用いた地層の浸透率の推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3837513B2 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1879022A1 (en) * | 2005-04-11 | 2008-01-16 | Keio University | Method for locally measuring mobility of protic solvent in sample, device for locally measuring mobility of protic solvent in sample, measuring device for locally measuring behavior of protic solvent in sample by using magnetic resonance method, measuring method, and program |
| JP4831407B2 (ja) * | 2006-01-12 | 2011-12-07 | 独立行政法人産業技術総合研究所 | くりこみを用いた多孔質物体の浸透率の推定方法 |
| US7417426B2 (en) * | 2006-03-11 | 2008-08-26 | Xigo Nanotools Llc | Compact and portable low-field pulsed NMR dispersion analyzer |
| JP5588713B2 (ja) * | 2010-03-31 | 2014-09-10 | イビデン株式会社 | シミュレーション装置及びそのコンピュータプログラム |
-
2002
- 2002-07-30 JP JP2002220536A patent/JP3837513B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004061327A (ja) | 2004-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mohnke et al. | Pore size distributions and hydraulic conductivities of rocks derived from magnetic resonance sounding relaxation data using multi-exponential decay time inversion | |
| Al-Mahrooqi et al. | An investigation of the effect of wettability on NMR characteristics of sandstone rock and fluid systems | |
| AU711508B2 (en) | Nuclear magnetic resonance borehole logging apparatus and method | |
| Luo et al. | Robust determination of surface relaxivity from nuclear magnetic resonance DT2 measurements | |
| Elsayed et al. | The effect of clay content on the spin–spin NMR relaxation time measured in porous media | |
| BRPI0410179B1 (pt) | método para registro de uma formação de terra e aparelho para utilização em uma formação de terra | |
| Costabel et al. | Noninvasive estimation of water retention parameters by observing the capillary fringe with magnetic resonance sounding | |
| US9541513B2 (en) | Method for nuclear magnetic resonance diffusion measurements | |
| Costabel et al. | Relative hydraulic conductivity in the vadose zone from magnetic resonance sounding—Brooks-Corey parameterization of the capillary fringe | |
| Müller-Petke et al. | Nuclear magnetic resonance average pore-size estimations outside the fast-diffusion regime | |
| Ronczka et al. | Optimization of CPMG sequences to measure NMR transverse relaxation time T 2 in borehole applications | |
| CN108779670B (zh) | 用于使地下地层的流体成像的化学选择性成像器及其使用方法 | |
| Liang et al. | Two-dimensional nuclear magnetic resonance method for wettability determination of tight sand | |
| Afrough et al. | Magnetic-resonance imaging of high-pressure carbon dioxide displacement: fluid/surface interaction and fluid behavior | |
| Romanenko et al. | An assessment of non-wetting phase relative permeability in water-wet sandstones based on quantitative MRI of capillary end effects | |
| Costabel et al. | Relative hydraulic conductivity and effective saturation from Earth’s field nuclear magnetic resonance–a method for assessing the vadose zone | |
| Grunewald et al. | Nonexponential decay of the surface-NMR signal and implications for water content estimation | |
| Zielinski et al. | Restricted diffusion effects in saturation estimates from 2D diffusion-relaxation NMR maps | |
| Afrough et al. | Magnetic-resonance imaging of fines migration in Berea sandstone | |
| Chen et al. | Nuclear magnetic resonance imaging for studies of flow and transport in porous media | |
| Tandon et al. | A new NMR-based method for wettability assessment in mixed-wet rocks | |
| Singh et al. | Comparison of saturated hydraulic conductivity estimated by surface NMR and empirical equations | |
| Coman et al. | Improved NMR Logging Approach to Simultaneously Determine Porosity, T 2 and T 1 | |
| Garcia et al. | Integrated workflow to estimate permeability through quantification of rock fabric using joint interpretation of nuclear magnetic resonance and electric measurements | |
| JP3837513B2 (ja) | 核磁気共鳴法を用いた地層の浸透率の推定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060214 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060403 |
|
| 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: 20060516 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 3837513 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |