Impact of Diagenetic Alterations on the Petrophysical and Multiphase Flow Properties of Carbonate Rocks Using a Reactive Pore Network Modeling Approach

Sedimentary reservoir rocks generally have complex and heterogeneous pore networks that are related to the original depositional rock texture and subsequent diagenetic alterations. Such alterations are in part controlled by the original mineralogy and sedimentological facies, the compaction history, the involved fluids (and rock/fluid interactions), the flow history and the related physico-chemical conditions. During the diagenetic evolution (paragenesis), cycles of alternating dissolution (porosity enhancement) and precipitation (porosity destruction) caused by changes in chemical and thermodynamic conditions may lead to heterogeneous rock structure at both local and reservoir scale. In the absence of cored plugs to measure the petrophysical properties (i.e. porosity, permeability and formation factor) and multiphase flow properties (i.e. capillary pressure, relative permeability and resistivity index), a numerical tool that calculates these properties from pore structure data by predicting its evolution during the diagenetic cycle is of great interest for the petroleum industry and reservoir characterization studies. A Pore Network Model (PNM) provides opportunities to study transport phenomena in fundamental ways because detailed information is available at the pore scale. It has been used over the last decades to understand basic phenomena such as capillarity, multiphase flow or coupled phenomena. In particular, this modeling approach is appropriate to study the rock/fluid interactions since the mass exchange at surfaces can be modeled explicitly. It can provide quantitative information both on the effective transport property modifications due to the reactions and on the structure evolution resulting from dissolution/precipitation mechanisms. In the present paper, this approach is used to study the effect of the diagenetic cycle on the petrophysical properties of carbonate rocks. It involves three discrete steps. The first step consists of replacing the original complex pore structure of real porous media by a conceptual network. The second step consists of resolving the governing equations of the precipitation and dissolution phenomena (i.e. reactive convection diffusion equation) in the conceptual 3D pore network and deducing the local reactive fluxes and the motion of the fluid-solid interface. The third step consists of updating the new pore structure and calculating the new petrophysical properties of the modified porous media. Those steps are re