Stockpilesizereduction
Stockpilesizereduction
Stockpilesizereduction
ABSTRACT
Comminution circuits generally have a coarse ore stockpile after the primary crusher to
decouple the relatively uneven mined ore supply from the downstream comminution and
concentrator circuits which require a stable feed. Whilst this assists with evening-out
throughput disturbances, it can influence the size distribution presented to the downstream
comminution circuits through size-segregation effects. This is particularly important in
Autogenous (AG) and Semi-autogenous (SAG) circuits. What is not often considered,
however, is that the stockpile may also cause size reduction. Size distribution data from
stockpile feeds and stockpile products were collected from 13 different plants around the
world. Analysis of these data supports the assertion that size reduction occurs in the
stockpile. A simple power-based model is described which accurately predicts the degree of
size reduction measured. Worked examples are provided which show how the model can
be used to predict the stockpile product size given data on stockpile geometry and ore
hardness.
INTRODUCTION
It is usual in comminution circuits to have a coarse ore stockpile after the primary crusher to
decouple the effects on the comminution/concentrator circuits of the relatively uneven
supply of material from the mine to the primary crusher. As primary crusher availability is
generally lower than the downstream comminution circuits the stockpile additionally
enables feed to continue to be delivered to the downstream comminution circuits when the
primary crusher is under maintenance. The coarse ore stockpile assists with evening-out
throughput disturbances; however, it can influence the size distribution presented to the
downstream comminution circuits through size-segregation effects. In the case of
Autogenous (AG) and Semi-autogenous (SAG) circuits this is particularly important (Morrell
and Valery, 2001). However, what is rarely considered is that the stockpile not only can
cause size-segregation it may also cause size reduction. This paper presents size distribution
data from a number of coarse ore stockpiles in an attempt to measure the extent to which
size reduction may occur. A mathematical model to describe this size reduction is also
presented together with worked examples.
MEASUREMENT APPROACH
The stockpile feed (primary crusher product) and stockpile product (typically AG or SAG mill
feed) size distributions used in this analysis were collected during plant surveys conducted
for optimization projects. The collected size distributions cover several different
commodities from operations globally.
The size distributions were measured by conducting sieving of belt cut samples. Due to the
coarseness of the primary crushed material, a very large sample size is required to provide a
representative sample (Barbery, 1972; Gy,1976). This presents a challenge with regards to
the material’s handling and sieving due to the quantities of material that must be handled.
To minimize the sample size but achieve a representative sample, two separate samples
were collected from the conveyor belts:
1. All the material collected from a measured length of belt (L1).
2. All material coarser than 75 mm collected from a second adjacent and much longer
measured length of belt (L2).
The necessary belt length for the two samples, L1 and L2 (as shown in Figure 1), depends on
the belt width, speed, and loading. In general, several full 200 litre drums of rock were
collected, the material from L1 being kept separate from the material from L2. When
collecting the +75 mm sample, all rocks that look big enough were initially checked using a
hand-held trial sieve and selected if they are of the right size. Later, full sieve analysis
removed any -75 mm particles that were accidentally collected. The data from sieve analysis
of the L1 and L2 samples were analysed in terms of mass/belt length then combined into a
single size distribution. An example of belt sampling is shown in Figure 2. A detailed
description of the procedure can be found in Napier-Munn et al (1996).
Figure 2 - Example of Belt Sampling - L1 (all material) on left and L2 (+ 75 mm) on right.
RESULTS
A summary of the results obtained is given in Table 1. In all cases the stockpile product 80%
passing size (D80 ) is finer than the stockpile feed D80. Differences between individual results
vary significantly and are believed to be due to differences in rock hardness and stockpile
height. This former effect can be seen in Table 1 in which the data has been grouped into
those with Mic (SMC Test® crusher hardness parameter) values < 5 kWh/t and those with
values > 5 kWh/t.It can be seen that the < 5 kWh/t group has a D80 difference between
stockpile feed and stockpile product almost double that of the > 5 kWh/t group. Details of
the individual size distributions are given in Appendix 1. To help visualise the consistent
trend seen at all sites the feed and product size distributions have all been combined and
are presented in Figure 3 in terms of the overall mean stockpile feed and product.
Table 1 – Results Summary
Where:
W = predicted circuit net specific energy (kWh/t)
Mi = hardness parameter (kWh/t)
𝑓𝑓 (𝑥𝑥 ) = −(0.295 + 𝑥𝑥 ⁄1000000)
x = 80% passing size in microns
P80 = 80% passing size of the product
F80 = 80% passing size of the feed
As applied to crushing circuits the form of equation 1 is as follows:
(𝑃𝑃80) (𝐹𝐹80)
𝑊𝑊𝑐𝑐 = 𝑀𝑀𝑖𝑖𝑖𝑖 × 𝑆𝑆𝑐𝑐 × 𝑘𝑘3 × 4 × �𝑃𝑃80 𝑓𝑓 − 𝐹𝐹80 𝑓𝑓 � (2)
Where:
Wc = predicted net specific energy of the crushing circuit (kWh/t)
Mic = SMC Test® crusher hardness parameter (kWh/t)
Sc = coarse feed factor
= 55 × (𝐹𝐹80 × 𝑃𝑃80)−0.2 : (0.5 < Sc < 1)
k3 = open/closed circuit factor
= 1.19 for open circuit
= 1.0 for closed circuit
As the rock falls it progressively converts this potential energy to kinetic energy, until on
impact some is absorbed by the rock. If what is absorbed exceeds the yield stress, breakage
will occur. How much kinetic energy is available on impact will depend on how much is lost
to drag resistance, whilst how much is subsequently absorbed by the rock depends on the
nature of the impact surface. If the impact surface is hard, rigid, flat and perpendicular to
the impact velocity (say a flat concrete floor), the proportion of energy absorbed will be
relatively large. However, if the impact surface comprises, say, an inclined loose bed of
material the proportion will be lower. This latter situation is more likely to be what rocks
falling onto a stockpile will encounter. To accommodate drag losses and the influence of
the impact surface conditions a fitted lumped “efficiency” factor, keff, is incorporated into
equation 3. Its value can vary in the range 0-1. Hence the potential energy absorbed by a
rock, PEabsorbed, is given by:
𝑃𝑃𝑃𝑃𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 = 𝑚𝑚 × 𝑔𝑔 × ℎ × 𝑘𝑘𝑒𝑒𝑒𝑒𝑒𝑒 (4)
There are also other factors which can influence the product size including segregation and
attrition and abrasion breakage which are not accounted for in the model; these are also
lumped into the keff efficiency factor. The proposed model is simple by necessity, based on
impact breakage only, using data that can practically and reliably be collected. Regardless of
the necessary simplification and assumptions the model achieves good accuracy as will be
demonstrated in the results presented later in Table 2.
The value of h (height of drop) in equation 4 will be related to the maximum stockpile height
and the operating level. Assuming that the operating level provides an impact surface
which has a height relative to an empty stockpile of hop and representing the maximum
stockpile height by hmax, the ratio hop/hmax can then be defined as a relative operating level
Lop where Lop varies in the range 0-1, 0 being the condition of an empty stockpile and 1
when it is full. Equation 4 can now be written:
𝑃𝑃𝑃𝑃𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 = 𝑚𝑚 × 𝑔𝑔 × ℎ𝑚𝑚𝑚𝑚𝑚𝑚 × (1 − 𝐿𝐿𝑜𝑜𝑜𝑜 ) × 𝑘𝑘𝑒𝑒𝑒𝑒𝑒𝑒 (5)
As specific energy is energy/mass then, using equation 5, the specific energy absorbed by
rocks falling on to a stockpile, Wsp, can be represented by:
𝑊𝑊𝑠𝑠𝑠𝑠 = 𝑔𝑔 × ℎ𝑚𝑚𝑚𝑚𝑚𝑚 × (1 − 𝐿𝐿𝑜𝑜𝑜𝑜 ) × 𝑘𝑘𝑒𝑒𝑒𝑒𝑒𝑒 (6)
The data in Table 1 were not originally collected with the view to developing a stockpile size
reduction model and hence there is limited information concerning stockpile heights (hmax)
and operating levels (Lop). To overcome this, stockpile heights were estimated based on
assuming a typical live capacity in tonnes, Cl. From a design perspective this is often
calculated using the number of hours that the live capacity is required to keep the
downstream comminution operating at its nominal throughput capacity during periods
where there is no material being fed to the stockpile. Hence, for example, in the case of
Newcrest Cadia operation (Dunne et al, 2001), it had a stockpile live capacity designed to
last 20 hours, which at its original nominal comminution circuit throughput of 2000 t/hr
gives a live capacity of 40,000 tonnes. In contrast, Newmont Boddington operation (Hart et
al, 2011) has a live capacity designed to last 8 hours, though as its nominal comminution
circuit throughput is 5000 t/hr it also has a live capacity of 40,000 tonnes. With reference to
Table 1, the mean live capacity in terms of hours, for those sites where it was known, was
14. For those sites where it was not known, 14 hours was assumed. Multiplying the live
capacity hours by the throughput of each comminution circuit, the live capacity of each
stockpile in terms of tonnes was estimated. From simple geometry it can be shown that for
a conical stockpile with a central offtake, the live capacity is 25% of the stockpile total
volume (see Appendix 2). This figure agrees with that quoted by Zamarano (2006). Hence
the total stockpile capacity is 4 times the live capacity:
𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇𝑇 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 (𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) = 4 × 𝐶𝐶𝑙𝑙 (7)
Assuming a void fraction of 0.3 and using the measured rock sg, the volume of the stockpile,
Vsp, is given by:
𝑉𝑉𝑠𝑠𝑠𝑠 (𝑐𝑐𝑐𝑐. 𝑚𝑚) = 4/(1 − 0.3) × 𝐶𝐶𝑙𝑙 /𝑠𝑠𝑠𝑠 (8)
Having defined hmax, then with reference to equation 6, this leaves Lop and keff as unknowns.
It is considered good practice to operate with high stockpile level (70% or higher) as much as
possible to avoid issues with segregation (Morrell and Valery, 2001). Therefore, it was
decided to assume that Lop was on average equal to 0.5 which is equivalent to assuming that
on average the stockpiles in Table 1 were operated with 87.5% of the live capacity available
(see Appendix 2). The parameter keff was initially left as a fitted parameter. Equation 6 can
now be reduced to:
1/3
𝐶𝐶
𝑊𝑊𝑠𝑠𝑠𝑠 (𝑘𝑘𝑘𝑘ℎ⁄𝑡𝑡) = 0.00204 × � 𝑙𝑙�𝑠𝑠𝑠𝑠� × 𝑘𝑘𝑒𝑒𝑒𝑒𝑒𝑒 (13)
With reference to equation 2, P80 is the 80% passing size of the stockpile product, F80 is the
80% passing size of the stockpile feed, Wsp replaces Wc, Mic is directly measured by
conducting SMC Tests on representative samples of stockpile feed, whilst k3 is set to the
open circuit value of 1.19 hence giving equation 14.
(𝑃𝑃80) (𝐹𝐹80)
𝑊𝑊𝑠𝑠𝑠𝑠 = 𝑀𝑀𝑖𝑖𝑖𝑖 × 𝑆𝑆𝑐𝑐 × 1.19 × 4 × �𝑃𝑃80 𝑓𝑓 − 𝐹𝐹80 𝑓𝑓 � (14)
Given that the objective is to rearrange equation 14 such that it can be solved for P80, the
fact that the exponent of the P80 term is a function of the P80 and the Sc term also contains
the P80, makes an analytical solution impossible. However, numerical approximation
techniques were found that provided very accurate estimates and gave rise to equations 15
and 16. Alternatively, the excel “solver” function can be used to search for a P80 value
which satisfies equation 14.
𝑊𝑊𝑠𝑠𝑠𝑠 𝐹𝐹80
−�0.295+ �
𝑃𝑃80 = −65500 × ��ln �(
4×𝑀𝑀𝑖𝑖𝑖𝑖×1.19×𝑆𝑆𝑐𝑐 )
+ 𝐹𝐹80 1000000 � � + 3� (15)
150
100
50
0
0 50 100 150 200 250
Observed P80 (mm)
Figure 4 – Stockpile Size Reduction Model Accuracy (Observed vs Predicted Stockpile P80)
Using 0.7 for keff in equation 13 and substituting for Wsp in equation 15, equation 15 can now
be reduced to:
𝐶𝐶 1/3
0.0003×� 𝑙𝑙�𝑠𝑠𝑠𝑠� −�0.295+
𝐹𝐹80
�
𝑃𝑃80 = −65500 × ��ln � (𝑀𝑀𝑖𝑖𝑖𝑖 ×𝑆𝑆𝑐𝑐)
+ 𝐹𝐹80 1000000 � � + 3� (17)
Equations 16 and 17 are best used in situations such as greenfield designs where the
stockpile height has yet to be finalised. However, in cases where the stockpile height is
known then by combining equation 6 with equation 15, the resultant equation 18 should be
used instead of equation 17.
𝐹𝐹80
0.0002 ×ℎ𝑚𝑚𝑚𝑚𝑚𝑚 −�0.295+ �
𝑃𝑃80 = −65500 × ��ln � (𝑀𝑀𝑖𝑖𝑖𝑖 ×𝑆𝑆𝑐𝑐 )
+ 𝐹𝐹80 1000000 � � + 3� (18)
CONCLUSIONS
Analysis of size distribution data from coarse ore stockpile feed and product size
distributions supports the assertion that size reduction occurs in stockpiles. A simple
power-based model assuming impact breakage to be the sole source of size reduction was
developed. It was found to predict the degree of size reduction, measured in terms of the
80% passing sizes, to within an accuracy of just over 10% at the 90% confidence level.
Worked examples which show how the model can be used to predict the stockpile product
size given data on stockpile geometry and ore hardness are provided in Appendix 3.
On average the size distribution data collected from the 13 plants included in this study
indicate that the stockpile products may contain slightly more fine material than a simple
impact comminution model would suggest. Whilst biases in the way that the data were
collected cannot be ruled out, it is possible that the apparent additional finer material may
have been the result of abrasion and attrition breakage not accounted for in the model.
DEM studies of stockpile dynamics may be able to shed light on this phenomenon.
REFERENCES
Barbery, G. 1972. Derivation of a formula to estimate the mass of a sample for size analysis,
Trans AusIMM, 81, (784), March, C49-C51.
Dunne, R., Morrell, S., Lane, G., Valery, W. and Hart, S. 2001. Design of the 40 foot diameter
SAG mill installed at the Cadia gold copper mine. Proc. Conf. SAG 2001, CIMM, Vancouver,
Canada
Froehlich, D.C. 2011. Mass Angle of Repose of Open-Graded Rock Riprap. Journal of
Irrigation and Drainage Engineering, July, pp 454-461. American Society of Civil Engineers
GMG Group. 2021. The Morrell Method to Determine the Efficiency of Industrial Grinding
Circuits. https://gmggroup.org/guidelines-and-publications/morrell-method-to-determine-
the-efficiency-of-industrial-grinding-circuits/
Gy, P.M. 1976, The sampling of particulate materials – a general theory. Symp of Sampling
Particles in the Min Ind, AusIMM, Melbourne, Australia.
Hart, S., Rees,T., Tavani, S., Valery, W. and Jankovic, A. 2011. Process Integration and
Optimisation of the Boddington HPGR Circuit. Proc. Conf. SAG 2011, CIMM, Vancouver,
Canada
Lewis-Gray, A. and Rasmussen, G. 2019. Design and Operation of the Vertical Shaft Impactor
in Fine Crushing Applications. Proc. SAG 2019 Conference, CIMM, Vancouver, Canada.
Morrell, S. 2004. An alternative energy-size relationship to that proposed by bond for the
design and optimisation of grinding circuits, Int. J. Miner. Process. 74:133–141.
Morrell, S. and Valery, W. 2001. Influence of feed size on AG/SAG mill performance. Proc
Conf. SAG 2001, CIMM, Vancouver, Canada.
Napier-Munn, T.J., Morrell, S., Morrison R.D. and Kojovic T.M. 1996. Mineral Comminution
Circuits Their Operation and Optimisation. JKMRC, University of Queensland
Zamarano, S. 2006. Coarse Ore Stockpiles. Bulk Solid Handling, Vol. 26, No. 4, pp 246-251.
APPENDIX 1 – Measured Size Distributions of the Stockpile Feed and Product
Input data:
Stockpile height (hmax) = 32m
Mic = 7.6 kWh/t
Expected primary crusher product (F80) = 152mm
Calculations:
0.0002 × 32 152000
𝑃𝑃80 = −65500 × ��ln � + 152000−�0.295+1000000� � � + 3�
(7.6 × 0.5)
= −65500 × (−5.03416 + 3)
= 133238 microns
Hence the stockpile product (AG/SAG mill feed) is expected to have on average a P80 of
133mm
B. Stockpile Height Unknown
Equations:
𝐶𝐶 1/3
0.0003×� 𝑙𝑙�𝑠𝑠𝑠𝑠� −�0.295+
𝐹𝐹80
�
𝑃𝑃80 = −65500 × ��ln � (𝑀𝑀𝑖𝑖𝑖𝑖 ×𝑆𝑆𝑐𝑐)
+ 𝐹𝐹80 1000000 � � + 3�
Input data:
Nominal comminution circuit throughput = 2500 tph
Required time for live stockpile to last = 12 hours
Mic = 4.2 kWh/t
Sg = 2.7
Expected primary crusher product (F80) = 120mm
Calculations:
𝐿𝐿𝐿𝐿𝐿𝐿𝐿𝐿 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑜𝑜𝑜𝑜 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑝𝑝𝑖𝑖𝑖𝑖𝑖𝑖 (𝐶𝐶𝑙𝑙 ) = 2500 × 12 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
Hence:
𝐶𝐶𝑙𝑙 = 30000 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡
1/3
0.0003�30000�2.7� 120000
𝑃𝑃80 = −65500 × ��ln � (4.2×0.534)
+ 120000−�0.295+1000000�� � + 3�
= −65500 × (−4.52935 + 3)
= 100172 microns
Hence the stockpile product (AG/SAG mill feed) is expected to have on average a P80 of
100mm