Secondary Porosity Index

The secondary porosity index (SPI) is a log-derived petrophysical parameter calculated as the difference between total porosity measured by the density-neutron combination (total porosity, reflecting all pore types including vugs and fractures) and the intergranular or matrix porosity measured by the sonic log (sonic porosity, sensitive primarily to the connected pore network that supports acoustic wave transmission), with the difference representing the fraction of total pore volume occupied by secondary pores — dissolution vugs, fractures, and moldic pores — that are large enough to be detected by density and neutron logs but are acoustically isolated from the primary pore network.

Key Takeaways

SPI = ΦDN − ΦSonic, where ΦDN is the density-neutron average porosity (reflecting total pore volume including secondary pores) and ΦSonic is the sonic-derived Wyllie time-average porosity (reflecting primarily the connected intergranular pore network), with positive SPI values indicating the presence of secondary porosity from dissolution, fracturing, or moldic pore development.
The physical basis for the SPI calculation is the different sensitivity of the sonic log versus the density and neutron logs to secondary pore types: sonic wave transit time is controlled by the primary intergranular pore network (the solid rock frame and its connected pores govern acoustic velocity), while density and neutron logs respond to total porosity including isolated vugs and fractures that don't contribute to acoustic wave travel time.
Secondary porosity is particularly important in carbonate reservoirs where diagenetic processes — dissolution during burial or subaerial exposure, stylolitization, dolomitization, and fracturing — create significant secondary pore volume that may provide excellent reservoir quality even when the original intergranular (intercrystalline in carbonates) porosity has been reduced by cementation.
The SPI is not a direct measure of fracture porosity or vug porosity individually — it represents the total secondary pore volume not sensed by the sonic log — and must be interpreted in the context of core observations, borehole image logs (which can distinguish fractures from vugs), and production test data to determine whether the secondary porosity is connected and contributes to permeability.
Limitations of the SPI include: the Wyllie time-average equation used to compute sonic porosity assumes fluid-saturated intergranular pores (gas in the secondary pores gives anomalously high transit times and may generate spurious positive SPI values unrelated to actual secondary porosity), the density and neutron logs are affected by hole size and mud invasion in fractured zones, and the SPI only detects secondary pores larger than the log tool resolution.

Fast Facts

The secondary porosity index concept was introduced in the early 1960s by Donald Pickett and subsequently refined by Schlumberger and other service companies as a quick-look carbonate reservoir quality indicator. In practice, SPI values above 0.02 (2 porosity units) in a carbonate zone are considered significant and warrant further investigation with borehole image logs or core to determine the type and connectivity of the secondary pore system. SPI values above 0.05 in a carbonate reservoir are typically associated with vuggy or fracture porosity that substantially enhances permeability and production rates relative to matrix-porosity-only zones. The SPI has a well-established place in carbonate petrophysics quick-look workflows for identifying zones with enhanced secondary porosity that may justify further evaluation even when matrix porosity appears modest.

What Is the Secondary Porosity Index?

Carbonate reservoirs — limestones and dolomites that host a large fraction of the world's oil and gas — commonly contain multiple coexisting pore types created by both original depositional processes and later diagenetic modification. Original intergranular or intercrystalline pores form during deposition or diagenetic recrystallization; secondary pores form afterward through dissolution of unstable minerals (creating molds, vugs, and caverns), fracturing from tectonic stress, or stylolitization from pressure solution. These different pore types contribute very differently to reservoir quality — connected secondary pores can provide the permeability pathways that make a low-matrix-permeability carbonate economically producible.

The challenge for log-based formation evaluation in carbonates is that conventional log interpretation methods developed for clastic (sandstone) reservoirs assume that all porosity is of similar type — intergranular. When secondary pores are present, the different logs respond differently to them, and this difference in response can be exploited to estimate secondary porosity volume without requiring a core sample.

The secondary porosity index exploits the fact that the sonic log, measuring how fast acoustic waves travel through the rock, is relatively insensitive to isolated large pores (vugs) and fractures because acoustic waves travel around them through the solid rock matrix rather than through them. The density and neutron logs, measuring electron density and hydrogen content respectively, respond to total pore volume including isolated vugs and fractures because they sample the bulk rock including the fluid in all pore types. The difference between total porosity (from density-neutron) and acoustic porosity (from sonic) therefore provides an estimate of the non-sonic-sensitive pore fraction — the secondary porosity.

SPI Calculation and Interpretation

The SPI calculation requires three log inputs: density log (RHOB, g/cm³), neutron porosity log (NPHI, decimal fraction, in limestone units), and sonic log (DT, microseconds per foot). The density-neutron porosity (ΦDN) is computed as the linear average: ΦDN = (ΦNPHI + ΦRHOB) / 2, where ΦRHOB = (ρmatrix − RHOB) / (ρmatrix − ρfluid). The sonic porosity (ΦSonic) uses the Wyllie time-average equation: ΦSonic = (DT − DTmatrix) / (DTfluid − DTmatrix), using appropriate matrix transit times (47.6 μs/ft for limestone, 43.5 μs/ft for dolomite, 55.5 μs/ft for sandstone) and fluid transit times (189 μs/ft for freshwater). Then SPI = ΦDN − ΦSonic.

