特許 6473524 高圧条件下の岩石熱物性測定システム及び方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6473524

(24)【登録日】2019年2月1日

(45)【発行日】2019年2月20日

(54)【発明の名称】高圧条件下の岩石熱物性測定システム及び方法

(51)【国際特許分類】

   G01N 25/18 20060101AFI20190207BHJP

   G01N 25/00 20060101ALI20190207BHJP

【ＦＩ】

   G01N25/18 E

   G01N25/00 K

【請求項の数】9

【全頁数】21

(21)【出願番号】特願2017-567374(P2017-567374)

(86)(22)【出願日】2016年4月19日

(65)【公表番号】特表2018-518686(P2018-518686A)

(43)【公表日】2018年7月12日

(86)【国際出願番号】CN2016079687

(87)【国際公開番号】WO2017152473

(87)【国際公開日】20170914

【審査請求日】2017年12月26日

(31)【優先権主張番号】201610130643.7

(32)【優先日】2016年3月8日

(33)【優先権主張国】CN

(73)【特許権者】

【識別番号】511282070

【氏名又は名称】中国科学院南海海洋研究所

(74)【代理人】

【識別番号】100095407

【弁理士】

【氏名又は名称】木村 満

(74)【代理人】

【識別番号】100109449

【弁理士】

【氏名又は名称】毛受 隆典

(74)【代理人】

【識別番号】100132883

【弁理士】

【氏名又は名称】森川 泰司

(74)【代理人】

【識別番号】100148633

【弁理士】

【氏名又は名称】桜田 圭

(74)【代理人】

【識別番号】100147924

【弁理士】

【氏名又は名称】美恵 英樹

(72)【発明者】

【氏名】楊 小秋

(72)【発明者】

【氏名】林 為人

(72)【発明者】

【氏名】多田井 修

(72)【発明者】

【氏名】徐 子英

(72)【発明者】

【氏名】施 小斌

【審査官】 野田 華代

(56)【参考文献】

【文献】 国際公開第２０１７／１４０００６（ＷＯ，Ａ１）

【文献】 特開平０８−２４７９７８（ＪＰ，Ａ）

【文献】 特開２０００−３２９７１９（ＪＰ，Ａ）

【文献】 特公昭４６−０３１３９４（ＪＰ，Ｂ１）

【文献】 中国特許出願公開第１０１９４９８０３（ＣＮ，Ａ）

【文献】 特開２０１５−２１４８０５（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｎ ２５／１８

Ｇ０１Ｎ ２５／００

(57)【特許請求の範囲】

【請求項1】

伝圧媒質が充填された第１キャビティが形成された第１耐圧タンクと、
伝圧媒質が充填された第２キャビティが形成された第２耐圧タンクと、
前記第１耐圧タンクに前記伝圧媒質を供給し、第１排出弁と第１圧力センサーが取り付けられた第１連通ダクトにより前記第１キャビティに接続された高圧ポンプと、
その入力端に、第１温度センサー、第２温度センサー、及び第３温度センサーの出力端が接続された温度観測モジュールと、
その入力端に、前記第１圧力センサー及び第２圧力センサーの出力端が接続された封圧観測モジュールと、
を含み、
前記第２キャビティ内に岩石サンプルが取り付けられ、前記岩石サンプルの中心に前記第１温度センサーが取り付けられ、前記岩石サンプルの外面に前記第２温度センサーが取り付けられ、前記伝圧媒質が充填された前記第２キャビティ内に前記第３温度センサーが取り付けられ、前記第１キャビティと前記第２キャビティとは第２排出弁と前記第２圧力センサーとが取り付けられた第２連通ダクトにより連通され、前記第２キャビティには第３排出弁が接続されていることを特徴とする高圧条件下の岩石熱物性測定システム。

【請求項2】

前記岩石サンプルの外面には、前記岩石サンプルを水密状態にするためのゴムカバーが設けられ、前記岩石サンプルの上端と下端とは硬質シリコンで密封されていることを特徴とする請求項１に記載の高圧条件下の岩石熱物性測定システム。

【請求項3】

前記岩石サンプルは円柱状であることを特徴とする請求項１に記載の高圧条件下の岩石熱物性測定システム。

【請求項4】

前記伝圧媒質はシリコンオイルであることを特徴とする請求項１に記載の高圧条件下の岩石熱物性測定システム。

【請求項5】

第１温度センサーを作製された円柱状の岩石サンプルの中心に取り付け、第２温度センサーを作製された円柱状の前記岩石サンプルの外面に取り付け、ゴムカバーにより前記岩石サンプルを水密状態にすると共に、硬質シリコンにより前記岩石サンプルの上端と下端を密封することにより、岩石サンプルモジュールを形成するステップ１と、
前記岩石サンプルモジュールと第３温度センサーを第２耐圧タンクに入れ、伝圧媒質が前記第２耐圧タンクに充填した後前記第２耐圧タンクを密封する同時に、第１排出弁と第１圧力センサーとが取り付けられた第１連通ダクトで高圧ポンプと第１耐圧タンクとを接続し、第２排出弁と第２圧力センサーが取り付けられた第２連通ダクトで前記第１耐圧タンクと前記第２耐圧タンクとを接続し、第３排出弁を前記第２耐圧タンクに取り付け、そして、前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーを温度観測モジュールに接続し、前記第１圧力センサーと前記第２圧力センサーとを封圧観測モジュールに接続することにより、岩石熱物性測定システムを形成し、前記温度観測モジュールと前記封圧観測モジュールとを作動させて、温度と封圧に対する観測を開始するステップ２と、
前記第１排出弁のみを開いて、前記第２排出弁と前記第３排出弁とを閉じ、前記高圧ポンプを作動させることにより、前記第１耐圧タンクの封圧を予定の圧力まで上昇させるステップ３と、
前記岩石熱物性測定システム全体のバランスを取ったときに、前記第１排出弁を閉じて、前記第３排出弁を閉じたまま、前記第２排出弁を開くことにより、前記第２耐圧タンクに瞬時に与圧するステップ４と、
前記温度観測モジュールによってリアルタイムで観測した前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーの温度変化と、前記封圧観測モジュールによってリアルタイムで観測した前記第２圧力センサーの封圧変化とに基づいて、有限元数値モデルを使って反転することにより、封圧状態における前記岩石サンプルの熱物性パラメータを得るステップ５と、
を含むことを特徴とする高圧条件下の岩石熱物性測定方法。

