土壌の熱特性

目的

地中温度（地温）は、１日および１年を周期として変動する。その振幅（温度の変化幅）は深くなるにつれて小さくなり、その位相（最大値・最小値を示す時間）も深いほど遅れる。こうした現象は、気温や日射により地表面に流入した熱が、熱伝導によって地中に移動することによっておこる。

本実験では、土壌中に熱電対温度計を挿入し、地表面から周期的に熱を与えたときの地温変化を計測する技術を学ぶ。また、得られたデータから熱伝導率を求める方法について学習する。

熱電対の原理

温度を測る簡単な道具の一つに熱電対（thermocouple）がある。熱電対は用いる種類に応じて-200^oC～2000^oC程度までの広い測定範囲をもち、比較的安価であること、遠隔測定ができること、測温部が小さいことなどの特徴がある。

熱電対は、その名の通り1対の金属線からなる。ある金属に温度差を与えると、金属内の自由電子が熱によって移動し低温側の密度が小さくなる(図-1)。このため、この金属の低温側は負に、高温側は正に帯電する。こうした自由電子の密度の変化は金属の種類によって異なる。そこで、2種類の金属線を図-2や図-3のようにつなぎ、接点A, Bに異なる温度を与えると、接点間に起電力（熱起電力）Eが発生することになる。この現象は発見者T. Seebeck (1821)にちなみ、ゼーベック効果と呼ばれる。

熱起電力の大きさは、金属がそれぞれ均質であれば金属線の組合せと接点間の温度差だけによって決まり、金属線の長さや太さ、接点以外の部分の温度などには無関係である。そこで、一方の接点を基準温度に保ち熱起電力を求めれば、もう一方の接点の温度が測定できる。銅－コンスタンタンの場合、１℃あたりの熱起電力はおよそ40μＶ（熱起電力1mVのときの温度差は25^oC）である。

図-1　ある金属内の温度勾配下の電子密度　　　　　　図-2　起電力の発生　（接点が２つの場合）

図-3　起電力の発生　（接点が１つの場合）

熱伝導の理論

１．基礎方程式

土壌中における熱移動は、フーリエの法則によって表される。

(1)

ここで、Jは熱フラックス(W/m²)、Tは温度(K)、xは距離(m)、λは熱伝導率(W/mK)である。

また、土壌中における熱の保存式は次式（２）で表される。

(2)

　ここで、Hは次式で与えられる体積熱容量(J/m³)、rは体積あたりの発熱量(W/m³)、tは時間(sec)である。

(3)

　ただし、Cは土壌の体積比熱(J/m³K)、T_refは基準温度である。r=0の場合、(2)式に(1)式と(3)式を代入し、Cを一定とすれば、次の熱伝導方程式を得る。

(4)

　さらに、λを一定と仮定すると、(4)式は次式となる。

(5)

ここで、κは次式で定義される熱拡散率(m²/s)である。

κ = λ/C (6)

２．解析解

(4)式と(5)式に適当な境界条件や初期条件を与えると、解析解が得られる。ここでは、地表面に周期的な温度変化が与えられたときの解析解を考える。

いま、境界条件として

at x = 0 (7)

at x = ∞ (8)

を与えると、(5)式の解析解は次式のようになる¹⁾。

(9)

ここで、

(10)

は、地表面温度の振幅が1/e倍（0.368倍）になる深さ（damping depth）、ωは角振動数、τは周期である。(10)式を用いると、土壌の熱拡散率κは次式で求められる。

(11)

３．ｄの求め方

3.1　振幅から求める方法（振幅法）

(9)式より、深さｘにおける地温の最大値、最小値、およびその差は次式となる。

(12)

(13)

(14)

ここで、２点の深さ（x₁, x₂）における地温の最大-最小の差を用いると、

(15)

(16)

3.2　位相から求める方法（位相法）

(9)式において、深さｘの地温は次式の関係を満たすとき最大となる。

(17)

ここで、２点の深さ（x₁, x₂）において地温の最大値を示す時刻を（t₁, t₂）とすると、

