Sampling
Sampling
Where: M= Minimum sample weight needed, grams
W= Weight of the entire lot being sampled, grams
C= Sampling constant for the material being sampled, g/cm
3
d
max
= Dimensions of the largest pieces in the sample, cm
s= value of the standard deviation that will be needed to give the
desired level of assurance (assay units, such as % wt, ounces
per ton, etc.)
When W is much larger than M, this equation can be simplified to:
2
3
max
s
Cd
M =
To calculate the desired value for s, first select the desired certainty range (for example, you
might want a copper assay to be accurate to within +/- 0.1% copper), then divide this value by
the mean probable assay value (so for a 5% copper sample, you get 0.1/5 = 0.02). Then, select
the number of standard deviations that will be needed to give the desired certainty. For a
probability of 95%, this is 1.96 standard deviations, so 1.96s=0.02, and so s=0.0102. If the
desired certainty is 99%, then 2.576 standard deviations are needed, and so for this case s =
0.02/2.576 = 0.0078
Calculating C is more complicated, as follows:
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
fglm C =
where: f= Shape factor
g= Size distribution factor
l= Liberation factor
m= Composition factor
The shape factor, f, is taken to be equal to 0.5 for most common natural materials, but in cases
where the material of interest is highly flattened and elongated, the shape factor is set at 0.2.
The size distribution factor, g, is determined by estimating the 95% passing size (d) and the 5%
passing size (d), and using the following guidelines:
If d/d > 4, then g=0.25;
If 2<d/d<4, then g=0.5;
If d/d < 2, then g=0.75;
If d/d = 1 (monosized particles), then g=1.
The liberation factor, l, is a measure of the degree of dispersion of the valuable material through
the bulk, and of the homogeneity of the material. It is calculated from the expression:
d
L
l =
where: L= the size where the values are essentially completely liberated
(grain size), cm.
The composition factor, m, is calculated from the formula:
( ) ( ) at r a
a
a
m +
= 1
1
where: r= specific gravity of the valuable component
t= specific gravity of the remainder of the material
a= fractional average assay of the valuable substance (for example, if an
iron-bearing slag is being analyzed, and the iron is present as Fe
2
O
3
, use the
fraction of the material that is Fe
2
O
3
, not the fraction that is iron, and it
would be expressed as a decimal fraction, not as a percent)
Calculation with Incomplete Information
For a truly unknown sample, most of the parameters that go into calculating C will be unknown,
and may even be impossible to guess with any accuracy. In this case, make the following
conservative assumptions:
f=0.5 (normal blocky particles);
g=0.75 (narrow size distribution. Use g=1 if the sample is obviously monosized);
l = 1 (grains are as large as the particles)
The value of m will still need to be calculated, based on your best estimate of the assay of the
sample and the densities of the components of interest. This initial sample will be larger than
necessary, but can be analyzed and used to refine the calculations for collection of any
subsequent samples.
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Example Calculation:
Bulk Material Parameters:
Material of Interest: CuFeS
2
in a silica matrix, 1.5% Cu (4.3318% CuFeS
2
); Top Size= 1.5 cm;
CuFeS
2
grain size= 0.01 cm
Desired Sampling Accuracy: +/- 0.02% Cu, certainty of 0.99 (2.576 standard deviations)
CuFeS
2
specific gravity= 4.2; Overall specific gravity: 2.8; Broad size distribution
Calculations:
(Calculations were done using a spreadsheet, and were carried out to more digits than are shown
here. Carrying out the calculations with the values as shown will give slightly different results
due to roundoff variations.)
f = 0.5 (regular shape)
g = 0.25 (d/d > 4 for a broad size distribution)
081650 . 0
5 . 1
01 . 0
= = =
d
L
l
a = 0.043318, r = 4.2, and t = 2.8, and so:
( ) ( ) ( ) 4177 . 91 8 . 2 043318 . 0 2 . 4 43318 . 0 1
43318 . 0
043318 . 0 1
= +
= m
( )( )( )( ) 933028 . 0 4177 . 91 081650 . 0 25 . 0 5 . 0 = = = fglm C
01333 . 0
5 . 1
02 . 0
576 . 2 = = s , so s = 0.0051758
( )
2
3
2
3
0051758 . 0
5 . 1 933028 . 0
= =
s
Cd
M =117,539 grams of sample needed for 99% certainty that the assay
will be within +/- 0.02% Cu of the true assay.
Since the original data was only carried to 2 significant digits, this should be rounded to 120,000
grams of sample, or 120 kg.
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Sampling Calculations using Gaudins Method
Reference:
A. M. Gaudin, Principles of Mineral Dressing, McGraw-Hill, 1939
Use:
This method is limited to materials where the species of interest is a small fraction (a few percent
or less) of the total volume. It does not allow for selecting the degree of certainty desired in the
results, and it is not as generally useful or as precise as Gys method. It would probably be
most useful for calculating the size of a preliminary sample where minimal information about the
material characteristics is available.
Basic Equations:
2
45 . 0
y
x
n =
3
6
a
n
s
=
n
n
S
=
where: x = volume fraction occupied by the species of interest
y = variation in volume fraction that would correspond to the
maximum allowable error
n = number of particles that will be needed in the final sample
n= estimated number of particles per gram
s
= density of the bulk sample, in grams/cm
3
a = top size of the material being sampled, in cm
S = total sample weight needed, in grams.
Example Calculation:
Bulk Material Parameters:
Same material as for the previous example:
CuFeS
2
, 1.5% Cu (4.3318% CuFeS
2
); Top Size= 1.5 cm;
Desired Sampling Accuracy: +/- 0.02% Cu (+/- 0.0578% CuFeS
2
)
CuFeS
2
specific gravity= 4.2; Overall specific gravity: 2.8
a = 1.5 cm
s
= 2.8 gm/cm
3
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
(Calculations were done using a spreadsheet, and were carried out to more digits than are shown
here. Carrying out the calculations with the values as shown will give slightly different results
due to roundoff variations.)
028879 . 0
8 . 2
100
2 . 4
3318 . 4
=
= x
00038505 . 0
8 . 2
100
2 . 4
0578 . 0
=
= y
87650
00038533 . 0
028879 . 0
45 . 0
2
=
= n
( )( )
63492 . 0
5 . 1 8 . 2
6
3
= = n
S = (87650)/(0.63492) = 138,050 grams
(compared to 117,599 grams using Gys method)
Again, this should be rounded to 2 significant digits, and so 140,000 grams, or 140 kg, of sample
is required.
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Comparison of Gaudins sampling equations and Gys
sampling equation
Gys equation is a general-purpose sampling equation,
2
3
max
s
Cd
M = , where M = minimum sample mass needed, d
max
= maximum particle size present,
and s
2
= allowable sampling variance. C is a sampling constant, which is calculated for any
given material based on particle shape, size distribution, degree of liberation, density, and
mineralogy.
Gaudins equation is specifically intended for use in calculating sample sizes for precious
metals, and makes a number of simplifying assumptions based on the fact that precious metals
make up a minute fraction of the mass of the ore (typically only a few grams per metric ton).
The equation is normally written as two parts: a calculation of n, the number of particles
needed in the sample, and of n, the expected number of particles per gram:
=
2
45 . 0
y
x
n , where x = ore grade expressed as volume fraction, and y = allowable sampling
error expressed as volume fraction.
3
6
a
n
s
= , where
s
= bulk ore density, and a = maximum particle size present.
The mass of the sample needed, S, is then ( )
2
3
3
2
)( 075 . 0
6
45 . 0
y
a
x
a
y
x
n
n
S
s
s
=
Comparison: We can put Gaudins equation into the same terms as Gys equation by using the
same symbols for equivalent values in each: a = d, S = M. The value of y (a volume fraction)
can be converted into (a weight fraction) by multiplying it by a constant, k, that depends on the
relative densities of the ore and the valuable minerals, and so y = /k. If we make these
substitutions, then Gaudins equation becomes ( )
2
3
2
)( 075 . 0
d
k x M
s
=
If we now take (0.075)(x
s
k
2
) = C, we see that Gaudins equation reduces to same form as Gys
equation. Gaudins equation is therefore simply a special case of Gys equation.
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Sampling Calculations Using Hassialis Method for
Increment Sampling
Reference:
Taggart, Handbook of Mineral Dressing, Wiley, 1945, Section 19
Use:
This method is intended for designing a mechanical sampling system that will account for
normal variations in the characteristics of the process stream being sampled. The sample cutter
is used to collect a series of test increments over a particular size interval, and then each
increment is assayed separately. The variations in the assays between increments are then used
to calculate the number of increments that the sample cutter needs to collect over that span of
time to give the desired level of certainty in the results.
Basic Equation:
2
=
z x
h
n
s
where n = number of increments that will be required to reach the
desired level of certainty
h = number of standard deviations that correspond to the desired
level of certainty (for example, 95% certainty corresponds to
1.96 standard deviations, see Table 2)
s
= experimental standard deviation of the assays of the test
increments collected by the mechanical sampler
x = experimental mean of the assays of the test increments
collected by the mechanical sampler
z = allowable sampling error, expressed as a fraction of the assay
value (for example, if the final sample is 15% Fe and z = .04,
then the allowable variation in analysis of the final sample
will be (0.04)(15) = 0.60% Fe. If the final sample is 30%Fe,
then the allowable variation would be (0.04)(30) = 1.2%Fe).
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Example Calculation
(Calculations were done using a spreadsheet, and were carried out to more digits than are shown
here. Carrying out the calculations with the values as shown will give slightly different results
due to roundoff variations.)
The following 21 test increment assays were obtained from increments collected at 15 minute
intervals, and took a total of 5 hours to collect. Values are weight % Fe.
10.6 10.9 10.2 10.8 9.4 9.6 10.0 10.1 9.9 9.8 9.5
9.2 9.1 9.7 8.8 10.2 10.4 10.2 9.8 9.2 9.3
Calculate the number of increments that must be collected over a 5 hour interval so that the
probability of the assay of the final composite sample being accurate to within 2% of the true
assay will be 99%
Standard Deviation:
s
= 0.57234
Mean: x = 9.8428 % Iron
Allowable variation: z = 0.02 (that is, the analysis of the sample will be accurate to within 0.02
of the mean, and so the analysis will be 9.8428 +/- 0.1968%Fe)
Number of standard deviations needed: h = 2.576
Table 2: Values of h corresponding to various levels
of probability that a particular sample of a lot will
have an assay value within the desired error limits for
the analysis.
Probability h Probability h
0.90 1.645 0.97 2.170
0.91 1.705 0.98 2.326
0.92 1.750 0.99 2.576
0.93 1.812 0.999 3.291
0.94 1.881 0.9999 3.890
0.95 1.960 0.99999 4.417
0.96 2.054 0.999999 4.892
CM3820.2009.Sampling.doc Last Updated: January 14, 2009 4:42 PM
Number of increments needed to make a composite sample of the desired degree of accuracy:
( )( )
( )( )
09128 . 56
02 . 0 8428 . 9
57234 . 0 576 . 2
2 2
=
=
z x
h
n
s
Since it is not possible to collect a fraction of an increment, this should be rounded up to 57
increments per composite sample, and these increments should be collected over a 5-hour
interval.