Porosité par compaction des sédiments

Lorsque les sédiments sont enfouis sous des dépôts plus jeunes, l'espace poral se réduit sous contrainte effective et diagénèse. Les modèles de bassin, la volumétrie de réservoir et le contrôle qualité des diagraphies ont besoin d'un moyen simple et reproductible d'estimer la porosité en profondeur à partir d'une valeur de surface et d'une tendance de compaction. La classique loi exponentielle de porosité–profondeur d'Athy exprime la porosité comme la porosité de surface multipliée par e à la puissance moins le coefficient de compaction fois la profondeur d'enfouissement. Ce calculateur de porosité par compaction des sédiments applique cette loi de façon déterministe : vous fournissez φ0 en fraction, la profondeur d'enfouissement z dans une unité de longueur cohérente et le coefficient de compaction c en inverse de longueur ; la réponse renvoie la porosité à cette profondeur, le pourcentage et le rapport résiduel φ/φ0. Un coefficient de compaction négatif est physiquement incompatible avec une compaction progressive et est rejeté comme entrée invalide.

Ce que signifie la loi exponentielle de porosité–profondeur

La relation exponentielle de porosité–profondeur de type Athy est le modèle de premier ordre de référence enseigné en géologie sédimentaire et utilisé comme tendance par défaut dans de nombreux workflows de bassins et de systèmes pétroliers. La porosité de surface φ0 est la fraction de vide dépositionnelle ou proche de la surface de la lithologie d'intérêt—typiquement de l'ordre de 0.4–0.5 pour des sables propres et plus élevée pour certaines roches argileuses avant une compaction significative.

Le coefficient de compaction c décrit la vitesse à laquelle cette porosité décroît avec la profondeur : un c plus grand signifie une réduction plus rapide de l'espace poral pour le même intervalle d'enfouissement. La profondeur d'enfouissement z doit partager l'unité de longueur implicite de c (par exemple des kilomètres lorsque c est exprimé par kilomètre, ou des mètres lorsque c est par mètre). En surface, z vaut zéro et l'exponentielle vaut un, de sorte que la porosité renvoie exactement φ0.

À mesure que la profondeur croît, exp(−c·z) tend vers zéro et la porosité asymptote vers zéro sous cette forme exponentielle pure. Le modèle est isotrope et unidimensionnel : il n'invente ni surpression, ni sauts de cimentation, ni troncatures d'inconformité. Ces processus comptent dans les bassins réels, mais se superposent à une tendance exponentielle de base ; ils ne justifient pas d'abandonner une formule de criblage transparente.

Ce calculateur n'applique que la ligne de base—φ(z) égale à φ0 fois exp(−c fois z)—afin que chaque client reproduise le même nombre à partir des mêmes trois scalaires. Lorsque vous avez besoin de valeurs de c spécifiques à la lithologie, prenez-les de votre calibration de puits ou de cartes régionales publiées et saisissez-les ici ; le point d'accès n'embarque pas de table secrète de coefficients.

La porosité de surface n'est acceptée que comme fraction sans dimension entre zéro et un inclus ; si vos diagraphies indiquent un pourcentage, divisez par cent avant d'appeler. La profondeur doit être non négative.

Quand estimer la porosité en profondeur

Les équipes d'exploration et de réservoir ont besoin d'estimations rapides de porosité en profondeur pour cribler des prospects, vérifier qu'une porosité dérivée des diagraphies est plausible pour un enfouissement donné, ou construire de simples profils unidimensionnels de compaction avant de lancer un modèle de bassin complet. Les pipelines de géomodélisation et les notebooks Jupyter peuvent appeler un petit point d'accès parse au lieu de réimplémenter exp et les cas limites dans chaque script, réduisant les erreurs silencieuses d'unités et les bogues de c négatif.

Les rapports publics et de conseil documentent souvent une tendance exponentielle à côté d'un puits ; automatiser ce calcul avec une chaîne de formule versionnée rend les annexes reproductibles. Une validation échouée ne consomme pas de solde. Gardez φ0 en fraction, figez les sorties arrondies dans les tests, et vous obtenez des nombres déterministes de porosité en profondeur pour l'enseignement, le QC et l'automatisation légère sans déployer une pile complète de géomécanique ou de modélisation de bassin.