Quartz Cementation In Sandstones Research Articles

Impact of grain-coating illite in preventing quartz cementation: Example from permo-carboniferous sandstone, Central Saudi Arabia

The origin and role of clay coatings (e.g., chlorite) in inhibiting quartz cementation in sandstones have been an area of research during the last three decades. However, the role of grain-coating illite in porosity preservation has received little attention. An integrated petrographic-mineralogical study is conducted to investigate the mode of occurrence of illite coatings, their coverage, thickness, and role in inhibiting the development of quartz overgrowth in a Permo-Carboniferous sandstone reservoir in Central Saudi Arabia. Results of the study reveal that illite coatings occur in two forms: tangential and radial to grain surfaces, both of which are effective in arresting quartz cementation. Illite-coating coverage exhibits and inverse relationship with quartz cement volume, and the coating coverage has been found to be grain-size dependent. Coarser-grained sandstones have higher coating coverage due to smaller surface area and smaller amount of illite required to form complete coatings. The study has found that thick illite coating (up to 27 μm) and high illite-coating coverage (≥80%) effectively arrest quartz cementation. Total pore-filling cements threshold of >20% has a negative impact on porosity. The findings of this study can be applied to other reservoirs of similar depositional environments to improve their reservoir quality prediction.

Rock Acoustics of Diagenesis and Cementation

We simulate the effects of diagenesis, cementation and compaction on the elastic properties of shales and sandstones with four different petro-elastic theories and a basin-evolution model, based on constant heating and sedimentation rates. We consider shales composed of clay minerals, mainly smectite and illite, depending on the burial depth, and the pore space is assumed to be saturated with water at hydrostatic conditions. Diagenesis in shale (smectite/illite transformation here) as a function of depth is described by a fifth-order kinetic equation, based on an Arrhenius reaction rate. On the other hand, quartz cementation in sandstones is based on a model that estimates the volume of precipitated quartz cement and the resulting porosity loss from the temperature history, using an equation relating the precipitation rate to temperature. Effective pressure effects (additional compaction) are accounted for using the Athy equation and the Hertz–Mindlin model. The petro-elastic models yield similar seismic velocities, despite the different levels of complexity and physics approaches, with increasing density and seismic velocities as a function of depth. The methodology provides a simple procedure to obtain the velocity of shales and sandstones versus temperature and pressure due to the diagenesis-cementation-compaction process.

The Importance of Effective Stress History on Quartz Cementation in Sandstones: Evidence from the Fulmar Formation, Central North Sea

The Importance of Effective Stress History on Quartz Cementation in Sandstones: Evidence from the Fulmar Formation, Central North Sea

The origin of silica cements revealed by spatially resolved oxygen isotope microanalysis and electron-beam microscopy; Heidelberg Formation, Germany

The origin and growth mechanisms of microcrystalline quartz cement in sandstones and its relationship to quartz overgrowths remain puzzling even though it has an important effect on rock properties in deeply buried sandstones; microcrystalline quartz cement inhibits the growth of quartz cement and so leaves sandstones with anomalously high porosity deep in sedimentary basins. Here we have examined the relationships between quartz cement, and duplicate layers of chalcedony/amorphous silica and microcrystalline quartz from outcrops of the Upper Cretaceous Heidelberg Formation, Germany. We have used an array of techniques including: NanoSIMS for spatially-resolved stable oxygen isotope analysis, and high-resolution chemical analysis of the silica cements using electron microprobe analysis, as well as crystallographic studies using TEM and EBSD. Based on burial history modeling and the low degree of compaction of the Heidelberg Formation sandstones, we conclude that maximum temperature due to conductive heating was about 42 °C; quartz overgrowths with δ18O of 18.8 ‰ developed from meteoric-derived formation water with δ18O of about −12 ‰ V-SMOW. Based on published precedents and an awareness that silica solubility decreases with decreasing temperature, we assume that chalcedony/amorphous silica and then microcrystalline quartz developed during cooling of flowing formation water. Chalcedony with δ18O of 27.4 ‰ developed from meteoric-derived formation water with δ18O of between −1 and −3 ‰ V-SMOW depending on whether the initial temperature of the cooling formation water was 45 or 55 °C. The interpreted higher δ18O of the warmer water responsible for chalcedony and microcrystalline quartz growth is typical of basin-center waters that have undergone considerable water–rock interaction. It is likely that the fluid movement was caused by thrust tectonics in the Coniacian to Campanian. The duplicate layers of chalcedony/amorphous silica and then microcrystalline quartz suggests that there were two thrust-induced fluid expulsion events.

