ビジネス

ダルシーの法則の導出方法は?成り立ちや証明も解説!(浸透流:ダルシー流速:地下水流動:比例関係など)

ダルシーの法則の基本関係
当サイトでは記事内に広告を含みます

ダルシーの法則の導出方法は?成り立ちや証明も解説!(浸透流:ダルシー流速:地下水流動:比例関係など)

地下水の流れや土中の浸透流を学ぶ際、基礎となるのがダルシーの法則です。

透水係数、動水勾配、流量といった用語が一度に登場するため、式だけを暗記してしまいやすい分野でもあります。

しかし、法則がどのような実験から見いだされ、なぜ流量が水頭差に比例するのかを順に確認すると、地下水流動の計算や透水試験の理解が深まります。

ここではダルシーの法則の導出、成立条件、ダルシー流速の意味、証明の考え方までを、数式と具体例を交えて解説します。

ダルシーの法則の基本関係

ダルシーの法則の基本関係

それではまずダルシーの法則の基本関係について解説していきます。

法則が示す浸透流の比例関係

ダルシーの法則は、土や砂、岩盤などの多孔質媒体の内部を水がゆっくり流れるとき、流量が動水勾配に比例することを示す法則です。

多孔質媒体とは、土粒子の間や岩石の割れ目に細かな空隙があり、その空隙を通って水が移動できる材料を指します。

水は水頭の高い側から低い側へ向かって流れますが、その流れの強さを表す指標が動水勾配です。

水頭差が大きいほど流れは強くなり、流れる距離が長いほど流れは弱くなるという関係を、実務で扱いやすい形にしたものがダルシーの法則といえるでしょう。

一般的な式は、次のように表されます。

Q=kAi

Qは流量、kは透水係数、Aは流れに直交する断面積、iは動水勾配を表します。

動水勾配iは、水頭差hを流路長Lで割った値として、i=h/Lで求めます。

流量Qの単位は通常m³/s、断面積Aの単位はm²、透水係数kの単位はm/sです。

動水勾配iは長さどうしの比であるため、単位を持たない量として扱われます。

式の右辺を確認すると、透水係数が大きい地盤ほど、また断面積が広いほど、単位時間あたりに流れる水の量が増える仕組みです。

ダルシー流速と実際の地下水流速

ダルシーの法則では、流量Qを断面積Aで割った値をダルシー流速と呼びます。

ダルシー流速vは、次の式で整理できます。

v=Q/A=ki

ダルシー流速は、断面全体を水が均一に通過したと仮定した見かけの流速です。

ここで注意したいのは、ダルシー流速が土粒子の間を実際に水が進む速さそのものではない点です。

水が通れるのは断面のすべてではなく、空隙部分に限られます。

そのため地下水の実流速を求める場合には、有効間隙率neを用いて、ダルシー流速を補正する必要があります。

実流速uは、概ねu=v/neとして求められます。

有効間隙率が小さい地盤では、同じダルシー流速であっても空隙内の水はより速く移動します。

地下水汚染の到達時間や浄化計画を検討する場面では、この違いが重要になります。

透水係数が表す地盤の通水性

透水係数kは、水が地盤や材料の中を通り抜けやすいかどうかを示す値です。

砂利や粗い砂では空隙が比較的大きいため透水係数が大きくなり、粘土のように細かな粒子が密に詰まった土では小さくなります。

ただし、透水係数は粒径だけで決まるわけではありません。

粒度分布、締固めの程度、間隙比、飽和状態、水温、土の構造なども結果に影響します。

地盤材料の例 透水性の傾向 地下水流動の特徴
砂利層 非常に高い 比較的速く水が移動しやすい
粗砂層 高い 井戸や集水施設の対象になりやすい
細砂層 中程度 動水勾配によって流量が変化しやすい
シルト層 低い 水の移動に時間がかかりやすい
粘土層 非常に低い 遮水層として利用されることがある

透水係数は現地透水試験や室内透水試験から求められます。

数値を用いる際は、試験条件と実際の地盤条件が近いかどうかまで確認することが大切です。

ダルシーの実験装置と成り立ち

続いてはダルシーの実験装置と成り立ちを確認していきます。

ヘンリー ダルシーによるろ過実験

