Carbonate Pore Types
Carbonate pore systems exhibit distinct velocity-porosity relationships based on pore geometry:
- Interparticle: Pores between grains, follows Wyllie time average, best permeability
- Moldic: Dissolved grains leaving molds, intermediate velocity-porosity
- Vuggy (separate): Isolated vugs, high velocity for porosity, poor permeability
- Vuggy (touching): Connected vugs, intermediate properties
- Fractured: Very high velocity, cycle skipping, excellent permeability
Lucia Classification
Lucia's petrophysical classification divides carbonates into three classes based on pore size:
- Class 1: Large (>20 μm) interparticle or intercrystalline pores. Highest permeability for given porosity.
- Class 2: Small (<20 μm) interparticle pores or separate-vug porosity. Moderate permeability.
- Class 3: Microporosity (<2 μm) in mud-dominated fabrics. Lowest permeability.
Each class has distinct porosity-permeability relationships
Velocity-Porosity Transform
The acoustic response discriminates pore types:
Vp = √(K/ρ)
Where K is bulk modulus, ρ is bulk density
- Interparticle: Follows Wyllie time average: 1/Vp = φ/Vfluid + (1-φ)/Vmatrix
- Vuggy: Higher Vp (lower ΔT) than expected - vugs don't reduce stiffness proportionally
- Moldic: Between interparticle and vuggy trends
- Fractured: Even lower ΔT, often with cycle skipping
Acoustic Formation Index (AFI)
Quantifies deviation from time-average relationship:
AFI = (ΔTma - ΔT) / (ΔTma - ΔTf)
Where ΔTma = matrix transit time, ΔTf = fluid transit time
- AFI > 0.8: Interparticle porosity (follows time average closely)
- 0.5 < AFI < 0.8: Moldic or touching-vug porosity
- AFI < 0.5: Separate-vug or fractured (high velocity)
Density-Neutron Separation
Helps identify vuggy porosity:
- Large separation: Density porosity >> neutron porosity indicates separate vugs
- Close agreement: Interparticle or well-connected porosity
- Mechanism: Neutron sees hydrogen in pore fluid, density sees bulk void space
- Vugs: Density detects void, neutron may miss if water-filled vugs are large
Anselmetti & Eberli Approach
Velocity-porosity-mineralogy relationships in carbonates:
- Demonstrates systematic velocity variations with pore type at constant porosity
- Separate trends for grain-dominated vs mud-dominated fabrics
- Links acoustic properties to depositional texture and diagenesis
- Provides framework for seismic reservoir characterization
Permeability Estimation
Lucia permeability-porosity transforms by class:
log(k) = a·φ + b
Where constants vary by pore class:
- Class 1: k ≈ 45 × 1015φ (best permeability)
- Class 2: k ≈ 2 × 1015φ (moderate)
- Class 3: k ≈ 0.1 × 1015φ (poorest)
- Fractured: Permeability dominated by fracture network, often 100-1000× higher
Development Implications
- Interparticle: Predictable porosity-permeability, good sweep, waterflood candidate
- Moldic: Moderate heterogeneity, fair sweep efficiency
- Separate-vug: Poor connectivity, low recovery factor, may need fracturing
- Touching-vug: Complex flow paths, dual-porosity behavior
- Fractured: Excellent productivity, high heterogeneity, early water breakthrough risk
Application Guidelines
- Combine indicators: Use velocity, density-neutron, and image logs together
- Calibrate locally: Core-calibrate transforms for your field
- Watch for: Mixed pore types in single interval
- Validate: Compare with core thin sections and permeability data
- Consider: Diagenetic history affects velocity-porosity relationships