Simple Kriging Demo

Download as xlsx, pdf, or txt
Download as xlsx, pdf, or txt
You are on page 1of 6

Simple Kriging Demonstration

Michael Pyrcz, Geostatistics at Petroleum and Geosystems Engineering, University of Texas at Austin (mpy

1. Data and Estimate Locations and Values 2. Distance Matrix


Point x y value residual 0.00
1 60 80 0.1 -0.040 46.10
2 25 50 0.12 -0.020 72.80
3 80 10 0.2 0.060
unknown 50 50
mean 0.140

5. Covariance Matrix 6.Inverse Left Side


1.000 0.771 0.643 0.842 2.667 -1.644
0.771 1.000 0.666 0.875 -1.644 2.810
0.643 0.666 1.000 0.752 -0.621 -0.813

1: Data and Estimate Locations


100
90
80
70
Y Coordinate (m)

60
50 Data
Unknown
40
30

20
10
0
0 10 20 30 40 50 60 70 80 90 100
X Coordinate (m)

Description
This sheet provides an illustration of Simple Kriging at a single estimated location.

Step 1: Input the data locations and values, the unknown simulated location. At any point thes
may be changed to observed their influence on the simulation.
Step 2: The distance matrix is automatically calculated, that is the distance between the data a
Step 3: Enter the model of spatial continuity in the form of an isotropic spherical variogram and
should sum to one). This model may be changed at any time to observed sensitivities to spatia
Step 4: Variogram matrix is calculated by applying the distance matrix to the isotropic variogram
Step 5: Covaraince matrix is calculated by subtracing the variogram from the variance (1 for sta
This is applied to improve numerical stability as a diagonally dominant matrix is more readily in
Step 6: The left hand side of the covariance matrix is inverted.
Step 7: The inverted left handside matrix is multiplied by the right hand side matrix to calculate
Step 8: The kriging estimate and kriging variance are calculated with the weights and covarianc
Step 9: With the Gaussian assumption the complete local conditional cumulative distribution fu
niversity of Texas at Austin ([email protected])

2. Distance Matrix 3. Variogram Model


46.10 72.80 31.62 Nugget 0
0.00 68.01 25.00 Spherical 1
68.01 0.00 50.00 Range 300

4. Variogram Matrix
0.000 0.229 0.357 0.158
0.229 0.000 0.334 0.125
7. Weights 0.357 0.334 0.000 0.248
-0.621 0.341
-0.813 0.462 8. Kriging Results
1.941 0.225 Kriging Estimate 0.131
Sum Weights 1.029 Kriging Variance 0.139
Mean Weight -0.029
p-value f-1(p)
5% -0.48
10% -0.35 9. Local CCDF
15% -0.26 100%
Legend 20% -0.18
90%
Information 25% -0.12
User Input 30% -0.06 80%

Data Calculation 35% -0.01 70%


Unknown Redundancy Measures 40% 0.04
60%
Closeness Measures 45% 0.08
Percentile

50% 0.13 50%

55% 0.18 40%


60% 0.22
30%
65% 0.27
70% 0.33 20%

75% 0.38 10%


80% 0.44
0%
85% 0.52 -3.00 -2.00 -1.00 0.00 1.00 2.00
90% 0.61 Value
estimated location. 95% 0.74

lated location. At any point these locations and values

the distance between the data and the unknown locations.


sotropic spherical variogram and nugget effect (contributions
o observed sensitivities to spatial continuity.
matrix to the isotropic variogram model.
gram from the variance (1 for standard normal distribution).
ominant matrix is more readily invertable.

ght hand side matrix to calculate the simple kriging weights.


d with the weights and covariances.
itional cumulative distribution function is available.
cal CCDF

0 0.00 1.00 2.00 3.00

Value

You might also like