【請求項6】

前記ステップ５は、
前記円柱状の岩石サンプルの中心を原点とし、前記円柱状の岩石サンプルの径方向と軸方向で形成されたカラム座標係において、熱伝導性微分方程式に基づいて、有限元数値モデルを確立するステップ５１と、
前記岩石サンプルの熱伝導率をλ、体積熱容量を（ρｃ）とし、地殻によくある岩石の熱伝導率の範囲は０．５〜６．０Ｗ・ｍ^−１・Ｋ^−１、体積熱容量の範囲は０．５×１０^６〜５．０×１０^６Ｊ・ｍ^−３・Ｋ^−１であり、計算領域である

【数1】

における二つのパラメータλと（ρｃ）をｍ等分し、初期の（ｍ＋１）×（ｍ＋１）個のグリッドノード（λ_ｉ、（ρｃ）_ｊ）を獲得する（ただし、ｉ、ｊ＝１、２、３、・・・、ｍ）ステップ５２と、
全てのグリッドノード（λ_ｉ、（ρｃ）_ｊ）を確立した有限元数値モデルに入力し、瞬時に与圧する過程においてリアルタイムで観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）を境界条件として、（λ、（ρｃ））＝（λ_ｉ、（ρｃ）_ｊ）をアナログ計算するときに、前記岩石サンプルの中心の温度変化をＴ_ｉ、ｊ^ｍоｄと示すステップ５３と、
最小二乗法を利用し、前記有限元数値モデルを計算することにより得たＴ_ｉ、ｊ^ｍоｄとリアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）を線形フィッティングし、即ち、

【数2】

当該フィッティング直線傾斜Ｋ_ｉ、ｊと相関係数Ｒ_ｉ、ｊを計算し、ただし、相関係数の計算式は、

【数3】

であり、ただし、ｎはサンプリング合計回数であり、ｔ_ｋは第ｋ回のサンプリングの時刻であり、Ｔ０１（ｔ_ｋ）は瞬時に与圧した後ｔ_ｋ時刻で第１温度センサーが獲得した温度変化であり、１≦ｋ≦ｎとするステップ５４と、
目標関数は

【数4】

と定義され、各グリッドポイントの目標関数値Ｆ（λ_ｉ、（ρｃ）_ｊ）を計算するステップ５５と、
目標関数値が最小のグリッドポイントを探し出し、即ち、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）＝ｍｉｎ｛Ｆ（λ_ｉ、（ρｃ）_ｊ）｝であり、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦ε、ここでεは計算請求に満足するかどうかを判定するために設定された閾値であるのであれば、（λ_ｉ^０、（ρｃ）_ｊ^０）は計算する必要がある岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））であり、そうではないと、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦εを満足するまで、（λ_ｉ^０、（ρｃ）_ｊ^０）を中心とする領域を計算領域とし、グリッドの密度を増加し、ステップ５３に戻すことにより計算して、岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））＝（λ_ｉ^０、（ρｃ）_ｊ^０）を得るステップ５６と、
最後に、熱伝導率λ、体積熱容量（ρｃ）及び熱拡散率κの関係式κ＝λ／（ρｃ）＝λ_ｉ^０／（ρｃ）_ｊ^０から計算して、岩石サンプルの熱拡散率を得るステップ５７と、
を含むことを特徴とする請求項５に記載の高圧条件下の岩石熱物性測定方法。

【請求項7】

前記ステップ５１において、カラム座標係における熱伝導性微分方程式は

【数5】

【数6】

と示され、
その初期条件は

【数7】

であり、
境界条件は岩石熱物性測定システムが観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）により制約され、即ち、

【数8】

であり、
ただし、γは各媒質の断熱応力変化の温度応答係数であり、Ａは封圧変化∂ｐ／∂ｔによる温度変化に対応する熱源であり、ｒ_０は岩石サンプルの半径であり、ｚ_０１は円柱状の岩石サンプルの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離であり、ｚ_０２は密封後の円柱状の岩石サンプルモジュールの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離であることを特徴とする請求項６に記載の高圧条件下の岩石熱物性測定方法。

【請求項8】

リアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）、岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）は、各時刻におけるリアルタイムで観測した温度と瞬時に与圧したときの温度との差であり、即ち、

【数9】

【数10】

【数11】

であり、
ただし、Ｔ_１（ｔ）Ｔ、Ｔ_２（ｔ）Ｔ、Ｔ_３（ｔ）Ｔはそれぞれに瞬時に与圧した後時刻ｔにおける第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ５により得られ、Ｔ_１（０）、Ｔ_２（０）、Ｔ_３（０）はそれぞれに瞬時に与圧した第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ２により得られることを特徴とする請求項７に記載の高圧条件下の岩石熱物性測定方法。

【請求項9】

第１温度センサー、第２温度センサーと第３温度センサーが観測した温度はいずれも安定したときに、岩石熱物性測定システム全体の温度バランスを取ることを特徴とする請求項５に記載の高圧条件下の岩石熱物性測定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、岩石熱物性測定の技術分野に関し、更に具体的には、高圧条件下の岩石熱物性測定システムに関する。

【背景技術】

【0002】

地球の内部の岩石熱物性は、地球の内部の熱構造、熱進化及び地球動力学の研究の中で最も基礎的な物性パラメータである。温度、圧力によって、岩石熱物性には違いがある。したがって、異なる封圧条件下における岩石熱物性測定に対する研究は重要な意義を持つ。

【0003】

従来の高圧条件下の岩石熱物性測定システム及び方法において、熱物性試験プローブ（熱源と温度センサーを含む）と岩石サンプルをあらかじめ組み立て、耐圧タンク内に配置する。高圧ポンプを作動させることにより耐圧タンクに作用する封圧を予定の圧力まで上昇させて一定時間維持し、そして、システム全体の温度が平衡に達した後、熱物性測定システムを作動させ電気加熱する同時に、岩石内部の温度変化を観測する。これにより、異なる封圧条件下で熱物性測定を実行することができる。

【0004】

