达西定律多孔介质流量计算器

使用达西定律计算流体在多孔介质中的流量、达西速度和渗流速度——地下水水文学和油藏工程必备。

输入渗透率、横截面积、压差、黏度、流动长度和孔隙度,即可计算体积流量和各类速度。

达西定律多孔介质流量计算器
使用达西定律计算流体在多孔介质中的流量、达西速度和渗流速度——地下水水文学和油藏工程必备。

关于达西定律计算器

达西定律是流体力学中最基础的方程之一,用来描述流体如何在多孔介质中流动。该定律由亨利·达西于 1856 年在法国第戎进行砂床水过滤实验后首次提出,建立了体积流量与施加压强梯度之间的线性关系。其表达式为 Q = kA·ΔP / (μL),其中 Q 为体积流量 (m³/s),k 为介质的本征渗透率 (m²),A 为垂直于流向的横截面积 (m²),ΔP 为驱动流动的压差 (Pa),μ 为流体的动力黏度 (Pa·s),L 为流动路径长度 (m)。 渗透率是达西定律中最重要的单一参数。它只是多孔介质自身的性质——与流体无关——反映介质根据孔隙结构、孔隙连通性和曲折度传输流体的能力。渗透率的数量级跨度很大:黏土约为 10⁻²⁰ 到 10⁻¹⁸ m²,细砂约为 10⁻¹⁶ 到 10⁻¹⁴ m²,粗砂和砾石约为 10⁻¹⁴ 到 10⁻¹⁰ m²,而高度裂隙化的岩石或多孔混凝土可达 10⁻¹⁰ m² 及以上。在石油工程中,渗透率通常以毫达西表示(1 mD = 9.869×10⁻¹⁶ m²)。 达西定律会导出两种速度。达西速度(或表观速度)为 v = Q/A,表示假设流体占据整个横截面(包括固体骨架)时的表观速度。渗流速度(或孔隙速度)是流体在孔隙空间中的实际平均速度:v_seepage = v/φ,其中 φ 为孔隙度。由于只有孔隙承担流动,渗流速度始终比达西速度高 1/φ 倍。对于孔隙度为 25% 的介质,流体在孔隙中的运动速度是达西速度所示的 4 倍。 达西定律是水文地质学(地下含水层流动和污染物迁移建模)、石油工程(油藏模拟和产量预测)、土壤科学(灌溉与排水设计)、化学工程(填充床反应器和膜过滤)以及土木工程(大坝渗流分析和基础排水)的基础。它适用于均质、各向同性、由牛顿流体饱和的多孔介质中的层流、稳态、不可压缩流动。在高流速下,当惯性效应变得显著时,Forchheimer 方程会加入一个二次速度项;在极小尺度下,滑移流(Knudsen 扩散)可能需要 Klinkenberg 修正。 本计算器使用压差的绝对值,因此无论符号约定如何,请输入样品两端的绝对压降。结果给出的流量和速度大小沿压强梯度方向。

算例

四个具有代表性的多孔介质流动场景,展示达西定律在不同工程应用中的表现。

场景结果说明
砂岩油藏:k=1×10⁻¹² m²,A=0.01 m²,ΔP=10⁶ Pa,μ=0.001 Pa·s,L=0.1 m,φ=0.25Q = 1×10⁻⁴ m³/s;v_darcy = 1×10⁻² m/s;v_seepage = 4×10⁻² m/s典型的石油油藏流动。高压差推动流体在该岩样中形成显著流动。
砂质土壤:k=1×10⁻¹⁰ m²,A=0.1 m²,ΔP=1000 Pa,μ=0.001 Pa·s,L=1.0 m,φ=0.35Q = 1×10⁻⁵ m³/s;v_darcy = 1×10⁻⁴ m/s;v_seepage ≈ 2.86×10⁻⁴ m/s砂质含水层中的地下水流动。较低的压强梯度带来缓慢但稳定的渗流。
工业陶瓷过滤器:k=1×10⁻¹⁴ m²,A=0.001 m²,ΔP=50,000 Pa,μ=0.001 Pa·s,L=0.05 m,φ=0.15Q = 1×10⁻⁸ m³/s;v_darcy = 1×10⁻⁵ m/s;v_seepage ≈ 6.67×10⁻⁵ m/s极致密的过滤介质需要较高压力才能获得可测流量。

如何使用达西定律计算器

  1. 输入介质的本征渗透率 k,单位为 m²。可查阅对应介质类型的公开表格,或使用实验测得的数值。若以毫达西表示,可用 1 mD = 9.869×10⁻¹⁶ m² 进行换算。
  2. 输入垂直于流向的横截面积 A,单位为 m²。对于圆柱岩心,A = π·r²;对于矩形板,A = 宽度 × 高度。
  3. 输入压差 ΔP,单位为帕——即驱动流动的入口到出口压降。请填写正值。
  4. 输入工况温度下流体的动力黏度 μ,单位为 Pa·s。20°C 时水的黏度为 0.001 Pa·s;油类黏度更高,且会随温度升高而降低。
  5. 输入流动长度 L(单位:m)以及孔隙度 φ(0 到 1 之间的小数,例如 0.30 表示 30% 孔隙度)。点击“计算”即可查看流量、达西速度和渗流速度。

常见问题

什么是渗透率?如何获取?
渗透率(k)是多孔介质的本征属性,描述流体通过它的难易程度。它只取决于孔隙结构,与流体本身无关。你可以通过对岩心样品进行实验室渗透率测试、查阅对应材料的公开数据表,或在已知几何条件下测得流量和压降后反推达西定律求得。
达西速度和渗流速度有什么区别?
达西速度(q = Q/A)是单位总横截面积上的体积流量——它把多孔介质视为一根实心管道。渗流速度是流体穿过连通孔隙的实际平均速度:v_seepage = q/φ。由于只有空隙部分(孔隙度)承担流动,因此它总是高于达西速度。
达西定律什么时候不适用?
在高流速(孔隙中的雷诺数较高)时,惯性力开始显著作用,达西定律会失效——通常当基于孔径的 Re > 1–10 时就是如此。对于这类情况,Forchheimer 方程会加入二次阻力项。它在极低压力下的气体流动(Klinkenberg 滑移流)以及高度非均质或裂隙化介质中也不适用,因为流道会绕过大部分基质。
如何将毫达西换算为 m²?
1 达西 = 9.869233×10⁻¹³ m²,因此 1 毫达西(mD)= 9.869233×10⁻¹⁶ m²。将 mD 数值乘以 9.869×10⁻¹⁶ 即可得到 m²。许多石油油藏的渗透率为 1–1000 mD,对应 10⁻¹⁵ 到 10⁻¹² m²。
温度会如何影响计算?
温度主要影响流体黏度。对于水而言,黏度会从 0°C 时的 0.00179 Pa·s 降至 100°C 时的 0.000283 Pa·s,减少约 6 倍。温度越高,黏度越低,在相同压强梯度下流量就越大。为了得到准确结果,请始终使用实际工况温度下的黏度。
什么是水力传导率?它与渗透率有什么关系?
水力传导率 K (m/s) 将渗透率与流体性质结合在一起:K = k·ρ·g/μ,其中 ρ 是流体密度,g 是重力加速度。它常用于地下水水文学,通常假设流体是已知温度下的水。渗透率 k (m²) 只是介质的纯物理性质;水力传导率 K 已经包含了流体因素。本计算器为了适用于所有流体,采用渗透率。