High-resolution truncated plurigaussian simulations for the characterization of heterogeneous formations