Vertical effective stress and temperature as controls of quartz cementation in sandstones: Evidence from North Sea Fulmar and Gulf of Mexico Wilcox sandstones

Abstract We present quantitative petrographic data, high spatial resolution oxygen isotope analyses of quartz cement, basin modelling and a kinetic model for quartz precipitation for two Paleocene-Eocene Wilcox Group sandstones from Texas and two Jurassic Fulmar Formation sandstones from the Central North Sea. At each location, one sandstone has been buried to ca. 145 °C and one to ca. 185 °C. A key difference between the Wilcox and Fulmar burial histories is that the Wilcox sandstones are currently at higher vertical effective stresses and, from basin modelling studies, have been subjected to generally higher vertical effective stresses through their burial history. The amounts of quartz cement in the Wilcox sandstones are between 12 and 18%, and between 2 and 6% in the Fulmar sandstones. High-spatial-resolution oxygen isotope data obtained from the quartz cements suggest temperature ranges for quartz precipitation from 60 to 80 °C to values approaching maximum burial temperature. Factors such as grain coatings or the timing of petroleum emplacement cannot explain the differences in the amounts of quartz cement. Petrographic data show that most of the silica for quartz cement can be derived from intergranular pressure dissolution. Although the sample set is small, we interpret the results to suggest that the differences in quartz cementation in Fulmar and Wilcox sandstones can be explained better by differences in their vertical effective stress history than their temperature history; in this case, the supply of silica rather than the precipitation of quartz becomes an important control on the rate and extent of cementation.

Clay‐coat diversity in marginal marine sediments

AbstractThe specific mineralogy of clay grain coats controls the ability of the coat to inhibit quartz cementation in sandstones during prolonged burial and heating. How and why clay‐coat mineralogy varies across marginal marine systems is poorly understood, even though these eogenetic phenomena strongly influence subsequent mesodiagenesis and reservoir quality. The novel development of the ability to predict the distribution of clay‐coat mineralogy would represent an important development for sandstone reservoir quality prediction. In marginal marine sediments, clay minerals occur as grain‐coats, floccules, mud intraclasts, clay‐rich rock fragments or as dispersed material. However, the relationships between clay mineralogy, the amount of clay, and its distribution is poorly understood. This study focused on the Ravenglass Estuary, UK. The key aim was to develop and apply a novel methodology utilising scanning electron microscope – energy dispersive spectrometry, for the first time, on grain coats in modern sediments, to differentiate the clay‐coat mineral signature from that of the bulk sediment, and reveal the distribution of clay minerals across marginal marine sediments. The study showed that marginal marine sediments principally have their clay mineral assemblage present as clay‐coats on sand grains. These clay‐coats have a mixed clay mineralogy and are spatially heterogeneous across the range of marginal‐marine depositional environments. The study further showed that clay‐coat mineralogy is governed initially by the hydrologically‐controlled segregation of the clay minerals within inner estuarine depositional environments, and subsequently by the selective abrasive removal of specific clay mineral types during reworking and transport into the outer estuary and the marine environment. The highest relative abundance of grain‐coating chlorite was in sand‐flat and tidal‐bar depositional environments. The availability of an analogue data set, and an understanding of the controlling processes of clay‐coat mineralogy, offer crucial steps in building a predictive capability for clay‐coat derived elevated reservoir quality in deeply buried sandstones.

Vertical effective stress as a control on quartz cementation in sandstones

Abstract Temperature-controlled precipitation kinetics has become the overwhelmingly dominant hypothesis for the control of quartz cementation in sandstones. Here, we integrate quantitative petrographic data, high spatial resolution oxygen isotope analyses of quartz cement, basin modelling and a kinetic model for quartz precipitation to suggest that the supply of silica from stress-sensitive intergranular pressure dissolution at grain contacts is in fact a key control on quartz cementation in sandstones. We present data from highly overpressured sandstones in which, despite the current burial temperature of 190 °C, quartz cement occurs in low amounts (4.6 ± 1.2% of bulk volume). In situ oxygen isotope data across quartz overgrowths suggest that cementation occurred over 100 Ma and a temperature range of 80–150 °C, during which time high fluid overpressures resulted in consistently low vertical effective stress. We argue that the very low amounts of quartz cement can only be explained by the low vertical effective stress which occurred throughout the burial history and which restricted silica supply as a result of a low rate of intergranular pressure dissolution at grain contacts.

