Analysis of different gridding
methods using “Surfer7”
Michigan
Technological
University
presentation by
Fehmi Kamcili
February 10, 2001
Fig.: Shaded Relief and Contour map
Analysis of different gridding
methods using “Surfer7”
Table of content
• Introduction
• Grid Files
• Grid Methods
1. Inverse Distance to a Power; 2. Kriging;
3. Minimum Curvature; 4. Modified Shepard’s Method;
5. Natural Neighbor; 6. Nearest Neighbor;
7. Polynomial Regression; 8. Radial Basis Function;
9. Triangulation w/ linear Interpolation
• Best Results
• Conclusion
Analysis of different gridding
methods using “Surfer7”
Introduction
• Surfer 7; Feb. 2000
Surface mapping system
Golden Software, Inc.
• Differences in created
grid-files visualized by
overlaying wireframe
Fig.: Wireframe and Contour Map and contour maps
• Vernon/Rosebuch oilfields in Isabella county, Michigan
with 245 well Data-Points (Lat-Long of Top Dundee)
Analysis of different gridding
methods using “Surfer7”
Grid Files
• gridding by specifing
source file
-Spreadsheets (Excel)
-manually in Surfer
(Worksheet)
• gridding method
• accuracy of grid
• faults and breaklines
Fig.: Contour and Post Map • creates “file.grd”
Analysis of different gridding
methods using “Surfer7”
Grid Methods
• Inv. Distance to power
• Kriging
• Minimum curvature
• Modified shepard’s
• Natural neighbor
• Nearest neighbor
• Polynomial regression
• Radial basis function
Fig.: Wireframe Map
• Triangulation w/ linear
interpolation
Analysis of different gridding
methods using “Surfer7”
Fig.: Power:
0.00001
-2905.1933
-2905.1933
-2905.1933
Inverse Distance to a
-2905.1933
-2905.1934
-2905.1934
Power
• weighting average
-2905.1934
-2905.1934
-2905.1934
-2905.1935
-2905.1935
-2905.1935
-2905.1935
interpolator
• Power parameter
-2905.1935
-2905.1936
-2860
between 1E-38 & 38
-2870
-2880 “0” = planar surface;
-2890
-2900
-2910
great weighting power
-2920
-2930
= less effect on points
-2940
-2950 far from the grid node
-2960
-2970
during interpolation
Fig.: Power: 2
Analysis of different gridding
methods using “Surfer7”
Fig.:
Power: 8
-2860
-2870
Inverse Distance to a
-2880
-2890
-2900
Power
-2910
-2920
-2930
• exact or smoothing
-2940
-2950
-2960
interpolator
• generate "bull's-eyes”
-2970
-2980
-2875
-2880
-smoothing reduce
-2885
-2890
-2895
this effect
-2900
• very fast method for
-2905
-2910
-2915
-2920
-2925
-2930
-2935
gridding
-2940
-2945
-2950
-2955
-till 500 Datapoints
Fig.: with -2960
-2965
smoothing: 0.01
Analysis of different gridding
methods using “Surfer7”
Kriging
-2860
-2870
-2880
-2890
• express trends
-2900
-2910
-2920
suggested in your data
-2930
-2940 • Point or Block
-2950
-2960
-2970
-Block is using
-2980
average values-
Fig.: Point Kriging smoother; not perfect
• by entering path & file • specify & add as
name production estimated many variogram-
standard deviation grid components as wished
Analysis of different gridding
methods using “Surfer7”
Minimum Curvature
• smooth but not exact
• recalculation of grid
-2860
-2880
-2900
node values until
-2920
-2940
-2960
-2980
-3000
-3020
reached less of max.
-3040
-3060
-3080
Residual value, or
-3100
-3120
-3140
max. Iteration
• Set Internal and
-3160
Fig.: default Boundary Tension
• Relaxiation Factor
Analysis of different gridding
methods using “Surfer7”
Fig.: 1 -2800 Modified Shepard’s
-2850
-2900
-2950
Method
-3000
-3050
• Like IDP Method
• exact or smoothing
-3100
-3150
-3200
-3250
-3300
• Weighting and
-2900
Quadratic Neighbors
-2950
-3000
parameters specifies
-3050
-3100
size (number) of local
-3150
-3200
neighborhood
-3250
-3300
• Fig.1: Q13/W19;
Fig.: 2
-3350
-3400
Fig.2: Q40/W60
Analysis of different gridding
methods using “Surfer7”
Natural Neighbor
• Natural Neighbor
-2860
-2870
interpolation
-2880
-2890
algorithm uses a
-2900
-2910 weighted average of
-2920
-2930
-2940
neighboring
-2950
-2960
observations, where
the weights are
-2970
-2980
proportional to new
Fig.: default
added polygons
Analysis of different gridding
methods using “Surfer7”
Nearest Neighbor
• assigns value of the
-2860
-2870
nearest point to each
-2880
-2890
grid node
• useful when data are
-2900
-2910
-2920
-2930
-2940 already evenly
-2950
-2960
-2970
spaced
• this method is
-2980
Fig.: default effective for filling
missing values
Analysis of different gridding
methods using “Surfer7”
Polynomial
Fig.: 1 -2875
-2880
-2885
-2890
Regression
-2895
-2900 • used to define large-
-2905
-2910
-2915
scale trends &
-2920
-2925 patterns
-2930
-2935
• not real interpolator
-2840
-2850 (does not predict
-2860
-2870
-2880
unknown Z values)
• Fig.1: Simple planar
-2890
-2900
-2910
-2920
-2930 surface; Fig.2: Bi-
-2940
-2950
linear saddle
Fig.: 2 -2960
Analysis of different gridding
methods using “Surfer7”
Radial Basis
Function
1. Inverse Multiquadric
-2870
-2880
-2890
-2900
2. Multilog
-2910
-2920
-2930
-2940
3. Multiquadric
-2950
-2960
-2970
-2980
4. Natural Cubic Spline
-2990
-3000
-3010
-3020
-3030
5. Thin plate Spline
• 1.; 2.; 4.- error
-3040
-3050
Fig.: 5. Thin plate spline • 5.- good; 3.- best
Analysis of different gridding
methods using “Surfer7”
Radial Basis
Function
-2870
• all exact
-2875
-2880
-2885
-2890
interpolators
-2895
-2900
-2905
-2910
+ smoothing factor
• = variogram in K.
-2915
-2920
-2925
-2930
-2935
-2940
-2945
-2950
(mathematically
-2955
-2960
-2965
-2970
specifies spatial
-2975
variability of data set
Fig.: 3. Multiquadric & resulting grid file
Analysis of different gridding
methods using “Surfer7”
Triangulation w/
linear Interpolation
• creates triangles by
-2860
-2870
-2880
drawing lines
-2890
-2900
between data points
• exact interpolator
-2910
-2920
-2930
• for evenly distributed
-2940
-2950
-2960
-2970
-2980
data over grid area
-sparse areas result in
Fig.: default distinct triangular
facets on map
Analysis of different gridding
methods using “Surfer7”
Best Results
Kriging is exact, has many options & modifications,
but needs knowledge -2860
-2870
-2880
-2890
-2900
-2910
-2920
-2930
-2940
-2950
-2960
-2970
Fig.: Point Kriging -2980
Analysis of different gridding
methods using “Surfer7”
Best Results
Radial Basis Function is exact, shows nice views
& is uncomplicated
-2870
-2875
-2880
-2885
-2890
-2895
-2900
-2905
-2910
-2915
-2920
-2925
-2930
-2935
-2940
-2945
-2950
-2955
Fig.: Radial Basic -2960
-2965
Function; Multiquadric -2970
-2975
Analysis of different gridding
methods using “Surfer7”
Conclusion
• I did not consider all aspects & details
• Surfer is very powerful for the 3D-Visualization
• Surfer 7 works:
-fast & without consuming much disk space
-uncomplicated with Object manager
• All processes (gridding, mapping) can be automated
with writing programs in Visual Basic
• Help content is very useful & describes also background
information
