Experimental Variography and Variogram Models

After data analysis and domaining comes the key part of the contaminant estimation workflow: variography. A theoretical variogram model is proposed and tested using experimental variography. Kriging, nearest neighbour, inverse distance weighting or RBF estimators are defined and potentially combined.

Variography is the analysis of spatial variability of values within a region. Understanding how sample values relate to each other in space is a vital step in informing estimates in a contaminant model. A variogram is used to quantify this spatial variability between samples.

In estimation, the variogram is used for:

• Selecting appropriate sample weighting in Kriging and RBF estimators to produce the best possible estimate at a given location
• Calculating the estimators’ associated quality and diagnostic statistics

In Leapfrog Works contaminant estimations, variograms are created and modified using the Spatial Models folder. To create a new variogram model, right-click on the Spatial Models folder and select New Variogram Model.

You can define as many variograms as you wish; when you define an estimator that uses a variogram model, you can select from those available in the contaminant estimation’s Spatial Models folder:

Working with variograms is an iterative process. The rest of this topic provides an overview of the Variogram Model window and describes how to use the different tools available for working with variograms. This topic is divided into:

• The Variogram Model Window
• Variogram Model Controls
• Experimental Controls

See The Ellipsoid Widget for information that is useful in working in the Variogram Model window.

The Variogram Model Window

The Variogram Model window is divided into three parts:

	Variogram model controls for adjusting the variogram model type, trend and orientation. See Variogram Model Controls below.
	Graphs plotting the selected variogram
	Experimental controls for verifying the theoretical variogram model. See Experimental Controls below.

When you edit a variogram model in the Variogram Model window, an ellipsoid widget is automatically added to the scene. The ellipsoid widget helps you to visualise the variogram in 3D, which is useful in setting variogram rotation and ranges and in defining search neighbourhoods.

If the Variogram Model window is docked as a tab, you can tear the window off. As a separate window, you can move and resize the window so you can see the ellipsoid change in the 3D scene while you make adjustments to the model settings. The detached window can be docked again by dragging the tab back alongside the other tabs, as described in Organising Your Workspace.

There are two ways to save the graph for use in another application:

• Click the Export button () and select Export graph image from the options to export the graph as a PDF, PNG or SVG file.
• Click the Copy graph button () to copy the graph to the clipboard. You can then paste it into another application.

Additionally, you can save variogram parameter data by selecting a theoretical variogram from the list, clicking the Export button () and selecting Export data from the options. Note that Export data is not available for the experimental variogram or the Radial Plot.

There are two buttons for refreshing the graphs when changes are made:

• When Auto Refresh is enabled, recalculations will be carried out each time you change a variogram value. This can produce a brief lag.
• When auto refresh is disabled, you can click the Refresh button () whenever you want the graphs to be updated. This is the best option to use when working with a large dataset.

If auto refresh is disabled and values have been changed without the graphs being updated, the chart will turn grey and a reminder will be displayed over it:

Variogram Model Controls

The Variogram Model controls adjust the variogram model type, trend and orientation. The graph will change as model parameters are adjusted. The ellipsoid in the scene will also reflect the changes you make to the Variogram Model.

Multi-structured variogram models are supported, with provision for nugget plus two additional structures.

The Nugget represents a local anomaly in values, one that is substantially different from what would be predicted at that point based on the surrounding data. Increasing the value of Nugget effectively places more emphasis on the average values of surrounding samples and less on the actual data point and can be used to reduce noise caused by inaccurately measured samples.

Each additional structure has settings for the component Sill and the normalised sill, labelled Norm. sill, model Type, Alpha (if the model type is spheroidal), and the Major, Semi-Major and Minor ellipsoid ranges.

The Sill defines the upper limit of the model, the distance where there ceases to be any correlation between values. The Sill can be set for the Nugget, Structure 1 and Structure 2. A spherical variogram reaches the sill at the range and stays there for increasing distances beyond the range. A spheroidal variogram approaches the sill near the range, and approaches it asymptotically for increasing distances beyond the range. The distinction is insignificant. A linear model has no sill in the traditional sense, but along with the ellipsoid ranges, the sill sets the slope of the model. The two parameters sill and range are used instead of a single gradient parameter to permit switching between linear, spherical and spheroidal interpolant functions without also manipulating these settings.

The Norm. sill represents the same information as the Sill, but proportionally scaled to a range between 0 and 1, where 1 represents the data Variance. As you select the radio buttons for Sill and Norm. sill, the Y-axis scale on the displayed chart will change to correspond to your selection.

