参数设置
多孔介质类型
砾石（k=1×10⁻⁹ m²）
砂（k=1×10⁻¹² m²）
砂岩（k=1×10⁻¹⁴ m²）
粘土（k=1×10⁻¹⁷ m²）
流体种类
水（μ=1×10⁻³ Pa·s, ρ=1000）
油（μ=5×10⁻² Pa·s, ρ=850）
空气（μ=1.8×10⁻⁵ Pa·s, ρ=1.2）
计算结果
多孔介质截面概念图（代表性体积元 REV）与渗流矢量示意
达西定律：$q = \dfrac{k}{\mu}\dfrac{\Delta P}{L}$
体积流量：$Q = q \cdot A$
水力传导系数：$K = \dfrac{k\,\rho g}{\mu}$
渗流雷诺数：$Re_D = \dfrac{q\sqrt{k}}{\nu}$
Forchheimer方程：
$\dfrac{\Delta P}{L}= \dfrac{\mu}{k}q + \beta\rho q^2$
$\beta \approx 1.75/(\sqrt{k}\,\phi^{1.5})$
什么是达西定律渗流
🎓
简单来说，它就是描述水（或油、气）如何在沙子、岩石这类“多孔”材料里流动的规则。你可以想象成水在咖啡粉里渗透的过程。在实际工程中，比如我们想知道地下水能流多快，或者石油能从岩层里采出多少，都得靠它。你试着在模拟器里把“多孔介质类型”从“砂”换成“砾石”，看看“达西流速”会怎么变化？
🙋
诶，真的吗？我换成砾石后流速果然变大了！那这个“渗透率k”是啥？为什么它对流速影响这么大？
🎓
渗透率k就像是多孔介质的“通畅程度”。砾石的孔隙大、连通好，所以k值大，水就流得快；粘土的孔隙极小，k值就非常小，水几乎流不动。工程现场常见的是，为了预测油井产量，我们必须精确知道地下岩层的渗透率。你再用那个对数滑块调整一下k值，观察“体积流量”的变化，是不是非常敏感？
🙋
原来如此！那公式里还有个“渗流雷诺数”，它有什么用？旁边那个叫Forchheimer的方程又是干嘛的？
🎓
问得好！渗流雷诺数 $Re_D$ 是个“裁判”，它告诉我们达西定律这个简单规则什么时候还管用。当流速很低，$Re_D$ 小于1到10时，达西定律是准确的。但比如在油井附近，流体被快速抽吸，流速很高，$Re_D$ 变大，流动就变得“湍急”起来，惯性力不能忽略了。这时就需要请出Forchheimer方程，它在达西定律后面加了一个 $q^2$ 项来描述惯性效应。你可以在模拟器里把压力梯度调得特别高，看看流速和压降的关系是不是不再是一条直线了？
物理模型与关键公式
达西定律是描述多孔介质中低速渗流现象的本构方程，其核心是流速与压力梯度成正比。
$$q = -\dfrac{k}{\mu}\dfrac{\Delta P}{L}$$
其中，$q$ 是达西流速（m/s），$k$ 是固有渗透率（m²），$\mu$ 是流体动力粘度（Pa·s），$\Delta P/L$ 是沿流动方向的压力梯度（Pa/m）。负号表示流动方向指向压力降低的方向。
当流速较高，渗流雷诺数 $Re_D \gt 1\text{~}10$ 时，需使用Forchheimer方程来考虑惯性项的影响。
$$-\dfrac{\Delta P}{L}= \dfrac{\mu}{k}q + \beta \rho q^2$$
这里，$\beta$ 是惯性系数（1/m），$\rho$ 是流体密度（kg/m³）。方程右边第一项是粘性阻力（达西项），第二项是惯性阻力。$Re_D = \dfrac{q\sqrt{k}}{\nu}$，其中 $\nu$ 是运动粘度。
现实世界中的应用
地下水与环境保护： 用于模拟污染物在地下水中的迁移扩散，评估污染范围与修复方案。例如，预测化工厂泄漏的化学物质在多长时间内会到达下游的饮用水井。
石油与天然气开采： 是油藏数值模拟的核心基础。工程师用它来预测单井产量、优化布井方案以及设计注水驱油策略，直接影响油田的经济效益。
土木与地质工程： 分析大坝坝基、边坡和隧道围岩中的渗流稳定性，防止管涌、滑坡等灾害。计算基坑开挖时的降水井布置和抽水量。
材料与化工过程： 应用于过滤、催化反应器设计、燃料电池气体扩散层分析以及复合材料树脂浸渍成型（LCM）工艺的模拟，优化产品性能与生产流程。
常见误解与注意事项
开始使用此工具时，尤其是一些被现场委以计算任务的新人容易陷入几个误区。首先要明白“达西流速 q 并非流体的实际平均流速” 。达西流速是用多孔介质整体截面积折算的表观流速。实际水流是在砂粒间的复杂路径中穿行的，因此其实质流速要快得多。例如，对于孔隙度为0.3（30%）的砂岩，实际平均流速等于达西流速除以0.3。
其次，注意渗透率 k 的单位与数量级 。教科书通常以 m² 为单位，但现场常用“达西（D）”或“毫达西（mD）”。1达西约等于 9.87×10⁻¹³ m²，是一个非常小的数值。优质砂岩的渗透率约为100-1000 mD，致密泥岩可能低于0.001 mD。看到滑块显示“1e-12”这类指数形式时不必慌张。
最后，需警惕“福希海默修正系数 β 并非独立于 k” 这一陷阱。工具中虽可独立设置这两个参数，但实际上它们通过多孔结构（颗粒大小与形状）相互关联。若无实验数据，有时会采用“β ≈ 1.8 / (k^0.5)”等经验公式进行估算。引入修正项会显著改变计算结果，务必核实系数的来源与依据。