上記従来の測定システム及び方法において、サージ法熱物性測定に必要な「熱源」として、自発的に電気加熱しなければならない（例えば定常電流で加熱糸を通過する）。従って、熱源と温度センサーを同時に岩石の内部に配置しなければならないので、測定システムは複雑になる。そして、熱物性パラメータ測定の環境温度の恒温性に対する要求は非常に高く、実験室で熱物性測定を行う過程で測定システムは常に空気と接触して、恒温の環境下で測定を行うことは難しい。環境の温度の変動は制御しにくいので、測定結果に大きな影響を与えてしまう。

【0005】

このような測定方法と技術において、サージ法熱物性測定に必要な「熱源」として、自発的に電気加熱しなければならない（例えば定常電流で加熱糸を通過する）。

【発明の概要】

【発明が解決しようとする課題】

【0006】

ところで、実験の結果から見ると、地殻によくある岩石の応力−温度応答係数（ΔＴ／Δσ）は小さくて（２〜６ｍＫ／ＭＰａ）である。一方、伝圧媒質（例えばシリコンオイル）の応力−温度応答係数は地殻によくある岩石の応力−温度応答係数の２００倍以上であり、１３８．７４ ｍＫ／Ｍｐａに達する。従って、封圧は瞬時に上昇した後、岩石サンプルと伝圧媒質との間に、温度の差が生じる。本発明は、耐圧タンク内に作用する封圧が瞬時に上昇する過程における岩石サンプルの表面温度及び中心温度と伝圧媒質の温度変化とをリアルタイムで観測する。そして、有限元数値反転方法を利用することにより、高圧条件下の岩石サンプルの熱物性パラメータ（熱伝導率：ｔｈｅｒｍａｌ ｃｏｎｄｕｃｔｉｖｉｔｙ、熱拡散率：ｔｈｅｒｍａｌ ｄｉｆｆｕｓｉｖｉｔｉｙ、及び体積熱容量：ｖｏｌｕｍｅｔｒｉｃ ｈｅａｔ ｃａｐａｃｉｔｙ）を得る。

【課題を解決するための手段】

【0007】

以上の課題を解決するために、本発明の目的は、岩石サンプルの中心及び表面と、伝圧媒質が充填されたキャビティ内には温度センサーが取り付けられ、急速に排出弁を開くことにより、岩石サンプルに対し瞬時に与圧し、そして、封圧が瞬時に上昇する過程における岩石サンプルの中心及び表面と及び伝圧媒質が充填されたキャビティ内の温度変化を観測し、確立した有限元数値モデルを利用し、最適化の方法を採用することにより、高圧条件下の岩石サンプルの熱物性パラメータを得て、「熱源」を必要としない高圧条件下の岩石熱物性測定システムを提供することである。それによって、「熱源」を必要としないサージ熱物性測定を実現し、岩石熱物性測定システム及び操作手順を大幅に簡略化する。

【0008】

上記目的を達成するために、本発明が採用する技術手段は以下である。

【0009】

高圧条件下の岩石熱物性測定システムは、伝圧媒質が充填された第１キャビティが形成された第１耐圧タンクと、伝圧媒質が充填された第２キャビティが形成された第２耐圧タンクと、前記第１耐圧タンクに前記伝圧媒質を供給し、第１排出弁と第１圧力センサーが取り付けられた第１連通ダクトにより前記第１キャビティに接続された高圧ポンプと、その入力端に第１温度センサー、第２温度センサー、及び第３温度センサーの出力端が接続された温度観測モジュールと、その入力端に、前記第１圧力センサーと前記第２圧力センサーの出力端が接続された封圧観測モジュールと、を含み、前記第２キャビティ内に岩石サンプルが取り付けられ、前記岩石サンプルの中心に前記第１温度センサーが取り付けられ、前記岩石サンプルの外面に前記第２温度センサーが取り付けられ、前記伝圧媒質が充填された前記第２キャビティ内に前記第３温度センサーが取り付けられ、前記第１キャビティと前記第２キャビティとは第２排出弁と第２圧力センサーとが取り付けられた第２連通ダクトにより連通され、前記第２キャビティには第３排出弁が接続されている。

【0010】

前記岩石サンプルの外面には、前記岩石サンプルを水密状態にするためのゴムカバーが設けられ、前記岩石サンプルの上端と下端とは硬質シリコンで密封されている。

【0011】

前記岩石サンプルは円柱状である。

【0012】

前記伝圧媒質はシリコンオイルであり、もちろん、植物油や脱イオン水であってもよい。

【0013】

本発明の他の目的は、岩石サンプルの中心及び表面と伝圧媒質が充填された第２耐圧タンクとにはそれぞれ一つの温度センサーが取り付けられ、急速に排出弁を開くことにより、岩石サンプルに対し瞬時に与圧し、そして、封圧が瞬時に上昇する過程における岩石サンプルの中心及び表面と伝圧媒質の温度変化を観測し、確立した有限元数値モデルを利用し、最適化の方法を採用することにより、高圧条件下の岩石サンプルの熱物性パラメータを得て、「熱源」を必要としない高圧条件下の岩石熱物性測定方法を提供することである。よって、「熱源」を必要としないサージ熱物性測定を実現し、岩石熱物性測定システム及び操作手順を大幅に簡略化する。

【0014】

上記目的を達成するために、本発明が採用する技術手段は以下である。

【0015】

高圧条件下の岩石熱物性測定方法は、
第１温度センサーを作製された円柱状の岩石サンプルの中心に取り付け、第２温度センサーを作製された円柱状の前記岩石サンプルの外面に取り付け、ゴムカバーにより前記岩石サンプルを水密状態にすると共に、硬質シリコンにより前記岩石サンプルの上端と下端を密封することにより、岩石サンプルモジュールを形成するステップ１と、
前記岩石サンプルモジュールと第３温度センサーを第２耐圧タンクに入れ、伝圧媒質が前記第２耐圧タンクに充填した後前記第２耐圧タンクを密封する同時に、第１排出弁と第１圧力センサーとが取り付けられた第１連通ダクトで高圧ポンプと第１耐圧タンクとを接続し、第２排出弁と第２圧力センサーが取り付けられた第２連通ダクトで前記第１耐圧タンクと前記第２耐圧タンクとを接続し、第３排出弁を前記第２耐圧タンクに取り付け、そして、前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーを温度観測モジュールに接続し、前記第１圧力センサーと前記第２圧力センサーとを封圧観測モジュールに接続することにより、岩石熱物性測定システムを形成し、前記温度観測モジュールと前記封圧観測モジュールとを作動させて、温度と封圧に対する観測を開始するステップ２と、
前記第１排出弁のみを開いて、前記第２排出弁と前記第３排出弁とを閉じ、前記高圧ポンプを作動させることにより、前記第１耐圧タンクの封圧を予定の圧力まで上昇させるステップ３と、
前記岩石熱物性測定システム全体のバランスを取ったときに、前記第１排出弁を閉じて、前記第３排出弁を閉じたまま、前記第２排出弁を開くことにより、前記第２耐圧タンクに瞬時に与圧するステップ４と、
前記温度観測モジュールによってリアルタイムで観測した前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーの温度変化と、前記封圧観測モジュールによってリアルタイムで観測した前記第２圧力センサーの封圧変化とに基づいて、有限元数値モデルを使って反転することにより、封圧状態における前記岩石サンプルの熱物性パラメータを得るステップ５と、
を含む。

