(57)【特許請求の範囲】
【請求項１】
様々なサイズの空隙を内包する不均一な多孔質物体の浸透率の推定方法において、推定対象となる前記多孔質物体の空隙率φと空隙サイズｄの分布データを計測し、空隙サイズｄの分布からランダムにデータを多数選び、一辺の長さがＬの立方体の３組の対向する２面間に直径ｄの２本のまっすぐなパイプが当該２面に直交して埋め込まれて成るくりこみのための空隙構造モデルを用いることによって得られたＬに依存しない以下の（４）式から浸透率ｋを計算してオリジナルの浸透率分布Ｋ１を作成し、該オリジナルの浸透率分布Ｋ１を、前記多孔質物体を前記立方体をランダムに３次元に積み上げたものと想定してくりこみ操作を行うことによって次第に小さな標準偏差の分布に収束させ、その平均値を推定対象となる前記多孔質物体全体の浸透率とすることを特徴とするくりこみを用いた多孔質物体の浸透率の推定方法。
【数４】
0004831407.tif 000010
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、様々なサイズの空隙を内包する不均一な多孔質物体、例えば地層のシステム全体としての実効的な浸透率を、くりこみ（ｒｅｎｏｒｍａｌｉｚａｔｉｏｎ）という手法を用いて、空隙サイズ分布データから推定する発明に関するものである。
【背景技術】
【０００２】
地下水や石油などの地下流体は、多孔質地層の空隙に存在し、地層に井戸を掘って吸引すると、流体は圧力勾配に駆動されて空隙内部を流れて井戸に向かって移動する。多孔質地層の浸透率（単位はm2）は、一定の圧力勾配をかけたときの単位面積当たりの流速を決める量であり、井戸の生産速度を規定する経済的に重要な水理学的特性である。多孔質地層の空隙サイズデータは、浸透率を計算するために必要な量であり、核磁気共鳴法（ＮＭＲ）法や、回収したコアの水銀ポロシメーター室内分析法などで計測している。
【０００３】
この場合、問題となるのは天然の岩石の空隙サイズの著しい不均一性である（例えば、図３参照。）。図３は、天然の秩父産多孔質砂岩の２次元Ｘ線ＣＴ画像で、直径４ｍｍの円柱形試料の断面図を示しており、黒い部分が空隙、灰色・白色部分が造岩鉱物である。わずか４ｍｍの中に大小さまざまなサイズの空隙が混在しているのがわかる。
一般に、数ｎｍサイズの空隙から数ｍｍオーダーの空隙が、一見均一に見える数グラムの岩石試料の中に混在している。カルマン・コゼニーモデル（円柱パイプモデル）によれば、浸透率は空隙サイズの２乗できいてくる（たとえば、非特許文献１参照。）。したがって、たとえば10ｎｍの空隙サイズから計算する浸透率と、10ｍｍの空隙サイズから計算するそれとでは、サイズが100万倍異なるので、浸透率では1兆倍ものひらきがある。井戸の生産速度を支配する浸透率を正しく評価するためには、図３のように微視的に桁違いに広い空隙サイズをもつ多孔質岩石から、たとえば数cm〜数100mサイズでの巨視的な（実効的な）地層の浸透率を求める方法が要求される。
【非特許文献１】Y. Nakashima and T. Yamaguchi (2004年) X-ray CT for  Geomaterials; Soils, Concrete, Rocks. (Otani, J. and Obara, Y.編集) A. A. Balkema Publishers, pp. 103-110
【発明の開示】
【発明が解決しようとする課題】
【０００４】
頻繁に使われている従来技術の一つとして、図１０に示すように、核磁気共鳴センサーを井戸に降ろして坑壁を連続スキャンして、水や石油で満たされた多孔質地層のプロトン横緩和時間分布データ（空隙サイズ分布データと等価）を原位置で得て、それから地層の浸透率を推定する物理探査方法がある（たとえばK.-J. Dunn, D.J. Bergman, G.A. Latorraca著、 Nuclear Magnetic Resonance Petrophysical and Logging Applications, Pergamon, Amsterdam, 2002年）。
図１０は、NMR検層の模式図であり、NMRゾンデ中に、プロトンを歳差運動させる永久磁石と、スピンを励起・検出するためのコイルが搭載されている。ゾンデを坑壁に押しつけて、ワイヤーを引き上げながら、地下水や石油で満たされた多孔質地層の坑壁表面から１．９ｃｍ奥にある感度領域（斜線部、概形は直径1.3cm, 長さ12cmの円柱）の間隙流体のプロトン横緩和時間をスキャンするものである。
【０００５】
この、いわゆる核磁気共鳴検層では、下記の（１）式または（２）式で地層の巨視的な浸透率kを推定している。
〔Schlumberger-Doll-Research法〕
【数１】
0004831407.tif 000002
φは、多孔質岩石の空隙率つまり空隙の体積分率（０ ≦φ≦ １）。T2LMは、プロトン横緩和時間(T2)の分布の対数平均値。この対数平均値がｍｓ単位ならば、（１）式のｋはミリダルシー単位である。
〔Timur-Coates法〕
【数２】
0004831407.tif 000003
FFVは自由流体体積分率(Free Fluid Volume)とよび、サイズの大きな空隙のしめる体積分率のことである。BFVは拘束流体体積分率(Bound Fluid Volume)とよび、サイズの小さな空隙のしめる体積分率である。φ= FFV + BFVの関係がある。NMR緩和時間の観点で言えば、FFVはある閾値より大きなT2をもつ空隙の体積分率であり、BFVはある閾値より小さなT2をもつ空隙の体積分率である。本明細書では、(２)式において、FFVとBFVの閾値（横緩和時間のカットオフ値）は、よく使われる値33msにした。なお、（２）式のｋはミリダルシー単位である。
【０００６】
上記の両方式に共通する点は、数桁にわたって広く連続的に分布している空隙サイズの分布から、分布の面積総和（空隙率φ）と平均値（対数平均値あるいは、FFV/BFV比）という２つの統計量だけを使って地層の浸透率を推定している点である。この方式は計算量が少ないので、プロトン横緩和時間分布から浸透率を迅速に知りたい場合は役にたつ。しかし、２つの統計量だけではサイズ分布の形状をユニークに決めることはできない。図１１は、２つの仮想的な空隙サイズ分布を斜線と黒で示したものである。斜線はバイモーダル分布で２つのピーク間隔が狭い。黒もバイモーダルだが、ピーク間隔は広い。斜線と黒の分布は、ともに総面積（φ）と対数平均値（T2LM）は同じであるが、あきらかに分布の形状はことなる。したがって、巨視的な浸透率が同じ値を示す保証はない。
このように、２つの統計量だけを使った（１）（２）式は、原油を貯留する砂岩や石灰岩などの特殊な地層には使えるかもしれないが、いかなる形状のサイズ分布（横緩和時間分布）にも適用できるかといえば、それは間違いである。石油の貯留岩以外の、様々な空隙構造特性を持つ自然界の地層のすべてに適用できる保証はない。実際、仮想的なデータについてではあるが、幾何平均値（対数平均値と同義）を使うと真の浸透率値を正しく計算できない例が既に報告されている(Transport in Porous Media, 4, 37-58, 1989)。
【０００７】
本発明は、統計物理学で使われているくりこみ(renormalization)理論を空隙構造モデルに適用し、図３のような様々なサイズの空隙を内包する不均一な多孔質物体、たとえば地層の、システム全体としての実効的な浸透率（「透水係数」ともいう。）を推定するものであって、空隙サイズ分布のすべてのデータポイントを使うことにより、面積総和と平均値の２つしか使わない従来の方法にくらべて精度良く浸透率を推定可能とすることを目的とする。
【課題を解決するための手段】
【０００８】
「くりこみ」とは，もともとは多体系の統計力学において粗視化の概念を定式化して巨視的な性質を導くための方法で、とくに臨界現象に適用されて成功をおさめている。くりこみを地下水理学における巨視的な浸透率の計算に初めて応用したのは、P.R. Kingである(Transport in Porous Media, 4, 37-58, 1989年)。
【０００９】
図１は、「くりこみ」の基本概念を示したもので、本発明は３次元くりこみであるが、わかりやすいように２次元の場合を示している。
今、４ｘ４のセルからなる不均一な多孔質物体である多孔質の岩石があるとする。これが、オリジナルの岩石組織である。１６個のセルの空隙サイズは異なっている。その結果、浸透率もある分布Ｋ１に従い、セルごとに異なる値(K1A, K1B,,,,,,K1P)をとる。くりこみの第１ステップは、４個のセルをサブシステムとして仮想的に考え、その４個の系の全体的な浸透率（いわば電気回路網の合成抵抗）をキルヒホッフとオームの法則を用いて厳密に解くことである。その操作を４回くりかえすと、サブシステム４個の浸透率組み合わせを得る。その結果、浸透率分布として、Ｋ２を得る（成分はK２A, K２B, K２C, K２D）。
くりこみの第２ステップでは、さらに上位の２ｘ２スーパーサブシステムを想定して、４個のサブシステム（成分はK２A, K２B, K２C, K２D）を統合した全体的な浸透率Ｋ３を、キルヒホッフとオームの法則を用いて厳密にとく。このようなセルの粗視化による、システム全体の浸透率の推定がくりこみ操作である。実際の岩石は４ｘ４ではなく無数のセルから構成されているので、くりこみ操作も多数回行うのが通常であり、得られた浸透率分布の列Ｋ２，Ｋ３、Ｋ４、、、が十分小さい標準偏差の分布に収束していくことを確認して、収束したその平均値をシステム全体の実効的な浸透率とみなすものである。
【００１０】
図２に基づいて、くりこみによる具体的な浸透率の計算方法を説明する。
図１の４ｘ４の２次元のオリジナルシステムを例にとると、２ｘ２=４個のセルごとにくりこみ操作を行う。たとえば、K１A, K１B, K１E, K１Fの浸透率をもつ４個のセルについていえば、各セルを抵抗値1/K１A, 1/K１B, 1/K１E, 1/K１Fの電気抵抗値をもつ素子と見なす。なお、抵抗値は流れにくさの指標であり、浸透率は流れやすさの指標なのでお互いに逆数の関係にある点に注意する。さらに、セル間の接続を可能にするため、たとえばK１Aのセルは、抵抗値1/2K１Aの素子４個に細分し、図２（ｂ）のように組み上げて、４種類の電気抵抗値をもつ合計１６個の抵抗素子からなる電気回路ネットワークを形成する。図２（ｃ）のように、この回路に仮想的に電流を流すことを考え、回路全体の合成抵抗値, 1/K2A, をキルヒホッフとオームの法則を用いて厳密に解く。その結果、
【数３】
0004831407.tif 000004
を得る。K2Aは、K１A, K１B, K１E, K１Fの浸透率値を持つ2x2個のセルをくりこみ操作で１つのセルにまとめた（粗視化した）ときの浸透率値である。
【００１１】
上記した精度良く浸透率を推定可能とすることを達成するため、本発明のくりこみを用いた様々なサイズの空隙を内包する不均一な多孔質物体の浸透率の推定方法において、推定対象となる前記多孔質物体の空隙率φと空隙サイズｄの分布データを計測し、空隙サイズｄの分布からランダムにデータを多数選び、以下の（４）式から浸透率ｋを計算してオリジナルの浸透率分布Ｋ１を作成し、該オリジナルの浸透率分布Ｋ１をくりこみ操作によって次第に小さな標準偏差の分布に収束させ、その平均値を推定対象となる前記多孔質物体全体の浸透率とすることを特徴としている。
【数４】
0004831407.tif 000005
【発明の効果】
【００１２】
本発明は、上記の手段を採用することにより、以下のような優れた効果を奏する。
（１）数桁にわたって広く分布している空隙サイズの分布から、たった２つの統計量（面積と平均値）しか抽出していない従来の核磁気共鳴法では石油の貯留岩（砂岩及び石灰岩など）の特殊な地層にしか使えないのに比較して、本発明は、様々な空隙構造をもつ自然界の多孔質物体のすべてに適用可能である。
（２）従来の核磁気共鳴法に比べて、格段に真の浸透率の値（実測値）に近い値を得ることができる。
【発明を実施するための最良の形態】
【００１３】
本発明に係るくりこみを用いた多孔質物体の浸透率の推定方法を実施するための最良の形態を実施例に基づいて図面を参照して以下に説明する。  
【実施例１】
【００１４】
くりこみ理論を使うために、まず、図４に示すようにくりこみのための空隙構造モデル（パイプモデル）を作った。多孔質セル１は、一辺の長さがＬの立方体とし、これが、図１の４ｘ４個のオリジナルセルの一個一個に相当する。各セル１には、この図のように３方向に直交して走るまっすぐなパイプ２の群が埋め込まれている。各方向に走っているパイプ２の数は同じとする。パイプ２の中を流体が低いレイノルズ数で流れる。パイプ２の直径ｄは、セル１ごとには異なるが一つのセル内では共通な値をとる。セル１の空隙率φはｄによらずすべてのセルで共通である。すなわち、セルとして十分大きなＬ値のモデルを想定した場合、すべてのセルで空隙率が共通な値φ（中心極限定理で収束した値）をとることができる。
【００１５】
0004831407.tif 000006
【数４】
0004831407.tif 000007
(４)式はＬには依存しないので、仮想的なセルサイズＬの恣意性はない。空隙サイズデータの分布の広さはｄの値に取り込まれる。すなわち、各セルのｄ値は実測した空隙サイズデータからランダムにサンプリングした値をとり、そのセルを図５に示すような３次元にランダムに積み上げた大きなシステムを仮想的に考える。図５は、３次元ランダム多孔質物体３の概念図（このケースでは、５０ｘ５０ｘ５０のセルから構成）である。各セル１は図４の多孔質立方体からなり、全体として空間相関のない等方でランダムな浸透率分布になっている。  画素数が無限大のシステムについて、そのシステムに流体を流したときの流れやすさ、すなわち、システム全体としての実効的な浸透率を図１のくりこみ法で推定するわけである。
【００１６】
ちなみに、K. Xu, J-F Daian, and D. Quenard の提出した空隙構造モデル（Transport in Porous Media, 26, 51-73, 1997年）は、数学的な処理の効率性を追求するあまり、現実的でない空隙構造モデルを採用したので不自然である。たとえば、彼らはセルサイズLの分布がべき乗であることを仮定し、さらにLとｄの比率が一定であるという前提条件を設けている。しかし、これは、現実の多孔質岩石の３次元空隙構造に関する最近の研究(Lindquist他、Journal of Geophysical Research, 105, 21509-21527, 2000年)から明らかなように、現実とは矛盾する内容である。一方当該特許は、このような非現実的な仮定は設けていないので、より現実に則したモデルといえる。
【００１７】
くりこみ法による浸透率の推定の手順は以下のとおりである。
【００１８】
手順（１）：
空隙率と空隙サイズ分布データを計測する。たとえば回収したコアを水銀ポロシメーターで分析する。あるいは、NMR検層を行って横緩和時間の分布を得る。ただしNMR検層の生データは横緩和時間分布であってサイズ分布ではないので、たとえばコア試料を使って岩石のプロトン表面緩和率 (surface relaxivity）,ρ2, を計測して横緩和時間, T2, をサイズ分布に換算する必要がある（たとえばK.-J. Dunn, D.J. Bergman, G.A. Latorraca著、 Nuclear Magnetic Resonance Petrophysical and Logging Applications, Pergamon, Amsterdam, 2002年）。ここでは、図４に示す円柱パイプを仮定しているので、d = 4T2ρ2という式で換算した。
【００１９】
手順（２）：
空隙サイズｄの分布からランダムにデータを多数選び、（４）式で与えられる浸透率ｋを計算し、浸透率の分布をつくる。これが、図１のオリジナル浸透率データセットK1に相当する。
【００２０】
手順（３）：
図１の概念で、３次元セルのくりこみシミュレーションを行う。２ｘ２ｘ２個からなる３次元立方体サブシステムについて、この８個からなる系の全体的な浸透率（いわば電気回路網の合成抵抗）をキルヒホッフとオームの法則を用いて厳密に解く。詳細な計算方法は、Transport in Porous Media, 4, 37-58, 1989年)に従う。広い標準偏差を持つオリジナルの浸透率分布K1がくりこみ操作によって次第に小さな標準偏差の分布に収束していく様子を計算機上でモニタリングしながら、十分な精度まで収束したと判断したらそこでくりこみ操作を終了して、その平均値をシステム全体の実効的な浸透率とみなす。
【００２１】
〔実験例〕
２つの多孔質砂岩と１つの多孔質溶岩について、水銀ポロシメーターによる空隙サイズ分布データから本発明のくりこみ手法で浸透率を推定した。また比較のため、従来の核磁気共鳴法、つまり（１）式および（２）式、による浸透率推定も行った。結果は表１に示す。
【表１】
0004831407.tif 000008
［実験に用いた岩石試料の説明］
（１）埼玉県秩父産の砂岩（図３）は、空隙率φ=0.14であり、プロトン表面緩和率（ρ2= 5.6x10-5 m/s）や浸透率の実測値などのデータは「Y. Nakashima, T. Nakano, K. Nakamura, K. Uesugi, A. Tsuchiyama, and S. Ikeda 著(2004年) Journal of Contaminant Hydrology, 74, 253-264」で開示されている。水銀ポロシメーター実験による空隙サイズｄの分布データを図６にしめす。縦軸は岩石試料１ｇあたりの空隙体積である。図６のｄをd = 4T2ρ2という式で換算した横緩和時間分布について解析した結果、FFV = 0.0155, BFV = 0.1245, T2LM = 2.61 msであった。
（２）秋田県澄川産の安山岩質溶岩は、澄川地熱発電所の構内でボーリングして地下57.3mから採取したコア由来のものである。空隙率φ=0.0286であり、コアの透水試験の実測値は3x10-19 m2である。水銀ポロシメーター実験による空隙サイズｄの分布データを図７にしめす。縦軸は岩石試料１ｇあたりの空隙体積である。水銀ポロシメーターで計測した空隙サイズ分布をNMR横緩和時間に換算する時に必要なプロトン表面緩和率ρ2は9.5x10-4 m/sである。その結果、図７のｄをd = 4T2ρ2という式で換算した横緩和時間分布について解析した結果、FFV = 0, BFV = 0.0286, T2LM = 0.0143 msであった。
0004831407.tif 000009
【００２２】
秩父産砂岩に関するシミュレーション結果を図９の（ａ）〜（ｅ）にしめす。ヒストグラムの縦軸は浸透率データの数。横軸は、浸透率の値（単位は10-15 m2）。Ｋ２，Ｋ３、Ｋ４，Ｋ８の各ヒストグラムのデータ点の総数は７０００個。ヒストグラムの算術平均値と標準偏差も書き込んである（単位は10-15 m2）。非常に広いオリジナル浸透率分布（Ｋ１）が、くりこみ操作を繰り返すごとに標準偏差が小さくなっていく様子がわかる。この場合では、７回くりこみを行うと十分小さな標準偏差になることがわかり、その分布Ｋ８の平均値1.32x10-15 m2 ≒1x10-15 m2をもって浸透率の推定値とした。
【００２３】
［実験結果の考察］
上記３つの岩石試料について、本発明のくりこみ手法と従来のNMR法（(2)式のSDR法と(3)式のTimur-Coates法）による推定結果の比較を表１に総括した。表１から明らかなように、本発明の推定値は、埼玉県秩父産の砂岩と秋田県澄川産の安山岩質溶岩について、従来の核磁気共鳴法による推定値よりも格段に真の値（実測した浸透率値）に近くでた。NMR法をくりこみ法と比べてみると、ベレア産の砂岩についてはSDR法はくりこみ法とかろうじて同じ推定値を出したものの、それ以外のケースでは、桁違いに推定精度が悪いといえる。秩父産砂岩の推定精度は、真の値（実測値）より１桁食い違っている。しかし、天然の地層は、本来約10-21 m2から約10-10m2まで10桁以上の広い不確定性をもっているものである。したがって、くりこみ法による秩父産砂岩の推定が真の値からのズレをわずか１桁に押さえることができたのは、許容できる誤差といえる。このように、３つの試料をもちいた実例によって、当該特許の有効性が確認された。
【図面の簡単な説明】
【００２４】
【図１】本発明において採用している「くりこみ」を説明するための概念図である。
【図２】「くりこみ」による具体的な浸透率の計算方法を示す図である。
【図３】天然の秩父産多孔質砂岩の２次元X線CT画像である。
【図４】繰りこみのためのパイプモデルを示した斜視図である。
【図５】図５の多孔質立方体からなる各セルの５０ｘ５０ｘ５０から構成した３次元ランダム多孔質物体の概念図である。
【図６】秩父産砂岩の水銀ポロシメーター実験による空隙サイズｄの分布データを示す図である。
【図７】澄川産溶岩の水銀ポロシメーター実験による空隙サイズｄの分布データを示す図である。
【図８】(a)は、ベレア産砂岩の水銀ポロシメーター実験による空隙サイズｄの分布データ、(b)は、ベレア産砂岩のＮＭＲ実験(2.2MHzのＮＭＲ検層ゾンデを使用)によるプロトン横緩和時間の分布データを示す図である。
【図９】秩父産砂岩のくりこみシミュレーション結果を示す図である。
【図１０】NMR検層の模式図である。
【図１１】２つの仮想的な空隙サイズ分布を示す図である。
【符号の説明】
【００２５】
１    セル
２    パイプ
３    ３次元ランダム多孔質物体
【図１】
0004831407.tif 000011
【図２】
0004831407.tif 000012
【図４】
0004831407.tif 000013
【図６】
0004831407.tif 000014
【図７】
0004831407.tif 000015
【図８】
0004831407.tif 000016
【図９】
0004831407.tif 000017
【図１０】
0004831407.tif 000018
【図１１】
0004831407.tif 000019
【図３】
0004831407.tif 000020
【図５】
0004831407.tif 000021