Case Studies
In order to demonstrate the FPG techniques, four datasets with high-skewed distributions had been chosen. They are: the synthetic data of GSLIB, the Amapari gold deposit in Brazil and the U and V Walker Lake datasets.
For a FPG application, a previous analysis of frequencies is carried out, in order to determine if there is any inflexion in the declustered cumulative frequency curve. In the first case study, the V variable of Walker Lake, there is no inflexion point and the FPG transformation is performed to all data and the estimated values are back transformed to grades.
In the second case study, the gold deposit of Amapari, there are evidences of two populations, and an inflexion point could be assigned. This point is treated as limiting threshold, determining that only the grades above it should be estimated by the FPG technique.
In the Amapari gold deposit, the first population, with low grades, is estimated by ordinary kriging of the original grades. The second population, with high grades, has their values processed by FPG, estimated by ordinary kriging and back transformed to the original grades. The final grade of each block is obtained by weighting, using the proportions given by indicator kriging.
For the sake of comparison, the grade estimation will also be made by conventional techniques, including ordinary kriging of the original values and ordinary kriging after capping the extreme values. Results obtained by FPG and other techniques help in defining the impact on overall statistical parameters of the deposit. For Walker Lake the estimated values are compared against the reference values, derived from the exhaustive dataset.
FPG was used before for the synthetic data of GSLIB (Armony, 2000) and for the U variable of Walker Lake (Machado et al, 2011), both case studies provided satisfactory results.
A two dimensional set of 2500 data are uniformly distributed in a square, and only 140 highly skewed data are available for estimation. The goal is to compare statistical results of FPG using only the 140 data, with the 2500 reference statistical results as shown in Table 1. For FPG estimation, conditional simulation was used to fulfil the square. Only one population was found.
| Variable | Data | Mean | Std. Dev. | Variance |
|---|---|---|---|---|
| GSLIB (Reference) | 2500 | 2.58 | 5.15 | 26.53 |
| FPG simulation | 140 | 2.87 | 5.16 | 26.65 |
The variograms using the samples and FPG are shown in figures 1 and 2.

Fig. 1
Variogram shape of the 140 samples of GSLIB synthetic data.

Fig. 2
Variogram shape of the FPG transform of the 140 GSLIB synthetic data.
Conclusion
At this stage FPG provided better results on Walker Lake both globally and locally.
Presenting the dataset
The Walker Lake dataset (Isaaks and Srivastava, 1989) is derived from a digital elevation model. It consists of the categorical variable T and the numerical variables U and V. The U and V variables will be used in this study, which are preferentially sampled at high-grade zones (Figures 1 and 2) and both present only one population. The exhaustive set has 78,000 measured values.

Fig. 1
Sample locations for U variable.

