GSLIB Sequential indicator simulation
To run this function: Choose Block model > Estimation > GSLIB > Sequential indicator simulation, or...
Files Tab
Input
Location, ID number, String range, D-field
You must complete each of these inputs to describe the string file or files from which the sample data is to be obtained. Enter the Location and Id range of the required string file. Enter the String range and the description field(s) to define the data from the string file which are to be used for estimation. The specified data will be converted from a string file to the geoEAS format input file that GSLIB requires.
Variable type
Select whether the data is categorical (eg. representing rock types) or continuous (eg. representing gold grades).
Minimum
All samples strictly less than this value are ignored. An example value is -1E+21.
Maximum
All samples strictly greater than this value are ignored. An example value is 1E+21.
Output
Estimate
Specify a block model attribute prefix. This will be combined with each realization number to create the attribute names that store the results of each realization. For example, if you specify a prefix of gold and four realizations, the attributes used will be gold1, gold2, gold3 and gold4.
Number of realizations
Specify the number of realizations. Each realization creates a unique set of estimates. An example value is 3.
Seed
A random number seed (a large, odd integer).
Constrain interpolation
If you want to estimate a subset of the block model blocks by use of a block model constraint, check this box.
Debug file
The name of the GSLIB debug file to create.
Debug level
The higher the debugging level, the more output. The normal levels (none, low) summarize the results. Medium and high provide all the kriging matrices and data used for the estimation of every point/block. It is recommended that a high debugging level not be used with a large block model.
Clean up
When the form is Applied, a number of (temporary) files are created, either for input to GSLIB or as output from GSLIB. If you check this box, the temporary files will be cleaned up (deleted).
Cutoffs Tab
Soft data
If you wish to use already transformed indicator values in addition to the sample data, check this box.
Location, ID Number, String range
You must complete each of these inputs to describe the string file or files from which the sample data is to be obtained. Enter the Location and Id range of the required string file. Enter the String range to define the data from the string file which are to be used for estimation. The specified data will be converted from a string file to the geoEAS format input file that GSLIB requires. A D-field containing soft data for each Threshold / Category must be specified in the table below.
Markov Bayes simulation
Check this box to consider Markov Bayes simulation for cokriging with soft indicator data.
Threshold / Category
The cutoff (for continuous data) or category (for categorical data). Example values are 1,2,3,4,6,10.
cdf / pdf
The global cdf value (for continuous data) or global pdf value (for categorical data). This is only used in Simple kriging - see the Simulation Options tab, or if there is no sample data for this cutoff). Example values are 0.32, 0.48, 0.62, 0.73, 0.85, 0.98.
D-field
If you are using already transformed indicator values (soft data) in addition to the sample data, specify which D-field of the soft data string file contains data for this Threshold / Category.
B(z)
If you chose the Markov Bayes simulation option, B(z) calibration values are required. An example value is 0.
Simulation Options Tab
Full IK
If Full IK is selected, then a Full IK is performed. If Median IK is selected, then a Median IK approximation is performed where the kriging weights from the category or cutoff closest to the specified IK cutoff are used for all cutoffs.
IK cutoff
The kriging weights from the category or cutoff closest to this value will be used for kriging. An example value is 2.5.
Kriging Type
Specify which type of kriging to use.
- Simple kriging
- Ordinary kriging
Maximum samples per octant
The maximum number of samples per octant (octant search is not used if this value is left at 0). The octant search ensures that data are taken on all sides of the point being estimated. This is particularly important with drillhole data. An octant search ensures that data is sourced from more than one drillhole.
Maximum original data
The maximum number of original data to use for the simulation of a node. An example value is 4.
Number of simulated nodes to use
The maximum number of previously simulated nodes to use for the simulation of another node. An example value is 8.
Maximum soft indicator nodes
The maximum number of soft data (at node locations) to use for the simulation of a node. An example value is 12.
Multiple grid search
Check this box to perform a multiple grid simulation. Leave the box unchecked to perform a standard spiral search of previously simulated nodes.
Number of multiple grid refinements
The number of multiple grid refinements to consider. An example value is 3.
Assign data to nodes
If this box is unchecked, the data and previously simulated grid nodes are searched separately: the data are searched with a super block search and the previously simulated nodes are searched with a spiral search. If the box is checked, the data are relocated to grid nodes and a spiral search is used.
Search radius
The search radius in the maximum horizontal direction (hMax), the minimum horizontal direction (hMin) and the vertical direction (Vert). See Ellipsoid Definition for more detail on specifying these radii. An example value is 20.
Search angles
The angle parameters that describe the orientation of the search ellipsoid. See Ellipsoid Definition for more detail on specifying these angles. An example value is 0.
Size of covariance lookup table
Enter the X, Y and Z dimensions of the covariance lookup table. Example values are X 8 Y 8 Z 11.
Variograms Tab
You must specify a variogram for each Threshold / Category.
An acceptable variogram model for indicator kriging in GSLIB consists of a nugget effect and any positive linear combination of standard variogram models: spherical, exponential, power or hole effect.
Nugget
Specify the nugget constant for the variogram. An example value is 0.
Type
Select the type of structure represented by this line of the table. The structure types supported by GSLIB for Indicator Kriging are:
- Spherical
- Exponential
- Power
- Hole effect
Cc
Specify the c parameter for this structure. An example value is 0.97.
hMax
An example value is 42.779.
hMin
An example value is 42.779.
Vert
The maximum horizontal range, minimum horizontal range and vertical range. See Ellipsoid Definition for more detail on specifying these values. An example value is 42.779.
Angle1
An example value is 0.
Angle2
An example value is 0.
Angle3
Specify the angles defining the geometric anisotropy. See Ellipsoid Definition for more detail on specifying these angles. An example value is 0.
Post Processing Tab
Post processing is possible for simulated data.
Output type
Select whether to compute:
E-type estimate and conditional variance, i.e., the mean value of the conditional distribution and conditional variance of the conditional distribution.
Estimate
Specify a block model attribute to store the estimate.
Conditional variance
Specify a block model attribute to store the variance.
Probability and means above and below ... a fixed threshold,
Threshold
The Threshold of interest.
probability > Threshold
Specify a block model attribute to store the probability of exceeding the threshold.
mean > Threshold
Specify a block model attribute to store the mean value above the threshold.
mean < Threshold
Specify a block model attribute to store the mean value below the threshold.
Z percentile corresponding to ... a fixed conditional cumulative distribution function (cdf) value,
cdf
The cdf value of interest.
z value
Specify a block model attribute to store the z value.
Symmetric probability interval with width...
Total width
The total width of the interval.
Lower, Upper
Specify block model attributes to store the symmetrical probability interval.
Indicator extrapolation
Min data value
The minimum estimated value for any realization. This is most commonly used for linear interpolation, but can be used for other interpolation types. An example value is 0.
Lower tail
Select the method of extrapolation to the lower limit.
w
Specify the power model parameter.An example value is 0.
Middle
Select the method of interpolation within the distribution.
w
Specify the power model parameter. An example value is 1.
Max data value
The maximum estimated value for any realization. This is most commonly used for linear interpolation, but can be used for other interpolation types. An example value is 0.
Upper tail
Select the method of extrapolation to the upper limit.
w
Specify the power model parameter. An example value is 1.
Tabulated values
If table lookup values are required by one of the interpolation / extrapolation options, a file of the tabulated values is required.
File
The name of the geoEAS formatted file of z data providing details between the IK thresholds.
Variable
The column for the values that will be used for the global distribution. An example value is 3.
Weight
The column for the (declustering) weights that will be used for the global distribution. An example value is 0.
Results
Press Cancel to cancel the function or Apply to invoke GSLIB.
The flow of data is:
Sequential indicator simulation with GSLIB is a two step process:
GSLIB (sisim) requires a sample (.dat) file, a parameter (.par) file and (optionally) a (.dat) file of "soft" sample data.
The sample and soft (.dat) file(s) are created from the specified string file(s) and the form input is used to create the parameter (.par) file.
The file names are created from the location of the sample file and the relevant file name extension.
GSLIB (sisim) creates two files, a (.out) file of simulated results and a (.dbg) file with information about the simulated results. The simulated results are automatically loaded into the specified block model attribute(s). The GSLIB screen output is written to the message window prefixed with:
GSLIB>
After GSLIB (sisim) has finished kriging, further post processing can be performed by GSLIB (postsim).
GSLIB (postsim) requires a GSLIB (sisim) output (.out) file and a parameter (.par) file.
The GSLIB (sisim) output (.out) file is created by GSLIB (sisim) and the form input is used to create the parameter (.par) file.
GSLIB (postsim) creates one file, a (.psm) file of results. The results are automatically loaded into the specified block model attribute(s). The GSLIB screen output is written to the message window prefixed with:
GSLIB>
If you selected Clean up on the Files tab, the temporary files will be deleted after the results have been loaded into the block model.
Messages
Removing GSLIB intermediate files.
You selected Clean up on the Files tab, so the temporary files created for or by GSLIB are being removed.
Error loading file filename0.str
There was an error loading the specified string file.
You must select a model first.
You must open a block model before you can use this function.
File c:\Program Files\Gslib90\sisim.exe not found
File c:\Program Files\Gslib90\postsim.exe not found
sisim.exe or postsim.exe could not be found in the specified directory. Make sure the gslib90 user option is pointing to the directory containing sisim.exe.
Invoking "c:\Program Files\Gslib90\sisim.exe" sample.par ...
Invoking "c:\Program Files\Gslib90\postsim.exe" sample2.par ...
sisim.exe or postsim.exe is being invoked with the specified parameter file.
You have terminated "c:\Program Files\Gslib90\sisim.exe" sample.par. Results were not loaded.
You have terminated "c:\Program Files\Gslib90\postsim.exe" sample2.par. Results were not loaded.
You chose to abort the simulation process. This terminated the executable using the progress feedback cancel button before the results were loaded.
"c:\Program Files\Gslib90\sisim.exe" sample.par terminated prematurely. Results were not loaded.
"c:\Program Files\Gslib90\postsim.exe" sample2.par terminated prematurely. Results were not loaded.
sisim or postsim crashed or otherwise ended prematurely (for example, it was terminated outside of Surpac by using Task Manager) before the results were loaded.
Error writing file "sample.dat". Error writing file "sample.par". Error writing file "sample2.par".
An error occurred while writing the specified file. Check that there is enough room to write the file.
Error opening file sample.dat Error opening file sample.par Error opening file sample2.par
An error occurred while trying to open the specified file for writing. Check that the file is not read-only or in use by another process.
Error opening file sample.out Error opening file sample.psm
An error occurred while trying to open the specified file for reading. Check that sisim.exe did not end prematurely and that the file is not in use by another process.
Created sample.dat. Created sample.par. Created sample2.par.
The specified file was successfully created.
You have terminated the loading of results. Results were only partially loaded.
You chose to abort the block model loading process.
Number of blocks estimated 4096
The number of blocks actually informed is reported.