(18)

(19)

４．土壌の熱物性

(4)式と(5)式中のC、λ、κは一般的な物質であれば、熱的な性質を表す物質固有の定数である。しかし、土壌は土粒子、水、空気で構成されるので、その構成割合や配列などによってこれらの物性値が変化する。特に、水分量の影響を大きく受けることが土壌の熱物性の特徴である。

4.1土壌の体積比熱

土壌の体積比熱は構成割合に基づいた次式により求められることが多い。

(20)

ここで、θは体積分率、ρは密度、cは質量比熱、添字a, w, siはそれぞれ空気、水、土粒子成分を表す。通常、右辺第一項は他に比べて無視できるので、(20)式は次式のように簡略化できる。

C = ρ_d ( c_s + c_ww )　 (21)

　　ただし、ρ_d は土壌の乾燥密度、w は土壌の含水比、c_s, c_wはそれぞれ土粒子と水の質量比熱で、それぞれ0.84(J/g K)、4.2 (J/g K)である。

4.2 土壌の熱拡散率

実測した地温変化から(11)式により決定する。

4.3 土壌の熱伝導率

　体積比熱と熱拡散率から(6)式を用いて算出する。

実験方法

準備するもの

１．土砂、クロボク土、立川ローム

２．熱電対４本　（長さ1m程度のもの３本）＋素線（１本）

３．試料容器　　　　　　　1組　（土壌試料を入れる底付の円筒容器）

４．ガラス板（金属板） 1枚

５．赤外線ランプ 1個

６．スタンド１脚　（赤外線ランプ固定用）

７．ストップウォッチ 1個

８．天秤１個

９．秤量缶３個

１０．ノギス１個

１１．薬さじ１本

１２．データロガー１組　（ロガー本体・電池・読取器）

実験手順

１．熱電対の素線を使って熱電対温度計を自作する。

２．熱電対４本をデータロガーに接続し、データロガーのスイッチを入れる。

３．データロガーの読取器の操作法を確認する。

４．熱電対の先端を指先で握ったり、水に浸けたりしながら、温度計として動作していることを確認する。

①　*0 でLOG1を確認

②　*6A　データ表示

③　03:　～　06:　チャンネルに温度が表示される

④　*1A1　と入力後、01:　30Aと入力（30秒ごとに温度測定）

⑤　*0 でLOG1を確認

⑥　*6A　でデータ表示

５．試料容器の空質量と内容積を測定する。

６．一様に湿らせた土を試料容器一杯に詰めて、質量を量る。試料は数回にわけて入れ、そのつど突き固め棒を使って均一に詰める。

７．土の表面から1cm, 2cm, 5cmの位置に”つまようじ”で穴をあけ、熱電対の先端を埋設する。熱電対の挿入口と先端から10cmくらいの２箇所をビニールテープで固定する。４本目の熱電対でランプから離れたところの気温を測定っする。

８．ガラス板と土の間に隙間ができないように、試料容器の上にガラス板を密着させる。

９．試料容器の側面と底面を断熱材で覆う。

１０．試料容器の真上30cmのところに赤外線ランプをつるす。

１１．ランプを6分間つけ6分間消すサイクルを２回繰り返し、この間の地温を30秒間隔で測定する。

１２．測定終了後、熱電対先端の位置（深さ）を正確に測定する。

１３．熱電対先端付近の土を３つの秤量缶に採り、含水比wを測定する。　

実験結果の整理

１．各深さの温度と時間の関係をグラフに表す。

２．グラフより適当な区間（12分間）を選ぶ。この間、実測点は25点含まれる。

３．選んだ区間の1番目と25番目の点を直線で結ぶ。

４．「３．」で引いた直線とグラフとの差を求め、次の表を作成する。

時間 t (min)	0	0.5	1	1.5	・・・	11.5	12
温度差(1cm深) T₁	y₀	y₁	y₂	y₃	・・・	y₂₄	y₂₅
温度差(2cm深) T₂	z₀	z₁	z₂	z₃	・・・	z₂₄	z₂₅