Fig. 2
Sample locations for V variable.
The statistical parameters for U and V are presented in Table 1.
Table 1 – Statistical parameters for variables U and V.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| U | 275 | 0 | 5190.1 | 604.08 | 766.01 | 586769.89 | 1.27 |
| U (declustered) | 275 | 0 | 5190.1 | 555.00 | 684.79 | 468932.61 | 1.23 |
| V | 470 | 0 | 1528.1 | 435.30 | 299.56 | 89738.06 | 0.69 |
| V (declustered) | 470 | 0 | 1528.1 | 408.53 | 293.23 | 85985.06 | 0.72 |
Results and discussions
Different methods have been tested and their results were compared with a reference model – U (Reference) and V (Reference) – using block size of 5m 5m. Tables 2 and 3 show the statistical parameters for U and V estimates respectively. U (OK) and V (OK) are the results of ordinary kriging of the original values; U (capping Q97) and V (capping Q90) are the results of ordinary kriging after capping the extreme values above the 90th percentile of the distribution (V = 818.9). The results of median indicator kriging were post processed with postik from GSLIB (Deutsch and Journel, 1998): U (E-type 1.5) and V (E-type 1.5) are the results with a hyperbolic upper tail extrapolation of parameter ω=1.5; and U (E-type 3.0) and V (E-type 3.0) are the results with a hyperbolic upper tail extrapolation of parameter ω=3.0.
Table 2 – Statistical parameters for the U estimates.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| U (Reference) | 3120 | 0.00 | 2844.55 | 266.04 | 358.81 | 128744.62 | 1.35 |
| U (OK) | 3101 | 52.51 | 2500.58 | 542.88 | 314.26 | 98759.35 | 0.58 |
| U (capping Q97) | 3120 | 61.86 | 1468.65 | 543.20 | 260.68 | 67954.06 | 0.48 |
| U (E-type 1.5) | 3120 | 7.20 | 3058.02 | 464.56 | 418.14 | 174841.06 | 0.90 |
| U (E-type 3.0) | 3120 | 7.20 | 2316.51 | 428.82 | 343.55 | 118026.60 | 0.80 |
| U (FPG) | 3042 | 4.84 | 1590.44 | 255.22 | 225.40 | 50805.16 | 0.88 |
Table 3 – Statistical parameters for the V estimates.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| U (Reference) | 3120 | 0.00 | 1378.12 | 277.98 | 228.66 | 52285.40 | 0.82 |
| U (OK) | 3120 | 3.81 | 1052.34 | 300.44 | 171.73 | 29491.19 | 0.57 |
| U (capping Q97) | 3120 | 8.70 | 814.56 | 304.92 | 153.93 | 23694.44 | 0.50 |
| U (E-type 1.5) | 3120 | 20.74 | 1295.85 | 302.10 | 217.77 | 47423.77 | 0.72 |
| U (E-type 3.0) | 3120 | 20.74 | 1125.01 | 296.64 | 202.58 | 41038.66 | 0.68 |
| U (FPG) | 3120 | 2.26 | 922.68 | 261.84 | 181.45 | 32924.10 | 0.69 |
For both variables the global average of deposit was best approximated by the FPG technique, especially for the U variable. The correlation coefficients are used as a comparative parameter, as it can be seen in Tables 4 and 5. FPG produces better results, and this improvement is evident when comparing the results with ordinary kriging of the original values.
Table 4 – Correlation coefficients for U.
| U (FPG) | U (OK) | U (capping Q97) | U (E-type 1.5) | U (E-type 3.0) | |
|---|---|---|---|---|---|
| Reference | 0.510 | 0.383 | 0.375 | 0.501 | 0.502 |
Table 5 – Correlation coefficients for V.
| U (FPG) | U (OK) | U (capping Q97) | U (E-type 1.5) | U (E-type 3.0) | |
|---|---|---|---|---|---|
| Reference | 0.867 | 0.832 | 0.808 | 0.855 | 0.859 |
Conclusion
At this stage FPG provided better results on Walker Lake both globally and locally.
Presenting the dataset
The area containing the deposit is located in the state of Amapá, Brazil. There are three types of mineralisation: sulphides (primary mineralisation), which occur at depth, in fresh rock; saprolite, made of its in situ weathered and oxidised portions; and the colluvial ores, overlaying the saprolite ore as a blanket. The colluvia domain of the deposit was chosen to develop this study, in which the grade distribution reflects the distribution observed in the primary mineralisation. A division of colluvia domain into northern and southern portions was performed, whose statistical parameters are presented in Table 1 and Table 2, respectively.
Table 1 – Statistical parameters of AU in the NORTHERN portion of the deposit.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| AU, ppm | 1720 | 0.00 | 24.04 | 1.47 | 2.32 | 5.38 | 1.58 |
| AU (declus), ppm | 1720 | 0.00 | 24.04 | 0.90 | 1.51 | 2.28 | 1.67 |
Table 2 – Statistical parameters of AU in the SOUTHERN portion of the deposit.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| AU, ppm | 4127 | 0.00 | 145.30 | 1.50 | 4.26 | 18.15 | 2.84 |
| AU (declus), ppm | 4127 | 0.00 | 145.30 | 0.97 | 2.64 | 6.97 | 2.71 |
Results and discussion
An analysis of the distribution and spatial arrangement shows that two different populations can be recognized: the first ranging from 0.00 ppm to 5.00 ppm and the second, equal or greater than 5.00 ppm, to be treated by FPG using cumulative extensions.
Different methods have been tested and their results were compared. Table 3 and table 4 show the overall statistics of block estimates in northern and southern portions respectively. AU (OK) is the result of ordinary kriging of the original values; AU (FPG) is the result obtained by FPG technique; AU (capping Q95) and AU (capping Q98) are the results of ordinary kriging after capping the extreme values above the 95th and the 98th percentile of the distribution (AU = 5.03 ppm and AU = 9.20 ppm, respectively).
Table 3 – Statistical parameters for AU estimates in the NORTHERN portion of the deposit.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| AU (OK), ppm | 213979 | 0.11 | 9.88 | 0.72 | 0.56 | 0.31 | 0.78 |
| AU (capping Q98), ppm | 213979 | 0.10 | 6.19 | 0.72 | 0.53 | 0.28 | 0.73 |
| AU (capping Q95), ppm | 213979 | 0.10 | 4.24 | 0.70 | 0.46 | 0.21 | 0.66 |
| AU (FPG), ppm | 213979 | 0.08 | 6.61 | 0.59 | 0.45 | 0.20 | 0.76 |
Table 4 – Statistical parameters for AU estimates in the SOUTHERN portion of the deposit.
| Variable | Count | Min. | Max. | Mean | Std. Dev. | Variance | Coef. of Variation |
|---|---|---|---|---|---|---|---|
| AU (OK), ppm | 267061 | 0.10 | 64.44 | 0.81 | 1.11 | 1.23 | 1.37 |
| AU (capping Q98), ppm | 267061 | 0.12 | 7.66 | 0.75 | 0.67 | 0.45 | 0.89 |
| AU (capping Q95), ppm | 267061 | 0.12 | 4.95 | 0.69 | 0.48 | 0.23 | 0.69 |
| AU (FPG), ppm | 267061 | 0.07 | 8.65 | 0.62 | 0.56 | 0.31 | 0.90 |
The method yields similar results to classic kriging in the northern part of the deposit. The results show a greater variation in the south, where the variance is greater.
The correlation coefficient is used as a comparative parameter as it can be seen in Table 5 and Table 6. In this case FPG produces results similar to ordinary kriging with a capping in the extreme values, but there is no way to know which top cut is the best.
Table 5 – Correlation coefficients for Amapari NORTH.
| AU (OK) | AU (capping Q98) | AU (capping Q95) | |
|---|---|---|---|
| AU (FPG) | 0.877 | 0.864 | 0.846 |
Table 6 – Correlation coefficients for Amapari SOUTH.
| AU (OK) | AU (capping Q98) | AU (capping Q95) | |
|---|---|---|---|
| AU (FPG) | 0.716 | 0.861 | 0.865 |
Conclusion
At this stage FPG provided better results on Walker Lake both globally and locally.