The Total sill is the sum of the component sills for both the data sills and the normalised sills.

The Variance is calculated automatically from the data and shows the magnitude of the variance for the dataset.

Linear, Spherical and Spheroidal Model Options

The model Type provides three options: Linear, Spherical and Spheroidal.

Linear provides a general purpose multi-scale option, but with the assumption that samples can have an influence at any range. Spherical and Spheroidal options are suitable for modelling most contaminants, as there is a finite range beyond which the influence of the data falls to zero. A Spherical model has a range beyond which the value is the constant sill. A Spheroidal model flattens out when the distance from the sample data point is greater than the range. At the range, the function value is 96% of the sill with no nugget, and beyond the range the function asymptotically approaches the sill. Each axis of the variogram ellipsoid has its own range, adjusted using Major, Semi-Major and Minor. The ranges are colour-coded in the variogram model plot.

Alpha is only available when the model Type is Spheroidal. The Alpha constant determines how steeply the interpolant rises towards the Sill. A low Alpha value will produce a variogram that rises more steeply than a high Alpha value. A high Alpha value gives points at intermediate distances more weighting, compared to lower Alpha values. An Alpha of 9 provides the curve that is closest in shape to a spherical variogram. In ideal situations, it would probably be the first choice; however, high Alpha values require more computation and processing time, as more complex approximation calculations are required. A smaller value for Alpha will result in shorter times to evaluate the variogram.

Normalised Y Axis

A key use of copying a variogram is to apply it to another correlated contaminant. This would typically be accomplished, having determined the appropriate variogram using the values for one contaminant, by copying the domained estimation and changing the Numeric values field to a different contaminant. However, the absolute values for the nugget and sills for each structure would be completely inappropriate for the new contaminant; while we desire the shape to be the same, the values for different compounds will inevitably be different. To make this work correctly, the variogram is normalised or standardised, rescaling the range of Y-axis values to between 0 and 1, where 1 is equivalent to the data Variance. This makes the variogram information portable between domained estimations. This is performed automatically, requiring no intervention on your part.

You can freely switch between the Sill and Norm. sill options. The selection only changes the Y-axis scale on the displayed chart.

As you adjust the Nugget or the Structure 1 or Structure 2 sill values, the Total sill will change both for the absolute Sill values and the Norm. sill values. The Norm. sill total may end up being something other than 1.0. This is expected, as the value reference for the normalised scale uses the data Variance for 1.0, not the total sill.

Note that the data Variance is recorded in the Y axis label so the chart scale is always meaningful, including when the chart is exported.

Direction

The trend Direction fields set the orientation of the variogram ellipsoid. Adjust the ellipsoid axes orientation using the Dip, Dip Azimuth and Pitch fields.

The Set From Plane button sets the trend orientation of the model ellipsoid based upon the current settings of the moving plane.

The View Ellipsoid button adds a 3D ellipsoid widget visualisation to the scene, in case it has been deleted from the scene since the variogram model window was opened for editing.

Experimental Controls

A variogram model can be verified through the use of the experimental variography tools that use sample data. Use these to find the directions of maximum, intermediate and minimum continuity.

The variogram displayed in the chart is selected from the variograms listed in the panel in the top left corner of the window.

Besides the theoretical Variogram Model, a set of experimental Axis Aligned Variograms are available, including a Radial Plot variogram and a variogram for each axis of the variogram ellipsoid.

Click on one of these experimental variograms to select it and display its parameters. The displayed graph will change to match this selection. For example, here the graph and settings for the combined Axis Aligned Variograms are displayed:

Selecting the Semi-major Axis variogram changes the chart and settings displayed:

Adjust the model variogram parameters to see the effect different parameters have when applied to the actual data.

Experimental Variogram Parameters

The experimental variogram controls along the left side of the window define the search space, define the orientation for custom variograms and change how the variograms are displayed.

Capping

Data Capping fields limit the values of the Lower bound and Upper bound for the data as specified. This is not a filter that discards these points, but values below or above the caps are treated as if the value was the lower or upper bound.

Defining the Search Space

The first set of parameters controls the search space.

• Lag distance controls the size of the lag bins. The first bin will be one quarter of the size of the Lag distance. Experimental variograms generally measure lags as distances along a direction vector, though downhole variograms measure lags as distances along the boreholes.
• Lag tolerance allows for the reality that data pairs are rarely the same distance apart. The data is scanned and pairs are assembled after applying a Lag tolerance to the Lag distance. If the Lag tolerance Auto box is ticked, Leapfrog Works defaults to using a Lag tolerance of half the Lag distance. Controlling the Lag tolerance explicitly allows you to test the sensitivity of the variogram. Typically, a larger value will be used for sparse datasets and a smaller value for dense datasets.