Cementation and structural diagenesis of fluvio-aeolian Rotliegend sandstones, northern England

Quartz cementation in sandstones is closely linked to grain coating phases and diagenetic alteration. Grain coatings consisting of illite–smectite stained with iron oxides and hydroxides are able to preserve large amounts of porosity by preventing the formation of syntaxial quartz overgrowth cement. The Penrith Sandstone Formation was chosen as an analogue for Rotliegend reservoirs to test the impact of grain coatings on quartz cementation. This adds to an existing model of cementation. Differences of grain coating coverage can be linked to grain size. Extensive grain coatings are present in finer grained laminae in some samples. Coarser grained laminae contain less extensive grain coatings. The analysis of grain coatings based on standard petrographic analyses is combined with high-resolution QEMSCAN® analyses. Structural features include deformation bands of different ages. Diagenetic alterations around faults, recorded by grain coatings, allow the delineation of relative temporal relations, revealing at least two generations of deformation band formation associated with normal faulting. In the Vale of Eden succession one normal faulting event postdates burial diagenetic quartz cementation as is evident by fault focused fluid flow and associated bleaching of iron and absence of quartz overgrowth.

Overpressure preventing quartz cementation? – Comment to Stricker et al. (2016)

There is increasing evidence that quartz cementation can be viewed as a process controlled by temperature and insensitive to effective stress. This view of quartz cementation in sandstones is often referred to as the illite-mica induced dissolution model (IMID), which assumes quartz dissolution to occur along stylolites and clay laminae rather than at quartz-quartz grain contacts. In the present comment it is argued that the exceptional reservoir quality in the Skagerrak Formation of the Heron Cluster, North Sea, is due to grain coatings and rapid Cenozoic burial limiting the exposure to quartz cementation. This line of reasoning implies overpressure has had neglectable porosity preserving effect in the Heron Cluster.

Insights into rates of fracture growth and sealing from a model for quartz cementation in fractured sandstones

A new model accounts for crystal growth patterns and internal textures in quartz cement in sandstone fractures, including massive sealing deposits, thin rinds or veneers that line open fracture surfaces, and bridge structures that span otherwise open fractures. High-resolution cathodoluminescence imaging of bridge structures and massive sealing deposits indicates that they form in association with repeated micron-scale fracturing of growing quartz crystals, whereas thin rinds do not. Model results indicate that the three morphology types develop in response to (1) the ratio of the rates of quartz growth to fracture opening and (2) the substantially faster growth rate that occurs on noneuhedral surfaces in certain crystallographic orientations compared to euhedral crystal faces. Rind morphologies develop when the fracture opening rate exceeds two times the fastest rate of quartz growth (along the c axis on noneuhedral surfaces) because growing crystals develop slow-growing euhedral faces. Massive sealing, on the other hand, develops where the net rate of fracture opening is less than twice the rate of quartz growth on euhedral faces because all quartz growth surfaces along the fracture wall seal the fracture between fracturing events. Bridge structures form at fracture opening rates that are intermediate between the massive sealing and rind cases and are associated with crystallographic orientations that allow growth to span the fracture between fracturing events. Subsequent fractures break the spanned crystal, introducing new, fast-growing noneuhedral growth surfaces where quartz grows more rapidly compared to the euhedral faces of nonspanning crystals. As the ratio of fracture opening to quartz growth rate increases, the proportion of overgrowths that span the fracture decreases, and the range in c -axis orientations for these crystals comes progressively closer to perpendicular to the fracture wall until the maximum spanning limit is reached. Simulation results also reproduce “stretched crystal,” “radiator structure,” and “elongate blocky” textures in metamorphic quartz veins. The model replicates a well-characterized quartz bridge from the Cretaceous Travis Peak Formation as well as quartz cement abundances, internal textures, and morphologies in the sandstone host rock and fracture zone using the same kinetic parameters while honoring fluid-inclusion and thermal-history constraints. The same fundamental driving forces, in both in the host rock and fracture system, are responsible for quartz cementation, with the only significant difference within the fracture zone being the creation of new pore space as well as new noneuhedral surfaces for cases where overgrowths span fractures between fracturing events. Rates of fracture growth and sealing may be inferred from fracture cement textures using model results.