【0016】

前記ステップ５は、
前記円柱状の岩石サンプルの中心を原点とし、前記円柱状の岩石サンプルの径方向と軸方向で形成されたカラム座標係において、熱伝導性微分方程式に基づいて有限元数値モデルを確立するステップ５１と、
前記岩石サンプルの熱伝導率をλ、体積熱容量を（ρｃ）とし、地殻によくある岩石の熱伝導率の範囲は０．５〜６．０Ｗ・ｍ^−１・Ｋ^−１、体積熱容量の範囲は０．５×１０^６〜５．０×１０^６Ｊ・ｍ^−３・Ｋ^−１であり、計算領域である

【数1】

における二つのパラメータλと（ρｃ）をｍ等分し、初期の（ｍ＋１）×（ｍ＋１）個のグリッドノード（λ_ｉ、（ρｃ）_ｊ）を獲得する（ただし、ｉ、ｊ＝１、２、３、・・・、ｍ）ステップ５２と、
全てのグリッドノード（λ_ｉ、（ρｃ）_ｊ）を確立した有限元数値モデルに入力し、瞬時に与圧する過程においてリアルタイムで観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）を境界条件として（λ、（ρｃ））＝（λ_ｉ、（ρｃ）_ｊ）をアナログ計算するときに、前記岩石サンプルの中心の温度変化をＴ_ｉ、ｊ^ｍоｄと示すステップ５３と、
最小二乗法を利用し、前記有限元数値モデルを計算することにより得たＴ_ｉ、ｊ^ｍоｄとリアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）を線形フィッティングし、即ち、

【数2】

当該フィッティング直線傾斜Ｋ_ｉ、ｊと相関係数Ｒ_ｉ、ｊを計算し、ただし、相関係数の計算式は、

【数3】

である（ただし、ｎはサンプリング合計回数であり、ｔ_ｋは第ｋ回のサンプリングの時刻であり、Ｔ０１（ｔ_ｋ）は瞬時に与圧した後ｔ_ｋ時刻で第１温度センサーが獲得した温度変化であり、１≦ｋ≦ｎ）ステップ５４と、
目標関数は

【数4】

と定義され、各グリッドポイントの目標関数値Ｆ（λ_ｉ、（ρｃ）_ｊ）を計算するステップ５５と、
目標関数値が最小のグリッドポイントを探し出し、即ち、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）＝ｍｉｎ｛Ｆ（λ_ｉ、（ρｃ）_ｊ）｝であり、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦ε（εは計算請求に満足するかどうかを判定するために設定された閾値である）であれば、（λ_ｉ^０、（ρｃ）_ｊ^０）は計算する必要がある岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））であり、そうではないと、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦εを満足するまで、（λ_ｉ^０、（ρｃ）_ｊ^０）を中心とする領域を計算領域とし、グリッドの密度を増加し、ステップ５３に戻すことにより、計算して、岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））＝（λ_ｉ^０、（ρｃ）_ｊ^０）を得るステップ５６と、
最後に、熱伝導率λ、体積熱容量（ρｃ）及び熱拡散率κの関係式κ＝λ／（ρｃ）＝λ_ｉ^０／（ρｃ）_ｊ^０から、計算して、岩石サンプルの熱拡散率を得るステップ５７と、
を含む。

【0017】

前記ステップ５１において、カラム座標係における熱伝導性微分方程式は

【数5】

【数6】

と示され、
その初期条件は

【数7】

であり、境界条件は岩石熱物性測定システムが観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）により制約され、即ち、

【数8】

であり、
ただし、γは各媒質の断熱応力変化の温度応答係数であり、Ａは封圧変化による温度変化∂ｐ／∂ｔに対応する熱源であり、ｒ_０は岩石サンプルの半径であり、ｚ_０１は円柱状の岩石サンプルの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離であり、ｚ_０２は密封後の円柱状の岩石サンプルモジュールの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離である。

【0018】

リアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）、岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）は、各時刻におけるリアルタイムで観測した温度と瞬時に与圧したときの温度との差であり、即ち、

【数9】

【数10】

【数11】

であり、
ただし、Ｔ_１（ｔ）Ｔ、Ｔ_２（ｔ）Ｔ、Ｔ_３（ｔ）Ｔはそれぞれに瞬時に与圧した後時刻ｔにおける第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ５により得られ、Ｔ_１（０）、Ｔ_２（０）、Ｔ_３（０）はそれぞれに瞬時に与圧した第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ２により得られる。

【0019】

第１温度センサー、第２温度センサーと第３温度センサーが観測した温度はいずれも安定したときに、岩石熱物性測定システム全体の温度バランスが取れる。

【発明の効果】

【0020】

