多孔質媒体断面の概念図(代表体積要素 REV)と浸透流矢印
体積流量:$Q = q \cdot A$
水理伝導率:$K = \dfrac{k\,\rho g}{\mu}$
Re数:$Re_D = \dfrac{q\sqrt{k}}{\nu}$
フォルヒハイマー:
$\dfrac{\Delta P}{L}= \dfrac{\mu}{k}q + \beta\rho q^2$
$\beta \approx 1.75/(\sqrt{k}\,\phi^{1.5})$
多孔質媒体と流体を選択し、圧力勾配・断面積を設定するだけで、ダルシー流速・体積流量・水理伝導率・浸透 Reynolds 数をリアルタイム算出。3 媒体の流量比較チャートで地下水と石油工学のスケール感が直感的に掴めます。
多孔質媒体断面の概念図(代表体積要素 REV)と浸透流矢印
地下水・地盤工学:井戸からの揚水量の予測、汚染物質の地下拡散シミュレーション、ダムや堤防の基礎からの漏水評価に利用されます。地盤の透過率を正確に把握することで、適切な対策を講じることが可能です。
石油・ガス資源開発:貯留岩(砂岩や石灰岩)の中を原油や天然ガスがどのように流れるかをモデル化し、効率的な回収方法(二次・三次回収法)の計画立案に不可欠です。坑井配置の最適化にも用いられます。
化学プロセス・濾過:触媒充填層を通る反応流体の流れ、またはフィルターによる濾過プロセスの設計と解析に応用されます。圧力損失と処理流量のバランスを計算することで、装置のサイズと運転条件を決定します。
建材・断熱材の性能評価:コンクリートの透水性評価や、断熱材などの多孔質材料内での湿気(水蒸気)の移動解析に使用されます。耐久性や断熱性能の予測に貢献します。
このツールを使い始める際、特に現場から計算を任された新人さんが陥りがちなポイントがいくつかあるよ。まず「ダルシー流速 q は、実際の流体の平均流速ではない」ということ。ダルシー流速は、多孔質媒体全体の断面積で割った見かけの流速だ。実際の水は砂の粒の間の複雑な経路を縫うように流れるから、その実速度はもっと速い。例えば、空隙率が0.3(30%)の砂岩なら、実平均流速はダルシー流速を0.3で割った値になるんだ。
次に、透過率 k の単位とオーダーを間違えないこと。教科書では m² が基本だけど、現場では「ダルシー(D)」や「ミリダルシー(mD)」がよく使われる。1ダルシーは約 9.87×10⁻¹³ m² という非常に小さな値だ。良い砂岩で100-1000 mD、緻密な泥岩だと0.001 mD以下なんてこともある。スライダーで「1e-12」みたいな指数表示が出ても慌てないようにね。
最後に、「フォルヒハイマー補正の係数 β は k と独立ではない」という落とし穴。ツール上では独立に設定できるけど、実際には多孔質の構造(粒子の大きさや形状)によって k と β は関連している。実験データがない場合、例えば「β ≈ 1.8 / (k^0.5)」のような経験則で推定することもある。補正を入れると劇的に結果が変わるから、その係数がどこから来た値なのか、必ず裏付けを取ろう。
珪砂(均一砂)を用いた透水試験:K=5×10⁻⁴m/s、ΔP=2kPa、断面積A=50cm²(直径80mm円形供試体)、長さL=20cmの場合、ダルシー流速q≈0.005m/sとなり、体積流量Q=2.5×10⁻⁵m³/s≈1.5L/minです。同じ圧力勾配で粘土(K=1×10⁻⁷m/s)では流速は1/5000に低下し、Q≈5×10⁻⁹m³/sとなります。油層の貯留岩評価では、砂岩(K=1×10⁻⁶m/s)と泥岩(K=1×10⁻¹⁰m/s)で流動性が大きく異なります。