Geostatistical Interpolation Constrained by Lithofacies

Summary Geostatistical technologies, which include two-point geostatistics (TPS) and multiple-point geostatistics (MPS), are significant in both geological modelling and geophysical inversion. However, TPS is incapable to characterise complex geological structures, due to the dependency on the variogram. When using MPS in the simulation of the continuous variables such as velocity, the calculation and memory burdens are heavy. We integrate TPS with MPS to release these issues. Based on kriging theory, we present a geostatistical interpolation strategy constrained by the lithofacies probability distribution, which is obtained by MPS. The proposed approach utilizes the advantages of TPS and MPS simultaneously. The model test illustrates the effect of the method on depicting complicated geological structure.