本発明において、二つの耐圧タンク（内部の封圧は異なって、岩石サンプルを封圧が低い第２耐圧タンクに配置する）の間に位置するバルブを開くことにより、封圧が低い第２耐圧タンクの封圧を瞬時に上昇させる。実験の結果から見ると、地殻によくある岩石の応力−温度応答係数（ΔＴ／Δσ）は小さくて（２〜６ｍＫ／ＭＰａ）、伝圧媒質（例えばシリコンオイル、植物油や脱イオン水などの伝圧媒質）の応力−温度応答係数は地殻によくある岩石の応力−温度応答係数の２００倍以上である（例えばシリコンオイルの応力−温度応答係数は１３８．７４ ｍＫ／Ｍｐａに達する）。従って、封圧が瞬時に上昇した後、岩石サンプルと伝圧媒質との間に、温度の差が生じる。耐圧タンク内の封圧（Ｃｏｎｆｉｎｉｎｇ ｐｒｅｓｓｕｒｅ）、岩石サンプルの中心、表面と第２耐圧タンクの伝圧媒質の温度変化をリアルタイムで観測し、有限元数値反転方法を利用することにより、高圧条件下の岩石サンプルの熱物性パラメータ（熱伝導率：ｔｈｅｒｍａｌ ｃｏｎｄｕｃｔｉｖｉｔｙ、熱拡散率：ｔｈｅｒｍａｌ ｄｉｆｆｕｓｉｖｉｔｉｙ、及び体積熱容量：ｖｏｌｕｍｅｔｒｉｃ ｈｅａｔ ｃａｐａｃｉｔｙ）を得た。本発明は、「熱源」を必要としなくて、岩石サンプルの中心、表面及び第２耐圧タンクの伝圧媒質にはそれぞれに一つの温度センサーが取り付けられ、封圧が瞬時に上昇する過程における岩石サンプルの中心、表面及び第２耐圧タンクの伝圧媒質の温度変化を観測し、確立した有限元数値モデルを利用し、最適化の方法を採用することにより、高圧条件下の岩石サンプルの熱物性パラメータを得るというメリットがある。高圧条件下の岩石熱物性測定システムを提供することである。それによって、「熱源」を必要としないサージ熱物性測定を実現し、高圧条件下の岩石熱物性測定システム及び操作手順を大幅に簡略化する。

【図面の簡単な説明】

【0021】

【図1】本発明に係る水中岩の断熱応力変化の温度応答測定システムの構造概略図である。

【図2】二次元カラム座標係下的有限元数値モデルである。

【図3】龍門山断層帯の砂岩Ｌ２８を瞬時に与圧する過程における温度応答曲線である。

【図4】インドＲａｊａｓｔｈａｎ砂岩ＲＪＳを瞬時に与圧する過程における温度応答曲線である。

【図5】龍門山断層帯の砂岩Ｌ２８の中心温度に対する測定結果と有限元数値モデルに対するシミュレーション結果との比較図。

【図6】インドＲａｊａｓｔｈａｎ砂岩ＲＪＳの中心温度に対する測定結果と有限元数値モデルに対するシミュレーション結果との比較図。

【発明を実施するための形態】

【0022】

以下、図面及び実施例を参照して、本発明の内容をさらに詳しく説明する。

【0023】

（実施例）
図１に示すように、本発明にかかる高圧条件下の岩石熱物性測定システムは、伝圧媒質（例えばシリコンオイル、植物油や脱イオン水など）が充填された第１キャビティ１１が形成された第１耐圧タンク１と、伝圧媒質が充填された第２キャビティ２１が形成された第２耐圧タンク２と、第１耐圧タンク１に伝圧媒質を供給し、第１排出弁５１と第１圧力センサー５２が取り付けられた第１連通ダクト５により第１キャビティ１に接続された高圧ポンプ３と、温度観測モジュール８と、封圧観測モジュール９と、を含む。前記第２キャビティ２１内には、円柱状の岩石サンプルモジュールが取り付けられている。円柱状の岩石サンプルモジュールは、円柱状の岩石サンプル４を含み、その中心には第１温度センサー６１が取り付けられ、その表面には第２温度センサー６２が取り付けられている。そして、円柱状の上硬質シリコン４１と下硬質シリコン４２とをそれぞれ円柱状の岩石サンプル４の上側と下側とに押し付ける。、上硬質シリコン４１、円柱状の岩石サンプル４、及び下硬質シリコン４２は、ゴムカバー４３により覆われ、水密状態にされた後、第２耐圧タンク２の第２キャビティ２１内に配置される。第２キャビティ内には第３温度センサー６３が取り付けられている。また、前記第１キャビティ１１と第２キャビティ２１とは、第２排出弁７１と第２圧力センサー７２とが取り付けられた第２連通ダクト７により連通されている。前記第２キャビティ２１には、第３排出弁２２が接続されている。第１温度センサー６１、第２温度センサー６２、及び第３温度センサー６３の出力端は、温度観測モジュール８の入力端に接続されており、温度観測モジュール８により三つの温度センサーの温度変化をリアルタイムで観測することができる。また、第１圧力センサー５２の出力端と第２圧力センサー７２の出力端とは、封圧観測モジュール９の入力端に接続されている。これにより、封圧観測モジュール９により、第２耐圧タンク２内の封圧変化をリアルタイムで観測することができる。温度観測モジュール８の出力端と封圧観測モジュール９の出力端とは、処理モジュール１０に接続されている。よって、処理モジュール１０により岩石サンプルの断熱応力変化の温度応答係数を計算することができる。なお、処理モジュール１０は高圧ポンプ３の動作を制御できる。これにより封圧観測モジュール９が第１圧力センサー５２により第１耐圧タンク１の封圧が予定の圧力に達したことを検測したときに、処理モジュール１０は高圧ポンプ３の動作を停止する。

【0024】

本発明にかかる高圧条件下の岩石熱物性測定システム及び方法において、まず、高圧ポンプ３を使って第１耐圧タンク１内の封圧を予定の圧力（例えば１３０ＭＰａ）に上昇させ、システム全体の温度を平衡にした後、手動で急速に第１耐圧タンク１と第２耐圧タンク２との間に配置された第２排出弁７１を開く。これにより、第２耐圧タンク２の封圧は１〜２ｓ以内で瞬時に上昇する。伝圧媒質（例えばシリコンオイル）の応力−温度応答係数は岩石サンプルの応力−温度応答係数の２００倍以上であるので、岩石サンプルと伝圧媒質との間に、温度差が生じる。封圧が瞬時に上昇した過程において、岩石サンプルの中心及び表面と、第２耐圧タンクの伝圧媒質との温度変化をリアルタイムで観測する。そして、確立した有限元数値反転モデルを利用し、最適化の方法を採用することにより、高圧条件下の岩石サンプルの熱物性パラメータを得た。これにより、「熱源」を必要としないサージ熱物性測定を実現し、高圧条件下の岩石熱物性測定システム及び方法を大幅に簡略化することができる。

【0025】

有限元数値モデルと方法は以下のように示される。

【0026】