ダルシーの法則は、フランスの技術者アンリ ダルシーが砂を用いたろ過実験を行った結果として知られています。

ダルシーは、水を通す砂層の上下に水位差を与え、そのときに流れる水量を測定しました。

実験では、砂層の長さ、断面積、水位差を変えながら、流量との関係を比較しています。

その結果、流量は断面積と水位差に比例し、砂層の長さには反比例する傾向が得られました。

この観測結果が、現在の透水計算で用いられる基本式につながっています。

ダルシーの法則は机上で仮定された式ではなく、ろ過実験から得られた経験的な関係です。

水頭差と流路長の意味

水が流れるためには、水頭に差が必要です。

水頭は位置水頭、圧力水頭、速度水頭などで構成されますが、地下水流動では位置と圧力に関する成分が特に重要です。

たとえば、水平な砂層の左右に異なる水位を与えると、水位の高い側にはより大きなエネルギーがあり、水は低い側へ移動します。

両端の水頭差をh、砂層の長さをLとすると、h/Lが動水勾配です。

同じ水頭差であれば、砂層が短いほど単位長さあたりのエネルギー差は大きくなります。

反対に、同じ長さで水頭差を大きくすれば、水を押し出す力が増えるため流量も増加します。

動水勾配は、単なる水面の傾きではありません。

地下水では見えない位置にある圧力水頭も含めて評価するため、観測井の水位やピエゾ水頭を使って求めます。

実験結果から得られる流量式

実験で得られた傾向を式にすると、流量Qは断面積Aと水頭差hに比例し、流路長Lに反比例します。

比例定数をkとすれば、Q=kAh/Lと書けます。

ここでh/Lを動水勾配iに置き換えることで、Q=kAiというよく知られた式になります。

この導出は、実験結果を比例式として整理する手順です。

厳密な意味での理論証明とは異なりますが、現実の地盤で有効な近似関係を明快に示しています。

土質力学や地下水学では、この経験式を出発点として、水みち、揚水、浸透破壊、遮水設計などを考えます。

ダルシーの法則の導出過程

続いてはダルシーの法則の導出過程を確認していきます。

比例式による段階的な導出

ダルシーの法則を理解するには、複数の比例関係に分けて考えると分かりやすくなります。

まず、流れる断面積Aが2倍になれば、同じ条件で通れる水の通路もおおむね2倍になるため、流量Qも2倍になります。

したがって、QはAに比例すると考えられます。

次に、水頭差hを大きくすると、水を移動させるエネルギー差が増えます。

低流速の範囲では、このエネルギー差の増加に合わせて流量も増えるため、Qはhに比例します。

一方で、流路長Lが長くなると、水はより長い経路を通る必要があります。

流動への抵抗が増えるため、QはLに反比例します。

QはAに比例します。

Qはhに比例します。

QはLに反比例します。

これらをまとめると、Q=C Ah/Lとなります。

Cを透水係数kと置き、i=h/Lを用いると、Q=kAiが得られます。

このように、ダルシーの法則の導出は、流量に影響する要素を整理した比例式から始まります。

式を暗記するよりも、断面積は増える方向、水頭差は増える方向、流路長は減る方向に働くと押さえると理解しやすくなります。

圧力損失との関係

水が土の空隙を進むとき、土粒子の表面との摩擦や流路の曲がりによってエネルギーが失われます。

このエネルギー損失が、水頭の低下として表れます。

ダルシー流れの範囲では、圧力損失と流速の関係がほぼ線形になります。

すなわち、流速を2倍にすると、必要な水頭差も概ね2倍になる関係です。

管路流れでは条件によって損失が流速の二乗に近づくことがありますが、多孔質媒体内の遅い流れでは線形関係が基本になります。

この違いは、流れの状態と抵抗の生じ方が異なるためです。

土中の水の流れを一般的な配管の流れと同じ感覚で扱わないことが、理解のポイントになります。

連続の式との組み合わせ

地下水流動を二次元や三次元で扱う場合、ダルシーの法則だけでは足りません。

ある場所に流れ込む水量と流れ出る水量、さらに貯留量の変化を考える連続の式と組み合わせます。

定常状態では、ある領域に入る水量と出る水量が釣り合います。

この条件にダルシーの法則を代入すると、地下水面や水頭分布を求めるための基礎式が導かれます。

