Petrophysical Rock Type
What Is a Petrophysical Rock Type?
Petrophysical rock type (also called hydraulic flow unit, reservoir rock type, or PRT) is a classification of reservoir rock intervals that groups core or log data into discrete populations with consistent and predictable relationships between porosity, permeability, capillary pressure, and fluid saturation. Unlike lithofacies, which are defined by depositional environment and grain composition, petrophysical rock types are defined by their flow behavior: two rocks with identical mineralogy but different pore geometry can belong to different rock types, while two rocks of different mineralogy but similar pore structure and flow capacity belong to the same rock type. PRTs are the fundamental building blocks used to assign permeability-porosity transforms, capillary pressure curves, and relative permeability curves to the cells of a reservoir simulation grid.
Key Takeaways
- Petrophysical rock types are defined by flow behavior, specifically the pore geometry that controls the porosity-permeability relationship, capillary pressure curve shape, and irreducible water saturation, not by mineralogy or depositional environment alone.
- The flow zone indicator (FZI) is the most widely used quantitative method to identify rock types from core data, calculated from the reservoir quality index (RQI) and normalized porosity using the Kozeny-Carman equation framework.
- Log-based rock typing using electrofacies clustering, NMR T2 distribution analysis, or self-organizing maps extends rock type assignments from cored wells into the uncored well population across a field.
- Each petrophysical rock type receives its own permeability-porosity transform, capillary pressure curve, and relative permeability curve, which are used to initialize fluid saturations and define flow behavior in reservoir simulation.
- Poor rock type definition is one of the most common reasons static reservoir models fail to match production history during simulation, because permeability and saturation heterogeneity are under-represented.
How Petrophysical Rock Types Are Defined
The theoretical basis for petrophysical rock typing rests on the Kozeny-Carman equation, which relates permeability to porosity through a pore geometry term that reflects pore throat size and tortuosity. Amaefule and Altunbay (1993) reformulated this equation into two dimensionless groups: the reservoir quality index (RQI = 0.0314 x sqrt(k/phi), where k is permeability in millidarcies and phi is porosity as a fraction) and the normalized porosity index (phi_z = phi / (1 - phi)). The ratio FZI = RQI / phi_z, the flow zone indicator, is a single number that characterizes the pore geometry of a rock sample. When FZI values are plotted on a log-log plot of RQI versus phi_z, samples from the same hydraulic unit fall on a straight line with unit slope. Distinct populations of samples that cluster around different parallel lines represent different hydraulic flow units, or petrophysical rock types, each with its own characteristic pore geometry and flow capacity.
In practice, a petrophysicist uses routine core analysis data (porosity and permeability measured on plug samples at ambient conditions) from one or more cored wells to calculate FZI for each sample. Statistical clustering methods (k-means clustering, hierarchical clustering, or probability distribution analysis) are applied to identify the natural breaks in the FZI distribution. The number of distinct rock types selected is a balance between capturing the real heterogeneity in the reservoir and keeping the model tractable: too few rock types over-smooth the permeability distribution and miss baffles; too many create an under-populated model where each type has insufficient data to define a reliable capillary pressure or relative permeability curve. Typical field studies define between 3 and 8 petrophysical rock types.
Special core analysis (SCAL) data are then grouped by rock type to build the input functions for simulation. Mercury injection capillary pressure (MICP) curves for samples within one rock type share a characteristic entry pressure and pore throat size distribution; these are averaged or representative curves are selected for use in fluid contact initialization. Similarly, unsteady-state or steady-state relative permeability measurements on plugs from each rock type define the saturation-dependent flow behavior for history matching.
- Abbreviations: PRT, HFU (hydraulic flow unit), RRT (reservoir rock type)
- Foundational paper: Amaefule and Altunbay, 1993, SPE-26436
- Key discriminator: Flow zone indicator (FZI), units of micrometers
- FZI formula: FZI = RQI / phi_z = 0.0314 x sqrt(k/phi) / (phi / (1 - phi))
- Typical number of rock types per field: 3 to 8
- Core data needed: Routine core analysis (porosity, permeability); SCAL (Pc, kr) for simulation
- Log-based extension methods: Electrofacies clustering, NMR T2, self-organizing maps (SOM)
- Primary simulation use: Permeability-porosity transform, capillary pressure, and relative permeability assignment per grid cell
When a new well is drilled without a core, assign petrophysical rock types from wireline logs using the electrofacies model calibrated in cored wells. Then cross-check the log-derived rock type assignment against the NMR T2 distribution if an NMR log was acquired: each rock type should have a characteristic T2 peak position reflecting its dominant pore throat size. If the two methods agree, you have high confidence in the assignment. Where they disagree, revisit the electrofacies model calibration because NMR responds to pore geometry directly and is more sensitive to rock type boundaries than resistivity or density-neutron crossplots alone.
Log-Based Rock Typing and Extension to Uncored Wells
Core data exist in only a fraction of the wells in a typical field. Extending rock type assignments to uncored wells requires identifying wireline log signatures that correlate reliably with core-derived FZI values. The classic approach is electrofacies analysis: log curves (gamma ray, density, neutron, photoelectric factor, deep resistivity) are input into a clustering algorithm that groups depth intervals with similar log responses. Each electrofacies cluster is then calibrated against core FZI values from the cored wells to assign a petrophysical rock type label. The calibrated model is applied to all uncored wells to produce a continuous rock type log.
NMR logging provides a more physically direct route to rock typing because the T2 relaxation time distribution measured by NMR reflects the pore size distribution of the rock. Smaller pores, with more surface-area-to-volume contact, relax faster (short T2), while larger pores relax slowly (long T2). A rock type with good pore connectivity and large pore throats shows a T2 distribution skewed toward longer times, while a tight, fine-pored rock type shows a distribution skewed toward short times. Clustering on the T2 distribution shape directly groups rock into hydraulic units with minimal reliance on the lithology-sensitive conventional logs. Self-organizing maps (SOM), a neural network approach, can simultaneously use conventional logs, NMR, and image log data to define rock type boundaries that would not be apparent from any single log in isolation.
Petrophysical Rock Type Synonyms and Related Terminology
- hydraulic flow unit (HFU) - the specific term used in the Amaefule-Altunbay framework, emphasizing the flow-capacity definition; often used interchangeably with petrophysical rock type in reservoir engineering literature
- reservoir rock type (RRT) - common industry abbreviation, used most frequently in carbonate reservoir studies where pore system complexity makes FZI alone insufficient
- electrofacies - a log-based rock classification defined purely by wireline log response clustering, without direct calibration to flow properties; must be calibrated to core to become a petrophysical rock type
- flow unit - broader term encompassing any subdivision of the reservoir based on flow capacity; can be defined at scales from pore (pore type) to field (reservoir zone)
Related terms: flow zone indicator, capillary pressure, relative permeability, reservoir simulation, electrofacies, NMR log
Frequently Asked Questions About Petrophysical Rock Types
What is the difference between a petrophysical rock type and a lithofacies?
A lithofacies is a classification based on sedimentological and mineralogical characteristics: grain size, sorting, cement type, depositional texture, and composition. A petrophysical rock type is a classification based on flow behavior: pore throat size distribution, permeability-porosity relationship, capillary entry pressure, and irreducible water saturation. The two systems often correlate but are not equivalent. A single lithofacies (for example, fine-grained bioturbated sandstone) can contain multiple petrophysical rock types if the degree of bioturbation, cementation, or diagenetic alteration varies significantly across the interval. Conversely, two different lithofacies deposited in different environments but cemented to the same pore geometry can fall into the same petrophysical rock type. For reservoir simulation, petrophysical rock types are the operationally meaningful classification because they directly define which permeability and saturation function to use.
How many core plugs are needed to define a reliable petrophysical rock type?
Industry practice generally requires a minimum of 20 to 30 routine core analysis plugs per rock type to define a statistically meaningful permeability-porosity transform, and a minimum of 3 to 5 SCAL measurements (capillary pressure or relative permeability) per rock type for simulation input. Rock types that are defined on fewer than 10 plugs should be flagged as poorly constrained in the model documentation, and the simulation team should perform sensitivity runs varying the properties of those rock types to understand the uncertainty. In fields with limited coring, this consideration can drive decisions to core additional wells specifically to fill data gaps in under-sampled rock types.
How are petrophysical rock types used to initialize fluid saturations in a reservoir model?
Each petrophysical rock type has an assigned capillary pressure curve that describes the relationship between water saturation and height above the free water level (where capillary pressure equals zero). In the static model, each grid cell is assigned a rock type based on its log-derived or interpolated electrofacies. The capillary pressure function for that rock type, combined with the cell's height above the computed free water level and the fluid density contrast, gives the initial water saturation for that cell. This process, called capillary pressure initialization or J-function initialization, is far more accurate than assigning a single average saturation-height function to the entire reservoir, because different rock types have different pore throat sizes and therefore different capillary entry pressures and transition zone thicknesses.
Why Petrophysical Rock Types Matter in Oil and Gas
Reservoir simulation models that skip rock typing and use a single porosity-permeability transform and capillary pressure curve for the entire reservoir routinely fail history matching because they cannot represent the flow heterogeneity that real reservoirs exhibit. Thin, high-permeability streaks within a lower-quality matrix, for example, can dominate water or gas breakthrough timing by orders of magnitude, and only a model that identifies these as a separate rock type with its own high-permeability transform will replicate the observed production behavior. The commercial consequence is significant: a simulation model that does not history-match well will generate unreliable production forecasts, leading to misplaced infill wells, incorrect waterflood pattern design, and poor EOR project economics. Investing in rigorous petrophysical rock typing during the static model build is one of the highest-return technical activities in field development planning.