（１．熱伝導性微分方程式）
当該測定システムの岩石サンプルを円柱状とし、簡単に計算するために、ここでは、第１温度センサー６１を岩石サンプルの中心に配置し、且つ岩石サンプルにおいて第２温度センサー６２と第１温度センサー６１は径方向に同一の円にある。従って、対応するカラム座標係（２ｄｒｚ）における熱伝導性微分方程式は、

【数12】

【数13】

と示される。
その初期条件は、

【数14】

である。
境界条件は岩石熱物性測定システムが観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）により制約される。即ち、

【数15】

ただし、λと、ρｃはそれぞれ各媒質の熱伝導率（ｔｈｅｒｍａｌ ｃｏｎｄｕｃｔｉｖｉｔｙ）と体積熱容量（ｖｏｌｕｍｅｔｒｉｃ ｈｅａｔ ｃａｐａｃｉｔｙ）であり、γは各媒質の断熱応力変化の温度応答係数（ａｄｉａｂａｔｉｃ ｐｒｅｓｓｕｒｅ ｄｅｒｉｖａｔｉｖｅ ｏｆ ｔｅｍｐｅｒａｔｕｒｅ）である。また、Ａは封圧変化による温度変化に対応する「熱源」（ｈｅａｔ ｓｏｕｒｃｅ ｔｅｒｍ ｄｒｉｖｅｎ ｂｙ ｃｈａｎｇｅ ｒａｔｅ ｏｆ ｃｏｎｆｉｎｉｎｇ ｐｒｅｓｓｕｒｅ ∂ｐ／∂ｔ）である。

【0027】

図２に示されるように、熱伝導性微分方程式に基づいて、カラム座標係（２ｄｒｚ）において有限元数値モデルを確立する。

【0028】

ステップ２において、前記岩石サンプルの熱伝導率をλ、体積熱容量を（ρｃ）とし、地殻によくある岩石の熱伝導率の範囲は０．５〜６．０Ｗ・ｍ^−１・Ｋ^−１、体積熱容量の範囲は０．５×１０^６〜５．０×１０^６Ｊ・ｍ^−３・Ｋ^−１である。反転の方法の適応性を広げるために、再び岩石の熱物性パラメータの計算領域である、

【数16】

を適当に増大できる。
ただし、二つのパラメータλと（ρｃ）をｍ等分し、初期の（ｍ＋１）×（ｍ＋１）個のグリッドノード（λ_ｉ、（ρｃ）_ｊ）を獲得する（ただし、ｉ、ｊ＝１、２、３、・・・、ｍ）。

【0029】

ステップ３において、全てのグリッドノード（λ_ｉ、（ρｃ）_ｊ）を確立したＰＴ−ＦＥ有限元数値モデルに入力し、瞬時に与圧する過程においてリアルタイムで観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質（例えばシリコンオイル）の温度変化Ｔ０３（ｔ）を境界条件として（図２参照）、（λ、（ρｃ））＝（λ_ｉ、（ρｃ）_ｊ）をアナログ計算するときに、岩石サンプルの中心の温度変化をＴ_ｉ、ｊ^ｍоｄと示す。

【0030】

ステップ４において、最小二乗法を利用し、前記有限元数値モデルを計算することにより得たＴ_ｉ、ｊ^ｍоｄとリアルタイムで観測した岩石サンプルの中心の温度変化Ｔ０１（ｔ）とを線形フィッティングする。即ち、

【数17】

当該フィッティング直線傾斜Ｋ_ｉ、ｊと相関係数Ｒ_ｉ、ｊとを計算する。ただし、相関係数の計算式は、

【数18】

である。

【0031】

ステップ５において、目標関数は、

【数19】

と定義され、各グリッドポイントの目標関数値Ｆ（λ_ｉ、（ρｃ）_ｊ）（ただし、ｉ、ｊ＝１、２、３、・・・、ｍ）を計算する。

【0032】

ステップ６において、目標関数値が最小のグリッドポイントを探し出し、即ち、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）＝ｍｉｎ｛Ｆ（λ_ｉ、（ρｃ）_ｊ）｝を探し出す。Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦ε（εは計算請求に満たすかどうかを判定するために設定された閾値である）であれば、（λ_ｉ^０、（ρｃ）_ｊ^０）は計算する必要がある岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））であり、そうではないと、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦εを満足するまで、（λ_ｉ^０、（ρｃ）_ｊ^０）を中心とする領域を計算領域とし、グリッドの密度を増加し、ステップ３に戻すことにより、計算して、岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））＝（λ_ｉ^０、（ρｃ）_ｊ^０）を得る。

【0033】

ステップ７において、最後に、熱伝導率、体積熱容量及び熱拡散率の関係式κ＝λ／（ρｃ）から、計算して、岩石サンプルの熱拡散率を得る。これにより、ある封圧状態における岩石の熱物性パラメータに対する計算を完成する。

【0034】

本発明にかかる高圧条件下の岩石熱物性測定方法のステップは以下のように示される。

【0035】

ステップ１：第１温度センサー６１を作製された円柱状の岩石サンプル４の中心に取り付け、第２温度センサー６２を岩石サンプル４の外面に取り付ける。そして、ゴムカバーにより岩石サンプル４を水密状態とすることにより、岩石サンプルモジュールを形成する。

【0036】

ステップ２：岩石サンプルモジュールと第３温度センサー６３とを第２耐圧タンク２に入れ、密封する。同時に、第１排出弁５１と第１圧力センサー５２とが取り付けられた第１連通ダクト５で、高圧ポンプ３と第１耐圧タンク１とを接続する。また、第２排出弁７１と第２圧力センサー７２とが取り付けられた第２連通ダクト７で第１耐圧タンク１と第２耐圧タンク２とを接続する。また、第３排出弁２２を第２耐圧タンク２に取り付け、そして、第１温度センサー６１、第２温度センサー６２、及び第３温度センサー６３を温度観測モジュール８に接続する。また、第１圧力センサー５２と第２圧力センサー７２とを封圧観測モジュール９に接続することにより、岩石熱物性測定システムを形成する。そして、温度観測モジュール８と封圧観測モジュール９とを作動させて、温度と封圧に対する観測を開始する。

【0037】

ステップ３：第１排出弁５１のみを開いて、第２排出弁７１と第３排出弁２２とを閉じ、高圧ポンプ３を作動させることにより、第１耐圧タンク１の封圧を予定の圧力まで上昇させる。

【0038】

