地盤に水が浸透するとき、その流れを支配する基本法則がダルシー則です。
飽和透水係数の計算は、このダルシー則を正確に理解し適用することから始まります。
実際の設計現場では、透水量や動水勾配、断面積といった変数を組み合わせて透水係数を求めたり、逆に透水係数から浸透流量を計算したりする作業が頻繁に発生します。
本記事では、ダルシー則の基本式から出発し、透水量・動水勾配・透水断面積を使った計算方法を丁寧に解説します。
さらに、フローネットや浸透流解析への展開についても取り上げ、実践的な知識の習得をサポートします。
目次
飽和透水係数の計算方法：ダルシー則が基本の出発点
それではまず、飽和透水係数の計算方法とダルシー則の基本について解説していきます。
ダルシー則（Darcy’s law）は、飽和した多孔質媒体（土や岩盤など）を通る流れを記述する基本法則です。
1856年にフランスの技術者アンリ・ダルシーが砂層を通る水の流れを実験的に研究して発見したことから、その名が付けられました。
ダルシー則の基本式：Q = k × i × A
ここで Q：透水量（m³/s）、k：透水係数（m/s）、i：動水勾配（無次元）、A：透水断面積（m²）
この式が飽和透水係数の計算すべての基礎となります。
ダルシー則が成立するためには、いくつかの前提条件があります。
まず、流れはなめらかで乱流でないこと（層流条件）、次に土が均質・等方性であること、そして完全飽和であることが必要です。
礫のような高透水性材料では乱流が生じることがあり、その場合はダルシー則が厳密には成立しないため注意が必要です。
レイノルズ数が1～10程度以下であればダルシー則が適用可能とされています。
動水勾配と透水量の計算
続いては、動水勾配と透水量の計算について確認していきます。
ダルシー則を使いこなすためには、動水勾配（i）の正確な理解が不可欠です。
動水勾配の定義と計算方法
動水勾配（i）とは、流れの方向に沿った単位長さあたりの全水頭の損失量を表す無次元の量です。
動水勾配の定義式：i = Δh / L
ここで Δh：水頭差（m）、L：流れの経路長（m）
例：2点間の水頭差が0.5 m、流れの距離が2.5 m の場合、i = 0.5 / 2.5 = 0.2
全水頭（全ポテンシャルヘッド）は、位置水頭・圧力水頭・速度水頭の和として定義されますが、多孔質体内の流れでは速度が小さいため速度水頭は無視されることが多いです。
動水勾配の方向は、全水頭が高い方から低い方へ向かいます。
これは川の水が高いところから低いところへ流れるのと同じ原理です。
透水量の計算手順
透水量（Q）はダルシー則の式から直接計算できます。
実際の計算では、k・i・A の値をそれぞれ確定し、単位を統一したうえで式に代入します。
【計算例】
条件：k = 1.0 × 10⁻⁴ m/s、Δh = 1.0 m、L = 5.0 m、断面積 A = 0.1 m²
動水勾配：i = 1.0 / 5.0 = 0.2
透水量：Q = 1.0 × 10⁻⁴ × 0.2 × 0.1 = 2.0 × 10⁻⁶ m³/s
1日あたりの透水量：Q × 86400 s = 0.173 m³/日
このように、透水係数・水頭差・断面積・距離という4つの情報から透水量を計算できます。
逆に、試験で流量・水頭差・断面積・距離を計測することで透水係数を求めることも可能です。
浸透速度と実際の流速の違い
ダルシー則で求められる流速 v = k × i は、ダルシー流速（見かけの流速）と呼ばれます。
これは断面積全体（土粒子の部分も含む）に対する平均的な流速であり、実際に水が間隙内を移動する速度（実流速）とは異なります。
実流速は、間隙を通過する経路の長さと間隙率によって補正が必要です。
実流速（真の流速）：v_r = v / n = (k × i) / n
ここで n：間隙率（土全体の体積に対する間隙体積の割合）
砂の場合 n ≈ 0.3～0.4 程度であるため、実流速はダルシー流速の2.5～3.3倍になります。
汚染物質の移動や地下水年代の計算では、実流速の使用が不可欠です。
ダルシー流速だけを使うと移動速度を過小評価してしまう可能性があるため注意が必要です。
フローネットを用いた浸透流解析
続いては、フローネットを用いた浸透流解析について確認していきます。
フローネット（流線網）は、二次元浸透流問題を図的に解く手法です。
ダルシー則と連続の方程式を組み合わせたラプラス方程式の解をグラフィカルに表現したものです。
フローネットの構成要素
フローネットは、流線（streamlines）と等ポテンシャル線（equipotential lines）が直交するメッシュ状の網で構成されます。
流線は水が移動する経路を示し、等ポテンシャル線は全水頭が等しい点を結んだ線です。
この2組の曲線群が直交して形成する各要素（フローエレメント）は、理想的な正方形に近い形になるように描くことが基本原則です。
手描きでフローネットを描く際には、流線の本数（Nf）と等ポテンシャル線の区間数（Nd）の比（Nf/Nd）が透水量計算のポイントになります。
フローネットからの透水量計算
フローネットを描いた後、単位幅あたりの全透水量 Q は次の式で求められます。
Q = k × Δh × (Nf / Nd)
ここで、k：透水係数（m/s）、Δh：全水頭差（m）、Nf：流線管の数、Nd：等ポテンシャル線の区間数
例：k = 5 × 10⁻⁵ m/s、Δh = 3 m、Nf = 4、Nd = 12 の場合
Q = 5 × 10⁻⁵ × 3 × (4/12) = 5.0 × 10⁻⁵ m³/(s・m)
フローネット法はダム・堰・矢板など二次元的な浸透流問題に広く適用されています。
手書きでの解析は経験と試行錯誤が必要ですが、近年ではFEM（有限要素法）や有限差分法を用いた数値解析が実務でも普及しています。
有効応力と浸透力の関係
浸透流が生じると、土粒子には水から浸透力（seepage force）が作用します。
上向き方向の浸透力は有効応力を減少させ、限界動水勾配（i_c）に達すると地盤が液状化（クイックサンド）する危険性があります。
限界動水勾配：i_c = (Gs – 1) / (1 + e)
ここで Gs：土粒子の比重（≈2.65）、e：間隙比
砂質土では i_c ≈ 0.8～1.2 程度が目安です。動水勾配がこの値を超えると地盤が不安定になります。
フローネットで求めた局所的な動水勾配と限界動水勾配を比較することで、パイピングや地盤の破壊リスクを事前に評価することができます。
数値解析による浸透流計算の実際
続いては、数値解析による浸透流計算の実際について確認していきます。
複雑な地盤条件や三次元的な浸透流問題には、有限要素法（FEM）や有限差分法（FDM）を用いた数値解析が有効です。
FEMによる浸透流解析の手順
有限要素法では、解析領域をメッシュ（要素）に分割し、各要素に土質定数（透水係数など）を割り当てます。
境界条件として、水頭が既知の境界（固定水頭境界）と流量が既知の境界（流量境界）を設定します。
異方性透水係数（水平方向と鉛直方向で透水係数が異なる場合）にも対応できることが、FEMの大きな強みです。
水平方向の透水係数（kh）と鉛直方向の透水係数（kv）の比が大きい場合（例：河川堤防の自然堆積土）は、等方性を仮定すると誤差が大きくなることがあります。
不飽和・飽和連成解析への展開
実際の地盤では、地下水位面より上は不飽和状態となります。
降雨浸透の解析や地すべりの安定評価では、不飽和・飽和連成解析が必要になるケースが増えています。
不飽和透水係数は飽和度（Sr）の関数として表され、飽和透水係数よりも非常に小さな値をとります。
リチャーズ方程式（Richards equation）が不飽和浸透流を記述する基本式として用いられます。
高度な解析ソフトウェアでは、飽和・不飽和の両領域を統一的に扱えるものが多く、実務での活用が広がっています。
解析結果の検証と感度分析
数値解析の結果は、常に現場計測データとの照合によって検証（キャリブレーション）することが重要です。
透水係数の不確かさは大きく、解析結果に直接影響するため、感度分析（透水係数を変化させたときの解析結果の変化を評価）を実施することが推奨されます。
感度分析によって、設計に最も影響する地層の透水係数が特定でき、調査・試験の優先順位付けにも役立ちます。
現場計測と数値解析を組み合わせたアプローチが、精度の高い浸透流評価につながるでしょう。
まとめ
本記事では、飽和透水係数の計算方法をダルシー則から始まり、動水勾配・透水量・透水断面積の計算、フローネット法、そして数値解析まで幅広く解説しました。
ダルシー則 Q = k × i × A という基本式の理解が、すべての計算の出発点です。
フローネット法は手法としての理解が設計思想の習得に役立ち、FEM解析は複雑な現実問題への対応力を高めます。
透水係数は測定精度に限界があることを念頭に置きつつ、感度分析や現場計測との照合によって解析の信頼性を高めていくことが実務での大切な姿勢です。
本記事が地盤工学・土質工学の理解促進に貢献できれば幸いです。
