Rock Types: Petrophysical Rock Typing, Flow Zone Indicators, and Permeability-Porosity Mapping in WCSB Reservoirs

In reservoir characterization, rock types are groupings of rock that share a common set of petrophysical characteristics, chiefly those governing how fluids move through the rock and how much fluid the rock can store. The two properties of greatest interest are permeability, which controls flow capacity, and porosity, which controls storage capacity. A petrophysical rock type, often abbreviated PRT, bundles rocks that exhibit a similar relationship between these two properties along with a similar pore-throat geometry, so that a single set of saturation-height and relative-permeability functions can represent the whole group. This matters because two samples can have identical porosity yet differ in permeability by two or three orders of magnitude depending on whether the pore space is connected through wide, well-sorted throats or choked by clay, cement, and microporosity. Geologists and petrophysicists separate rock types using a blend of core data and log response: routine and special core analysis supplies porosity, permeability, grain density, and mercury-injection capillary pressure, while wireline logs supply continuous porosity, resistivity, and lithology indicators that let the discrete core-defined types be propagated across uncored intervals and into the full three-dimensional reservoir model. Several quantitative schemes formalize the grouping. The Winland r35 method computes the pore-throat radius at 35 percent mercury saturation and bins rocks into petrophysical units called port types. The flow zone indicator, or FZI, derived from the Kozeny-Carman relationship, collapses porosity and permeability onto a single hydraulic-unit parameter so that rocks falling on one straight line in a normalized-porosity plot belong to one hydraulic flow unit. The Lucia classification, built for carbonates, sorts rock by particle size and interparticle porosity. In the Western Canadian Sedimentary Basin these methods are applied daily: the siltstone-dominated Montney is subdivided into rock types by the proportion of dolomite cement and the balance of intergranular versus intragranular porosity, the Cardium is split between clean conglomerate-rich pay and tight muddy facies, and Devonian carbonate plays such as Leduc and Nisku are typed by depositional fabric and diagenetic overprint. Permeability is reported in millidarcies, porosity as a decimal fraction or percent, and pore-throat radius in microns, with the same data often expressed in both metric and field units across a basin that straddles Alberta, British Columbia, and Saskatchewan. Accurate rock typing underpins net-pay cutoffs, volumetric reserves bookings under AER and BCER reporting standards, and the placement of horizontal wells and hydraulic fracture stages within the most permeable, most connected intervals.

Key Takeaways

  • Flow and storage define the class: A rock type groups samples by shared permeability-porosity behaviour and pore-throat geometry, not just lithology. Two rocks at 8 percent porosity can differ from 0.01 mD to 10 mD depending on throat connectivity, which is why a Montney siltstone at one porosity can be pay in one facies and non-reservoir in another.
  • Winland r35 and FZI quantify it: The Winland equation estimates the pore-throat radius at 35 percent mercury saturation to assign port types, while the flow zone indicator from the Kozeny-Carman relation collapses porosity and permeability onto a single hydraulic-unit value. Rocks on one FZI trend share one set of saturation-height and relative-permeability curves.
  • Core anchors, logs propagate: Discrete rock types are first cut on core porosity, permeability, and mercury-injection capillary pressure, then tied to wireline porosity, resistivity, and lithology logs so the classification can be carried through hundreds of metres of uncored WCSB section and into the full geocellular model.
  • WCSB plays are routinely typed: The Montney is divided by dolomite cement and pore class, the Cardium between conglomeratic pay and tight mud facies, and Leduc and Nisku carbonates by Lucia depositional and diagenetic fabric. Each rock type carries its own permeability distribution that drives completion design.
  • Reserves and well placement depend on it: Net-pay cutoffs, volumetric estimates reported to the AER and BCER, and the geosteering of horizontal laterals all rely on correct rock typing. Misassigning a tight facies as pay inflates booked reserves; missing a connected facies leaves recoverable barrels behind pipe.

Building Rock Types from Mercury-Injection Capillary Pressure

Mercury-injection capillary pressure, or MICP, is the workhorse measurement for rock typing because it directly probes pore-throat size distribution. Mercury is forced into a cleaned core plug at rising pressure; each pressure step corresponds to the throat radius being entered, so the saturation-versus-pressure curve is effectively a map of how the pore network is plumbed. A Montney sample with a sharp entry pressure near 1,000 psia and a steep plateau indicates well-sorted throats around 0.1 micron, a distinct rock type from a sample with a gradual curve and dominant microporosity. Petrophysicists fit the Thomeer or Brooks-Corey parameters to each curve, cluster the results, and assign each cluster a rock type with its own saturation-height function for the static model.

Carrying Rock Types into the Dynamic Reservoir Model

Once defined, rock types must reach the flow simulator, because each one needs its own relative-permeability and capillary-pressure table. In a Cardium waterflood near Pembina, assigning a single averaged relative-permeability curve to conglomerate and muddy facies alike would smear the water front and mispredict breakthrough by months. Modellers instead populate the geocellular grid with rock-type codes, usually through a porosity-permeability transform and indicator geostatistics, then attach the matching SCAL curves. The result feeds material-balance checks, history matching, and the recovery-factor estimates that justify infill drilling. Getting the rock-type proportions right is often the single largest control on forecast ultimate recovery in tight WCSB plays.

Fast Facts

The Kozeny-Carman equation that underlies the flow zone indicator dates to work by Josef Kozeny in 1927 and Philip Carman in 1937, decades before reservoir engineering existed as a discipline. Amaefule and co-authors repackaged it as the FZI hydraulic-unit method in a 1993 SPE paper, and it remains one of the most cited rock-typing techniques in the industry. A single tight-gas rock type in the Montney can span four orders of magnitude in permeability while holding porosity within a two-percent band, which is precisely why porosity alone never defines a rock type.

Rock typing sits at the intersection of several core concepts. It is built directly on porosity and permeability, the two properties that define storage and flow capacity for every group. Capillary pressure supplies the pore-throat geometry through mercury-injection curves that distinguish one type from another and set the saturation-height function. The grouped types ultimately feed the reservoir model, where each carries its own relative-permeability behaviour and controls how recoverable hydrocarbons are calculated and produced.

Rock Typing a Montney Horizontal Program near Dawson Creek

An operator planning a multi-well Montney pad near Dawson Creek in northeast British Columbia cut 36 metres of core across the upper and middle Montney, ran MICP and routine core analysis on 120 plugs, and resolved four petrophysical rock types ranging from a dolomite-cemented silt with 0.08 mD permeability down to a clay-rich silt below 0.001 mD. Calibrating the porosity and resistivity logs to these core types, the team mapped a continuous 14 metre interval of the two best rock types and landed the 2,800 metre lateral within it, spacing fracture stages at 60 metres through the highest-FZI rock. Core and analysis carried roughly CAD 450,000 of the well budget.

The rock-type-guided landing delivered a thirty-day initial gas rate about 35 percent above the offset wells that had been geosteered on gamma ray alone, and the BCER reserves filing booked a tighter, better-supported recovery factor. The operator carried the same four-type scheme into the next eight wells on the pad, eliminating further coring cost while preserving the completion design that the rock typing had justified.