Positive SPI values (ΦDN greater than ΦSonic) indicate secondary pores. The magnitude of SPI provides a quick-look estimate of secondary pore volume, but interpretation must account for the limitations: gas effect on the sonic log (gas increases DT, reducing ΦSonic artificially and creating positive SPI that is not secondary porosity), compaction effects on the sonic log (overpressured formations have faster sonic velocities, leading to negative ΦSonic errors in deep wells), and borehole size effects on density and neutron logs (washouts increase apparent density porosity and can inflate ΦDN above true value).

In practice, SPI is used as a qualitative indicator rather than a quantitative measurement of secondary porosity volume. Zones with elevated SPI relative to the local background are flagged for further investigation, and the SPI trend is compared to borehole image log fracture intensity, core vug count, and production test flow rates to calibrate the SPI response in the specific geological setting.

Secondary Porosity Index Across International Jurisdictions

Canada (AER / WCSB): SPI analysis is applied in WCSB carbonate reservoirs including Devonian Leduc, Nisku, and Swan Hills Formation reefs and platform carbonates, where vuggy and fracture porosity are important components of reservoir quality. AER well log data from WCSB carbonate pool wells includes the density, neutron, and sonic logs needed for SPI calculation, and SPI logs have been computed as standard outputs in WCSB carbonate formation evaluation programs since the 1960s. Alberta carbonate well evaluations submitted to the AER for pool delineation and reserve determination use SPI in conjunction with porosity-permeability cross-plots and core data to characterize secondary pore contribution to reservoir quality.

United States (API / SPE): SPI is a standard carbonate quick-look petrophysical tool in the Permian Basin (San Andres, Yates, Grayburg formations), the Anadarko Basin (Hunton, Arbuckle limestones), and the Rocky Mountain provinces (Madison Formation). SPE formation evaluation literature extensively documents SPI applications in US carbonate plays, including calibration of SPI against core thin section point counts and production test kh values. API RP 40 provides the core analysis framework (thin section petrography, mercury injection capillary pressure) that calibrates SPI-derived secondary porosity estimates in US carbonate well studies.

Norway (Sodir / NORSOK): Chalk reservoirs on the NCS (Ekofisk, Valhall) are characterized by very low matrix permeability (nanoDarcy to microDarcy range) but significant fracture permeability from natural fracture networks. SPI analysis in chalk wells can identify fracture-dominated intervals where the sonic-density separation is elevated, though chalk's unusual acoustic properties (highly compressible matrix) require special treatment of the Wyllie time-average equation. Equinor's petrophysical analysis of chalk wells uses SPI as one indicator of natural fracture intensity, calibrated against borehole image logs and production logging data.

Middle East (Saudi Aramco): Arab Formation carbonates in Saudi Arabia show complex secondary porosity development from both vuggy dissolution (in the Arab-D reservoir) and fracturing (in the tighter Arab-A, B, and C members). Saudi Aramco's formation evaluation program for Arab Formation wells routinely includes SPI calculations calibrated against extensive core thin section and SEM porosity data. The correlation between SPI and production performance in Arab Formation carbonates has been documented in Aramco SPE publications, confirming that elevated SPI zones show higher productivity indices and better well test performance than low-SPI zones with similar total porosity.

The secondary porosity index is also called the secondary porosity indicator or the sonic-density porosity separation. Related terms include vuggy porosity, fracture porosity, sonic log, density log, neutron log, carbonate reservoir, Wyllie time-average equation, and borehole image log. The total porosity measured by density-neutron is sometimes written as ΦT (total porosity) versus ΦP (primary or intergranular porosity from sonic), with SPI = ΦT − ΦP explicitly defining the secondary pore fraction.

Tip: Before trusting elevated SPI values as indicators of secondary porosity, cross-check with the caliper log to rule out borehole washout as the cause. In naturally fractured or vuggy carbonate zones, the borehole wall is often mechanically weak and the drilling fluid can erode the borehole, creating a larger-than-bit diameter hole. This washout affects the density and neutron logs (reading too much porosity due to the caved material or mud filling the washout) but does not proportionally affect the sonic log, creating an artificial positive SPI that reflects borehole deterioration rather than actual secondary pores. The rule of thumb: if the caliper shows a borehole diameter more than 10 to 15 percent larger than the bit size, treat the SPI value with caution and check whether a borehole size correction or pad tool measurement (microlog, ultrasonic) is available to verify the secondary porosity interpretation.

FAQ

Does a high SPI guarantee good permeability?
No — the SPI indicates the presence of secondary pore volume but says nothing about whether those secondary pores are connected to each other or to the primary pore network in a way that provides flow paths. Large isolated vugs detected by density and neutron logs but not acoustically connected (the sonic wave travels around them) contribute to total porosity without contributing to permeability. Fractures, by contrast, typically provide excellent permeability even at low fracture porosity because the fracture aperture is much larger than the matrix pore throat — a fracture with 10 micron aperture can have permeability thousands of times higher than the surrounding matrix. Combining SPI with borehole image logs (to identify fractures specifically), core measurements (to calibrate SPI against measured vug porosity and permeability), and production testing (which directly measures the effective kh) is necessary to determine whether elevated SPI translates to commercially significant permeability enhancement.