透水係数がどの方向にも同じ等方性地盤では、比較的単純な形で計算できます。

地層の水平と鉛直で透水係数が異なる異方性地盤では、方向ごとの透水性を考慮しなければなりません。

実際の地盤では層構造や局所的な割れ目もあるため、モデル化の前提を明確にする必要があります。

ダルシーの法則が成り立つ条件

続いてはダルシーの法則が成り立つ条件を確認していきます。

層流領域にある浸透流

ダルシーの法則は、主に水がゆっくり流れる層流領域で成り立ちます。

層流とは、水の流れが比較的整っており、乱れが小さい状態です。

土の空隙は複雑な形をしていますが、流速が小さい範囲では、流量と動水勾配の比例関係が保たれやすくなります。

流速が大きくなると、慣性力の影響が無視しにくくなり、流量と動水勾配の関係が直線から外れる場合があります。

特に粗い砂利、破砕帯、岩盤の大きな割れ目などでは、非ダルシー流れへの注意が必要です。

ダルシーの法則を使う前に、対象となる流れが低流速で安定しているかを確認することが重要です。

飽和地盤と不飽和地盤の違い

飽和地盤では、土の空隙がほぼ水で満たされているため、透水係数を一定値として扱いやすい傾向があります。

地下水面より下の帯水層では、この考え方を基本として地下水流動を計算します。

一方、不飽和地盤では空隙の一部に空気が存在し、水の通り道が連続しない場合もあります。

不飽和状態では、含水比やサクションの変化によって透水性が大きく変わります。

このため、飽和透水係数だけを使って不飽和浸透を精密に表現することは難しいでしょう。

降雨浸透や表層土の排水を扱うときは、不飽和透水係数や土壌水分特性曲線を考慮する場合があります。

透水係数を一定とみなす前提

ダルシーの法則を使った単純な計算では、透水係数kが場所や時間によらず一定であると仮定することがあります。

ただし実地盤は均質ではなく、砂層と粘土層が交互に分布していることも珍しくありません。

透水係数が層ごとに異なる場合、水は流れやすい層へ集中しやすくなります。

また、細粒分の移動、目詰まり、化学反応、地盤の圧密などにより、透水係数が時間とともに変化するケースもあります。

計算式が簡潔でも、透水係数の設定は慎重に行う必要があります。

室内試験の値だけで判断せず、地層構成、現地透水試験、地下水位の観測結果を合わせて評価することが望まれます。

設計上の安全性が求められる場面では、ばらつきを考慮した代表値や保守的な値を選ぶこともあります。

ダルシーの法則の証明と検証方法

続いてはダルシーの法則の証明と検証方法を確認していきます。

経験則としての証明の考え方

ダルシーの法則は、もともと実験で確かめられた経験則です。

そのため数学の定理のように、限られた前提から完全に演繹する形式の証明とは少し性質が異なります。

法則の妥当性は、条件を管理した透水試験で水頭差と流量を測り、両者の関係が直線になることによって検証できます。

横軸に動水勾配i、縦軸にダルシー流速vを取ったグラフを考えてみましょう。

ダルシーの法則が成り立つ範囲では、原点付近を通る直線的な関係が得られます。

その傾きが透水係数kに対応します。

実験値が直線関係を示すことが、ダルシーの法則の最も基本的な検証になります。

毛細管束モデルによる説明

ダルシーの法則を理論的に理解する方法の一つに、土の空隙を多数の細い毛細管の束とみなす考え方があります。

一本の円管を通る粘性流れでは、圧力差と流量の関係を示すポアズイユの法則が知られています。

ポアズイユの法則では、低速の層流において流量は圧力差に比例します。

多孔質媒体の内部には複雑で曲がりくねった流路が多数存在しますが、それらを平均的に扱うと、全体の流量も水頭差に比例する形で表せます。

この平均化された比例係数が、地盤固有の透水係数として現れます。

もちろん実際の土の空隙は円管ではなく、太さも方向も一様ではありません。

それでも、低速域における粘性抵抗の支配という共通点が、ダルシーの法則の物理的な裏付けになります。

非ダルシー流れとの比較

流れが速くなると、粘性抵抗に加えて水の慣性の影響が大きくなります。