Global sensitivity analysis through polynomial chaos expansion of a basin-scale geochemical compaction model

We present a model-driven uncertainty quantification methodology based on sparse grid sampling techniques in the context of a generalized polynomial chaos expansion (GPCE) approximation of a basin-scale geochemical evolution scenario. The approach is illustrated through a one-dimensional example involving the process of quartz cementation in sandstones and the resulting effects on the dynamics of the vertical distribution of porosity, pressure, and temperature. The proposed theoretical framework and computational tools allow performing an efficient and accurate global sensitivity analysis (GSA) of the system states (i.e., porosity, temperature, pressure, and fluxes) in the presence of uncertain key mechanical and geochemical model parameters as well as boundary conditions. GSA is grounded on the use of the variance-based Sobol indices. These allow discriminating the relative weights of uncertain quantities on the global model variance and can be computed through the GPCE of the model response. Evaluation of the GPCE of the model response is performed through the implementation of a sparse grid approximation technique in the space of the selected uncertain quantities. GPCE is then be employed as a surrogate model of the system states to quantify uncertainty propagation through the model in terms of the probability distribution (and its statistical moments) of target system states.

Smectite-illite-muscovite transformations, quartz dissolution, and Silica release in shales

AbstractQuantitative analysis of the smectite-to-illite and illite-to-muscovite transformations indicates that 17–28 wt.% SiO2 and 17–23 wt.% SiO2, respectively, are liberated during these reactions, assuming that Al is conserved. Dissolution of quartz silt in shales yields up to 6–9% SiO2 in the range up to 200°C and a further 10–15% SiO2 in the 200–500°C range. For muds altered to shales at 200°C, 14–20 wt.% silica is evolved. From 200 to 500°C, a further 18–28 wt.% silica is evolved. Additional small amounts of silica may be released in the alteration of feldspar to clay and by stylolitization of quartz silt. Thus, in the burial and temperature range of diagenesis into the epizone, major quantities of silica are released from clays and by quartz dissolution in shales. Within this range of alteration, concomitant decline of whole-rock Si/Al (SiO2/Al2O3) in the transformation of smectite to illite to muscovite suggests the liberated silica migrates from the source shale. As a result, the metamorphosed shales are more micaceous and less quartzose than their progenitor muds. In the diagenetic zone and anchizone, the evolved silica is probably a major source of quartz cement in sandstones. In the epizone, evolved silica is commonly present in quartz veins in the parent rocks. Fluid-inclusion temperatures in quartz overgrowths and fracture fillings in some sandstones suggest that some cements may have been derived from downdip basinal shales and pressure solution in sandstones.

Compaction microstructures in quartz grains and quartz cement in deeply buried reservoir sandstones using combined petrography and EBSD analysis

Sandstone diagenesis involves porosity loss by mechanical and chemical compaction. In this study we have examined relations of diagenesis and microstructures in quartz grains and quartz cement in sandstones from relatively deep burial depths offshore mid Norway. Electron backscatter diffraction analysis (EBSD) has been combined with optical and cathodoluminescence petrography study of samples with different degrees of compaction and quartz cementation. Quartz cement is shown to be syntaxial to the nearest host domain (grain or subgrain), both for monogranular, polygranular, undeformed and deformed quartz. The quartz cement growth may have been initiated at different subgrain surfaces, and preserving pre-depositional deformation structures. EBSD inverse pole figure imaging shows dauphiné twins to be common in all samples, both in quartz grains and quartz cement. The dauphiné twins appear in grain-grain contacts and in cement-crystal boundaries, and commonly crossing grain–cement boundaries. On basis of twin distributions we suggest that both inherited twins from the source area and twins formed by compaction-induced grain boundary deformation are present. A general problem in sedimentary rocks is to distinguish inherited from sedimentary compaction-induced intra-grain deformation structures. This study shows that the role of ductile deformation mechanisms (e.g. twinning) in sediment compaction can be clarified by more comprehensive following up studies using EBSD analysis in combination with other techniques.

