Calculateur de la loi de Darcy pour milieu poreux
Calculez le débit des fluides, la vitesse de Darcy et la vitesse de filtration dans les matériaux poreux grâce à la loi de Darcy — indispensable en hydrogéologie et en ingénierie des réservoirs.
Saisissez la perméabilité, la section transversale, la différence de pression, la viscosité, la longueur d'écoulement et la porosité pour calculer le débit volumique et les vitesses.
Calculateur de la loi de Darcy pour milieu poreux
Calculez le débit des fluides, la vitesse de Darcy et la vitesse de filtration dans les matériaux poreux grâce à la loi de Darcy — indispensable en hydrogéologie et en ingénierie des réservoirs.
À propos du calculateur de la loi de Darcy
La loi de Darcy est l'une des équations les plus fondamentales de la mécanique des fluides, car elle décrit comment les fluides s'écoulent à travers les milieux poreux. Formulée pour la première fois par Henry Darcy en 1856 après des expériences de filtration de l'eau à travers des lits de sable à Dijon, en France, elle établit une relation linéaire entre le débit volumique et le gradient de pression appliqué. Elle s'écrit Q = kA·ΔP / (μL), où Q est le débit volumique (m³/s), k la perméabilité intrinsèque du milieu (m²), A la section transversale perpendiculaire à l'écoulement (m²), ΔP la différence de pression motrice (Pa), μ la viscosité dynamique du fluide (Pa·s) et L la longueur du trajet d'écoulement (m).
La perméabilité est le paramètre le plus important de la loi de Darcy. C'est une propriété propre au milieu poreux — indépendante du fluide — qui quantifie sa capacité à transmettre le fluide selon sa structure poreuse, la connectivité des pores et la tortuosité. La perméabilité couvre de nombreux ordres de grandeur : l'argile se situe entre 10⁻²⁰ et 10⁻¹⁸ m², le sable fin entre 10⁻¹⁶ et 10⁻¹⁴ m², le sable grossier et le gravier entre 10⁻¹⁴ et 10⁻¹⁰ m², et les roches très fracturées ou le béton poreux à partir de 10⁻¹⁰ m². En ingénierie pétrolière, la perméabilité est souvent exprimée en millidarcys (1 mD = 9.869×10⁻¹⁶ m²).
Deux vitesses émergent de la loi de Darcy. La vitesse de Darcy (ou vitesse superficielle) est v = Q/A, c'est-à-dire la vitesse apparente comme si le fluide occupait toute la section, matrice solide comprise. La vitesse de filtration (ou vitesse de pore) est la vitesse moyenne réelle du fluide à travers les vides connectés : v_seepage = v/φ, où φ est la porosité. Comme seuls les pores conduisent l'écoulement, la vitesse de filtration est toujours supérieure à la vitesse de Darcy d'un facteur 1/φ. Pour un milieu de porosité 25 %, le fluide se déplace dans les pores quatre fois plus vite que ne le suggère la vitesse de Darcy.
La loi de Darcy sous-tend l'hydrogéologie (modélisation des écoulements dans les aquifères et du transport des contaminants), l'ingénierie pétrolière (simulation des réservoirs et prévision de la production), la science du sol (conception de l'irrigation et du drainage), le génie chimique (réacteurs à lit fixe et filtration sur membrane) et le génie civil (analyse des suintements de barrages et drainage des fondations). Elle est valable pour un écoulement laminaire, stationnaire et incompressible dans un milieu poreux homogène et isotrope saturé par un fluide newtonien. À haut débit, lorsque les effets inertiels deviennent significatifs, l'équation de Forchheimer ajoute un terme quadratique de vitesse ; à très petite échelle, un écoulement avec glissement (diffusion de Knudsen) peut nécessiter la correction de Klinkenberg.
Le calculateur utilise la valeur absolue de la différence de pression ; saisissez donc la chute de pression absolue à travers l'échantillon, quel que soit le signe adopté. Les résultats donnent la magnitude du débit et des vitesses dans le sens du gradient de pression.
Exemples résolus
Quatre scénarios représentatifs d'écoulement en milieu poreux illustrant la loi de Darcy dans différentes applications d'ingénierie.
| Scénario | Résultat | Remarques |
|---|---|---|
| Réservoir de grès : k=1×10⁻¹² m², A=0.01 m², ΔP=10⁶ Pa, μ=0.001 Pa·s, L=0.1 m, φ=0.25 | Q = 1×10⁻⁴ m³/s; v_darcy = 1×10⁻² m/s; v_seepage = 4×10⁻² m/s | Écoulement typique d'un réservoir pétrolier. Un fort différentiel de pression entraîne un débit important dans cet échantillon de roche. |
| Sol sableux : k=1×10⁻¹⁰ m², A=0.1 m², ΔP=1000 Pa, μ=0.001 Pa·s, L=1.0 m, φ=0.35 | Q = 1×10⁻⁵ m³/s; v_darcy = 1×10⁻⁴ m/s; v_seepage ≈ 2.86×10⁻⁴ m/s | Écoulement d'eau souterraine dans un aquifère sableux. Un faible gradient de pression produit une infiltration lente mais régulière. |
| Filtre céramique industriel : k=1×10⁻¹⁴ m², A=0.001 m², ΔP=50,000 Pa, μ=0.001 Pa·s, L=0.05 m, φ=0.15 | Q = 1×10⁻⁸ m³/s; v_darcy = 1×10⁻⁵ m/s; v_seepage ≈ 6.67×10⁻⁵ m/s | Un milieu filtrant très serré nécessite une forte pression pour obtenir un débit mesurable. |
Comment utiliser le calculateur de la loi de Darcy
- Saisissez la perméabilité intrinsèque k en m². Consultez les tableaux publiés pour votre type de milieu ou utilisez des valeurs mesurées en laboratoire. Convertissez depuis les millidarcys avec 1 mD = 9.869×10⁻¹⁶ m².
- Saisissez la section transversale A en m² perpendiculaire à l'écoulement. Pour un carottage cylindrique, A = π·r² ; pour une plaque rectangulaire, A = largeur × hauteur.
- Saisissez la différence de pression ΔP en Pascals — la chute de pression de l'entrée à la sortie qui entraîne l'écoulement. Entrez une valeur positive.
- Saisissez la viscosité dynamique μ en Pa·s de votre fluide à la température de fonctionnement. L'eau à 20°C vaut 0.001 Pa·s ; la viscosité augmente pour les huiles et diminue avec la température.
- Saisissez la longueur d'écoulement L en mètres et la porosité φ sous forme décimale entre 0 et 1 (par ex. 0.30 pour 30 % de porosité). Cliquez sur Calculer pour voir le débit, la vitesse de Darcy et la vitesse de filtration.
Foire aux questions
Qu'est-ce que la perméabilité et comment la déterminer ?
La perméabilité (k) est une propriété intrinsèque du milieu poreux qui décrit la facilité avec laquelle un fluide peut le traverser. Elle dépend uniquement de la structure des pores, pas du fluide. Vous pouvez la déterminer à partir d'essais de perméamètre en laboratoire sur des carottes, de tableaux de données publiés pour votre matériau, ou en inversant la loi de Darcy si vous pouvez mesurer le débit et la chute de pression dans une géométrie connue.
Quelle est la différence entre vitesse de Darcy et vitesse de filtration ?
La vitesse de Darcy (q = Q/A) est le débit volumique par unité de section totale — elle traite le milieu poreux comme un tube plein. La vitesse de filtration est la vitesse moyenne réelle du fluide à travers les pores connectés : v_seepage = q/φ. Elle est toujours plus élevée que la vitesse de Darcy parce que seule la fraction vide (porosité) transporte l'écoulement.
Quand la loi de Darcy n'est-elle pas valable ?
La loi de Darcy cesse d'être valable à fort débit (nombre de Reynolds élevé dans les pores), lorsque les forces d'inertie deviennent importantes — en général lorsque Re > 1–10 selon la taille des pores. L'équation de Forchheimer ajoute alors un terme quadratique de traînée. Elle échoue aussi pour l'écoulement des gaz à très basse pression (glissement de Klinkenberg) et dans les milieux très hétérogènes ou fracturés où les canaux d'écoulement court-circuitent la majeure partie de la matrice.
Comment convertir la perméabilité des millidarcys en m² ?
1 darcy = 9.869233×10⁻¹³ m², donc 1 millidarcy (mD) = 9.869233×10⁻¹⁶ m². Multipliez votre perméabilité en mD par 9.869×10⁻¹⁶ pour obtenir des m². De nombreux réservoirs pétroliers ont des perméabilités de 1 à 1000 mD, soit 10⁻¹⁵ à 10⁻¹² m².
Comment la température affecte-t-elle le calcul ?
La température affecte principalement la viscosité du fluide. Pour l'eau, la viscosité passe de 0.00179 Pa·s à 0°C à 0.000283 Pa·s à 100°C, soit une réduction d'environ six fois. Une température plus élevée signifie une viscosité plus faible et donc un débit plus important pour le même gradient de pression. Utilisez toujours la viscosité à la température réelle de fonctionnement pour obtenir des résultats précis.
Qu'est-ce que la conductivité hydraulique et quel est son lien avec la perméabilité ?
La conductivité hydraulique K (m/s) combine la perméabilité avec les propriétés du fluide : K = k·ρ·g/μ, où ρ est la masse volumique du fluide et g l'accélération de la pesanteur. Elle est couramment utilisée en hydrogéologie, où le fluide est l'eau à une température connue. La perméabilité k (m²) est la propriété physique pure du milieu ; la conductivité hydraulique K intègre déjà le fluide. Ce calculateur utilise la perméabilité pour rester général pour tous les fluides.