• Number of lags constrains the number of lag bins in the search space.
• The In Plane and Off Plane Angle tolerance and Bandwidth settings define the search shape, and the effects of these settings are discussed in more detail below.

The search shape usually approximates a right rectangular pyramid on a rectangular parallelepiped. The pyramid and parallelepiped will be square if the In Plane Angle tolerance and Bandwidth are used without defining Off Plane values. Using the Off Plane Angle tolerance and Bandwidth fields will make the shape rectangular. The plane being referred to is the major-intermediate axis plane of the variogram ellipsoid, the same plane used for the radial plot. The Angle tolerance is the angle either side of a direction vector from the data point origin. Once the sides of the pyramid defined by the angle tolerances extend out to the limits specified by the Bandwidth, the search neighbourhood is constrained to the bandwidth dimensions.

The search shape becomes a more complex “carpenter’s pencil” shape when a wide In Plane angle and a narrow Bandwidth are defined along with a narrow Off Plane angle and a wide Bandwidth, or vice versa. This rendering should assist in visualising the shape; the major axis is shown in red, the semi-major axis is shown in green, and the orthogonal minor axis in a transparent blue:

Note that although not shown here, the outer ends of the search shapes are not flat, but rounded, being defined by the surface of a sphere with a radius of the maximum distance defined by the number of lags and their size.

Radial Plot Parameters

Radial Plot has parameters specifically for the radial plots in Axis Aligned Variograms.

Increasing the Radial divisions slices the space into a larger number of sectors, with each block in the radial plot covering a smaller arc of the compass. As a result, each block has a smaller volume; this also means that the amount of data in each block is reduced. Because the bandwidth angle above and below the major-intermediate plane matches the angle used to slice the plot into its sectors, increasing the Radial divisions also reduces the number of data points used above and below the plane. Using a smaller number of Radial divisions will be faster. If you increase the number of divisions, you may want to turn off Auto Refresh Graphs first and click Refresh graphs afterwards. Experiment with the number of divisions and choose the lowest number of Radial divisions that helps you gain the best understanding of continuity; this will maximise the data that falls in each division.

Mean Value will display radial plot bins coloured to indicate the mean value for each bin. Contour will display a plot showing lines of equal value. Fill shades the chart between the contour lines.

Changing Axis Limits

The Axis Limits settings control the chart scaling.

X axis limits and Y axis limits control the ranges for the X-axis and Y-axis and can effectively be used to zoom the chart. You can directly control these by manipulating the axes with your mouse. Click and drag an axis to increase or decrease the maximum limit of the axis. Right-click and drag an axis to reposition the axis so the minimum value on the axis is not zero. Double-click the axis to reset the axis minimum and maximum range to the default values. Fit Graph will auto-fit the graph to the available data.

Changing Variogram Display

The Display settings change how the variogram is displayed:

• Show model variogram plots the variogram on the chart.
• Show pair count annotates the chart with the count of data point pairs.
• Tick the Scale points to pair size checkbox and the chart will display the plotted points with dots that scale to the size of the pair count. Choose Actual so the dots represent the pair size proportionally and Inverse to see larger dots for lower pair counts.
• Tick Moving avg for a rolling mean variogram value, drawn on the variogram plots in orange. Next to it, select a value for the window size relative to the Lag distance.The window size will be the given number times the lag distance. This window, centred on x, will be used to plot y, the average of all the variogram values falling inside the window. The lag multiplier field has a useful tooltip reminder; hover your mouse pointer over the field to see the tooltip.

Axis Aligned Variograms

The Axis Aligned Variograms are useful for determining the direction of maximum, intermediate and minimum continuity. When you select the Axis Aligned Variograms option, all four variograms are displayed in the chart:

The radial plot and each of the axes variograms can be viewed in greater detail by clicking on them in the variogram tree:

In the axes variogram plots, you can click-and-drag the plus-shaped dragger handle to adjust the range and the sill. A triangular dragger handle adjusts the nugget. A solid line shows the model variogram, and dotted horizontal lines show the Total Sill and, if Sill is chosen instead of Norm. sill, the variance.

In the radial plot, you can click-and-drag the axes arrows to adjust the pitch setting, between the values 0 and 180 degrees. Each bin in the mean value radial plot shows the mean semi-variogram value for pairs of points binned by direction and distance. When the contour plot is selected, lines follow the points of equal value.

 

Got a question? Visit the Seequent forums or Seequent support