PCTO-SIM

Multiplepoint Geostatistics is a well-known general statistical framework by which complex geological phenomena have been modeled efciently. Pixel-based and patch-based are two major categories of these methods. In this paper, the optimization-based category is used which has a dual concept in texture synthesis as texture optimization. Our extended version of texture optimization uses the energy concept to model geological phenomena. While honoring the hard point, the minimization of our proposed cost function forces simulation grid pixels to be as similar as possible to training images. Our algorithm has a self-enrichment capability and creates a richer training database from a sparser one through mixing the information of all surrounding patches of the simulation nodes. Therefore, it preserves pattern continuity in both continuous and categorical variables very well. It also shows a fuzzy result in its every realization similar to the expected result of multi realizations of other statistical models. While the main core of most previous Multiplepoint Geostatistics methods is sequential, the parallel main core of our algorithm enabled it to use GPU efficiently to reduce the CPU time. One new validation method for MPS has also been proposed in this paper. The main core of the algorithm has a totally parallel structure and fully adapted for GPU.It has the fastest exhaustive search in the training database in MPS methods.Our method has a self-enrichment capability and creates new patches out of TIs.It outperforms DS, Bunch DS, Quilting in sparse and dense conditional simulation.Our algorithm is able to extend to higher dimensions (4) so simply.