Sandstone porosity related to vitrinite reflectance

Vitrinite reflectance (VR) and the volume fraction of quartz cement in sandstones are both measures of thermal exposure. This paper shows how porosity reduction by quartz cementation can be related to VR for heating at a constant rate, using simple models for quartz cementation and VR. This relationship reveals that the heating rate is an explicit parameter, which means that there is no pure relationship between porosity (volume fraction quartz cement) and VR. However, the rate of quartz cementation and the rate of change of TTI (time–temperature integral) have a close exponential dependence on temperature. Porosity reduction by quartz cementation as a function of VR therefore nearly follows the same curve for very different temperature histories. This property of porosity reduction by quartz cementation is used to calibrate average specific surface parameters as a function of porosity. The calibration procedure is demonstrated on two datasets, where the average specific surface depends on porosity to the power of an exponent that is somewhat greater than 3.

Quartz Cementation in Sandstones.

Quartz Cementation in Sandstones.

Making diagenesis obey thermodynamics and kinetics: the case of quartz cementation in sandstones from offshore mid-Norway

Calculation of the quantity and distribution of quartz cement as a function of time and temperature/depth in quartzose sandstones is performed using a coupled dissolution/diffusional–transport/precipitation model. This model is based on the assumptions that the source of the silica cement is quartz surfaces adjoining mica and/or clay grains at stylolite interfaces within the sandstones, and the quantity of silica transport into and out of the sandstone by advecting fluids is negligible. Integration of the coupled mass transfer/transport equations over geologically relevant time frames is performed using the quasi-stationary state approximation. Results of calculations performed using quartz dissolution rate constants and aqueous diffusion coefficients generated from laboratory data, are in close agreement with both the overall porosity and the distribution of quartz cement in the Middle Jurassic Garn Formation only after optimizing the product of the effective surface area and quartz precipitation rate constants with the field data. When quartz precipitation rate constants are fixed to equal corresponding dissolution rate constants, the effective surface area required to match field data depends on the choice of laboratory generated quartz rate constant algorithm and ranges from 0.008 cm −1 to 0.34 cm −1. In either case, these reactive surface areas are ∼2 to 4 orders of magnitude lower than that computed using geometric models.

Can oil emplacement prevent quartz cementation in sandstones?

In the current debate on the ability of oil emplacement to halt quartz cementation in sandstone oil reservoirs, this paper supports a yes, it can answer. Although we find some of the empirical data petroleum fluid inclusion abundances, quartz-depth trends) to be inconclusive, we have made progress in this debate by examining theoretical controls for three end-member quartz cementation scenarios. (1) Silica externally sourced via advection. Oil emplacement must halt quartz cementation because of the associated massive reduction in relative permeability to water. (2) Silica supplied internally, water-wet reservoir. Quartz cementation will be inhibited because the rate of diffusion will be reduced by orders of magnitude, making diffusion the rate controlling step. The degree of inhibition will depend on rock fabric and oil saturation; at very high oil saturations quartz cementation could effectively cease. (3) Silica supplied internally, oil-wet reservoir. Quartz cementation will be halted as water cannot access to grain surfaces.

TEMPERATURE CONTROLLED POROSITY/PERMEABILITY REDUCTION, FLUID MIGRATION, AND PETROLEUM EXPLORATION IN SEDIMENTARY BASINS

Findings mainly from the Norwegian Continental Shelf indicate that at temperatures greater than approximately 60°C, internally sourced quartz cementation and diagenetic clay become important porosity and permeability reduction factors, respectively. Porosity loss due to quartz cementation in sandstones and siltstones proceeds mainly independent of effective stress or fluid pressure. Porosity loss rates approach those required to generate high fluid overpressures at approximately 120°C, which can result in hydrofracturing of the overlying low permeability shales. In thick low permeable sediment sequences, the probability of km scale vertical fluid migration increases. These thermally driven processes are capable of generating and sustaining fluid overpressure and facilitate fluid migration for up to tens of millions of years. Unlike mechanical compaction models, this model predicts that porosity loss and therefore fluid migration will continue despite the buildup of overpressure, even during periods of no sedimentation. In sealed compartments, fluid flow from more deeply buried high permeable lithologies will occur by hydrofracturing of overlying low permeable lithologies, preferentially along near vertical faults if present. When hydrofractures are induced from high permeable sediments, they will propagate vertically through overlying low permeable sediments, unless they enter lithologies with sufficient permeability and volume to bleed off the fracture fluid propagation pressure. Quantitative analysis of the mineral reactions causing porosity loss cannot only identify sediments in thermal zones which are expelling fluids, but also lithologies in thermal zones which are likely to receive those fluids, in addition to the timing and rates of fluid migration. The model has important applications for evaluating petroleum exploration risks and the potential for remigrated hydrocarbon plays both at basin scale and prospect level.

