Stochastic simulation of patterns using distance based pattern modeling
This presentation is the property of its rightful owner.
Sponsored Links
1 / 37

Stochastic Simulation of Patterns using Distance-Based Pattern Modeling PowerPoint PPT Presentation


  • 90 Views
  • Uploaded on
  • Presentation posted in: General

Multiple-Point Geostatistics. Search Template 9x9 Inner Patch 5x5 Multiple-Grids 3 Number of Clusters 100. Proposed Method. Filtersim. Stochastic Simulation of Patterns using Distance-Based Pattern Modeling. Mehrdad Honarkhah. Motivation.

Download Presentation

Stochastic Simulation of Patterns using Distance-Based Pattern Modeling

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Stochastic simulation of patterns using distance based pattern modeling

Multiple-Point Geostatistics

Search Template 9x9

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100

Proposed Method

Filtersim

Stochastic Simulation of Patterns using Distance-Based Pattern Modeling

Mehrdad Honarkhah


Motivation

Motivation

  • Improving pattern reproduction

  • Less parameters – less user interaction

  • Reduce simulation time

  • Checking pattern/MPS reproduction


Stochastic simulation of patterns using distance based pattern modeling

Introduction

  • better, more realistic models should not require an increase in user-set parameters (Filtersim)

  • Why use a complexmethod when a much simpler one works just as well ?


Distance based methods what will be demonstrated

Distance-Based Methods:(what will be demonstrated)

In distance-based modeling, many of the tasks usually performed in multiple-point geostatistical algorithms can be carried out in a surprisingly simple yet powerful way

  • Second

  • First

  • Answer queries through several inference mechanisms

  • Store and organize a domain of knowledge about the TI

MPS

Realization

Distance-Based

Outline


Broad outlook on the workflow

Broad Outlook on the Workflow


Implementation details 1 pattern database

Implementation details(1) Pattern Database

Pattern Database

3×3 SearchTemplate

Training Image


Implementation details 2 mds mapping

Implementation details(2) MDS mapping


Implementation details 3 kernel mapping

Implementation details(3) Kernel Mapping

Cartesian Space

Feature Space Projection


Implementation details 4 kernel k means clustering

Implementation details(4) Kernel K-means Clustering

Cartesian Space

Feature Space Projection


Implementation details summary example

Implementation detailsSummary Example

δ12

Training Image

Multi Dimensional Scaling

δ13

δ14

δ23

δ24

δ34

Kernel Space Mapping

Kernel K-Means


Workflow summary

WorkFlow Summary

?


Implementation details 5 simulation

Implementation details(5) Simulation

Recall: Simulation Algorithm in Filtersim

  • Classify training image patterns into clusters using filter scores and partition

  • Loop through all nodes of the simulation grid

    • Retrieve the data event at that node

    • Find the most similar cluster prototype to that data event

    • Randomly pick a pattern from that cluster

    • Paste it on the simulation grid

  • end


Implementation details 5 simulation1

Implementation details(5) Simulation

Simulation Algorithm in the Proposed Method

  • Classify training image patterns into clusters using Kernel k-means

  • Loop through all nodes of the simulation grid

    • Retrieve that data event at that node

    • Find the most similar cluster prototype to Data Event

    • Randomly pick a pattern from that cluster

    • Paste it on the simulation grid

  • end

Note: The basic Filtersim algorithm is maintained.

only the modeling of patterns changes


Mps simulation examples 1 2

MPS Simulation Examples (1/2)

Search Template 9x9

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Mps simulation examples 2 2

MPS Simulation Examples (2/2)

Search Template 11x11

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Workflow summary1

WorkFlow Summary


Reducing user interaction 1

Reducing user interaction (1)

Automatic Template Size Selection:

  • Calculating mean Entropy at different template Dimensions

  • Calculating second derivative of the mean Entropy Curve

  • Calculating Profile Log-Likelihood of the resulting Curve

  • Template Size = Maximum in the Profile


Reducing user interaction 11

Reducing user interaction (1)

Automatic Template Size Selection:

5 x 5

Training image

Mean Entropy

Log-Likelihood

Template Size

Template Size

1

2

3

4

5

1

2

3

The actual smallest template size is 5 x 5

4

5


Reducing user interaction 12

