(57)【要約】地下の貯留層の切削片などの岩石の屑または破片についての測定値から、貯留層の岩石の全飽和範囲にわたる毛管圧曲線を、短時間で低コストで評価する方法。  本方法は、岩石の屑の浸透率ｋを測定するステップと、初期状態で液体（たとえば食塩水）により飽和しているこれらの破片を遠心分離することにより、飽和率の関数として毛管圧曲線Ｐｃを測定するステップと、毛管圧曲線の全体を求めるために、遠心分離により測定された毛管圧曲線の漸近部分に対して、及び、切削片について測定された浸透率ｋの値に対して適合するように制約されていて、調整可能なパラメータに依存している経験的関係を満足する毛管圧曲線Ｐｃをパラメータ化するステップと、を有している。  応用分野：たとえば、炭化水素用貯留層の評価
2005525574.tif 000002
【特許請求の範囲】
【請求項１】
地下の貯留層から採取された岩石の切削片についての測定値から前記岩石の毛管圧曲線を求める方法であって、前記切削片の浸透率ｋを測定するステップを有する方法において、
初期状態で液体により飽和している前記切削片を遠心分離することにより、飽和率の関数として毛管圧曲線Ｐｃを測定するステップと、
前記毛管圧曲線の全体を求めるために、遠心分離により測定された前記毛管圧曲線の漸近部分に対して、及び、前記切削片について測定された前記浸透率ｋの値に対して適合するように制約されていて、調整可能なパラメータに依存している経験的関係を満足する、毛管圧曲線Ｐｃをパラメータ化するステップと、を有することを特徴とする方法。
【請求項２】
前記パラメータ化するステップは、毛管圧Ｐｃの低い飽和率に対応している前記漸近部分の較正を可能にする一式の前記パラメータを既定値で選択するステップと、使用された前記経験的関係の１つにより得られた前記浸透率の予測が、前記切削片に対して実施された前記浸透率の測定値と前記漸近部分に最も良好に適合するように前記パラメータを段階的に修正するステップとを有することを特徴とする、請求項１に記載の方法。
【請求項３】
前記切削片の浸透率ｋは、前記切削片を含んだ液体が満たされている容器を、加圧されている同じ液体を収容しているタンクと所定の期間連通させた後の、前記容器内の圧力の変化の測定値と、前記切削片によって実際に吸収された体積と、前記容器内のモデル化された圧力の変化が前記切削片の物理パラメータの測定された変化に対して最も良好に適合するように反復して調整される、前記切削片の物理パラメータに対して選択された初期値から、前記容器内の圧力または体積の変化をモデル化することと、から測定されることを特徴とする、請求項１または請求項２に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地下の貯留層から採取された岩石の屑についての測定値から、その貯留層の岩石の毛管圧曲線を評価する方法に関する。
【背景技術】
【０００２】
コアや切削片の実験室での測定
地盤構造を貫通した坑井の掘削中に得られる岩石の破片の浸透率、間隙率、毛管特性などの貯留層物性パラメータの測定は、坑井の運営会社にとっては、坑井が横切る領域を形成する貯留層物性の特徴を最初に素早く得ることができる興味深い機会である。
【０００３】
本出願人により出願されたフランス特許第２，８０９，８２１号は、掘削による泥に含まれて地表に戻された切削片から、地下の油層領域の多孔質の岩石の絶対浸透率などの物理パラメータを評価するシステムを記述している。切削片が粘性の液体に浸されている封入容器内において、岩石の孔に捕らえられている気体が圧縮されるように、所定の圧力のしきい値に達するまで時間と共に増加する圧力を印加して、この液体がいくらか注入される。この注入段階の後に、注入が停止される緩和段階が続く。注入処理中の圧力の変化は、切削片の物理パラメータに対して選択されている初期値からモデル化されていて、モデル化された圧力曲線と実際に測定された圧力曲線とが可能な限り最も良好に一致するように、コンピュータはそれらのパラメータを反復して調整する。
【０００４】
出願人により出願されているフランス特許出願公開明細書第０２／００２３号は、同様に切削片から、地下の油層領域の岩石の絶対浸透率や間隙率などの物理パラメータを評価する別の方法について記述している。岩石の破片を収容していて、粘性の液体が満たされている封入容器を、岩石の孔に捕らえられている気体が圧縮されるように、所定の圧力を有していて、同じ液体を収容している容器と連通させる。この圧力を作用させる時間が短いか長いかにより、封入容器内の圧力の変動または岩石の破片により実際に吸収される体積の変動を測定することができる。それから、封入容器内の圧力または体積の変化が、破片の物理パラメータに対して選択された初期値からモデル化され、モデル化された変化が、封入容器内の物理パラメータの測定された変化に対して最も良好に適合するように、岩石の破片の物理パラメータの値が反復して調整される。
【０００５】
貯留層物性の特徴付けにおいては、毛管圧も、運営会社にとっては、非常に重要なデータであり、それは、毛管圧により、以下の項目が定まるためである。
−帯水層領域（当業者は、水と油の接触を表す、ＷＯＣと呼ぶ）から貯留層の上部（遷移領域）までの貯留層内の液体の初期の分布。貯留層の岩石と貯留層内の液体の性質に関連している毛管圧曲線により、この遷移領域は、数メートルあるいは数十メートルにわたって延びることがあり、これは、貯留層内の蓄積状態を求める場合に重要な影響を持つ。
−岩石の入力圧力、これは、帽岩について特に重要である。たとえば、ガスの貯蔵タンクに対して、帽岩の入力圧力により、漏洩のない貯留レベル内の許容可能な超過圧力が定まる。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
現状の技術によれば、毛管圧曲線は、貯留層のコアを実験室で測定することにより求めることができる。これらの方法は、コアについての測定以外に、コアを取得する作業が必要になるために、高価になり、また、掘削後数か月しか利用できないことが多い。
【０００７】
毛管圧曲線を迅速に求めるためのアプローチ
しかし、掘削中、または、そのわずか後に、毛管圧曲線を素早く評価する代替の方法が文献に記述されている。
【０００８】
最も一般的に使用されるアプローチは、切削片から空気／水銀の毛管圧曲線Ｐｃを直接測定するための水銀圧入技術を使用する方法である。しかし、求められた曲線は、湿潤用液体により高飽和状態にあるコアについて得られた基準曲線とは、著しく異なる。また、このアプローチは水銀の使用に基づいていて、水銀の使用は非常に汚染性が高く、多くの国で法律により次第に禁止されつつあり、そのため、この技術を適用することは近い将来大きな問題となる。
【０００９】
その他の既知の方法では、掘削直後に坑井で測定された層序データから毛管圧曲線を迅速に予測するために、核磁気共鳴（ＮＭＲ）技術を使用する。この技術は、以下の刊行物に特に記述されている。
−Ｂｏｗｅｒｓ，Ｍ．，Ａ．ｅｔ  ａｌ：＜＜Ｐｒｅｄｉｃｔｉｏｎ  ｏｆ  ｐｅｒｍｅａｂｉｌｉｔｙ  ｆｒｏｍ  ｃａｐｉｌｌａｒｙ  ｐｒｅｓｓｕｒｅ  ｃｕｒｖｅｓ  ｄｅｒｉｖｅｄ  ｗｉｔｈ  ＮＭＲ＞＞，１７  Ｓｅｐｔｅｍｂｅｒ  １９９８（Ｂｏｗｅｒｓ，Ｍ．，Ａ他：＜＜ＮＭＲにより得られた毛管圧曲線からの浸透率の予測＞＞、１９９８年９月１７日）
−Ｍａｒｓｈａｌｌ，Ｄ．，ｅｔ  ａｌ：＜＜Ｍｅｔｈｏｄ  ｆｏｒ  ｃｏｒｒｅｌａｔｉｎｇ  ＮＭＲ  ｒｅｌａｘｏｍｅｔｒｙ  ａｎｄ  ｍｅｒｃｕｒｙ  ｉｎｊｅｃｔｉｏｎ  ｄａｔａ＞＞、ＳＣＡ  Ｎｏ．  Ｓｏｃｉｅｔｙ  ｏｆ  Ｃｏｒｅ  Ａｎａｌｙｓtｓ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｓｙｍｐｏｓｉｕｍ  １９９５（Ｍａｒｓｈａｌｌ，Ｄ．他：＜＜ＮＭＲ緩和測定と水銀圧入データとを相関させる方法＞＞、ＳＣＡ番号：コア分析者学会国際シンポジウム１９９５）
−Ｖｏｌｏｋｉｔｉｎ，Ｙ．，Ｗ．Ｊ．  ｅｔ  ａｌ：＜＜Ａ  ｐｒａｃｔｉｃａｌ  ａｐｐｒｏａｃｈ  ｔｏ  ｏｂｔａｉｎ  １ｓｔ  ｄｒａｉｎａｇｅ  ｃａｐｉｌｌａｒｙ  ｐｒｅｓｓｕｒｅ  ｃｕｒｖｅｓ  ｆｒｏｍ  ＮＭＲ  ｃｏｒｅ  ａｎｄ  ｌｏｇ  ｄａｔａ＞＞、ＳＣＡ  Ｎｏ．  Ｓｏｃｉｅｔｙ  ｏｆ  Ｃｏｒｅ  Ａｎａｌｙｓｔｓ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｓｙｍｐｏｓｉｕｍ  １９９９（Ｖｏｌｏｋｉｔｉｎ，Ｙ．，Ｗ．Ｊ．他：＜＜ＮＭＲによるコアデータと記録データから最初の排液毛管圧曲線を求める実際的なアプローチ＞＞、ＳＣＡ番号：コア分析者学会国際シンポジウム１９９９）
ＮＭＲ緩和信号が、孔の大きさの分布に関してまず変換され、それから、しきい値の大きさの分布に関して変換され、これにより疑似毛管圧曲線を計算することができる。このアプローチは、既知の曲線Ｐｃのいくつかのサンプルに対して検討されている。その結果から、試験する岩石の性質により個別に実施される厳格な較正段階を実施した場合にだけ、基準曲線と良好に一致することがわかった。この較正段階が必要になるのは、以下の項目が不確定であるためである。
−岩石により変動する表面の緩和度の値に依存しているＮＭＲ信号−孔の大きさの分布の変換
−岩石の性質と続成作用に依存している孔の大きさの分布−しきい値の大きさの分布の変換
そのため、このアプローチは、予測的な調査には適していない。いずれにしても、このアプローチは切削片には適用できない。
【００１０】
曲線Ｐｃを得るための処理としては、画像分析もある。多孔質の媒体が、走査電子顕微鏡、つまりＳＥＭにより撮影された薄い断面の形態で最初に用意される。それから、取得された画像は、岩石に関する空隙の割合と形状を表すパラメータを得るために分析される。特に、疑似毛管圧曲線Ｐｃを再現するために、しきい値の大きさの分布を得ることができる。この方法の主な制限事項は、定義上、毛管圧は３次元（３Ｄ）の性質を有しているのに対して、薄い断面が２次元（２Ｄ）の性質を有していることである。さらに、この技術は、かなり大幅な条件付けが必要であり、これは、わずか後に得られる結果とは実際に見合ったものではない。画像分析は、切削片に適用可能であるが、良好な予測可能性を得るには、慎重な較正が必要になる。
【００１１】
最後に、現場で遠心分離技術が切削に適用されることがあるが、これは、岩石から切削液をできるだけ大量に抽出して、汚染物質の環境への排出を最低限にし、回収された切削液を再利用することでコストを抑えるためである。我々の知る限りでは、毛管特性を得るために、切削片に対する遠心分離が検討されたことはない。
【課題を解決するための手段】
【００１２】
本発明による方法では、地下の貯留層から採取された岩石の屑または破片（切削片など）についての測定値から、これらの岩石の全飽和範囲にわたって、貯蔵層の岩石の毛管圧曲線を、短時間に低コストで求めることができる。本方法は、
−岩石の屑の浸透率ｋを測定するステップと、
−初期状態で液体により飽和している切削片を遠心分離することにより、岩石の屑の飽和率の関数として毛管圧曲線Ｐｃを測定するステップと、
−毛管圧曲線の全体を求めるために、遠心分離により測定された毛管圧曲線の漸近部分に対して、そして、切削片について測定された浸透率ｋの値に対して適合するように制約されていて、調整可能なパラメータに依存している経験的関係を満足する、毛管圧曲線Ｐｃをパラメータ化するステップと、を有している。
【００１３】
曲線のパラメータ化は、毛管圧Ｐｃの低い飽和率に対応している漸近部分の較正を可能にする一式のパラメータを既定値で選択し、使用された経験的関係の１つにより得られた浸透率の予測が、岩石の屑に対して実施された浸透率の測定値と漸近部分に最も良好に適合するようにパラメータを段階的に修正することにより、有利に実施することができる。
【００１４】
切削片の浸透率ｋは、たとえば、切削片を含んだ液体が満たされている容器を、加圧されている同じ液体を収容しているタンクと所定の期間連通させた後の、容器内の圧力の変化の測定値と、切削片によって実際に吸収された体積と、容器内のモデル化された圧力の変化が切削片の物理パラメータの測定された変化に対して最も良好に適合するように反復して調整される、切削片の物理パラメータに対して選択された初期値から、容器内の圧力または体積の変化をモデル化することと、から測定される。
【００１５】
本方法は、入手がより容易で、より安価に取得できる、単なる切削片に基づいて岩石の毛管圧曲線が得られる点で、特に有利である。結果も、コアを使用した場合に比べて、はるかに迅速に得ることができる。
【００１６】
本発明による方法のその他の特徴と利点は、添付の図面を参照して、非限定的な例による実施形態についての以降の説明を読むことで、明らかになろう。
【発明を実施するための最良の形態】
【００１７】
本発明による岩石の破片や切削片からの毛管圧曲線Ｐｃの迅速な評価を目的とした方法を図７に示す。本方法は、全飽和領域にわたってＰｃを再現するための、２つの実験的な測定段階と、それに続く既知のカーブを参照することによるパラメータ化段階に基づいている。両測定段階の一方の間に、初期状態で水で飽和している岩石の破片を遠心分離することにより毛管圧データＰｃが取得される。他方の実験段階により、前述のフランス特許出願公開明細書第０２／００２３号に記述されている方法から、岩石の浸透率の値ｋを計算することができる。パラメータ化された形式を使用して、全飽和領域にわたって曲線Ｐｃが再現される。曲線のパラメータは、曲線で、実験的に得られた毛管圧データＰｃと、以降で再度記述しているＴｈｏｍｅｅｒの方法やＳｗａｎｓｏｎ−Ｋａｍａｔｈの方法などの既知の方法により得られた曲線Ｐｃから予測された浸透率の値を、岩石の破片について測定された値に一致させることにより求められる。一切の事前の較正手続きなしに、基準曲線と非常によく一致している、いくつかの適用事例を示す。
【００１８】
Ｉ）切削片の遠心分離によるＰｃの測定
本方法を実装するために、標準的な遠心分離手段や、たとえば、出願人により出願されている欧州特許第６０３，０４０号（米国特許第５，４６３，８９４号）、フランス特許第２，７６３，６９０号、フランス特許第２，７７２，４７７号（米国特許第６，１８５，９８５号）、またはフランス特許第２，７９８，７３４号に記述されているような、生成される液体の体積を自動的に監視する、より洗練された手段を使用することができる。
【００１９】
掘削作業中に地表に取り出された切削片は、まずソックスレー型の装置内で溶剤により洗浄され、それから乾燥され、３０ｇ／ｌの塩水で飽和される。それから、切削片は、様々な切削片ｄの間に捕らえられていた水を除去するために、湿った布で水を吸い取られ、それから、回転アームの端部に固定されているセルまたはカップ内に配置される。遠心分離により切削片から放出された水は、格子を通過して流れ、カップの基部で回収される。コアの遠心分離の状況と同じようにして、実験データが得られる。（所定の回転速度での）遠心分離の段階で、顕著な変動が観測されなくなるまで、水の生成の変化が測定され、それから、新しい段階を開始するために、回転速度が上昇する。
【００２０】
図１からわかるように、回転速度の上昇に伴う多孔質の媒体の脱飽和を示している、ＮＭＲ信号の逐次減少と、短い緩和時間Ｔ２への移行が観測される。この種類の測定により、切削片間に毛管状の接触が明らかに存在することが示され、それにより、遠心分離により切削片からＰｃが測定可能になる。
【００２１】
実験中に発生する水の体積は、切削片に初期状態で含まれていた水の体積から飽和データに変換される。後者は、重量の測定（飽和の前後の切削片の重量の差）または直接ＮＭＲ測定により求められる。
【００２２】
図２Ａないし５Ａは、コアの遠心分離による従来の方法で測定された曲線Ｐｃが得られている、既知の性質を有している岩石から実験室で作られた、大きさが１から２ｍｍの、モデル切削片に対して実施された実験の結果を示している。漸近部分のレベル（湿潤用液体による低飽和状態）で、基準曲線とよく一致していることがわかる。それに対して、湿潤用液体による飽和が高い状態については、大きな差が観測されている。したがって、切削片に対して水銀を使用した圧入測定で得られた結果と同等の結果が、汚染の危険なしに得られる。
【００２３】
測定された毛管圧曲線Ｐｃは、漸近部分でだけ正しい値を表している。そのため、全飽和領域にわたって曲線のふるまいを評価するには、再現手続きが必要である。
【００２４】
II)切削片の浸透率ｋの測定
前述の特許出願公開明細書第０２／０２，２４２号に記述されている方法が、切削片の浸透率の測定に適用される。そのため、切削片は、粘性の液体を収容している封入容器内でその液体に浸される。それから、岩石の孔に捕らえられている気体を圧縮するために、封入容器は、圧力が加えられている同じ液体を収容している容器と連通させる。第１の実施形態によると、この連通期間を、非常に短くし、待ち時間の後、封入容器の中の圧力の変化の測定を行うことができる。別の実施形態によると、切削片に実際に吸収された体積の変動を観測し、測定できるように、連通期間を十分に長くすることができる。
【００２５】
それから、封入容器内の圧力または体積の変化が、切削片の物理パラメータに対して選択された初期値からモデル化され、モデル化された変化が、封入容器内の物理パラメータの測定された変化に対して最も良好に適合するように、岩石の切削片の物理パラメータの値が反復して調整される。
【００２６】
この手続きにより、良好な結果が得られる。切削片の浸透率の値ｋは、コアから得られた基準測定値に完全に合致する。
【００２７】
ＩＩＩ）毛管圧曲線全体の再現
第３段階では、岩石の物理パラメータのモデル化に適していると評価されている経験的関係を使用して、毛管圧曲線Ｐｃの漸近部分を生成可能にした前述の測定値と、浸透率ｋの測定値とを合成する。以下の刊行物は、毛管圧曲線から浸透率を評価する方法を記述している。
−Ｔｈｏｍｅｅｒ、Ｊ．Ｈ．Ｍ．：“Ｉｎｔｒｏｄｕｃｔｉｏｎ  ｏｆ  ａ  ｐｏｒｅ  ｇｅｏｍｅｔｒｉｃａｌ  ｆａｃｔｏｒ  ｄｅｆｉｎｅｄ  ｂｙ  ｔｈｅ  ｃａｐｉｌｌａｒｙ  ｐｒｅｓｓｕｒｅ  ｃｕｒｖｅ”，Ｔｒａｎｓ  ＡＩＭＥ，ｖｏｌ．  Ｍａｒｃｈ，ｐｐ．７３−７７，１９６０（Ｔｈｏｍｅｅｒ、Ｊ．Ｈ．Ｍ．：“毛管圧曲線により定まる孔の幾何学的要素の採用”、ＡＩＭＥ論文集、３月巻、７３〜７７ページ、１９６０年
−Ｔｈｏｍｅｅｒ，Ｊ．Ｈ．Ｍ．：“Ａｉｒ  ｐｅｒｍｅａｂｉｌｉｔｙ  ａｓ  ａ  ｆｕｎｃｔｉｏｎ  ｏｆ  ｔｈｒｅｅ  ｐｏｒｅ  ｎｅｔｗｏｒｋ  ｐａｒａｍｅｔｅｒｓ”，Ｔｒａｎｓ  ＡＩＭＥ，ｖｏｌ．  Ａｐｒｉｌ，ｐｐ．８０９−８１４，１９８３（Ｔｈｏｍｅｅｒ，Ｊ．Ｈ．Ｍ．：“孔のネットワークの３つのパラメータの関数としての空気の浸透率”、ＡＩＭＥ論文集、４月巻、８０９〜８１４ページ、１９８３年）
毛管圧曲線Ｐｃは、以下の数式１の形態でモデル化されている。
【００２８】
【数１】
2005525574.tif 000003
【００２９】
ただし、
−Ｇは、毛管圧曲線の曲率を説明している（孔の大きさの分布の形状に関連している）、形状パラメータである。
−Ｐdは、０に等しいＳHgについて外挿された排除圧力。
−Ｖb∞は、実験の最後で、毛管圧が無限大（φ×ＳHgに等しい）の場合に、水銀により占められる体積のパーセント。
【００３０】
このモデルの３つのパラメータの、浸透率との関係は、以下の数式２により表されている。
【００３１】
【数２】
2005525574.tif 000004
【００３２】
以下の刊行物では、著者は、浸透率の値を、水銀圧入曲線上の比率（Ｖb／Ｐc）の最大値に、相関させることを提案している。
−Ｓｗａｎｓｏｎ，Ｂ．Ｆ．：“Ａ  ｓｉｍｐｌｅ  ｃｏｒｒｅｌａｔｉｏｎ  ｂｅｔｗｅｅｎ  ｐｅｒｍｅａｂｉｌｉｔｙ  ａｎｄ  ｍｅｒｃｕｒｙ  ｃａｐｉｌｌａｒｙ  ｐｒｅｓｓｕｒｅｓ”，ＪＰＴ，ｖｏｌ．Ｄｅｃｅｍｂｅｒ，ｐｐ．２４９８−２５０４，１９８１（Ｓｗａｎｓｏｎ，Ｂ．Ｆ．：「浸透率と水銀毛管圧との簡単な相関付け」、ＪＰＴ、１２月巻、２４９８〜２５０４ページ、１９８１年）
この特定の点は、著しい毛管圧の増大の直前の、浸透状態の最後に発生する、状態変化におおむね対応している。多孔質の媒体に流入する液体の量に関して、この点は、すべての孔のネットワークが接続される場合の孔の大きさ、したがって、流れを制御する孔の大きさを表すため、特に重要である。著者が示した、もっとも一般的な相関は次の数式３（Ｓｗａｎｓｏｎ  １９８１）で表される。
【００３３】
【数３】
2005525574.tif 000005
【００３４】
以下の刊行物も、水銀圧入曲線から浸透率の値を求めるために、経験的相関や物理モデルを使用しているアプローチの比較評価に関連している。
−Ｋａｍａｔｈ，Ｊ．：“Ｅｖａｌｕａｔｉｏｎ  ｏｆ  ａｃｃｕｒａｃｙ  ｏｆ  ｅｓｔｉｍａｔｉｎｇ  ａｉｒ  ｐｅｒｍｅａｂｉｌｉｔｙ  ｆｒｏｍ  ｍｅｒｃｕｒｙ  ｉｎｊｅｃｔｉｏｎ  ｄａｔａ”，ＳＰＥ  Ｆｏｒｍａｔｉｏｎ  ｅｖａｌｕａｔｉｏｎ，ｖｏｌ．７，ｐｐ．３０４−３１０，１９８２（Ｋａｍａｔｈ，Ｊ．：「水銀注入データからの空気の浸透率の予測精度の評価」、ＳＰＥ地層評価、第７巻、３０４〜３１０ページ、１９８２年）
Ｓｗａｎｓｏｎ  １９８１による前述の文献の特性長さに基づいた新しい相関により、最も良好な一致が得られる。
【００３５】
【数４】
2005525574.tif 000006
【００３６】
【数５】
2005525574.tif 000007
【００３７】
Ｌｍａｘは、次の数式６で定義される。
【００３８】
【数６】
2005525574.tif 000008
【００３９】
ただし、
−λ：毛管圧曲線の曲率（孔の大きさの分布の形状に関連している）を表している指数。
−Ｐｃ：０に等しいＳHgについて外挿された排除圧力。
−Ｓr：湿潤用液体により占められている残留飽和率（％）
Ｔｈｏｍｅｅｒ、Ｓｗａｎｓｏｎ、またはＫａｍａｔｈのアプローチから得られた経験的な毛管圧曲線Ｐｃをパラメータ化するために、毛管圧曲線Ｐｃは、本方法の第１段階の間に遠心分離により得られた漸近部分に適合するように制約されている。曲線全体も、本方法の第２段階の間に切削片について測定された浸透率ｋの値を使用して、制約されていて、測定された浸透率ｋは経験的な関係から得られた結果と比較される。それから、測定された漸近的なふるまいと浸透率の予測値が一致するまで、毛管圧Ｐｃのパラメータが修正され、それにより、全飽和領域Ｓｗにわたって、毛管圧曲線を制約することができる。
【００４０】
水による飽和率が低い毛管圧Ｐｃの漸近的なふるまいを較正可能な、一式の既定値のパラメータを使用して、反転過程が開始される。それから、水による飽和率が低い場合のＰｃの測定値との良好な一致を保ったまま、前述の関係の１つから得られた浸透率の予測値が、切削片について得られた浸透率の測定値と良好に一致するように、これらのパラメータ（主に入力圧力ＰｅまたはＰｄあるいは形状要素λまたはＧ）が段階的に修正される。
１．１得られた結果
図２Ｂ、３Ｂ、４Ｂ、および５Ｂは、前述の手続きにより再現された曲線Ｐｃと、コアについて得られた基準曲線Ｐｃとの比較を示している。どのような例でも、この再現方法により、岩石の浸透率に関わらず、飽和領域全体で、そして、特に水による飽和率が高い場合に、適切なＰｃの変化が得られることがわかる。
【００４１】
【表１】
2005525574.tif 000009
【図面の簡単な説明】
【００４２】
【図１】切削片に対して実施された遠心分離の様々な段階でのＮＭＲ信号の変化を示している図である。
【図２Ａ】１つの岩石（Ｂ７）について、切削片とコアから得られた曲線Ｐｃを示している図である。
【図２Ｂ】コアについて測定された基準曲線Ｐｃと、遠心分離による測定と浸透率の測定から飽和領域全体に渡って再現された曲線Ｐｃとの比較を示している図である。
【図３Ａ】図２Ａに示している曲線に相当する、別の岩石ＧＤＶ１についての結果を示している図である。
【図３Ｂ】図２Ｂに示している曲線に相当する、別の岩石ＧＤＶ１についての結果を示している図である。
【図４Ａ】図２Ａに示している曲線に相当する、別の岩石Ｒｏｔ１についての結果を示している図である。
【図４Ｂ】図２Ｂに示している曲線に相当する、別の岩石Ｒｏｔ１についての結果を示している図である。
【図５Ａ】図２Ａに示している曲線に相当する、別の岩石Ｓｔ  Ｍａｘについての結果を示している図である。
【図５Ｂ】図２Ｂに示している曲線に相当する、別の岩石Ｓｔ  Ｍａｘについての結果を示している図である。
【図６】基準の浸透率と、前述のフランス特許出願公開明細書第０２／００２３号で説明されている方法の範疇内で切削片について測定された浸透率との比較を示している図である。
【図７】本方法を実装した流れ図の概略を示している図である。
【図１】
2005525574.tif 000010
【図２Ａ】
2005525574.tif 000011
【図２Ｂ】
2005525574.tif 000012
【図３Ａ】
2005525574.tif 000013
【図３Ｂ】
2005525574.tif 000014
【図４Ａ】
2005525574.tif 000015
【図４Ｂ】
2005525574.tif 000016
【図５Ａ】
2005525574.tif 000017
【図５Ｂ】
2005525574.tif 000018
【図６】
2005525574.tif 000019
【図７】
2005525574.tif 000020
【国際調査報告】
2005525574.tif 000021
2005525574.tif 000022
2005525574.tif 000023
2005525574.tif 000024
2005525574.tif 000025
2005525574.tif 000026