Minerals 08 00232 PDF
Minerals 08 00232 PDF
Article
Integrating the Theory of Sampling into
Underground Mine Grade Control Strategies
Simon C. Dominy 1,2, * ID , Hylke J. Glass 1 ID , Louisa O’Connor 2 ID
, Chloe K. Lam 1 ID
,
Saranchimeg Purevgerel 3 and Richard C.A. Minnitt 4 ID
1 Minerals Engineering Research Group, Camborne School of Mines, University of Exeter, Penryn,
Cornwall TR10 9FE, UK; [email protected] (H.J.G.); [email protected] (C.K.L.)
2 Department of Mining and Metallurgical Engineering, Western Australian School of Mines,
Curtin University, Bentley, WA 6102, Australia; [email protected]
3 Department of Mineral and Energy Economics, Western Australian School of Mines, Curtin University,
Murray Street, Perth, WA 6000, Australia; [email protected]
4 School of Mining Engineering, University of the Witwatersrand, Johannesburg 2050, South Africa;
[email protected]
* Correspondence: [email protected]
Received: 7 May 2018; Accepted: 25 May 2018; Published: 29 May 2018
Abstract: Grade control in underground mines aims to deliver quality tonnes to the process plant
via the accurate definition of ore and waste. It comprises a decision-making process including data
collection and interpretation; local estimation; development and mining supervision; ore and waste
destination tracking; and stockpile management. The foundation of any grade control programme is
that of high-quality samples collected in a geological context. The requirement for quality samples has
long been recognised, where they should be representative and fit-for-purpose. Once a sampling error
is introduced, it propagates through all subsequent processes contributing to data uncertainty, which
leads to poor decisions and financial loss. Proper application of the Theory of Sampling reduces
errors during sample collection, preparation, and assaying. To achieve quality, sampling techniques
must minimise delimitation, extraction, and preparation errors. Underground sampling methods
include linear (chip and channel), grab (broken rock), and drill-based samples. Grade control staff
should be well-trained and motivated, and operating staff should understand the critical need for
grade control. Sampling must always be undertaken with a strong focus on safety and alternatives
sought if the risk to humans is high. A quality control/quality assurance programme must be
implemented, particularly when samples contribute to a reserve estimate. This paper assesses grade
control sampling with emphasis on underground gold operations and presents recommendations for
optimal practice through the application of the Theory of Sampling.
Keywords: underground mine grade control; Theory of Sampling; sampling errors; representative
sampling; sampling protocol optimisation; quality assurance/quality control
1. Introduction
1.1. Overview
Sampling is a critical component throughout the mine value chain; it includes the sampling of
both in-situ and broken material for geological (resource and grade control), geoenvironmental, and
geometallurgical purposes [1–7]. Sampling errors generate both monetary and intangible losses [3,8].
The data produced must be fit-for-purpose in order to contribute to Mineral Resources/Ore Reserves
reported in accordance with the 2012 Joint Ore Reserves Committee (JORC) [9] or other codes.
Quality assurance/quality control (QAQC) is critical to maintain data integrity through documented
procedures, sample security, and monitoring of precision, accuracy and contamination. Samples and
their associated assays are key inputs into important decisions throughout the mine value chain.
The Theory of Sampling (TOS) defines sampling errors, where actions may lead to uncertainty
and create a total measurement error [1,10]. Theory of Sampling attempts to break down this error
into a series of contributions along the sampling value chain (e.g., the planning to assay-measurement
process) (Table 1).
Errors are generally additive throughout the sampling value chain and generate both monetary
and intangible losses. The aim is to collect representative samples to accurately describe the material
in question. Sample collection is followed by reduction in both mass and fragment size to provide
a sub-sample for assay. An assay is the quantitative measurement of the concentration (e.g., mass
fraction such as g/t gold) of a metal by a given methodology, for example a 30 g fire assay followed by
measurement of gold using an instrumental method (e.g., inductively coupled plasma spectroscopy).
control process. A best practice grade control programme will provide information on grade and its
Minerals 2018, 8, x FOR PEER REVIEW 3 of 43
relationship to geology (Figure 1).
Grade control samples such as chip or channel samples are relatively cheap to collect, and apart
Grade control samples such as chip or channel samples are relatively cheap to collect, and apart
from the effort involved, easy to take. The potentially large numbers of samples used in estimation
from the effort involved, easy to take. The potentially large numbers of samples used in estimation
smooth out individual assay variability. However, the consequence is that the sampling distributions
smooth out individual assay variability. However, the consequence is that the sampling distributions
drawn from these data may show large numbers of low grade assays, positive disjointed tails, and
drawn from these data may show large numbers of low grade assays, positive disjointed tails, and
high skewness. This makes the estimation process more challenging and generally indicates that the
high skewness. This makes the estimation process more challenging and generally indicates that the
samples are not representative.
samples are not representative.
(A) (B)
Figure 1. Face photographs showing sample lines based on geology: (A) A single chip sample line
Figure 1. Face photographs showing sample lines based on geology: (A) A single chip sample line
divided into a central vein sample with associated hanging and footwall samples; (B) Dual face chip
divided into a central vein sample with associated hanging and footwall samples; (B) Dual face chip
sample lines across a sheeted vein zone.
sample lines across a sheeted vein zone.
with extreme cases of 90% above 100 µm and including 60% about 1000 µm. Coarse particles may
be more clustered and related to high-grade zones. In rare cases, fine gold particles may cluster to
form pseudo-coarse gold masses of up to 500 µm in size. Clustering can have a significant impact on
sampling protocols, particularly during the early stages of sampling broken rock, whether it is bulk
sampling of development rounds, grab sampling of stockpiles, field sampling of RC chippings, or
laboratory sampling of crushed drill core. Once the sample has been pulverised, the additional in-situ
heterogeneity related to clustering is removed, though it may still be relatively high if coarse gold
is present.
2. Theory of Sampling
Sampling Effect on
Acronym Error Type Source of Error Error Definition
Error Sampling
Results from grade heterogeneity of the
broken lot. Of all sampling errors, the FSE
does not cancel out and remains even after
Fundamental FSE Correct Random Characteristics a sampling operation is perfect. Experience
Sampling Error Errors—Precision of the ore type. shows that the total nugget effect can be
(CSE) Generator Relates to artificially high because sample masses are
constitution not optimal.
and
Relates to the error due to the combination
distribution
of grouping and segregation of rock
heterogeneity
Grouping and fragments in the lot. Once rock is broken,
GSE
Segregation there will be segregation of particles at any
scale (e.g., surface stockpile or
laboratory pulp).
Results from an incorrect shape of the
Delimitation DE Sampling
Incorrect volume delimiting a sample.
Systematic equipment and
Sampling Error
Errors—Bias materials Results from the incorrect extraction of a
(ISE)
Generator handling sample. Extraction is only correct when all
Extraction EE
fragments within the delimited volume are
taken into the sample.
Samples should represent a consistent mass
Weighting WE
per unit (e.g., kg/m).
Refers to issues during sample transport
and storage (e.g., mix-up, damage, loss and
alteration), preparation (contamination
Preparation PE and/or losses), and intentional (sabotage
and salting), and unintentional (careless
actions and non-adherence of
protocols) actions.
Relates to all errors during the assay and
analytical process, including issues related
Analytical
Analytical AE - to rock matrix effects, careless actions, and
process
analytical machine maintenance
and calibration.
Its application to underground grade control sampling has been variable, with leading use in the
South African gold sector [19–21]. Some of the resistance to the application of TOS relates to the use of
the Fundamental Sampling Error (FSE) equation and confusion as to its calibration and application.
Theory of Sampling has a wider usage than simply the FSE equation and includes a number of errors
that must be considered to achieve representative samples (Table 2).
Theory of Sampling provides an insight into the causes and magnitude of errors that may occur
during the sampling of particulate materials (e.g., broken rock). Theory of Sampling strictly does not
include sampling of in-situ material such as face samples, as the sample does not have an equiprobable
chance of being collected. However, its application is relevant and helps to understand some of the
challenges of field sampling through analysis of the so-called incorrect sampling errors during the
sample collection process [4,21]. The application of TOS includes the mandatory use of an appropriate
QAQC system, including the use of duplicates. Table 3 shows the cause and effect of the common
sampling errors and provides some solutions.
Minerals 2018, 8, 232 6 of 45
Nugget Effect
Sample Type Error Type Cause Error Effect Solution
Component
In-situ (e.g., linear In-situ Local Poor Larger samples
GNE
samples and core) heterogeneity representativity precision More samples
Optimised protocols
Constitution Poor
FSE SNE Larger samples
Random Heterogeneity precision
More samples
Broken rock (from
rock chips to blasted Constitution Optimised protocols
Poor
material and and Distribution GSE SNE More sample
precision
laboratory pulps) Heterogeneity increments
Poor quality Optimised protocols
DE, EE, WE,
Systematic sampling High bias SNE Efficient training
PE, and AE
protocols Strict QAQC
GNE: geological nugget effect; SNE: sampling nugget effect.
2.4.1. Background
The FSE can be modelled before material is sampled, provided certain characteristics are
determined [1,10,22]. The FSE is applicable to grade control samples once they are collected (e.g., a face
sample is broken rock in a sample bag). For any process where the FSE escalates, there is an associated
loss due to uncertainty, which in terms of grade control relates to ore/waste misclassification. The FSE
equation (“François-Bongarçon modified” version) is given as [22]:
where f = shape factor; g = granulometric factor; dn = nominal material size; c = mineralogical factor;
d95Au = liberation diameter; b = (3 − α), where α is determined experimentally from Duplicate Series
Analysis (DSA) tests [25], or a default value of α = 1.5 is generally applied [22]; MS = sample mass;
and ML = lot mass. The sampling constant (K) is a product of: f·g·c (d95Au )b . Using this formula, it is
Minerals 2018, 8, 232 7 of 45
possible to: (a) calculate the error for a given sample size split from the original, or (b) calculate what
sub-sample size should be used to obtain a specified variance at a given reliability.
• Realisation of gold deportment, in particular the partitioning of gold as free gold, gold in sulphides,
and refractory gold.
• Gold particle size curve(s), including effects of clustering and relationship between gold particle
size and grade (e.g., high grade versus background grade).
• Definition of key FSE equation inputs (d95Au versus dAuclus ) and the sampling constant K.
• Recommendations as to optimum in-situ sample mass requirements.
The direct approach to d95Au may include a combination of methods such as core logging,
underground face observation, optical and automated microscopy, X-ray tomography, and staged
liberation and concentration [7,28–30]. Characterisation studies allow the practitioner to set the ore
sampling expectations across a number of d95Au -grade scenarios. In the production environment,
optimisation at the breakeven cut-off grade is usually the critical point. Where this requires specialist
and potentially costly protocols, then a cost-benefit analysis is required to understand the theoretical
need versus practical reality.
effect has two principal components [13,14,36]: the geological nugget effect (GNE) and the sampling
nugget effect (SNE).
The GNE relates to the microscopic differences in composition throughout the
mineralised zone/domain. It principally refers to the distribution of single grains or clusters
of gold particles distributed through the ore zone [16]. Variability is greatest where there are
small-scale, low continuity structures such as high-grade gold carriers within the main structure or
veinlets within wall rocks [14].
The SNE is directly related to the sampling process: that is, to the size of the field samples
taken; the effectiveness of collection; sample preparation; the size of sub-samples after splitting; and
analysis [13,14]. Throughout the mine value chain, optimised sampling protocols aim to reduce
the total sampling variance thereby also reducing the total nugget variance, skewness of the data
distribution, and number of extreme data values.
Table 4. Distribution of errors across stages of a gold grade control sampling programme. Refer to
Table 1 for the sampling value chain.
1 Duplicate 3Other
Error SVC Stage Sample Type/Activity TOS Errors Component 2 FSE TOS
Error Range Errors
In-situ sampling (e.g., core GNE
GNE, DE, EE, WE, PE
Sampling 2&3 and linear samples) ±20–70%
±16%
Broken rock sampling
(e.g., grab, sludge, and RC FSE, GSE, DE, EE, WE, PE
samples) ±23%
Drying PE
Preparation 4 ±5–20% ±11%
Crushing/grinding PE
Splitting FSE, GSE, DE, EE, PE
Splitting FSE, GSE, DE, EE, PE
Analytical 5 ±1–15% ±8%
Analysis AE
Total error ±20–70% ±21% ±23%
(±20%)
Target error (fine-gold)/coarse-gold ±32%
±40%
1 Potential component error range as determined from duplicate sample (pair) analysis [34,41]; 2 Maximum
recommended FSE distribution across the sampling stages [37]; 3 Maximum recommended other TOS error
distribution across the sampling stages [37]; SVC: sampling value chain; RC: reverse circulation.
assumptions are (3) and (5) above. Assumption (3) is unlikely to be true since the population will have
variable particle sizes and shapes. For assumption (5), the larger nuggets will be rare, but the fine
background gold is likely to be more abundant and disseminated. The Poisson approach defines a
maximum theoretical sample mass based on the largest gold nuggets present.
For example, gold mineralisation at a potential underground breakeven cut-off grade of 3.5 g/t
Au and a d95Au of 1 mm, requires a 250 kg theoretical sample mass to yield ±20% precision at 90%
reliability. If each face is represented by a composite channel sample of 4 kg, then 63 samples are
required to evaluate the ore zone. Such a figure could be used to guide the sampling around a
stope panel. A 2800 t stope panel could be blocked out by raises and drives totalling 120 m. If each
advance was 1.8 m, then this is 70 samples around the panel. Seventy 4 kg samples yields a total mass
of 280 kg, which is above the theoretical optimum.
Given that theoretical sample mass, face separation, and stope panel size are relatively fixed,
the variables that can be changed are individual sample mass and number of samples per face.
An individual field sample may not be representative, particularly in coarse-gold mineralisation.
A theoretical field sample mass of many tonnes could be required to achieve an acceptable precision.
As noted previously, a reasonable strategy is to collect multiple samples across a given domain (e.g.,
stope panel). Each sample may be locally unrepresentative, though by taking numerous spatially
distributed samples the mean grade may be ultimately approached and variability reduced.
Table 5. Example from a South American project of sampling characteristics for vein-stockwork wall
rock mineralisation domains.
Table 6. Underground grade control sampling types and methods. Note, strictly drilling is a linear
sampling method, though it is classified separately here as it is different from the main grade
control methods.
Table 7. Sample type for different stoping methods during the mining production chain.
(A) (B)
Figure 2. (A) Chip-channel sampling. Samples collected from sidewall of shanty drive. Each sample
Figure 2. (A) Chip-channel sampling. Samples collected from sidewall of shanty drive. Each sample
comprises near continuous chips collected across the vein; (B) Hand channel sampling of side walls.
comprises near continuous chips collected across the vein; (B) Hand channel sampling of side walls.
Table 8. Likely sampling errors during the collection of chip or chip-channel samples.
Figure 3. Orientation of linear samples (e.g., chip-channel) to match geological requirements. Sub-vertical
quartz reef in foreground with area of flat veining on the footwall (blue line marks footwall contact).
Sub-vertical reef sampled horizontally (samples 1 and 7), whereas flat footwall veins sampled
vertically (samples 2–6).
Table 8. Likely sampling errors during the collection of chip or chip-channel samples.
Sampling 1 Error
Error Effect
Error Magnitude
(A) (B)
Minerals 2018, 8, 232 13 of 45
Figure 2. (A) Chip-channel sampling. Samples collected from sidewall of shanty drive. Each sample
comprises near continuous chips collected across the vein; (B) Hand channel sampling of side walls.
Figure 3. Orientation of linear samples (e.g., chip-channel) to match geological requirements. Sub-vertical
Figure 3. Orientation of linear samples (e.g., chip-channel) to match geological requirements.
quartz reef in foreground with area of flat veining on the footwall (blue line marks footwall contact).
Sub-vertical quartzreef
Sub-vertical reefsampled
in foreground with
horizontally area of
(samples flat 7),
1 and veining onflat
whereas thefootwall
footwall (blue
veins line marks
sampled
footwallvertically
contact).(samples
Sub-vertical
2–6). reef sampled horizontally (samples 1 and 7), whereas flat footwall veins
sampled vertically (samples 2–6).
Table 8. Likely sampling errors during the collection of chip or chip-channel samples.
(A) (B)
(C) (D)
Figure 4. Collection of saw-cut channel sample (A) and resulting channel (B); Cutting a sidewall
Figure 4. Collection
channel sample (C)of saw-cut
and channel
using an air pick sample (A)
to remove theand resulting
sample (D). channel (B); Cutting a sidewall
channel sample (C) and using an air pick to remove the sample (D).
A typical 1 m length saw-cut channel cut should yield a sample mass of 4 kg (2.5 cm by 5 cm
profile).
A typical It is1 unlikely
m length that perfect sample
saw-cut channelextraction will yield
cut should be achieved due mass
a sample to irregularities at the
of 4 kg (2.5 cmbase
by 5 cm
profile). It is unlikely that perfect sample extraction will be achieved due to irregularities at the base of
the cut (e.g., where the sample block breaks free, EE). Where a large sample is required, channel width
and/or depth can be increased or multiple channels cut. The principal barrier to low DE in saw-cut
channels is rock face roughness. Where the face is very uneven, then even a saw cut is difficult to place
and may result in high DE and EE across the sample zone.
Single- or double-blade saw units are generally powered by compressed air and require water
for dust suppression. Single-blade units are the most appropriate, from both a safety and efficiency
perspective. Double-blade units have the advantage that two parallel cuts are made at the same time
2 cm apart, though double-blade units are generally heavier and pose more of a manual handling and
safety issue. Multiple cuts, one below the other, can be made using the upper or lower original cut as
a guide.
Saw-cut channels are generally of higher quality, though bias may occur when fines are washed
away during cutting, particularly in soft sections (EE). With rough/uneven rock surfaces, sample
quality can be compromised since the extracted sample has a variable support along the channel.
A larger saw blade may help to ameliorate a rough surface. The cost and time for taking saw-cut channel
samples depends on rock hardness and the availability of services to run the saw, but experience
Minerals 2018, 8, 232 15 of 45
shows that a well-trained two-person team can cut and collect a metre of channel sample in 10–15 min.
Proper training and the use of personal protective equipment will reduce safety concerns with the use
of saws.
A reasonable quality channel sample can be extracted using an air pick chipping device
(Figure 4D), but it is prone to operator induced DE bias, either purposeful (e.g., a deep cut around a
high grade section) or accidental (e.g., a deeper cut around soft material). Air picks have the tendency
to produce fly-rock, which is both an EE and safety issue. Channel sample attributes and errors are
given in Tables 10 and 11.
Table 10. Likely sampling errors during the collection of channel samples.
1 Error Magnitude
Sampling Error Error Effect
Hand Cut Saw Cut
Generally good as the sample is similar to a
DE Moderate Moderate
drill hole.
High if different rock properties are
encountered across a sample. For hand
collected samples, hard rock leads to loss of
sample in hard rock as fly particles. In
EE High Moderate–low
soft/friable material will lead to bias due to
oversampling or loss. Saw-cut channels
generally better, though cooling water may
lead to loss in soft/friable areas.
Dependent upon EE and specifically rock
WE Moderate Low hardness. Can be variable. Reduced if saw-cut
channel.
1 Indicative error magnitude: Red: high (>±15%); Orange: moderate (±5–15%); Green: low/acceptable (<±5%).
Attribute Comment
Moderate–difficult: requires time and/or compressed air saw
Logistics and planning
to collect.
Moderate: any reasonable number of samples can be collected, but
Frequency of collection
will take more time than chip
Low–high: Depends upon sampling location. Local safety protocols
must be applied. Requires sampler to be at rock face. Risk of fly rock
Hazard exposure
if collected by hand or air pick. Care required if a diamond saw is
used due to rotating parts.
High: can be undertaken at short notice. Slow to cut by hand, 30 min
Flexibility of method or more for 1 m sample. Can be quite quick to collect a saw-cut
sample once the saw is set up.
Cost High: labour intensive and can require compressed air saw.
Sample quality High: potentially a high quality sample.
Moderate–high: whilst relatively time consuming and costly, the
Value proposition
higher quality sample provides increased value.
Dispersed clusters of coarse gold in a vein system at the Gwynfynydd Mine (UK) created a
negative bias in grade control samples, leading to the replacement of standard chip-channel samples
by panel-channel
Minerals 2018, 8, 232 samples [12]. Faces were divided into horizontal quarters and close-chipped to 16give
of 45
a 25 kg sample from across each of two panels. Processing of the combined 50 kg micro-bulk sample
was through a gravity separation unit located on surface. The resulting data were used to map grade
throughalong
trends a gravity separationand
development unitto
located on surface.
provide The resultinggrade.
a semi-quantitative data were
Theused to map grade
potentially larger trends
mass
along development and to provide a semi-quantitative grade. The potentially larger mass
(>100 kg) of panel samples compared to chip or channel samples, means they are typically used for (>100 kg) of
panel samplessample
metallurgical compared to chip (Figure
collection or channel
5). samples, means they
Panel sampling areare
errors typically
given used for metallurgical
in Table 12 and key
sample collection
attributes in Table 13.(Figure 5). Panel sampling errors are given in Table 12 and key attributes in Table 13.
(A) (B)
(C) (D)
Figure 5. (A) Panel sample preparation. Large (100 kg) samples were required for metallurgical testing
Figure 5. (A) Panel sample preparation. Large (100 kg) samples were required for metallurgical testing
and grade verification; (B) Collecting panel samples with hammer and chisel. Samples were chipped
and grade verification; (B) Collecting panel samples with hammer and chisel. Samples were chipped
from the drive sidewalls onto a plastic floor sheet; (C) Part of panel sample. Chip fragment size is
from the drive sidewalls onto a plastic floor sheet; (C) Part of panel sample. Chip fragment size is
reasonably well controlled, with a maximum fragment size of about 8 cm, and most below 5 cm. With
reasonably well controlled, with a maximum fragment size of about 8 cm, and most below 5 cm. With
samples
samples collected
collected from
from brittle
brittle material
material such
such as
as quartz
quartz there
there are
are aa substantial
substantial proportion
proportion of
of sub-1
sub-1 mm
mm
fragments present; (D) Alternate sample capture method using rubber
fragments present; (D) Alternate sample capture method using rubber mats. mats.
Table 12. Likely sampling errors during the collection of panel samples.
Table 12. Likely sampling errors during the collection of panel samples.
Sampling Error 1 Error Magnitude Error Effect
Sampling Error Low–moderate1 Error Magnitude Error on
Effect
DE Generally good as the sample is delineated the face.
DE Low–moderate Generally
Can be acceptable, but maygood as thehigh
be locally sample is delineated
as rock on the face.
chips are collected. May be
EE Moderate–high Can be acceptable,
very high if different but may
rock properties arebe encountered
locally high as rock chips
across a panel.
are collected. May be very high if different rock
Oversampling of “shiny” and/or soft/friable material.
EE Moderate–high
properties are encountered across a panel. Oversampling
WE Low–moderate Dependent upon EE and specifically rock hardness. Can be variable.
of “shiny” and/or soft/friable material.
1 Indicative error magnitude: Orange: moderate (±5–15%);upon
Dependent Green:
EE low/acceptable
and specifically(<±5%).
rock hardness. Can
WE Low–moderate
be variable.
1 Indicative error magnitude: Orange: moderate (±5–15%); Green: low/acceptable (<±5%).
Minerals 2018, 8, 232 17 of 45
Attribute Comment
Logistics and planning Moderate–difficult: requires time to collect by hand.
Low–moderate: any number of samples can be collected, but will take
Frequency of collection
time and depend upon sample size.
Low–high: Depends upon sampling location. Local safety protocols
Hazard exposure must be applied. Requires sampler to be at rock face. Risk of fly rock
during collection.
High: can be undertaken at short notice, but may take at least 30 min
Flexibility of method
per sample to collect.
Cost Moderate–high: potentially cost intensive for large samples.
Moderate–high: potentially a high quality sample with large area of
Sample quality
face coverage.
Moderate–high: whilst relatively time consuming and costly, will
Value proposition
provide a large sample mass.
Where multiple channels or chips are being collected across a face or sidewall, they should be
collected from the bottom sample upwards. This approach avoids contamination of lower samples
from those above from falling fragments. Where the sample band/channel is inclined to cut the ore
zone at right angles, sampling should commence at the lower end and work upwards.
The greatest EE during face sampling arises from differential, hard versus soft rock properties, as
a greater quantity of softer material usually enters the sampling bag (Fourie and Minnitt, 2016). If this
is high-grade, then the sample will be positively biased, and if low-grade, vice versa. In very hard
sections, there is a likelihood of sample loss, as rock chips from each hammer blow will fly away from
the sampling bag or a lack of sample collected.
In a study based on Witwatersrand (RSA) gold mines, Fourie and Minnitt [21] identified an
“in-stope sample discard process”, where samplers were prone to make a selection of what they believe
was better looking material. Such discard of waste and preferential selection of mineralised material
leads to a positive “waste discard bias” and hence higher sample assays [21].
In well-jointed and/or bedded rocks, collection by hammering presents difficulties in that adjacent
non-sample material may fall into the sample lot (e.g., contamination; PE). This alien material is
unwanted and needs to be either removed or avoided by screening or barring down the loose zone
prior to sampling.
The quality of chip and channel samples is highly dependent upon the sampling crew and the
underground environment. The best-trained crews are unlikely to collect a high-quality sample.
Quality will vary from crew to crew and from shift to shift. The mine environment will also effect
quality through local temperature/humidity (e.g., fatigue effects), sample location (e.g., drive vs. stope
vs. raise), and time of shift (e.g., fatigue/motivation effects). Additional bias (DE) also comes from
the human eye, where the geologist will be drawn to collecting a larger sample from visible gold or
high-grade areas (e.g., sulphide-rich). Collector bias can be reduced by ensuring that all samples are
collected along a marked grade line and staff training (Figure 1).
In many mines, face samples are taken at a height above the drive floor, often 1 m ± 0.5 m, though
this depends on the type of ore body. For example, in the shallow-dipping Witwatersrand reefs, face
samples are taken across the reef irrespective of how high or low it is from the floor of the stope.
In some operations, the face is divided into two horizontal halves, with a single channel taken along
the half line or one sample taken along the centre of each half (i.e., two face samples). For larger drive
profiles, 3 m by 3 m and greater, it is difficult to sample safely above 2 m. In this case, a service vehicle
or rugged step-ladder is required to provide a safe working platform. In a production environment
this may be difficult, and in some cases a sample may not be collected. If relevant, the geologist will
have to use their experience to evaluate the effect of not sampling a particular geological domain.
Rowland and Sichel [47] compared chip and channel sampling within selected Witwatersrand
mines, concluding that: cutting chip or channel samples by hand was difficult, sampling supervision
and training are critical, and when using compressed air saws or picks it was possible to collect better
quality samples quicker. Storrar [11] emphasises the need to minimise DE and EE during sample
collection, and the need for the sampler to be motivated, trained, and encouraged to maintain a high
standard in sampling practice.
few mines will neither realise or admit this. However, it can be very difficult to get the cooperation
of the drilling crews if their bonus depends on meters of advance or number of holes drilled. Grade
control activity should be factored into the mining contract through geological-related quality checks,
for example, stope and drive profile measurement and sampling.
Where shanty drives are cut in narrow low-angle ore zones, the exposure of mineralisation in
sidewalls allows sampling to be potentially undertaken behind the development face if immediate
ore/waste decisions are not required (e.g., Nalunaq [32]; Figure 2A).
There can be an inverse relationship between the size of a mine and sample mass collected at
face, generally where smaller mines collect larger samples. Similarly, there can also be a relationship
between mine size and sample quality, where smaller mines collect higher quality samples. However,
some large mines, notably those on the Witwatersrand, have put much time and effort into improving
sample quality [19–21]. Mass and quality issues relate to: the rate of face advance, time pressures, and
availability of the geologist. Some mines require 24/7 coverage by the geological team to ensure all
development faces are sampled.
Table 14. Sample quality rating as applied at the Gwynfynydd mine. Rating applied to both
chip-channel 1 and panel 2 samples.
An alternate method is to record the mass of each sample at the laboratory after drying.
This measures the quantity of recovered rock over the expected mass. A target acceptance of ±5%
is nominally taken around the expected. Above or below this range, the samples are either over- or
under-sampled. This method is not perfect, for example an under-extracted sample could yield a close
to normal mass based on contamination by fallen material out of the delimited area.
Example of a mass-recovery QI for linear samples is shown in Table 15, where the rating is based
on percentage of samples on target. The data was collected from the same mine by the same crew
from 2.5 m by 2.5 m faces, which comprised 1.0–1.3 m of quartz vein hosted in weak–strongly altered
dolerite wall rocks. The highest quality sample was the saw-cut channel, which achieved 75% of
samples on target. The lowest quality was the chip samples, achieving 0% with all under-sampled.
The air-pick cut channel yielded a low value, which was due to excessive loss (EE) by fly rock.
Table 15. Example of mass recovery (quality index (QI)) for linear samples across quartz-dominated
vein mineralisation from a North American mine [42].
4. Grab Sampling
4.1. Introduction
There are a number of mines using grab samples (or muck or broken rock samples) as a method of
grade control [51,52]. The goal of grab sampling is to predict the grade of ore from a particular source
(e.g., stope, stockpile, or truck) and/or reconcile with the predicted grade. It is often used because of
access issues (e.g., non-entry stopes—sub-level caving operations), safety (e.g., to avoid unsupported
backs), or lack of other sample data. On occasion, samples may also be taken from drive or stope floors
to monitor gold loss. Surface grab sampling is often used to monitor the grade of stockpiles prior to
blending or feeding to the plant or for reconciliation purposes. A fuller discussion and case studies are
given in Dominy [51], and Carswell et al. [39] present a case study based on the high-grade Gosowong
mine (Indonesia).
Figure 6. Grab sampling from a surface stockpile showing the challenges of collecting a representative
Figure 6. Grab sampling from a surface stockpile showing the challenges of collecting a representative
sample and oversampling fine material.
sample and oversampling fine material.
Figure 7. Drawpoint with sampling grid (Left) and sample collection (Right).
Figure 7. Drawpoint with sampling grid (Left) and sample collection (Right).
The use of grab sampling in block cave mines, whilst problematic, provides management with
The use of grab sampling in block cave mines, whilst problematic, provides management with
information which allows them to make decisions where the grade appears to be reducing. Future
information which allows them to make decisions where the grade appears to be reducing. Future very
very large-scale caving operations will require a new approach to production sampling to meet safety
large-scale caving operations will require a new approach to production sampling to meet safety and
and health expectations. The numbers of draw points (potentially up to 1000) combined with the use
of more automation mean that a faster and smarter sampling will be required to meet the demands
of bigger mines [52].
health expectations. The numbers of draw points (potentially up to 1000) combined with the use of
more automation mean that a faster and smarter sampling will be required to meet the demands of
bigger mines [52].
the fact that the samples were representing the <3 cm fraction well, but not the coarser >3 cm fractions.
In this case there was low heterogeneity in the high-grade fine fraction related to the consistent
liberation of gold-bearing sulphides. However, this fraction contained around 50–65% of the
total gold. Full characterisation, including a granulometric study is required to understand grab
sample representativity.
Table 16. Examples of grab sample performance from two contrasting mineralisation styles from South
America [42].
4.5. Summary
The key issue with grab sampling is the size of primary sample required to be representative.
The few kilogrammes that are usually collected are generally inadequate and leads to a large FSE,
potentially to ±500% or more. This error may be reduced by collecting multiple small samples over a
pile, for example fifty 5 kg samples giving a composite of 250 kg. Pitard [23] suggests that a minimum
of 32 individual samples should be taken from a pile to reduce the GSE. In the most challenging
of cases (e.g., some coarse gold systems), it is possible that hundreds of kilogrammes or tonnes are
required for each sample composite. Where large samples are collected, the best option is sample
size and fragment reduction via a sampling tower. The FSE equation can be used to optimise sample
reduction, where Equation (1) can be rearranged to solve for MS .
The method is prone to high sampling errors (Table 17); key attributes are given in Table 18.
Grab sampling may be effective where grades vary little between size fractions (e.g., fine disseminated
Minerals 2018, 8, 232 25 of 45
gold that is locked in sulphides; Table 16) and where both the ore and waste break into pieces of
approximately equal size. In many gold mines, despite best efforts, the geologist should expect
precision to be poor and ore/waste misclassification high.
There are situations when grab sampling is the only way to evaluate grade, for example where
there is no drilling and/or development cannot be safely accessed. Grab sampling is known to
work sometimes in low-nugget effect mineralisation styles (e.g., disseminated gold deposits), but in
heterogeneous coarse gold systems a positive bias is likely given that the fine fraction is generally
enriched in gold (Table 18) [39]. If grab sampling is the only option, then gold particle characterisation,
granulometric, and duplicate programmes must be undertaken to determine representativity.
Attribute Comment
Logistics and planning Easy: relatively easy to undertake.
Frequency of collection High: numerous small (<6 kg) samples can be collected.
Low–high: Depends upon sampling location; risk with multiple
small samples as composite high mass results in manual handling
Hazard exposure issues. Sampling of surface or underground stockpiles may pose
traffic hazards. May require sampler to be at the rock face if
collecting from a muck pile.
Flexibility of method High: can be undertaken quickly and at relatively short notice.
Low: relatively easy and quick to collect, cost principally sample
Cost
preparation and assaying.
Sample quality Poor: prone numerous sampling errors.
Value proposition Low: strong likelihood of bias, and thus ore/waste misclassification.
5. Drilling Methods
Diamond core drilling is used in a number of mines for grade control and short- to medium-term
planning purposes [7,40,53–55].
(A) (B)
Figure 8. Types of sludge sampling. (A) Collection of stope blasthole sludge via plastic boxes placed
Figure 8. Types of sludge sampling. (A) Collection of stope blasthole sludge via plastic boxes placed
below
below the
the hole
hole leads
leads to
to very
very high
high DE
DEand
andEE;
EE;(B)
(B)Collection
Collectionof
offace
faceblasthole
blastholesludge
sludgevia
viascooping
scoopingfrom
from
the drive floor leads to high FSE, GSE, DE, and EE.
the drive floor leads to high FSE, GSE, DE, and EE.
A small number of operations have practiced the collection of sludge samples from development
faces by scooping cuttings from the floor of the level (Figure 8). This method is prone to very high
sampling errors and best described as a gambling process. In all sludge sampling exercises, there are
usually problems with loss of fines and the sample being biased by contamination from the sides of
the hole. Care needs to be taken whilst splitting a sub-sample to minimise DE, EE, and PE. Like RC,
sludge drilling produces a large primary sample, though generally of low quality. A similarly large
Minerals 2018, 8, 232 27 of 45
A small number of operations have practiced the collection of sludge samples from development
faces by scooping cuttings from the floor of the level (Figure 8). This method is prone to very high
sampling errors and best described as a gambling process. In all sludge sampling exercises, there are
usually problems with loss of fines and the sample being biased by contamination from the sides of
the hole. Care needs to be taken whilst splitting a sub-sample to minimise DE, EE, and PE. Like RC,
sludge drilling produces a large primary sample, though generally of low quality. A similarly large
FSE can result, where a small sub-sample is taken from the primary sample.
Ovalle et al. [58] propose the “Sampdrill” method which utilises a hand-held drill to place multiple
small diameter holes into the drive face or walls. The cuttings from each hole are then taken as an
individual sample or composited to a larger sample. This methodology was applied to porphyry
copper mineralisation, but is less likely to be relevant to gold unless it is fine-grained and disseminated.
Sludge sampling is often used to check where the ore zone is greater than the drive width. In other
cases, it is used to check that blastholes are in ore via assay and/or visual inspection of channels.
Sludge-based assays are generally not used to support resource estimates.
5.4. Comment
A reduction in the need for linear samples results from the use of more pre-development drilling
at a spacing to allow local estimation [53,55]. McBeath and Whitworth [55] present a case study
where drilling was reduced to 5 m to improve ore/waste definition in a geometrically complex and
high-nugget effect ore body. Both diamond core and RC drilling can be used underground, with RC
drilling providing a relatively quick and cheap alternative.
independent and accredited third party. If the laboratory is run by the mining company, then it is
critical that it be independently audited on a regular basis and that umpire assays are undertaken as
part of its QAQC programme.
Sample mass and fragment size reduction must be undertaken within the framework of the TOS;
the FSE equation can be applied (Equation (1)). Particular care is required in the presence of coarse
gold, with possibilities for gold loss and/or contamination [24,30,32,54,59,62]. Broken rock samples
from the mine (e.g., a face channel sample) must be crushed prior to splitting and pulverisation to
avoid high DE, EE, and GSE. Large pulverisers such as the LM5 (maximum capacity 3 kg) are useful;
however, removal of the pulp needs to be maximised (EE) and contamination controlled (PE). It is
recommended that laboratories include a barren flush after visible gold-bearing samples. However,
if the ore bears coarse gold, then flushing after each sample is optimal. Assaying of the pulveriser
wash (e.g., sand) is recommended at a rate of 1 in 10–20 for coarse gold ores. DE and EE are potentially
high during sample mass reduction, and a Jones riffle splitter or rotary sample divider should be used.
Mat-rolling and grab sampling from the pulp must be avoided (DE, EE, and GSE) [1,4,24].
The traditional paradigm of crushing and pulverising the entire sample (e.g., chip or channel
sample, or half drill core) and taking a 30–60 g sub-sample for fire assay is potentially flawed [4,7,17].
The approach is prone to high CSE and ISE, particularly when the assay charge is scooped from
the pulp. The propensity of gold not to pulverise efficiently (exacerbated in the presence of coarse
gold), promotes FSE and GSE, and potentially loss or contamination (e.g., PE) during sub-sampling for
the fire assay charge.
The use of an entire sample for total assay effectively yields FSE and GSE values
of zero [7,12,54,59,63]. With good laboratory practice, the PE and AE can be minimised.
Arguments against whole core sampling in particular, revolve around no reference core remaining for
re-logging, check sampling, etc. However, the authors contend that with modern digital photography,
detailed logging and internal and/or external peer review that this should not be an issue. A reference
sample can be split by rotary sample divider after crushing.
There are a number of techniques employed to assay samples for gold (Table 19); these include
traditional fire assay, and large charge methods such as screen fire assay, LeachWELL, pulverise and
leach, and whole sample by gravity. The LeachWELL method is particularly suited to gold assaying,
as samples of up to 5 kg can be leached [64]. LeachWELL is a catalyst that increases the dissolution
rate of gold, where samples can potentially be leached in less than 24 h. The gold laden solution is
then analysed via an instrumental method such as atomic adsorption spectrometry. LeachWELL is not
a total assay, therefore evaluation of its performance on a given ore type is required. In most cases, the
leach tails are subjected to a duplicate fire assay to provide a total gold grade.
Sample protocols should not be left solely to the laboratory to decide, but established by the project
team based on characterisation. In the absence of duplicate samples, the FSE equation (Equation (1))
can be used to evaluate protocol precision across different ore types. Examples are given in Table 20
for a series of protocols based on a fine- and coarse-gold ore types. It must be remembered that the
FSE equation only evaluates the FSE; an acceptable FSE value can be swamped by high GSE, DE, EE,
etc., thus making it inappropriate for use.
Table 20. Sensitivity of the FSE and sampling errors based on ore type (grade and liberation diameter)
and sample protocol.
Table 21. Best practice QAQC for a grade control sampling programme that informs a resource estimate.
6 Operation and
Blanks 1 in 20 Less than 0.05 g/t Au
laboratory
7 Pulp quality 1 in 20 Laboratory 95% to be P95 −75 µm
8 Barren flush 1 in 20–50 Laboratory <<0.5% gold loss
9 Umpire 1 in 20 Operation 90% ±10% HARD 90% ±13% RSV
Laboratory audit Quarterly Operation Full adherence to agreed practices and performance levels
Minerals 2018, 8, 232 31 of 45
Test work from a narrow vein deposit exemplifies the impact of sampling error through the
comparison of chip and channel samples (Table 22). Seventy-five sample triplicates (chip, hand-cut
channel, and saw-cut channel) were collected from around a 40 m by 20 m stope block. The block
was sampled from faces located every 1.5 m along its upper and lower drives and two raises. After
cleaning, a reference line was drawn across each face centre and the samples collected from the bottom
up; chip sample (7 cm wide by 1.5 cm deep), hand-cut and then saw-cut channels (4.5 cm wide by
2.5 cm deep). The sample delimitation dimensions were designed to achieve a theoretical sample
support of 3 kg/m. The mineralisation was known to have a moderate variability, containing visible
gold to 1.5 mm in size. In the test area, the vein was near vertical and its width relatively consistent,
ranging from 0.8 m to 1.2 m. Only the vein was sampled. All samples were prepared and assayed in
the same way, via a total sample preparation and screen fire assay route. The FSE for the protocol was
effectively zero. A QAQC programme was applied, with all CRMs and blanks within expectation.
The results show a marked reduction in RSV and nugget effect between the three sample sets
(Table 22). The rigorous laboratory protocol and QAQC indicate that errors within the laboratory were
at a minimum. Therefore, the remaining variability relates to the geological nugget effect and sample
collection. The dominant error for the channel samples relates to the geological nugget effect, given
that sampling error was minimal. The dominant difference between the chip and saw-cut channel
samples relates to sampling error. The results corroborate previous findings, showing that saw-cut
channel samples provide the best sample quality.
Table 22. Example from an Australian mine showing comparison between chip and channel sample
replicates. Total number of samples in each type was 75. Mean grades cut and diluted to stope width.
Reconciled stope head grade 13.7 g/t Au.
(<±20%). 2 Indicative preparation and assay error: Red: high (>±20%); Orange: moderate (±10–20%); Green: low
(<±10%).
Duplicate sample analysis across the sampling and laboratory process allows errors to be
apportioned to each phase of activity [34,41]. Dominy [42] reports test work from a moderate coarse
gold-bearing (d95Au : 750 µm) vein system. A series of fifty non-contiguous underground faces were
sampled via chip, hand-cut channel, and saw-cut channel samples. A field duplicate was taken for
each, and at the laboratory coarse and pulp duplicates were taken (Table 23). The results demonstrate
the dominance of sample collection error and higher chip sample error compared to the channel
samples (e.g., reduced DE and EE). The high preparation error most likely represents liberated coarse
gold present in the crushed material. The analytical error is at the high end of acceptable, reflecting
some remnant heterogeneity from the coarse gold. The total error could be reduced to around 23% by
putting the entire saw-cut sample through a LeachWELL assay.
There a need to move towards the quantification of errors, where a first step is the application of
the RSV as defined in DS3077 [38]. Resolution of individual relative errors across sampling, preparation,
and analysis can be gained from duplicate sample pairs (Tables 23 and 24). From a spatial perspective,
variography provides a method to determine the relative proportion of in-situ/geological versus
sampling nugget effect based on duplicate field pairs.
Minerals 2018, 8, 232 33 of 45
Table 23. Example from a North American mine of RSV values for duplicates taken across each mine
face [42]. See also Table 24, Mine B example.
Table 24. Examples from South (Mine A) and North (Mine B) America showing error comparison
between two differing styles of gold mineralisation [42].
methods. Smaller, often narrow vein operations with a small resource base will be more reliant on
definition of stope blocks by linear sampling methods. For these mines, the decision as to what will be
mined has effectively already been made. A retrospective decision to process a development round
can only be made if it has been stockpiled separately. Small mines generally require selective mining,
therefore optimal sampling and integration with mapping becomes critical. A low-grade stockpile
is recommended to ensure that misclassified ore can be recovered, even if at the end of the mine
life [13,14,65].
An asset of the geologist is the ability to work and communicate with mining crew(s).
Communication needs to provide a rationale for good samples, why poor samples are problematic
and provide an explanation for the methods used. There are often excuses the geologist does not
have enough time at the face because the miners do not want them there. There is some truth in
this, though a good working relationship with crews will pay dividends. In particular, if a crew
understands why a geologist is mapping and sampling a face, then their attitude is more likely to be a
positive one. A progressive operation will include grade control familiarisation in its induction and
training programme.
can apply an appropriate resource classification to provide a measure of this risk in the estimate
(e.g., Inferred, Indicated and Measured Mineral Resources) [9]. For example, in a high-nugget
deposit where multiple high-quality samples have been collected around a series of stope blocks,
the classification may only be Inferred, whereas exactly the same scenario in a low-nugget deposit
may yield a Measured resource. As the nugget effect increases, local estimation becomes more
difficult, and thus reconciliation becomes potentially poorer on a short-term basis [14]. Many high-
to extreme-nugget effect ore bodies will operate solely on Inferred Mineral Resources, though this is
dependent upon the data spacing available [14,54]. Moderate- to low-nugget effect ore bodies will
generally define at least Indicated Mineral Resources with Probable Ore Reserves [66].
9.4. Reconciliation
The ultimate challenge of any mine is that of achieving acceptable reconciliation between the
long-term reserves and short-term grade control model with actual plant performance. In many gold
mines, the reserve estimate overstates the achieved grade. Whilst reconciliation problems can be
related to estimation approach and/or mining, there are often underlying reasons related to sampling
error particularly bias (e.g., in-stope sample discard or soft reef bias: refer Section 3.3) and incompatible
datasets (e.g., mixing drillhole and channel sample results).
In higher nugget mineralisation, reconciliation over short periods of time (e.g., month by month)
will generally be difficult as variability may be ±50% or more [14]. Thus, longer periods of time
are likely to show more acceptable results (e.g., ±20% over three months). Operators should avoid
knee-jerk reactions such as altering cut-off grades or blending strategies based on one or two months of
poor reconciliation. Viewing a resource/reserve estimate as an “estimate” is important and operators
should understand the likely tolerances of grade and tonnage figures. In particular, they should
understand the sampling process and its effect on estimates and decision making.
Figure 9. Backs channel sampling in historical workings. A chip-channel is produced, where sample
Figure 9. Backs channel sampling in historical workings. A chip-channel is produced, where sample
fragments are supposedly being collected in a plastic bowl. In this case the EE is high, as many extracted
fragments are supposedly being collected in a plastic bowl. In this case the EE is high, as many extracted
fragments are not being collected. In addition, DE on the shoulders of the drive is particularly high since
fragments are not being collected. In addition, DE on the shoulders of the drive is particularly high
the delimited samples are biased within the delimited zone on the curved drive shoulder.
since the delimited samples are biased within the delimited zone on the curved drive shoulder.
The digital mine is high on the agenda of many companies [67]. From a grade control perspective,
assay data can feed into the digital realm and can be interfaced with digital photographs, outputs from
automated mapping, and surveys. Once into the digital platform, artificial intelligence may then allow
direct feed into development control and design, and grade models. Data analytics will play a role in
analysing grade and other data to produce more effective and agile short-term mine plans, potentially
real-time information, and lead to automated reconciliation. Sensors in drill holes, at the mine face and
in the plant, will be critical for this development.
Table 25. Potential contribution of grade control samples to geometallurgical ore control. Assumes
best practice sample collection for all samples. Green: applicable; Orange: potentially; Red: no.
Table 26. Example of risk review for a grade control sampling programme used to support a resource
estimate [6].
Table 27. Summary of steps to achieve optimal grade control sampling programmes.
• Select a sample strategy based on geology, mining method/selectivity and safety requirements
• Consider sampling method(s)
• Ensure that sample collection, preparation, and assaying protocols are designed with reference to the ore type
2: Develop • Ensure TOS and QAQC is applied
• Document within the framework of DS3077 [38]
• Review cost-benefit of more grade control drilling to provide better grade control/local estimates and thus
negating the need for face sampling
• Ensure that staff are properly equipped and trained, with a focus on reducing sampling errors and increasing
3: Execute sample quality
• Ensure regular staff supervision and introduce a sampling quality index
5: Refine • Update where required and/or return to stages 2 and 3 for further work
Minerals 2018, 8, 232 39 of 45
A five 2018,
Minerals stage8, xapproach is recommended starting with (1) scoping; moving to (2) development
FOR PEER REVIEW 37 of 43
through to (3) execution; (4) validation; and (5) refinement. Beyond establishing the business case
and aimsAoffive thestage
gradeapproach
control is recommended starting with (1) scoping; moving to (2) development
programme, a major precursor is the gold characterisation study that
through to (3) execution; (4) validation; and (5) refinement. Beyond establishing the business case and
focusses on the evaluation of gold particle size and distribution. It is important to understand the
aims of the grade control programme, a major precursor is the gold characterisation study that
mineralisation and controls to gold distribution through evaluation of geology, ore body shape and
focusses on the evaluation of gold particle size and distribution. It is important to understand the
architecture, and grade
mineralisation distribution.
and controls to goldSome mineralisation
distribution styles areofmore
through evaluation forgiving
geology, than
ore body others
shape andgiven
their architecture,
sampling characteristics (Figure 10).
and grade distribution. Some mineralisation styles are more forgiving than others given
Characterisation leads to development
their sampling characteristics (Figure 10). of the sampling strategy and selection of sampling
methods.Characterisation
Each method has leadsrelative
to development
positivesofandthe negatives,
sampling strategy and selection
which should of sampling
be considered during
methods.
scoping (TableEach method
28). Early inhas
the relative
planningpositives
processand negatives,
(Stage which
2), there mayshould be considered
be minimal during data
actual sample
scoping
available, (Table 28). duplicates.
particularly Early in the planning process
In this case, (Stage 2),by
supported there may be minimal
characterisation actualthe
results, sample data
model-based
available, particularly duplicates. In this case, supported by characterisation results, the model-based
FSE equation (Equation (1)) can be applied to evaluate sampling protocols. As the programme develops,
FSE equation (Equation (1)) can be applied to evaluate sampling protocols. As the programme
then duplicate samples permit an empirical approach to the estimation of error.
develops, then duplicate samples permit an empirical approach to the estimation of error.
d
o
C 75-100%
m
o
i
a
n
r
a
s
t
e 50-75%
e
d
d
o
25-50%
m
F
i
i
n
n
a
e
t
0-25%
e
d
Figure 10. Generalised applicability of sampling methods based on mineralisation type (dominance
Figure 10. Generalised applicability of sampling methods based on mineralisation type (dominance of
of coarse vs. fine gold and nugget effect). Red: avoid method; Orange: avoid method where possible;
coarse Green:fine
vs. gold applicable,
method and nugget effect). good
assuming Red: avoid method;
practices, Orange:
optimised avoidand
protocols, method
QAQC. where possible;
[1] Coarse
Table 28. Summary of sampling errors and risk rating for grade control sampling. For linear/broken
rock sampling, it assumes a fully equipped two-person team to achieve best practice sample collection.
Errors are best-case scenarios based on careful sample collection/operation.
1 Sampling Errors
Sample Type Summary 2 Total Error Rating
DE EE WE
Linear
Easy to collect, relatively fast
Chip High number can be collected High
Cheap
Relatively easy to collect and fast
Chip-channel Moderately high number can be collected High
Moderately cheap
Less easy to collect, requires effort
Channel
Moderately high number can be collected Mod
(hand-cut)
Higher cost
Requires effort and specialist equipment
Channel
Moderate number can be collected Low
(saw-cut)
Higher cost
Medium-large sample size (>25 kg)
Large sample size
Panel Mod
Less easy to collect
Moderate number can be collected
Broken Rock
Relatively easy to collect
Grab Moderate number can be collected High
Cheap
Drilling
Well established method
Diamond core Good geological information Low
Relatively costly
FSE potentially mod-high when sub-sample split from
Reverse
original, though modern technology much improved Low
circulation
Cheaper than core drilling
FSE potentially mod-high when sub-sample split from
original
Sludge
Sample collection/splitting at the rig problematic, though Mod
(blasthole)
some mines have devised a reasonable catcher system
Cheap
Mini-bulk (<1 t)/Bulk (>1 t)
Comprise (i) panel or (ii) blasted/excavated blocks
Careful planning required
Bulk Low
Good geological and geotechnical information
Relatively costly depending upon mass
1 Indicative sampling error definition: FSE: Red: high (>±30%); Orange: moderate (±20–30%); Green: low (<±20%).
2 Indicative total error: Red: high (>±50%); Orange: moderate (±30–50%); Green: low (<±30%).
13. Conclusions
(1) Practitioners can be justifiably critical of linear and grab sampling methods, given their potential
for poor precision and high bias. However, they have a role to play in underground grade control,
depending on their ultimate use and effort put into making them representative. Maintaining
representative samples is advantageous; however, there will always have to be a balance between
what is theoretical, practical, and safe. It is important to understand what decisions the sample
data will be used for and the inherent risks involved both in collection and decision-making.
Critical end-member applications are (a) to support a reserve estimate for short-term mine
planning, or (b) as ”indicators” to map grade trends and/or deportment along development and
in specific geological domains.
(2) A range of sampling methods are available, all of which require evaluation before routine
acceptance. The greatest error is generally introduced during sample collection (e.g., DE and EE).
A cost-effective way to reduce sample collection error involves staff training in TOS and correct
collection methods. Whilst correctly collected channel samples are of a higher quality than chip
samples, chip samples will continue to have importance given their relative ease of collection.
A reduction in the need for linear samples will only come from innovation or the use of more
pre-development drilling at a spacing to allow local estimation.
Minerals 2018, 8, 232 41 of 45
(3) Where grade control samples support a reserve estimate, a QAQC programme must be introduced
to ensure on-going quality of sampling and assaying. Written protocols and procedures, periodic
auditing of protocols and people, and training are required. DS3077 provides a framework on
which to compile transparent protocols [38].
(4) Sample protocol optimisation is best undertaken at the breakeven cut-off grade given that it is
the ore/waste decision point. An optimised sampling programme will never be perfect and as a
result, ore/waste misclassification will occur particularly where the nugget effect is high; thus a
low-grade stockpile is recommended. Many small mines do not have drilling, therefore the most
practical and representative linear samples are required.
(5) Samples feed into reserve and grade control models that are reconciled with plant performance.
All errors along the sampling value chain are additive and impart data variability making
estimation less reliable. Errors should be evaluated by empirical means to quantify their impact
on the estimation and decision making process. Where sets of mixed data such as channel and
drill core samples are used for estimation, it is important to ensure that any bias between the
samples is identified and acted upon if appropriate.
(6) The human impact on sampling is great, given that the operating environment is generally
unforgiving and hazardous. The sampling process is a physical one and often carried out by
relatively junior and inexperienced staff. Team training is critical to improve sample quality
and representativity.
(7) Geologists must appreciate the rigors of the TOS and also understand the needs of production.
The “ideal” is advantageous, but not always practical. However, where the ideal becomes
impractical, the potential effects of the non-ideal on the decision making process must be
understood. Management and technical staff should be educated to ‘sample better’ throughout
the mine value chain.
(8) In some cases, the application of grade control samples is an exercise in futility. This usually
reflects extreme-nugget effect coarse-gold mineralisation. Challenges are particularly exacerbated
when grades are low (<5 g/t Au) and/or where the dominant gold particle size is >1 mm or
clusters. In these cases, the best option is likely to be bulk (or mini-bulk) sampling and strong
geological control during mining.
(9) In context of other mineral deposits, particularly those of base metals, the application of TOS
and its error evaluation is relevant. Additionally, use of the FSE equation is also relevant,
though the original FSE equation is usually applied in this case. The challenges of applying the
various sample types are the same. QAQC is required, particularly when samples feed into a
resource estimate.
(10) Company management in their pursuit to reduce costs, often reduce grade control expenditure
with no consideration of the impact on the quality of samples and decisions made thereon.
Persuading managers to spend the money to rectify the problems, even when they are clearly
identified, can be challenging.
Author Contributions: S.C.D. conceived the paper. S.C.D., H.J.G., and L.O. analysed the numerical data presented.
C.K.L. and S.P. reviewed publically available reports to assess industry practice. R.C.A.M. contributed his expertise
from the South African gold mines. All authors contributed to the writing and editing of the paper.
Acknowledgments: The authors acknowledge a number of current and former mining companies across Australia,
Africa, East Asia, Europe, and the Americas for the opportunity to input into grade control programmes. S.C.D.
acknowledges discussions with the late Allan (Bon) Royle and Alwyn Annels on many aspects of mine sampling.
R.C.A.M. acknowledges extensive support from the South African gold mining industry. Thanks are due to the
Minerals reviewers for helpful comments on the manuscript.
Conflicts of Interest: The authors declare no conflict of interest.
Minerals 2018, 8, 232 42 of 45
References
1. Pitard, F.F. Pierre Gy’s Sampling Theory and Sampling Practice; CRC Press: Boca Raton, FL, USA, 1993; p. 488.
2. Glass, H.J.; Zegers, T.W. Statistical aspects of sampling particulate matter. In Proceedings of the Surface
Mining Conference, Johannesburg, South Africa, 30 September–4 October 1996; Southern African Institute of
Mining and Metallurgy: Johannesburg, South Africa, 1996; pp. 157–164.
3. Minnitt, R.C.A. Sampling: The impact on costs and decision making. J. South. Afr. Inst. Min. Metall. 2007,
107, 451–462.
4. Dominy, S.C. Importance of good sampling practice throughout the gold mine value chain. Min. Technol.
2016, 125, 129–141. [CrossRef]
5. Holmes, R.J. Sampling procedures. In Advances in Gold Ore Processing; Elsevier: Amsterdam, The Netherlands,
2016; pp. 45–56.
6. Dominy, S.C.; O’Connor, L.; Glass, H.J.; Purevgerel, S.; Xie, Y. Towards representative metallurgical sampling
and recovery testwork programmes. Minerals 2018, 8, 193. [CrossRef]
7. Dominy, S.C.; O’Connor, L.; Glass, H.J.; Xie, Y. Geometallurgical study of a gravity recoverable gold orebody.
Minerals 2018, 8, 186. [CrossRef]
8. Carrasco, P.C.; Carrasco, P.; Jara, E. The economic impact of incorrect sampling and analysis practices in the
copper mining industry. Chemom. Intell. Lab. Syst. 2004, 74, 209–214. [CrossRef]
9. Joint Ore Reserves Committee (JORC). Australasian Code for Reporting of Exploration Results, Mineral Resources
and Ore Reserves—The JORC Code; Joint Ore Reserves Committee of the Australasian Institute of Mining and
Metallurgy, Australian Institute of Geoscientists and Minerals Council of Australia: Melbourne, Australia,
2012; p. 44.
10. Gy, P.M. Sampling of Particulate Materials: Theory and Practice; Elsevier: Amsterdam, The Netherlands, 1982; p. 431.
11. Storrar, C.D. Underground sampling with special reference to gold reefs. In South African Mine Valuation;
Southern African Institute of Mining and Metallurgy: Johannesburg, South Africa, 1981; pp. 36–63.
12. Dominy, S.C.; Platten, I.M.; Fraser, R.M.; Dahl, O.; Collier, J.B. Grade control in underground gold vein
operations: The role of geological mapping and sampling. In Proceedings of the International Mining
Geology Conference, Perth, Australia, 17–19 August 2009; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2009; pp. 291–307.
13. Pitard, F.F. The in-situ nugget effect: A major component of the random term of a variogram. In Proceedings
of the World Conference on Sampling and Blending, Porto Alegre, Brazil, 23–25 October 2007; Fundacao
Luiz Englert: Porto Alegre, Brazil, 2007; pp. 91–110.
14. Dominy, S.C. Predicting the unpredictable—Evaluating high-nugget effect gold deposits. In Mineral Resource
and Ore Reserve Estimation; Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2014;
pp. 659–678.
15. Dominy, S.C.; Xie, Y.; Platten, I.M. Gold particle characteristics in narrow vein deposits: Implications
for evaluation and metallurgy. In Proceedings of the Narrow Vein Mining Conference, Ballarat, Australia,
14–15 October 2008; Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2008; pp. 91–104.
16. Dominy, S.C.; Platten, I.M. Clustering of gold particles and implications for sampling. Appl. Earth Sci. 2007,
116, 130–142. [CrossRef]
17. Lyman, G.J.; Robertson, I.G.; Day, T. The simple facts about sampling gold ores. In Proceedings of the Mill
Operators Conference, Perth, Australia, 10–12 October 2016; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2016; pp. 159–172.
18. François-Bongarçon, D.M. A brief history of Pierre Gy’s sampling theory and current trends in its acceptance
in the mining world. In Proceedings of the Sampling Conference, Perth, Australia, 27–29 May 2008;
Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2008; pp. 13–14.
19. Freese, B.; Flitton, T.; Pillay, A. Quality control and assurance of underground chip sampling at Kopanang
Mine, South Africa. In Proceedings of the Sampling and Analysis Conference, Muldersdrift, South Africa,
4–6 June 2013; Southern African Institute of Mining and Metallurgy: Johannesburg, South Africa, 2013;
pp. 277–278.
20. Fourie, D.; Minnitt, R.C.A. The simulated chip-sample model as a method for qualifying error and bias in
sampling thin carboniferous reef types. J. South. Afr. Inst. Min. Metall. 2014, 114, 19–30.
Minerals 2018, 8, 232 43 of 45
21. Fourie, D.; Minnitt, R.C.A. Review of gold reef sampling and its impact on the mine call factor. J. South. Afr.
Inst. Min. Metall. 2016, 116, 1001–1009. [CrossRef]
22. François-Bongarçon, D.M.; Gy, P.M. The most common error in applying Gy’s formula in the theory of
mineral sampling and the history of the Liberation factor. J. South. Afr. Inst. Min. Metall. 2002, 102, 475–479.
23. Pitard, F.F. Pierre Gy’s Theory of Sampling and Ingamells Poisson Process Approach: Pathways to Representative
Sampling and Appropriate Industry Standards; Aalborg University: Esbjerg, Denmark, 2009; p. 309.
24. Royle, A.G. Splitting gold assay pulps containing coarse gold. J. Leeds Univ. Min. Assoc. 1989, 89, 63–68.
25. Minnitt, R.C.A.; Assibey-Bonsu, W. A comparison between the duplicate series method and the heterogeneity
test as methods for calculating the sampling constants. J. South. Afr. Inst. Min. Metall. 2010, 110, 251–268.
26. Minnitt, R.C.A. Calibrating K and Alpha in Gy’s formula: A new approach. Math. Geosci. 2016, 48, 211–232.
[CrossRef]
27. Dominy, S.C.; Xie, Y. Optimising sampling protocols via the heterogeneity test: Challenges in coarse gold
mineralisation. Min. Technol. 2016, 125, 103–113. [CrossRef]
28. Dominy, S.C.; Platten, I.M.; Xie, Y.; Minnitt, R.C.A. Underground grade control protocol design: Case study
from the Liphichi gold project, Larecaja, Bolivia. Appl. Earth Sci. 2011, 119, 205–219. [CrossRef]
29. Cintra, E.C.; Scabora, J.A.; Viegas, E.P.; Barata, R.; Maia, G.F. Coarse gold sampling at Sao Francisco
mine, Brazil. In Proceedings of the World Conference on Sampling and Blending, Porto Alegre, Brazil,
23–25 October 2007; Fundacao Luiz Englert: Porto Alegre, Brazil, 2007; pp. 187–197.
30. Gonzalez, P.; Cossio, S. A review of sampling protocol for a gold ore based on liberation study. In Proceedings
of the World Conference on Sampling and Blending, Porto Alegre, Brazil, 23–25 October 2007; Fundacao
Luiz Englert: Porto Alegre, Brazil, 2007; pp. 163–174.
31. François-Bongarçon, D.M. Extensions to the demonstration of Gy’s formula. Explor. Min. Geol. 1998, 7,
149–154.
32. Dominy, S.C.; Petersen, J.S. Sampling coarse gold-bearing mineralisation-developing effective protocols
and a case study from the Nalunaq mine, Southern Greenland. In Proceedings of the World Conference on
Sampling and Blending, Sunshine Coast, Australia, 10–12 May 2005; Australasian Institute of Mining and
Metallurgy: Melbourne, Australia, 2005; pp. 151–165.
33. Fallon, M.; Rowley, M.; Du Plessis, E.; Golenya, F.; Robertson, I.; Swanson, R. Confirming the sampling,
sample preparation and assaying protocols at two Yilgarn gold deposits. In Proceedings of the Sampling
Conference, Perth, Australia, 21–22 August 2012; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2012; pp. 47–62.
34. Stanley, C.R.; Smee, B.W. Strategies for reducing sampling errors in exploration and resource definition
drilling programmes for gold deposits. Geochem. Explor. Environ. Anal. 2007, 7, 329–340. [CrossRef]
35. Simon, A.; Gosson, G. Considerations on quality assurance/quality control and sample security.
In Proceedings of the Sampling Conference, Perth, Australia, 27–29 May 2008; Australasian Institute of
Mining and Metallurgy: Melbourne, Australia, 2008; pp. 135–140.
36. François-Bongarçon, D.M. Theory of sampling and geostatistics: An intimate link. Chemom. Intell. Lab. Syst.
2004, 74, 143–148. [CrossRef]
37. Pitard, F.F. Guidelines for acceptable allotted sampling uncertainty. In Proceedings of the World Conference
on Sampling and Blending, Lima, Peru, 19–22 November 2013; Gecamin: Santiago, Chile, 2013; pp. 89–98.
38. DS3077. Representative Sampling–Horizontal Standard; Danish Standards Foundation: Copenhagen, Denmark,
2013; p. 41.
39. Carswell, J.T.; Yulia, K.; Lesmana, D.; Steamy, K. Grade control sampling quality assurance/quality
control in a high-grade gold mine—Gosowong, Indonesia. In Proceedings of the International Mining
Geology Conference, Perth, Australia, 17–19 August 2009; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2009; pp. 285–290.
40. Sims, D.A.; Benton, R. Life on the Rollercoaster—Mining and milling the Kencana K1 orebody, Indonesia.
In Proceedings of the International Mining Geology Conference, Perth, Australia, 17–19 August 2009;
Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2009; pp. 319–327.
41. Stanley, C.R.; Lawrie, D. Average relative error in geochemical determinations: clarification, calculation and
a plea for consistency. Explor. Min. Geol. 2007, 16, 267–275. [CrossRef]
42. Dominy, S.C. Sampling and Evaluation of Vein Gold Deposits; Final Technical Report for DevR-MIDAS Project;
University of Greenwich: Chatham, UK; Welsh Gold PLC: Dolgellau, UK, September 1998; p. 230.
Minerals 2018, 8, 232 44 of 45
43. François-Bongarçon, D.M. Fishy samples: How big a sample to avoid the infamous Poisson effect.
In Proceedings of the World Conference on Sampling and Blending, Cape Town, South Africa, 21–23
October 2009; Southern African Institute of Mining and Metallurgy: Johannesburg, South Africa, 2009;
pp. 43–46.
44. Gouws, J. The use of diamond saws in underground sampling at Randfontein Estates Gold Mining Company,
Witwatersrand. J. Inst. Mine Surv. S. Afr. 1986, 23, 85–97.
45. Magri, E.J.; McKenna, A.P. A geostatistical study of diamond saw sampling versus chip sampling. J. South.
Afr. Inst. Min. Metall. 1986, 86, 335–347.
46. Potter, D.; Sheriff, C.; Collins, P. Significance of geological control on assay data in the narrow, high-grade
gold quartz veins at the Sand Queen mine, Comet Vale, Western Australia. In Proceedings of the Narrow
Vein Mining Conference, Ballarat, Australia, 14–15 October 2008; Australasian Institute of Mining and
Metallurgy: Melbourne, Australia, 2008; pp. 71–78.
47. Rowland, R.S.; Sichel, H.S. Statistical quality control of routine underground sampling. J. South. Afr. Inst.
Min. Metall. 1960, 60, 251–284.
48. Van Der Merwe, H.C. Sirovision: A proposed solution for the implementation of a digital geological mapping
system at AngloGold Ashanti’s Moab Khotsong mine. In Proceedings of the World Gold Conference, Misty
Hills, South Africa, 26–30 October 2009; Southern African Institute of Mining and Metallurgy: Johannesburg,
South Africa, 2009; pp. 183–191.
49. Lerm, C. Evaluation of chip sampling vs the gold analyser. J. Inst. Mine. Surv. S. Afr. 1994, 27, 77–87.
50. Cawood, F.T. Underground face sampling on narrow reefs—Sampling practice and its impact on the mine
call factor on Witwatersrand gold mines. J. Inst. Mine Surv. S. Afr. 2003, 31, 202–224.
51. Dominy, S.C. Grab sampling for underground gold mine grade control. J. South. Afr. Inst. Min. Metall. 2010,
110, 277–287.
52. Ross, I.T. Sampling in block cave mines. In Proceedings of the Sampling Conference, Perth, Australia, 21–22
August 2012; Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2012; pp. 81–85.
53. Arslan, N.; Ozkayhan, M. Diamond drilling at Efemcukuru gold mine, Turkey—A case study on the
benefits of timely and cost-effective drilling to short-term planning. In Proceedings of the International
Mining Geology Conference, Hobart, Australia, 20–22 September 2017; Australasian Institute of Mining and
Metallurgy: Melbourne, Australia, 2017; pp. 285–290.
54. Dominy, S.C. Sampling coarse gold-bearing mineralisation-developing effective protocols and a case study
from the Ballarat mine, Australia. In Proceedings of the World Conference on Sampling and Blending, Perth,
Australia, 9–11 May 2017; Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2017;
pp. 3–19.
55. McBeath, R.; Whitworth, D. What drilling density is optimal—A case study from the Golden Pig
mine. In Proceedings of the Underground Operators Conference, Townsville, Australia, 29–31 July 2002;
Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2002; pp. 33–37.
56. Clark, F.; Carswell, J.T. Using reverse circulation drilling to improve sampling in a complex underground
gold operation, Sunrise Dam mine, Western Australia. In Proceedings of the International Mining Geology
Conference, Adelaide, Australia, 18–20 August 2014; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2014; pp. 143–152.
57. Prouhet, J.-P. A quick method of sampling development headings. In Proceedings of the Ore Reserve
Estimation and Grade Control Conference, Ville d’Esterel, QC, Canada, 25–27 September 1967; CIM: Montreal,
QC, Canada, 1968; pp. 97–102.
58. Ovalle, A.; Magri, E.; Zuniga, P.; Pitard, F.F. Sampdrill. Quant. Geol. Geostat. 1997, 8, 767–779.
59. Johansen, G.F.; Dominy, S.C. Development of sampling and assaying protocols at the new Bendigo gold
project, Victoria, Australia. In Proceedings of the World Conference on Sampling and Blending, Sunshine
Coast, Australia, 10–12 May 2005; Australasian Institute of Mining and Metallurgy: Melbourne, Australia,
2005; pp. 175–183.
60. Roberts, L.S.; Dominy, S.C.; Nugus, M.J. Problems of sampling and assaying in lode-gold deposits—Case
studies from Australia and north America. In Proceedings of the International Mining Geology Conference,
Bendigo, Australia, 17–19 November 2003; Australasian Institute of Mining and Metallurgy: Melbourne,
Australia, 2004; pp. 387–400.
Minerals 2018, 8, 232 45 of 45
61. Dominy, S.C. Effects of sample mass on gravity recoverable gold test results in low grade ores. Appl. Earth Sci.
2014, 123, 234–242. [CrossRef]
62. Dominy, S.C.; Annels, A.E.; Johansen, G.F.; Cuffley, B.W. General considerations of sampling and assaying in
a coarse gold environment. Appl. Earth Sci. 2000, 109, 145–167. [CrossRef]
63. Hutchison, R. Geology and definition of the Wattle Dam coarse gold deposit. In Proceedings of the Gravity
Gold Conference, Ballarat, Australia, 21–22 September 2010; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2010; pp. 109–118.
64. Elder, D.J. Leachwell versus fire assay: Comments from assay quality assurance at Tarmoola gold mine.
In Proceedings of the International Mining Geology Conference, Coolum, Australia, 14–17 May 2000;
Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2009; pp. 43–52.
65. Pitard, F.F. Sampling correctness: A comprehensive guide. In Proceedings of the World Conference on
Sampling and Blending, Sunshine Coast, Australia, 10–12 May 2005; Australasian Institute of Mining and
Metallurgy: Melbourne, Australia, 2005; pp. 55–66.
66. Sims, D.A. Sampling practices at Vera Nancy gold mine. In Proceedings of the International Mining
Geology Conference, Coolum, Australia, 14–17 May 2000; Australasian Institute of Mining and Metallurgy:
Melbourne, Australia, 2009; pp. 35–42.
67. Benndorf, J.; Buxton, M.W.N. Sensor-based real-time resource model reconciliation for improved mine
production control—A conceptual framework. Min. Technol. 2016, 125, 54–64. [CrossRef]
68. Kurth, H. Measuring ore quality using real-time on-belt elemental analysis—More than grade control.
In Proceedings of the International Mining Geology Conference, Hobart, Australia, 20–22 September 2017;
Australasian Institute of Mining and Metallurgy: Melbourne, Australia, 2017; pp. 213–220.
69. Rifai, K.; Laflamme, M.; Constantin, M.; Vidal, F.; Sabsabi, M.; Blouin, A.; Bouchard, P.; Fytas, K.; Castello, M.;
Nguegang Kamwa, B. Analysis of gold in rock samples using laser-induced breakdown spectroscopy: Matrix
and heterogeneity effects. Spectrochim. Acta B 2017, 134, 33–41. [CrossRef]
70. Tickner, J.; Ganly, B.; Lovric, B.; O’Dwyer, J. Improving the sensitivity and accuracy of gamma activation
analysis for the rapid determination of gold in mineral ore. Appl. Radiat. Isot. 2017, 122, 28–36. [CrossRef]
[PubMed]
© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).