Conditional simulation of multi-type non stationary Markov object models respecting specified proportions (Q2454625)

Boolean model simulation has been widely used during the last two decades to simulate sedimentary bodies (especially in fluvio-deltaic environments). A key issue however subsisted: one of the most important parameters for object model simulation, namely the (non stationary) intensity of the underlying object process, is not provided by the end user, but must instead be computed from other input parameters, such as local proportions of lithofacies, erosion rule between objects of different types, and interaction between objects. The present paper revisits a birth and death algorithm for simulating multi-type object models based on point processes with interaction. It provides accurate approximations for computing the local intensity of the underlying point process taking into account (i) non stationary properties and non stationary object parameters; (ii) erosion rules in the case of multi-type objects; (iii) attraction or repulsion between objects. Simulated examples show that this algorithm is able to reproduce the desired proportions. Finally, important issues for implementing the algorithm are discussed.