ステップ４：３〜６時間を経てシステム全体のバランスを取ったときに、第１排出弁５１を閉じて、第３排出弁２２を閉じたまま、第２排出弁７１を素早く開くことにより、第２耐圧タンク２に瞬時に与圧する。

【0039】

上記の操作により、岩石サンプルに対して瞬時に与圧することができ、且つ、この過程において温度の変化と封圧の変化とをリアルタイムで検出し記録することができる。そして、確立した有限元数値モデルと方法を使って反転することにより、封圧状態における岩石サンプルの熱物性パラメータを得ることができる。

【0040】

本発明において、岩石サンプルに対して瞬時に与圧することができることにより、熱物性測定を行うが、ステップ４以降で、第１排出弁５１と第２排出弁７１とを閉じた後で、第３排出弁２２を開くことにより、第２耐圧タンク２の圧力を瞬時に下げるようにしてもよい。これにより、圧力を瞬時に下げることができる。このように、圧力を瞬時に下げる前後における温度の変化と封圧の変化とから、封圧状態における岩石サンプルの熱物性パラメータを得ることもできる。

【0041】

図３と図４は、龍門山断層帯の砂岩Ｌ２８とインドＲａｊａｓｔｈａｎ砂岩ＲＪＳを瞬時に与圧する過程における温度応答曲線である。テーブル１は、１５．３１ＭＰａ、１３．６１Ｍｐａの封圧状態において砂岩Ｌ２８と砂岩ＲＪＳの熱物性パラメータに対する測定結果である。

【0042】

表１は、龍門山断層帯の砂岩（Ｌ２８）とインドＲａｊａｓｔｈａｎ砂岩（ＲＪＳ）の熱物性パラメータを示している。

【表1】

【0043】

図５、図６に示すように、岩石サンプルＬ２８とＲＪＳの中心温度の測定結果と有限元数値のシミュレーション結果を比べると、本発明にかかる方法とシステムは高圧条件下の岩石熱物性パラメータ測定に適用可能であり、岩石熱物性測定システム及び操作手順を大幅に簡略化することができる。

【0044】

具体的な実施例に基づいて本発明を説明したが、本発明の範囲から逸脱しない限り、本発明に対して様々な変更及び等効な置換がなされ得ることは、当業者には容易に理解できるであろう。また、所定な場合または用途では、本発明の範囲から逸脱しない限り、本発明に対して様々な修正を行うことができる。よって、本発明は、記載された具体的な実施例に限定されず、本発明の特許請求の範囲に入る全ての実施形態を含むものとする。

【0045】

（付記）
（付記１）
伝圧媒質が充填された第１キャビティが形成された第１耐圧タンクと、
伝圧媒質が充填された第２キャビティが形成された第２耐圧タンクと、
前記第１耐圧タンクに前記伝圧媒質を供給し、第１排出弁と第１圧力センサーが取り付けられた第１連通ダクトにより前記第１キャビティに接続された高圧ポンプと、
その入力端に、第１温度センサー、第２温度センサー、及び第３温度センサーの出力端が接続された温度観測モジュールと、
その入力端に、前記第１圧力センサー及び前記第２圧力センサーの出力端が接続された封圧観測モジュールと、
を含み、
前記第２キャビティ内に岩石サンプルが取り付けられ、前記岩石サンプルの中心に前記第１温度センサーが取り付けられ、前記岩石サンプルの外面に前記第２温度センサーが取り付けられ、前記伝圧媒質が充填された前記第２キャビティ内に前記第３温度センサーが取り付けられ、前記第１キャビティと前記第２キャビティとは第２排出弁と第２圧力センサーとが取り付けられた第２連通ダクトにより連通され、前記第２キャビティには第３排出弁が接続されていることを特徴とする高圧条件下の岩石熱物性測定システム。

【0046】

（付記２）
前記岩石サンプルの外面には、前記岩石サンプルを水密状態にするためのゴムカバーが設けられ、前記岩石サンプルの上端と下端とは硬質シリコンで密封されていることを特徴とする付記１に記載の高圧条件下の岩石熱物性測定システム。

【0047】

（付記３）
前記岩石サンプルは円柱状であることを特徴とする付記１に記載の高圧条件下の岩石熱物性測定システム。

【0048】

（付記４）
前記伝圧媒質はシリコンオイルであることを特徴とする付記１に記載の高圧条件下の岩石熱物性測定システム。

【0049】

（付記５）
第１温度センサーを作製された円柱状の岩石サンプルの中心に取り付け、第２温度センサーを作製された円柱状の前記岩石サンプルの外面に取り付け、ゴムカバーにより前記岩石サンプルを水密状態にすると共に、硬質シリコンにより前記岩石サンプルの上端と下端を密封することにより、岩石サンプルモジュールを形成するステップ１と、
前記岩石サンプルモジュールと第３温度センサーを第２耐圧タンクに入れ、伝圧媒質が前記第２耐圧タンクに充填した後前記第２耐圧タンクを密封する同時に、第１排出弁と第１圧力センサーとが取り付けられた第１連通ダクトで高圧ポンプと第１耐圧タンクとを接続し、第２排出弁と第２圧力センサーが取り付けられた第２連通ダクトで前記第１耐圧タンクと前記第２耐圧タンクとを接続し、第３排出弁を前記第２耐圧タンクに取り付け、そして、前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーを温度観測モジュールに接続し、前記第１圧力センサーと前記第２圧力センサーとを封圧観測モジュールに接続することにより、岩石熱物性測定システムを形成し、前記温度観測モジュールと前記封圧観測モジュールとを作動させて、温度と封圧に対する観測を開始するステップ２と、
前記第１排出弁のみを開いて、前記第２排出弁と前記第３排出弁とを閉じ、前記高圧ポンプを作動させることにより、前記第１耐圧タンクの封圧を予定の圧力まで上昇させるステップ３と、
前記岩石熱物性測定システム全体のバランスを取ったときに、前記第１排出弁を閉じて、前記第３排出弁を閉じたまま、前記第２排出弁を開くことにより、前記第２耐圧タンクに瞬時に与圧するステップ４と、
前記温度観測モジュールによってリアルタイムで観測した前記第１温度センサー、前記第２温度センサー、及び前記第３温度センサーの温度変化と、前記封圧観測モジュールによってリアルタイムで観測した前記第２圧力センサーの封圧変化とに基づいて、有限元数値モデルを使って反転することにより、封圧状態における前記岩石サンプルの熱物性パラメータを得るステップ５と、
を含むことを特徴とする高圧条件下の岩石熱物性測定方法。