この場合、動水勾配と流速の関係は単純な比例では表しにくくなります。

代表的な拡張として、流速の二乗項を含めるフォルヒハイマー型の式が使われることがあります。

比較項目 ダルシー流れ 非ダルシー流れ
主な支配要因 粘性抵抗 粘性抵抗と慣性力
流速と動水勾配 おおむね比例 直線関係から外れる
想定されやすい地盤 砂、シルト、一般的な帯水層 砂利、破砕岩、大きな割れ目
計算上の扱い Q=kAi 追加の抵抗項を考慮する

地下水解析でダルシーの法則を用いるときは、法則そのものを疑うよりも、適用範囲を超えていないかを点検する姿勢が重要です。

地下水流動計算への活用

続いては地下水流動計算への活用を確認していきます。

一次元浸透流の計算例

ここでは、水平な砂層を通る一次元の浸透流を考えます。

透水係数kが0.001m/s、流れに直交する断面積Aが2m²、水頭差hが0.5m、流路長Lが10mとします。

動水勾配はi=0.5/10=0.05です。

これをダルシーの法則に代入すると、流量Qは0.001×2×0.05となります。

Q=kAi

Q=0.001×2×0.05

Q=0.0001m³/s

この値は毎秒0.0001m³、すなわち毎秒0.1Lの流量に相当します。

計算では単位のそろえ方が重要です。

透水係数をcm/sで与えられた場合や、断面積をcm²で測った場合には、単位を統一してから計算します。

数値そのものだけでなく、m、s、m²、m³の対応を確認すると計算ミスを防ぎやすくなります。

揚水と地下水位低下の検討

井戸から地下水をくみ上げると、周囲の地下水位は低下します。

地下水位が低くなった場所へ向かって水が流れるため、帯水層内には新たな動水勾配が生じます。

ダルシーの法則は、この流れの大きさを把握する基礎になります。

実際の揚水計算では、井戸に向かう流れが放射状になることや、帯水層の厚さ、境界条件、貯留係数などを考慮します。

単純な一次元式だけでは十分でない場合でも、水頭差が地下水を動かし、透水性が流量を左右するという考え方は変わりません。

井戸周辺の地盤沈下や周囲の井戸への影響を検討する際にも、地下水流動の基本原理として役立ちます。

土木工事と浸透破壊の評価

ダム、堤防、掘削工事、地下構造物では、浸透流が安全性に直接関係します。

水が地盤内を流れることで、土粒子を持ち去る現象が生じると、パイピングやボイリングにつながるおそれがあります。

特に下向きではなく上向きに浸透流が働く場合、土の有効応力が小さくなり、地盤が不安定になることがあります。

この評価では、動水勾配と限界動水勾配の比較が重要です。

動水勾配が大きい場所ほど、単に水量が増えるだけでなく、地盤の安定性にも影響する可能性があります。

浸透流の計算は、排水量を求める作業だけではありません。

水圧、揚圧力、地盤の有効応力、浸透破壊のリスクを一体で確認することが、土木設計では求められます。

遮水壁やフィルター層、排水工を設ける計画では、流れをどこへ導き、どの程度弱めるかを考える必要があります。

ダルシーの法則のまとめ

ダルシーの法則は、多孔質媒体を通る浸透流の流量が、透水係数、断面積、動水勾配によって決まることを表す基本法則です。

基本式はQ=kAiであり、ダルシー流速はv=kiとして整理できます。

導出の考え方は、流量が断面積と水頭差に比例し、流路長に反比例するという実験的な関係に基づいています。

水頭差を流路長で割った動水勾配が、地下水を動かす強さを表すと理解すると、式の意味がつかみやすくなるでしょう。

ただし、法則が有効なのは主に低流速の層流領域であり、粗い地盤や高速の流れでは非ダルシー流れを検討する必要があります。

また、ダルシー流速は見かけの流速であるため、実際の地下水移動速度を扱うときは有効間隙率による補正も欠かせません。

透水試験、地下水位観測、地層構成の把握を組み合わせながら、対象地盤に合った透水係数と境界条件を設定することが大切です。

ダルシーの法則を比例関係として理解できれば、地下水流動、排水設計、揚水、浸透破壊の検討まで、幅広い場面で応用しやすくなります。