Sedimentkompaktion Porosität

Wenn Sedimente unter jüngeren Ablagerungen versenkt werden, schrumpft der Porenraum unter effektiver Spannung und Diagenese. Beckenmodelle, Reservoir-Volumetrie und Bohrlochlog-QC brauchen eine einfache, reproduzierbare Möglichkeit, die Porosität in der Tiefe aus einem Oberflächenwert und einem Kompaktionstrend zu schätzen. Das klassische exponentielle Porositäts–Tiefen-Gesetz nach Athy schreibt die Porosität als Oberflächenporosität mal e hoch minus Kompaktionskoeffizient mal Versenkungstiefe. Dieser Rechner für Porosität durch Sedimentkompaktion wendet dieses Gesetz deterministisch an: Sie geben φ0 als Bruchteil, die Versenkungstiefe z in einer konsistenten Längeneinheit und den Kompaktionskoeffizienten c in inverser Länge an; die Antwort liefert die Porosität in der Tiefe, den Prozentwert und das Restverhältnis φ/φ0. Ein negativer Kompaktionskoeffizient ist mit fortschreitender Kompaktion physikalisch unvereinbar und wird als ungültige Eingabe abgelehnt.

Was das exponentielle Porositäts–Tiefen-Gesetz bedeutet

Die exponentielle Porositäts–Tiefen-Beziehung im Athy-Stil ist das bewährte Modell erster Ordnung, das in der Sedimentgeologie gelehrt und in vielen Becken- und Petroleum-Systems-Workflows als Standardtrend genutzt wird. Die Oberflächenporosität φ0 ist der deposizionale oder oberflächennahe Hohlraumanteil der betrachteten Lithologie—typischerweise in der Größenordnung 0.4–0.5 für reine Sande und höher für manche Tonsteine vor signifikanter Kompaktion. Der Kompaktionskoeffizient c beschreibt, wie schnell diese Porosität mit der Tiefe abnimmt: größeres c bedeutet schnellere Reduktion des Porenraums für dasselbe Versenkungsintervall.

Die Versenkungstiefe z muss die in c implizierte Längeneinheit teilen (zum Beispiel Kilometer, wenn c pro Kilometer angegeben ist, oder Meter, wenn c pro Meter gilt). An der Oberfläche ist z null und die Exponentialfunktion eins, sodass die Porosität genau φ0 zurückgibt. Mit wachsender Tiefe fällt exp(−c·z) gegen null und die Porosität asymptotiert in dieser reinen Exponentialform gegen null.

Das Modell ist isotrop und eindimensional: es erfindet weder Überdruck, Zementationssprünge noch Diskordanz-Abschneidungen. Diese Prozesse zählen in realen Becken, liegen aber auf einem exponentiellen Basistrend; sie sind kein Grund, eine transparente Screening-Formel zu verwerfen. Dieser Rechner wendet nur die Basislinie an—φ(z) gleich φ0 mal exp(−c mal z)—damit jeder Client dieselbe Zahl aus denselben drei Skalaren reproduziert.

Wenn Sie lithologiespezifische c-Werte brauchen, entnehmen Sie sie Ihrer Bohrlochkalibrierung oder veröffentlichten Regionalkarten und speisen Sie sie hier ein; der Endpunkt enthält keine geheime Koeffiziententabelle. Oberflächenporosität wird nur als dimensionsloser Bruchteil zwischen null und eins einschließlich akzeptiert; wenn Ihre Logs Prozent melden, teilen Sie vor dem Aufruf durch hundert. Die Tiefe muss nichtnegativ sein.

Wann sich eine Porositätsabschätzung in der Tiefe lohnt

Explorations- und Reservoirteams brauchen schnelle Porositätsschätzungen in der Tiefe beim Screening von Prospects, beim Prüfen, ob eine log-abgeleitete Porosität für eine gegebene Versenkung plausibel ist, oder beim Aufbau einfacher eindimensionaler Kompaktionsprofile vor einem vollen Beckenmodell. Geomodellierungs-Pipelines und Jupyter-Notebooks können einen kleinen Parse-Endpunkt aufrufen, statt exp und Randfälle in jedem Skript neu zu implementieren, und so stille Einheitfehler und Bugs durch negatives c reduzieren. Behörden- und Beratungsberichte dokumentieren oft einen exponentiellen Trend neben einer Bohrung; die Automatisierung dieser Rechnung mit einer versionierten Formelzeichenkette hält Anhänge reproduzierbar.

Fehlgeschlagene Validierung wird nicht berechnet. Halten Sie φ0 als Bruchteil, pinnen Sie gerundete Ausgaben in Tests, und Sie erhalten deterministische Porositätszahlen in der Tiefe für Lehre, QC und leichte Automatisierung, ohne einen vollen geomechanischen oder Beckenmodellierungs-Stack aufzubauen.