Geostatistical modeling of clay spatial distribution in siliciclastic rock samples using the plurigaussian simulation method

In order to implement secondary and enhanced oil recovery processes in complex terrigenous formations as is usual in turbidite deposits, a precise knowledge of the spatial distribution of shale grains is a crucial element for the fluid flow prediction. The reason of this is that the interaction of water with shale grains can significantly modify their size and/or shape, which in turn would cause porous space sealing with the subsequent impact in the flow. In this work, a methodology for stochastic simulations of spatial grains distributions obtained from scanning electron microscopy images of siliciclastic rock samples is proposed. The aim of the methodology is to obtain stochastic models would let us investigate the shale grain behavior under various physicochemical interactions and flux regimes, which in turn, will help us get effective petrophysical properties (porosity and permeability) at core scale. For stochastic spatial grains simulations a plurigaussian method is applied, which is based on the truncation of several standard Gaussian random functions. This approach is very flexible, since it allows to simultaneously manage the proportions of each grain category in a very general manner and to rigorously handle their spatial dependency relationships in the case of two or more grain categories. The obtained results show that the stochastically simulated porous media using the plurigaussian method adequately reproduces the proportions, basic statistics and sizes of the pore structures present in the studied reference images.