Dolomitization: from conceptual to numerical models
Authors: Whitaker Fiona F. , Smart Peter L. , Jones Gareth D. 
Published: April 2008

Abstract
Dolomitization requires not only favourable thermodynamic and kinetic conditions, but also a fluid-flow mechanism to transport reactants to and products from the site of dolomitization. This paper reviews work that seeks to provide a quantitative framework for conceptual models of dolomitization, using analytical and, particularly, numerical simulation models of fluid flow and rock-water interaction. This approach is starting to yield new insights into the major controls on the rate and pattern of fluid flux, and the resultant dolomitization.
Three sets of forces can drive the fluid flow required for dolomitization: elevation (topographic) head of meteoric water and/or seawater; gradients in fluid density due to variation in salinity and/or temperature; and pressure due to sedimentological and/or tectonic compaction. However, in many situations individual flow mechanisms may not operate in isolation. Rather fluid flow will commonly be a product of a number of different drives acting simultaneously. The balance between drives will change over time with variations in relative sea-level, climate, platform geometry and palaeogeography (which collectively comprise the critical boundary conditions). The simplistic prediction of dolomite body geometry from a single driving force may be misleading, as fluid flow will critically depend both on the boundary conditions and the distribution of permeability. Indeed, even for single driving forces, model predictions change significantly as simplistic assumptions are relaxed and these key parameters are specified with increasing realism.
The coupled modelling of dolomitization reactions within the flow field is less tractable than that of groundwater circulation because the kinetics of dolomitization are less well understood, particularly at lower temperatures. Dolomitization is likely to occur along a reaction front, where a favourable balance is struck between mass transport and reaction kinetics. For instance, in simulations of geothermal convection dolomitization focuses along the 50–60 °C isotherm. Dolomitization reactions are favoured by higher temperatures in deeper zones, but rates are limited by low flow because of lower permeability. Although flow rates are higher in shallow more permeable carbonates, lower temperatures limit reactions. High flow rates during reflux of platform-top brines give rapid dolomitization. This is associated with porosity occlusion in front of and behind the broad zone of replacement dolomitization driven by anhydrite cementation and overdolomitization, respectively. Lithological heterogeneities strongly affect the pattern of dolomitization, which is highly focused within more permeable beds and those with a higher reactive surface area.
While we focus here on dolomitization, models can also provide insights into diagenetic processes such as marine calcite cementation and aragonite, calcite and evaporite dissolution by refluxing brines, and by seawater circulation below the aragonite and calcite compensation depths. However, it is important to be aware of the assumptions and limitations of the numerical model(s) used. Particular attention must be paid to specification of boundary conditions, permeability and reactive surface area. The uncritical application of numerical techniques to particular cases of dolomitization is at best uninformative and at worst misleading. Careful application of these techniques offers great promise for well-constrained field problems, with greater inclusion of natural heterogeneity and time-variant boundary conditions. We also need to model feedbacks between diagenesis and porosity-permeability, and to include platform growth in simulations of slower diagenetic processes.