【0050】

（付記６）
前記ステップ５は、
前記円柱状の岩石サンプルの中心を原点とし、前記円柱状の岩石サンプルの径方向と軸方向で形成されたカラム座標係において、熱伝導性微分方程式に基づいて、有限元数値モデルを確立するステップ５１と、
前記岩石サンプルの熱伝導率をλ、体積熱容量を（ρｃ）とし、地殻によくある岩石の熱伝導率の範囲は０．５〜６．０Ｗ・ｍ^−１・Ｋ^−１、体積熱容量の範囲は０．５×１０^６〜５．０×１０^６Ｊ・ｍ^−３・Ｋ^−１であり、計算領域である

【数20】

における二つのパラメータλと（ρｃ）をｍ等分し、初期の（ｍ＋１）×（ｍ＋１）個のグリッドノード（λ_ｉ、（ρｃ）_ｊ）を獲得する（ただし、ｉ、ｊ＝１、２、３、・・・、ｍ）ステップ５２と、
全てのグリッドノード（λ_ｉ、（ρｃ）_ｊ）を確立した有限元数値モデルに入力し、瞬時に与圧する過程においてリアルタイムで観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）を境界条件として、（λ、（ρｃ））＝（λ_ｉ、（ρｃ）_ｊ）をアナログ計算するときに、前記岩石サンプルの中心の温度変化をＴ_ｉ、ｊ^ｍоｄと示すステップ５３と、
最小二乗法を利用し、前記有限元数値モデルを計算することにより得たＴ_ｉ、ｊ^ｍоｄとリアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）を線形フィッティングし、即ち、

【数21】

当該フィッティング直線傾斜Ｋ_ｉ、ｊと相関係数Ｒ_ｉ、ｊを計算し、ただし、相関係数の計算式は、

【数22】

であり、ただし、ｎはサンプリング合計回数であり、ｔ_ｋは第ｋ回のサンプリングの時刻であり、Ｔ０１（ｔ_ｋ）は瞬時に与圧した後ｔ_ｋ時刻で第１温度センサーが獲得した温度変化であり、１≦ｋ≦ｎとするステップ５４と、
目標関数は

【数23】

と定義され、各グリッドポイントの目標関数値Ｆ（λ_ｉ、（ρｃ）_ｊ）を計算するステップ５５と、
目標関数値が最小のグリッドポイントを探し出し、即ち、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）＝ｍｉｎ｛Ｆ（λ_ｉ、（ρｃ）_ｊ）｝であり、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦ε、ここでεは計算請求に満足するかどうかを判定するために設定された閾値であるのであれば、（λ_ｉ^０、（ρｃ）_ｊ^０）は計算する必要がある岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））であり、そうではないと、Ｆ（λ_ｉ^０、（ρｃ）_ｊ^０）≦εを満足するまで、（λ_ｉ^０、（ρｃ）_ｊ^０）を中心とする領域を計算領域とし、グリッドの密度を増加し、ステップ５３に戻すことにより計算して、岩石サンプルの熱伝導率と体積熱容量（λ、（ρｃ））＝（λ_ｉ^０、（ρｃ）_ｊ^０）を得るステップ５６と、
最後に、熱伝導率λ、体積熱容量（ρｃ）及び熱拡散率κの関係式κ＝λ／（ρｃ）＝λ_ｉ^０／（ρｃ）_ｊ^０から計算して、岩石サンプルの熱拡散率を得るステップ５７と、
を含むことを特徴とする付記５に記載の高圧条件下の岩石熱物性測定方法。

【0051】

（付記７）
前記ステップ５１において、カラム座標係における熱伝導性微分方程式は

【数24】

【数25】

と示され、
その初期条件は

【数26】

であり、
境界条件は岩石熱物性測定システムが観測した岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）により制約され、即ち、

【数27】

であり、
ただし、γは各媒質の断熱応力変化の温度応答係数であり、Ａは封圧変化∂ｐ／∂ｔによる温度変化に対応する熱源であり、ｒ_０は岩石サンプルの半径であり、ｚ_０１は円柱状の岩石サンプルの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離であり、ｚ_０２は密封後の円柱状の岩石サンプルモジュールの頂面又は底面から岩石サンプルの中心円柱面までの垂直距離であることを特徴とする付記６に記載の高圧条件下の岩石熱物性測定方法。

【0052】

（付記８）
リアルタイムで観測した前記岩石サンプルの中心の温度変化Ｔ０１（ｔ）、岩石サンプル表面の温度変化Ｔ０２（ｔ）と伝圧媒質の温度変化Ｔ０３（ｔ）は、各時刻におけるリアルタイムで観測した温度と瞬時に与圧したときの温度との差であり、即ち、

【数28】

【数29】

【数30】

であり、
ただし、Ｔ_１（ｔ）Ｔ、Ｔ_２（ｔ）Ｔ、Ｔ_３（ｔ）Ｔはそれぞれに瞬時に与圧した後時刻ｔにおける第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ５により得られ、Ｔ_１（０）、Ｔ_２（０）、Ｔ_３（０）はそれぞれに瞬時に与圧した第１温度センサー、第２温度センサーと第３温度センサーがリアルタイムで観測した温度であり、且つステップ２により得られることを特徴とする付記７に記載の高圧条件下の岩石熱物性測定方法。

【0053】

（付記９）
第１温度センサー、第２温度センサーと第３温度センサーが観測した温度はいずれも安定したときに、岩石熱物性測定システム全体の温度バランスを取ることを特徴とする付記５に記載の高圧条件下の岩石熱物性測定方法。

【符号の説明】

【0054】

１…第１耐圧タンク、１１…第１キャビティ、２…第２耐圧タンク、２１…第２キャビティ、２２…第３排出弁、３…高圧ポンプ、４…岩石サンプル、４１…上硬質シリコン、４２…下硬質シリコン、４３…ゴムカバー、５…第１連通ダクト、５１…第１排出弁、５２…第１圧力センサー、６１…第１温度センサー、６２…第２温度センサー、６３…第３温度センサー、７…第２連通ダクト、７１…第２排出弁、７２…第２圧力センサー、８…温度観測モジュール、９…封圧観測モジュール、１０…処理モジュール。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】