Flow simulation with reactive transport applied to carbonate rock diagenesis

Significant volumes of the known hydrocarbon reserves are found in carbonate rocks, many of these dolomitized. The spatial distribution of diagenesis on these rocks is one of the main challenges in oil reservoir modeling. Reactive transport models can be a powerful tool to understand the active diagenetic processes and their effects on the quality of these reservoirs. In this study it was used, for the first time, the CMG-GEM simulator to model diagenetic evolution of a carbonate sequence, subjected to compaction-driven and geothermal flow in a simulated period of 200 thousand years. It was simulated carbonate cementation, dolomitization and dissolution, with and without presence of faults. Among the analyzed variables, the volume of circulating fluid was the most important factor. For both mechanisms, flow simulated velocities obtained had magnitudes smaller than 10−6 m/day. Diagenesis was insignificant for these low speeds at simulated time interval. Only dolomitized facies presented relevant diagenesis in form of calcite dissolution and dolomite precipitation. Simulations with flow rates of 1 m/day revealed a considerable increase in observed diagenesis, especially in carbonate cementation and in porosity enhancement. Diagenesis was more pronounced in more permeable sediments, highlighting the role of fluid flow in diagenetic reactions. Relative dissolution was greatly reduced during simulations performed with absence of dolomite and dolomitization reactions. The presence of faults strongly influences spatial distribution of diagenesis, especially relatively to dissolution. More permeable facies were more dissolved near fault, decreasing with increasing distance. Low permeability facies, as mudstones, are not dissolved, even near fault. Spatial distribution of diagenesis would then be dependent mainly on the quality of original pore structure, of fault presence and mineral composition of rock.