Reducing user interaction (1)

Automatic Template Size Selection:

13 x 13

Training image

Mean Entropy

Log-Likelihood

Template Size

Template Size


Reducing user interaction 2

Reducing user interaction (2)

Find Number of Clusters:

  • Eigenvalue decomposition of K =

  • Plot

  • Calculate Profile Log-Likelihood

    Number of Clusters = Maximum in the profile


Comparison with filtersim

Comparison with Filtersim

Find Number of Clusters:

  • Eigenvalue decomposition of K =

  • Plot

  • Calculate Profile Log-Likelihood

    Number of Clusters = Maximum in the profile

SIMULATION EXAMPLES


Comparison with filtersim 1 6

Comparison with Filtersim (1/6)

101 x 101 Training Image

Proposed Methodology Realizations

FilterSim Realizations

Search Template 9x9

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Comparison with filtersim 1 61

Comparison with Filtersim (1/6)

101 x 101 Training Image

Proposed Methodology Realizations

FilterSim Realizations

Search Template 13x13

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Comparison with filtersim 2 6

Comparison with Filtersim (2/6)

101 x 101 Training Image

Proposed Methodology Realizations

FilterSim Realizations

Search Template 9x9

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Comparison with filtersim 2 61

Comparison with Filtersim (2/6)

101 x 101 Training Image

Proposed Methodology Realizations

FilterSim Realizations

K = 100

K = 100

Search Template 13x13

Inner Patch 7x7

Multiple-Grids 3

Number of Clusters -----

K = 400


Comparison with filtersim 3 6

Comparison with Filtersim (3/6)

111 x 111 Training Image

FilterSim Realizations

Proposed Methodology Realizations

Search Template 15x15

Inner Patch 9x9

Multiple-Grids 3

Number of Clusters 1000


Comparison with filtersim 4 6

Comparison with Filtersim (4/6)

101 x 101 Training Image

FilterSim Realizations

Proposed Methodology Realizations

Search Template 11x11

Inner Patch 5x5

Multiple-Grids 3

Number of Clusters 100


Comparison with filtersim 5 6

Comparison with Filtersim (5/6)

169 x 169 Training Image

FilterSim Realizations

Proposed Methodology Realizations

Search Template 15x15

Inner Patch 11x11

Multiple-Grids 3

Number of Clusters 200


Comparison with filtersim 6 6

Comparison with Filtersim (6/6)

69 x 69x39 Training Image

Proposed Methodology Realizations

FilterSim Realizations

Search Template 15x15x9

Inner Patch 9 x 9 x 5

Multiple-Grids 3

Number of Clusters 200


Studying variability between realizations

Studying variability between realizations

How can we check the pattern/multiple point statisticalreproduction of the realizations with respect to the training image ?


Studying variability between realizations1

Studying variability between realizations

Variability between realizations correspond to the variability obtained while honoring the multiple-point information of the training image.

Variability Quantification

  • Generate N realizations

  • Find a distance function that characterizes the variability between any two realizations

  • Map the N realizations into a Cartesian space using MDS

  • The cloud of points in that space represents variability


Studying variability between realizations2

Studying variability between realizations

Measure of Variability between two realizations : l1and l2

f ( l1 ,l2 ) =JSD{ MPH(l1),MPH ( l2) }

Multiple-Point Histogram (MPH)

Using a multiple-point template, scan the realizations

Store the frequency of a specific configuration of outcomes

Results in the multiple-point histogram

Jensen Shannon Divergence (JSD)

A measure of similarity between two probability distributions: pandq


Studying variability between realizations3

Studying variability between realizations

Red points : Filtersim realizations Black points : Proposed Method realizations

Training Image Used

Filtersim Method

Proposed Method


Studying variability between realizations4

Studying variability between realizations

Training Image Used

CONCLUSION

Filtersim has less pattern reproduction than our method, and therefore increased artificial variability.

Great way to diagnose pattern reproduction (i.e. by looking at ALL realizations at the same time).


Conclusion

Conclusion

  • Distance method: easy to implement, very few user-set parameters

  • Distance methods also allow easy evaluation of MPS reproduction

  • Future work: use model expansion techniques to generate additional patterns to further increase geological realism and data conditioning


Conclusion1

Conclusion


  • Login