Ion microprobe analysis of 18O/16O in authigenic and detrital quartz in the St. Peter Sandstone, Michigan Basin and Wisconsin Arch, USA: Contrasting diagenetic histories

The oxygen isotopic compositions of authigenic quartz cements in sandstones provide a monitor of the temperatures, compositions, and origins of pore-occluding fluids during diagenesis, but quartz overgrowths are too fine-grained to be amenable to conventional isotopic analysis. We have used a Cameca ims-4f ion microprobe to determine oxygen isotopic variations in authigenic and detrital quartz in four samples of the Ordovician St. Peter Sandstone from the Michigan Basin and Wisconsin Arch, midwestern USA. Ion microprobe isotopic analyses have been successfully accomplished with an internal precision of ±1‰ (1 σ) and a spatial resolution of 20–30 μm at low mass resolution using a high voltage offset technique. Repeated analyses of the quartz standard demonstrate a reproducibility of close to ±1‰ (1 sd) in good agreement with that expected from counting statistics. The four ion probe samples were chosen as representative of thirteen samples that were mechanically and chemically disaggregated, sieved, and analysed as a function of sieve size by conventional bulk-mineral techniques. Conventional and ion microprobe analyses are mutually consistent, supporting the accuracy of the ion microprobe analyses. Within-sample isotopic variations of up to 13%. and micro-scale isotopic variations of at least 4‰ over a distance of 100 μm have been measured within quartz overgrowths in a sandstone from the Wisconsin Arch. Overgrowths are uniformly higher in δ18O than detrital grains, and gradients of up to 25% exist across a few microns.18O-enriched quartz overgrowths in sandstones from the Wisconsin Arch show complex CL zonation and reflect one of two possible processes: (1) low-temperature quartz precipitation during mixing of meteoric waters with upwelling basinal fluids; (2) higher temperature quartz precipitation during episodic gravity-driven upwelling of warm basinal fluids (of comparable isotopic composition to Michigan Basin fluids) from the Illinois Basin, related to evolution of Mississippi Valley type PbZn ore-forming fluids. Quartz overgrowths in Michigan Basin sandstones, which derived their silica locally by pressure solution, are thought to have precipitated continuously over a range of temperatures from hot basinal fluids of restricted isotopic composition and circulation. Detrital quartz shows significant intea-grain and intergrain isotopic variation (5–13‰) but no significant inter-sample variation, consistent with a uniform source of granitoid and metamorphic quartz in space and time.

Porosity/depth trends in reservoir sandstones: assessing the quantitative effects of varying pore-pressure, temperature history and mineralogy, Norwegian Shelf data

AbstractEmpirical porosity/depth trends from the Norwegian Continental Shelf are discussed in the light of models on compactional processes during burial. Linear best-fit lines are useful as a first approximation, illustrating the general decline in porosity with depth. Porosity increases non-linearly with pore-pressure, and notably higher than normal porosity is found in very high overpressure reservoirs. Hydrostatic pressure sandstones from the Haltenbanken area are on average more porous than equivalents from the northern Viking Graben. This may be related to the burial history, as generally more subsidence occurred in Late Pliocene and Pleistocene in the Haltenbanken basins. The porosity/depth data show considerable scatter, and this is to a large degree related to the initial mineralogical composition. Correlated with mineralogical analyses, the reservoir sandstones show a trend towards lower porosity with increasing clay content at shallow burial, probably due to more extensive mechanical compaction. At depths >3.0–3.5 km, the porosity loss is highest in sandstones with low clay contents. This may be explained by a lower rate of quartz cementation in sandstones with abundant intergranular clays or clay coatings around quartz grains. These data show that pre-burial mineralogy is a very important factor in porosity prediction of sandstones, in many cases more important than the degree of overpressure.