The present disclosure describes a method that improves the long-range geobody continuity in Multiple Point Statistical methods, wherein the coarsest multi-grid level cells are simulated in a regular path, and the subsequent level cells are simulated in a random path as usual. The method is general and is applicable to different cases: such as hard data conditioning, soft data conditioning, non-stationarity modeling, 2 or more than 2 types of facies modeling, and 2D and 3D modeling. The method is particularly useful in reservoir modeling, especially for the channelized systems, but can be generally applied to other geological environments.