|Development Geology Reference Manual
|Methods in Exploration
|Carbonate reservoir models: facies, diagenesis, and flow characterization
|F. Jerry Lucia
The flow characteristics of carbonate reservoirs are controlled by a combination of depositional and diagenetic processes. Depositional processes control the initial pore size distribution and the geometry of the individual depositional facies. The diagenetic overprint modifies the pore size distribution and controls the productivity of depositional facies. In some cases, reservoir quality and flow characteristics are totally controlled by diagenesis, as in karsted reservoirs.
Carbonate reservoir descriptions are based on observations of depositional and diagenetic fabrics and pore space from core and cuttings samples. The descriptions are correlated with wireline log responses and incorporated into geological facies and/or diagenetic models in order to map porosity, water saturation, and permeability.
Carbonate sediments and environments
The majority of carbonate sediments are produced in shallow, warm ocean waters by extraction of calcium carbonate from seawater by organisms to form their shells or skeletal material. The sediments are composed of a spectrum of sizes and pore geometries. The Dunham classification (Figure 1) describes depositional textures in a manner that can be related to pore geometries. The grain-supported textures tend to have larger pore sizes than do mud-supported textures. The textures have different geometries in different depositional environments.
There are five basic carbonate depositional environments. From shore to basin, they are peritidal (tidal flat), shallow shelf interior, shelf margin complex, slope, and basin (Figure 2). (For more information on carbonate depositional environments, see Scholle et al.)
The peritidal depositional environment is complex (Figure 2). Sediments deposited between mean high and mean low tide are called intertidal sediments, sediments deposited above mean high tide are called supratidal sediments, and sediments deposited below mean low tide are called subtidal sediments. In arid and semi-arid climates, evaporite flats (sabkhas) are present from which gypsum and halite are deposited. Eolian sand dunes composed of siliciclastic or carbonate grains may form on the supratidal surface.
The shallow shelf interior environment (Figure 2) is dominated by low-energy waters that allow lime mud to accumulate. Storms, however, churn the sediment into suspension, winnowing out the fine-sized material and concentrating the coarse material. Near shorelines, the shelf environment may be composed of offshore bars and spits oriented parallel to shoreline. Shorelines that face heavy wave action accumulate carbonate sand or gravel. Tidal currents are concentrated in channels between islands and produce tidal deltas on the lee side of the island.
The shelf margin complex (Figure 2) is characterized by the presence of carbonate sands and reefs. Reefs are commonly found at the shelf edge where their rigid framework can withstand strong wave action and they can take advantage of the nutrients upwelling from the deeper waters. Carbonate sands derived from a reef or from plants and animals inhabiting the shelf edge accumulate along a wide belt that follows the break between the shelf edge and the slope. Tidal and storm currents mold the sand belt into tidal channels and bars.
The slope (Figure 2) is dominated by sediment transport seaward from the shelf margin. Fine-grained sediment settles to the bottom forming thin-bedded mudstones, while slumps, debris flows, and turbidity currents form coarse-grained bodies of breccia, conglomerate, and carbonate sand. The resulting facies patterns depend upon the relief of the shelf margin and the nature of the shallow water portion of the margin.
Sediments of the basin environment (Figure 2) are dominated by very fine-grained skeletons of planktonic microorganisms, which when lithified, become chalks. They are found in deep basins as well as on drowned shelves.
The reservoir quality of sediments is modified by diagenetic processes. Some diagenetic processes are directly related to the physical nature of the sediment, and the resulting geometries can be directly related to sedimentation patterns. Other diagenetic processes involve the interaction of sediment with flowing interstitial water. Understanding hydrodynamics is very important in predicting the geometry of the diagenetic effects.
Calcium carbonate cementation and compaction
Calcium carbonate cementation and compaction are diagenetic processes that are initiated soon after deposition and by which intergranular pore space is progressively reduced producing systematic changes in petrophysical properties (see Evaluating diagenetically complex reservoirs and Reservoir quality). Compaction is a physical/chemical process, while cementation requires fluid flow. Microporosity within grains or between lime mud particles may be retained even though intergranular pore space is lost.
Dolomitization is a diagenetic process that converts limestones to dolostones through a microchemical process of calcium carbonate dissolution and dolomite precipitation. Dolomitization can change the rock fabric and the petrophysical properties significantly because the dolomite crystals are commonly larger than the replaced limestone particles. Dolomite cement systematically grows on the dolomite crystal faces, reducing reservoir quality. Dolomitization requires the addition of large quantities of magnesium through fluid flow.
The most common evaporite mineral found with carbonate rocks is anhydrite and its hydrous form, gypsum. Gypsum is the common form at shallow depths, but it converts to anhydrite at depth in response to higher temperatures. Bedded anhydrite is commonly found in tidal flat environments and is an effective reservoir seal. Diagenetic anhydrite is found in reservoir rocks as nodules and poikilotopic crystals, which have little effect on reservoir properties, and as pore-filling crystals, which reduce reservoir quality.
Dissolution and associated processes
Dissolution is the diagenetic process by which carbonate and evaporite minerals are dissolved and removed, thus creating and modifying pore space in reservoir rocks (see Reservoir quality). The effect of this process on permeability depends upon the geometry and location of the resulting voids relative to the rock fabric. In some cases, dissolution is fabric selective and results in formation of isolated vugs. In other cases, dissolution enlarges fractures and interparticle pores resulting in large, connected vugs. If the vugs are large enough, the roof may collapse, forming a floor breccia and fractured roof.
Carbonate rock fabric and petrophysical relationships
Petrophysical classification of carbonate pore space
The Archie classification was the first attempt at relating rock fabrics to petrophysical rock properties. Lucia improved on the Archie classification by defining the pore size distribution in relationship to the particle size and the spatial relationship of the pore space to the particles.
Interparticle porosity is defined as that pore space located between grains (intergranular porosity) or between crystals (intercrystalline porosity), but which is not significantly larger than the grains or crystals. Vuggy porosity is defined as pore space larger than or within rock particles. Isolated vugs are called separate vugs, while vugs that form a connected pore system on a reservoir scale are called touching vugs.
Rock fabric and petrophysical property relationships
The rock fabric elements important to petrophysical properties are particle size, interparticle porosity, separate vug porosity, and the presence or absence of touching vugs. These elements can be identified and quantified by microscopic examination of core material and related to petrophysical properties measured on the same material.
In nonvuggy carbonates, permeability and capillary properties are a function of interparticle porosity and particle size. The addition of separate vug porosity to a carbonate reservoir rock increases storage capacity, but has little effect on permeability. Separate vug porosity, however, has a large effect on the Archie m factor used in water saturation calculations from wireline logs and on the acoustic transit time. Touching vug pore systems dominate fluid flow and override the effects of interparticle porosity, particle size, and separate vug porosity.
Relationships between rock fabric and wireline log response
Since most fields have few cores available, wireline logs must be used to identify rock fabric elements. This requires the calibration of wireline log responses with core data.
Particle size can be determined from gamma ray, porosity, and resistivity logs. Grain-supported rocks commonly have lower gamma ray activity then do mud-supported rocks. Other fine-grained rocks, such as shaley and organic-rich carbonates, commonly have the highest gamma ray activity. However, the level of gamma ray activity in some carbonates (dolostones in particular) is not related to particle size because of the presence of anomalous concentrations of uranium. Water saturation is a function of particle size and interparticle porosity, and crossplots of porosity, water saturation, and reservoir height can be used to determine particle size.
Interparticle porosity can be determined by subtracting separate vug porosity from total porosity. Total porosity can be calculated from porosity logs, while separate vug porosity can be estimated from crossplots of acoustic transit time versus crossplot porosity. Touching vug pore systems can be identified using borehole televiewer and resistivity scanner logs. (For more details, see Borehole imaging devices.)
Lithology can have a major effect on porosity logs. Gypsum is of particular concern in carbonate reservoirs because it contains large quantities of bound water. The bound water is seen by neutron logs as porosity, resulting in erroneously high porosity indications.
Carbonate reservoir models
The key to describing the geometry and petrophysical properties of carbonate reservoirs lies in combining sedimentary and diagenetic descriptions, geological models, rock fabric analysis, and petrophysical rock properties. While each hydrocarbon accumulation has a unique geometry, certain combinations of sedimentation, diagenesis, and pore fabric seem to result in a similar distribution of petrophysical rock properties. Four such reservoir models are (1) the upward-shoaling cementation and compaction model, (2) the subtidal-supratidal dolomitization and sulfate emplacement model, (3) the karst-collapse model, and (4) the geological reef model.
Upward-shoaling cementation and compaction reservoir model
The upward-shoaling model is based on a depositional model of sediment aggrading to sea level. As the water shallows, the energy conditions increase, resulting in a vertically stacked sequence from mudstones and wackestones at the bottom to packstones and grainstones at the top (Figure 3).
Cementation and compaction occur with burial, reducing the reservoir quality of the mud-supported mudstones and wackestones more than the grain-supported packestones and grainstones. The result is a vertical sequence of lower porosity and permeability mud-supported rocks at the base and higher porosity and permeability grain-supported rocks at the top. Consequently, the best quality flow unit occurs at the top of the sequence. The permeable units are confined to the grainstone bars.
If the grainstone bars are exposed to meteoric diagenesis, significant separate vug porosity can develop, causing a loss of permeability while retaining porosity. In the dolomite environment, grainstone bars are commonly cemented with anhydrite, and intercrystalline porosity in mud-supported sediments forms the permeable facies.
Subtidal-supratidal dolomitization and sulfate emplacement reservoir model
The subtidal-supratidal model (Figure 3) is based on the transport of carbonate sediment onto the shore by storm and tidal currents resulting in the progradation of the tidal flat environment over the subtidal environment. Subtidal intervals are commonly composed of mudstones, wackestones, packstones, and grainstones in no predictable order. When present, grain-supported sediments may be concentrated in the upper part of the subtidal section in the form of offshore bars and high-energy shoreface deposits. Intertidal and supratidal sediments are typically muddy except in association with high-energy subtidal sediments. A typical vertical sequence would show intercalated mud- and grain-supported sediments in the subtidal interval overlain by algal mats in the intertidal interval, and mud-cracked and desiccated wackestones and mudstones in the supratidal interval.
The subtidal-supratidal sequence is commonly dolomitized and contains anhydrite and gypsum. In the subtidal interval, dolomitized grainstones retain their intergranular pore space, except where cemented by anhydrite, and form permeable units. Dolomitization of the subtidal mud-supported sediments converts the tight, mud-supported limestones to permeable units because of the larger dolomite crystals and intercrystalline pore space. This produces two types of flow units in the subtidal interval: a dolomud-supported flow unit and a dolograin-supported flow unit. Each will have a unique porosity-permeability transform.
Karst-collapse reservoir model
The karst-collapse model describes a touching vug pore system that is formed by massive dissolution of carbonate resulting from meteoric groundwater flow and subsequent collapse and filling of caverns (Figure 4). This process is independent of the sediment's original environment of deposition. Most karsted carbonates are thought to be related to major dissolution concentrated in the vadose and upper phreatic zones producing a horizontality to the caverns. Cavern geometry is also controlled by fracture orientation, often resulting in caverns with linear trends.
Three karst facies can be described that relate to touching vug pore geometries. The roof facies is characterized by dissolution-enhanced fractures formed by cavern collapse. The cave facies is characterized by infill sediment. The floor facies is characterized by collapse breccia from roof collapse and small caverns. Karsted reservoirs are generally highly compartmentalized and very complex.
Geological reef reservoir model
The geological reef model is a composite of the upward-shoaling subtidal-supratidal and karst-collapse reservoir models. The difference is that the facies tracts are compressed onto a carbonate shelf of limited areal extent with high relief above the seafloor and with steeply sloping sides. The interior shelf or lagoon facies (Figure 2) located landward of the shelf edge normally contains a high percentage of mud. Grainstones, packstones, and boundstones associated with the reef facies are typically found along the shelf edge.
Compaction and cementation typically destroy the permeability of the lagoonal muds, leaving the grain-dominated sediments and boundstones of the reef edge as reservoir rocks. However, selective leaching, dolomitization, and karsting can significantly alter the permeability patterns, as discussed in previous sections. The reservoir flow units can be very complex due to the numerous possible combinations of depositional and diagenetic events.
- Introduction to geological methods
- Evaluating tight gas reservoirs
- Reservoir quality
- Fluid contacts
- Geological heterogeneities
- Evaluating fractured reservoirs
- Evaluating stratigraphically complex fields
- Evaluating diagenetically complex reservoirs
- Dunham, R. J., 1962, Classification of carbonate rocks according to depositional texture, in, Classifications of Carbonate Rocks—A Symposium: AAPG Memoir 1, p. 108–121.
- Swanson, R. G., 1981, Sample examination manual: AAPG Methods in Exploration 1, 65 p.
- Scholle, P. A., D. G. Bebout, and C. H. Moore, eds., 1983, Carbonate depositional environments: AAPG Memoir 33, 708 p.
- Mcllreath, I. A., James, N. P., 1984, Carbonate slopes, in Walker, R. G., ed., Facies Models: Geoscience Canada Reprint Series 1, p. 245–258.
- Archie, G. E., 1952, Classification of carbonate reservoir rocks and petrophysical considerations: AAPG Bulletin, v. 36, p. 278–298.
- Lucia, F. J., 1995, Rock-Fabric/Petrophysical Classification of Carbonate Pore Space for Reservoir Characterization: AAPG Bulletin, v. 79, p. 1275-1300.
- Lucia, F. J., 1983, Petrophysical parameters estimated from visual descriptions of carbonate rocks—a field classification of carbonate pore space: Journal of Petroleum Technology, March, p. 629–637.
- Kerans, C. 1989, Karst-controlled reservoir heterogeneity and an example from the Ellenburger Group (Lower Ordovician) of west Texas: The Univ. of Texas Bureau of Economic Geology Report of Invetigations, n. 186, 40 p.