Geostatic Pressure: Overburden Gradient, Lithostatic Load, and WCSB Frac and Pore-Pressure Design

Geostatic pressure is the vertical pressure exerted at a given depth by the combined weight of the overlying rock column and the fluids contained in that rock, and it is used interchangeably with the term lithostatic pressure and closely related to overburden pressure. It is the total downward load that the earth's mass above a formation presses onto that formation, and it sets the upper bound against which pore pressure, fracture pressure, and effective stress are all measured. The value is calculated by integrating rock bulk density over depth: geostatic pressure equals the integral of bulk density times gravitational acceleration times depth, which in practice reduces to a gradient multiplied by true vertical depth. For a broad average of sedimentary rock with a bulk density near 2.3 g/cm3, the overburden gradient works out to roughly 1.0 psi/ft, or in metric terms about 22.6 kPa/m, and for quick estimates the industry rounds the overburden gradient up to 1 psi/ft. That average conceals real variation, because geostatic gradient changes with depth and lithology as sediments compact and porosity falls. Shallow, poorly consolidated, high-porosity sediments carry more pore fluid and less rock mass per unit volume, so their gradient is lower, while deep, compacted carbonates and tight sandstones approach and exceed the 1 psi/ft benchmark. Published Gulf Coast examples show a geostatic gradient of 0.85 psi/ft at 1,000 ft rising to 0.95 psi/ft at 14,000 ft, illustrating that the number is not a single constant but a depth-dependent curve that must be built from a density log or a regional compaction trend. In the Western Canadian Sedimentary Basin the overburden above a Montney or Duvernay target is dominated by Cretaceous and Triassic clastics and Devonian carbonates, and operators construct a basin-specific overburden curve from bulk-density logs run on offset wells rather than defaulting to 1 psi/ft. Geostatic pressure matters because of the fundamental effective-stress relationship, where the total geostatic (overburden) stress equals the effective grain-to-grain stress plus the pore-fluid pressure. That single equation governs abnormal pressure detection, since it is very unlikely for formation pressure to exceed the overburden gradient, so the geostatic curve defines the ceiling that a pore-pressure prediction cannot physically cross. It also anchors fracture gradient estimation through poroelastic models such as Eaton's equation, where minimum horizontal stress and therefore the pressure required to part the rock during a hydraulic frac are derived from overburden stress, pore pressure, and a Poisson's ratio term. Reservoir compaction, subsidence over depleting pools, casing-design load cases, and wellbore-stability mud-weight windows all trace back to the geostatic load, making it one of the most foundational quantities in geomechanics and drilling engineering across the WCSB.

Key Takeaways

  • Total weight of the overburden: Geostatic pressure is the vertical stress from the full column of rock plus contained fluids above a depth, identical in meaning to lithostatic pressure. It is computed by integrating bulk density with depth, and it represents the maximum load a formation carries, which is why pore pressure predictions are capped by the overburden curve and any prediction exceeding it is treated as a data or model error.
  • The 1 psi/ft rule of thumb: An average sedimentary bulk density near 2.3 g/cm3 gives an overburden gradient of about 1.0 psi/ft (roughly 22.6 kPa/m), the number engineers use for rapid checks. Real gradients run lower in shallow uncompacted section (near 0.85 psi/ft) and rise toward and past 1 psi/ft in deep compacted carbonates, so a proper curve comes from density logs, not the single constant.
  • Effective stress ties it together: The core relationship is overburden stress equals effective stress plus pore pressure. As a WCSB pool depletes and pore pressure falls, effective grain-to-grain stress rises under a nearly fixed overburden, driving reservoir compaction and, in extreme cases, surface subsidence, a concern tracked under AER reservoir-management and measurement requirements for large pools.
  • Sets the fracture ceiling: Fracture-gradient models (Eaton, Matthews-Kelly) derive minimum horizontal stress from overburden stress, pore pressure, and Poisson's ratio, so an accurate geostatic curve is a prerequisite for designing a Montney or Duvernay hydraulic frac. Underestimating overburden underpredicts breakdown pressure and treating pressure, risking screenouts or off-target frac geometry.
  • Depth and lithology dependent: Because compaction expels pore fluid and increases bulk density with burial, the geostatic gradient is a curve, not a point. WCSB operators build overburden profiles from offset bulk-density logs spanning the Cretaceous clastic overburden and Devonian carbonate section, since a single assumed gradient can misstate stress by hundreds of kPa over a 3,000 m column.

Building a WCSB Overburden Curve from Density Logs

A geostatic curve for a Montney well near Dawson Creek is assembled by integrating the bulk-density log from surface to total depth. Where near-surface log data is missing, the top interval is filled with an assumed shallow density of about 2.0 g/cm3, then the compensated density log carries the integration through the Cretaceous shales and sands into the Triassic Montney at roughly 2,500 m. At that depth the accumulated overburden stress typically reaches 55,000 to 60,000 kPa (about 8,000 to 8,700 psi), giving an average gradient close to 22.5 kPa/m. That curve then feeds the minimum-horizontal-stress model, so any density-log error near surface propagates through the entire stress profile used for casing and frac design.

Overburden Stress and Reservoir Compaction

In a depleting WCSB pool the overburden stress stays essentially constant while pore pressure declines with production, so effective stress on the rock frame increases and the reservoir compacts. In high-porosity, weakly cemented formations such as some Clearwater or oil-sands intervals, this compaction can transmit to surface as subsidence and can shear casing in extreme cases. Operators monitor it with repeat formation pressure surveys and, for the largest pools, geomechanical modeling that couples the fixed geostatic load to falling pore pressure, informing AER Directive 051 injection and pressure-maintenance decisions that keep effective stress within safe casing and caprock limits.

Fast Facts

The near-universal 1 psi/ft overburden benchmark is really a statement about rock density: it corresponds almost exactly to an average bulk density of 2.31 g/cm3, close to the density of quartz sandstone at moderate porosity. The coincidence that sedimentary basins worldwide average near this value is why a number first popularized in 1950s Gulf Coast drilling still travels as a global rule of thumb, even though a deep WCSB carbonate at 2.7 g/cm3 actually loads the formation at closer to 1.17 psi/ft.

Geostatic pressure is one term in the effective-stress balance, so it connects directly to Pore Pressure, the fluid pressure that, subtracted from overburden, leaves the effective stress on the rock frame. It sets the ceiling for Formation Pressure because pore pressure cannot physically exceed the overburden load. It is the primary input to Fracture Gradient models that predict breakdown pressure, and it is the foundational stress in Overburden analysis used throughout drilling and completions design.

Real-World WCSB Scenario: Underestimated Overburden on a Duvernay Pad Near Kaybob

A Duvernay operator near Kaybob planned a multi-well frac program using a generic 1 psi/ft overburden assumption instead of a log-derived curve. Because the thick Cretaceous shale overburden in that area carries a bulk density slightly below the assumed value in its upper section but higher Devonian carbonate loading below, the simplified curve underpredicted overburden stress at the 3,400 m target by roughly 4,000 kPa. The Eaton fracture-gradient model fed by that low overburden underestimated closure stress, so the treating pressure schedule and proppant ramp were set too low, and the first two stages screened out. Remedial coiled-tubing cleanout and lost rig time added about CAD 620,000 to the program.

The operator rebuilt the overburden profile from three offset bulk-density logs, corrected the minimum-horizontal-stress model, and re-tuned the treatment schedule. The remaining stages placed on design with no further screenouts, and the corrected geomechanical model became the pad-development standard, avoiding an estimated CAD 300,000 per well in future intervention costs.