Difference between revisions of "Model:Bond/Barratt SABC Models"

From SAGMILLING.COM
Jump to: navigation, search
(Formulae)
(Optional parameters)
 
(26 intermediate revisions by the same user not shown)
Line 2: Line 2:
 
[[Category:Specific Energy Models]]
 
[[Category:Specific Energy Models]]
 
[[Category: Bond/Barratt Model]]
 
[[Category: Bond/Barratt Model]]
  +
[[Category: P80 adjustment]]
  +
[[Category: Monte Carlo]]
 
== Bond/Barratt Specific Energy Consumption SABC Model ==
 
== Bond/Barratt Specific Energy Consumption SABC Model ==
   
Line 12: Line 14:
 
* [[Testwork: Bond rod mill work index|Bond rod mill work index]]
 
* [[Testwork: Bond rod mill work index|Bond rod mill work index]]
 
* [[Testwork: Bond crushing work index|Bond crushing (low-energy impact) work index]]
 
* [[Testwork: Bond crushing work index|Bond crushing (low-energy impact) work index]]
  +
  +
===Required parameters===
  +
* '''F<sub>80</sub>, µm''' is the 80% passing size of the fresh feed to the circuit (expected to be a Bond-compatible size distribution).
  +
* '''P<sub>80</sub>, µm''' is the 80% passing size of the circuit product (expected to be a Bond-compatible size distribution).
  +
* '''Availability''', expressed as a decimal (0.90 = 90% availability) is used to convert t/h to t/d.
  +
* Ball mill '''Operating strategy''' when circuit is SAG-limited (note, fixed speed can only be ''vary P<sub>80</sub>'')
  +
  +
===Optional parameters===
  +
* '''Description''' and '''Comment''' are optional text fields
  +
* '''Maximum t/h limit''' is a t/h throughput limit above which the mills will turn down by, for example, grinding out.
  +
* '''T<sub>80</sub> min''' and '''T<sub>80</sub> min''' override the transfer size restrictions built into the model
  +
* [[Ball mill work index adjustment]] used to adjust Wi<sub>BM</sub> for different P<sub>80</sub> sizes.
  +
** '''exponent (a)''', the fitted Hukki exponent to the adjustment equation (enter as a positive number)
  +
* '''Essbm calibration factor''' (CF), model tuning factor for E<sub>total</sub>
  +
* '''Easag crushing circuit product''', model tuning factor that dictates where the Barratt Easag equation switches from using the crushing Wi to the rod mill Wi
  +
* '''Easag RM circuit product''', model tuning factor that dictates where the Barratt Easag equation switches from using the rod mill Wi to the ball mill Wi
  +
* '''Essbm crushing circuit product''', model tuning factor that dictates where the Bond Essbm equation (the basis for the Barratt Etotal equation) switches from using the crushing Wi to the rod mill Wi
   
 
===Formulae===
 
===Formulae===
  +
<math>
<math> E_{ssbm} = {{BondWiEquation|Wi=Wi_{C}|F80=F_{80}|P80=9400 }} + {{BondWiEquation|Wi=Wi_{RM}|F80=9400|P80=2100 }}\times EF_4^{RM} + {{BondWiEquation|Wi=Wi_{BM}|F80=2100|P80=P_{80} }} \times EF_4^{BM} \times EF_5</math> <sup>[[Bibliography:_Specific_energy_consumption_models|Rowland, 2006]]</sup>
 
  +
E_{ssbm} = Wi_{C} \times \left ( \tfrac {10}{\sqrt{ Pc_{ssbm} }} - \tfrac {10}{\sqrt{ F_{80} }} \right ) + Wi_{RM} \times \left ( \tfrac {10}{\sqrt{ Prm }} - \tfrac {10}{\sqrt{ Pc_{ssbm} }} \right )\times EF_4^{RM} + Wi_{BM} \times \left ( \tfrac {10}{\sqrt{ P_{80} }} - \tfrac {10}{\sqrt{ Prm }} \right ) \times EF_4^{BM} \times EF_5</math>
  +
<sup>[[Bibliography:_Specific_energy_consumption_models|Rowland, 2006]]</sup>
   
  +
where
<math> E_{total} = E_{ssbm} \times (1+CF) </math>
 
  +
* Pc<sub>ssbm</sub> is the calibration value entered for "Essbm crushing circuit product" (default: 9400)
  +
* Prm is the calibration value entered for "Easag RM circuit product" (default: 2100)
   
<math> {{BarrattEsagEquation}} </math>
+
<math> E_{total} = E_{ssbm} \times (CF) </math>
  +
  +
  +
<math>\begin{align}
  +
E_{SAG} = \Big[ & Wi_{C} \times \left ( \tfrac {10}{\sqrt{ P_{C} }} - \tfrac {10}{\sqrt{ F_{80} }} \right ) + Wi_{RM} \times \left ( \tfrac {10}{\sqrt{ P_{R} }} - \tfrac {10}{\sqrt{ P_{C} }} \right ) \times EF_4^{RM} \\
  +
& + Wi_{BM} \times \left ( \tfrac {10}{\sqrt{110} } - \tfrac {10}{\sqrt{ P_{R} }} \right ) \times EF_4^{BM} \Big] \times 1.25 - Wi_{BM} \times \left ( \tfrac {10}{\sqrt{ {110} }} - \tfrac {10}{\sqrt{ T_{80} }} \right )
  +
\end{align}
  +
</math>
 
<sup>[[Bibliography:_Specific_energy_consumption_models|Barratt, 1989]]</sup>
 
<sup>[[Bibliography:_Specific_energy_consumption_models|Barratt, 1989]]</sup>
  +
  +
where
  +
* P<sub>C</sub> is the calibration value entered for "Easag crushing circuit product" (default: 18850)
   
 
<math> E_{bm}= E_{total} - E_{sag}</math>
 
<math> E_{bm}= E_{total} - E_{sag}</math>
Line 31: Line 64:
 
<math>F_{O}= 16000 \times \sqrt{\frac{13}{Max(Wi_{RM}, Wi_{C}) \div 1.1023} } </math>
 
<math>F_{O}= 16000 \times \sqrt{\frac{13}{Max(Wi_{RM}, Wi_{C}) \div 1.1023} } </math>
   
The EF<sub>4</sub><sup>RM</sup> is calculated using only the rod mill work index.
+
The EF<sub>4</sub><sup>RM</sup> is calculated greater of the crushing work index or the rod mill work index.
   
<math>EF_{4}^{RM} = \frac{RR+(Wi_{RM} \div 1.1023 -7) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}</math>
+
<math>EF_{4}^{RM} = \frac{RR+ \bigl(\frac{Max(Wi_C, Wi_{RM})}{1.1023} -7 \bigr) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}</math>
   
 
The ball mill oversize feed factor (EF<sub>4</sub><sup>BM</sup>) is always calculated with the rod mill work index for the optimum feed size and the ball mill work index for EF<sub>4</sub><sup>BM</sup>.
 
The ball mill oversize feed factor (EF<sub>4</sub><sup>BM</sup>) is always calculated with the rod mill work index for the optimum feed size and the ball mill work index for EF<sub>4</sub><sup>BM</sup>.
Line 39: Line 72:
 
<math>F_{O}= 4000 \times \sqrt{\frac{13}{Wi_{RM} \div 1.1023} } </math>
 
<math>F_{O}= 4000 \times \sqrt{\frac{13}{Wi_{RM} \div 1.1023} } </math>
   
<math>EF_{4}^{RM} = \frac{RR+(Wi_{BM} \div 1.1023 -7) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}</math>
+
<math>EF_{4}^{BM} = \frac{RR+ \bigl( \frac{Wi_{BM}}{1.1023} -7 \bigr) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}</math>
  +
  +
The model assumes that the calculated value of E<sub>ssbm</sub> + the specified contingency is greater than the sum of E<sub>sag</sub> + E<sub>peb</sub>; otherwise E<sub>bm</sub> will be negative. The model does a check to see if E<sub>bm</sub> calculated is negative (or close to zero) and, if so, will substitute the [[Model:RawBondModel|Raw Bond Model]] equation estimate for E<sub>bm</sub>.
  +
  +
== Recommended Usage ==
  +
  +
This model should be the lower end-member of the range of Bond/Barratt style models. It is useful when interpreting pilot plant results to determine where the pilot plant results fall relative to the spectrum of models. It can be used as a realistic model for design work in softer ore types like copper porphyries. The Optimized Bond/Barratt model is best used for '''low competency ores''', such as those found in the Western Cordillera of North America or the Andes, where the impact crushing work index is similar to, or less than the rod mill and ball mill work index values. These ores tend to demonstrate a phantom cyclone effect.
  +
  +
[[file:WorkIndex-CompentIncompetent.png|alt=Diagram of work index by size|Work index by size examples of competent and incompetent ores]]
   
The model assumes that the calculated value of E<sub>ssbm</sub> + the specified contingency is greater than the sum of E<sub>sag</sub> + E<sub>peb</sub>; otherwise E<sub>bm</sub> will be negative. The model does a check to see if E<sub>bm</sub> calculated is negative and, if so, will substitute the [[Raw Bond Model]] equation estimate for E<sub>bm</sub>.
 
   
 
=== Discussion, pebble crushing ===
 
=== Discussion, pebble crushing ===
 
The pebble crusher is assumed to be part of the SAG mill for the purposes of this calculation. E<sub>peb</sub> is calculated as the power drawn by the pebble crusher divided by the circuit throughput. This E<sub>peb</sub> is then deducted from E<sub>sag</sub>.
 
The pebble crusher is assumed to be part of the SAG mill for the purposes of this calculation. E<sub>peb</sub> is calculated as the power drawn by the pebble crusher divided by the circuit throughput. This E<sub>peb</sub> is then deducted from E<sub>sag</sub>.
   
The total circuit specific energy consumption is generally higher when no pebble crusher is present in a circuit. The rule of thumb is apply the CF value of 0.15 for SAB circuits (5% higher overall specific energy consumption).
+
The total circuit specific energy consumption is generally higher when no pebble crusher is present in a circuit. The rule of thumb is apply the CF value of 1.15 for SAB circuits (5% higher overall specific energy consumption) <sup>[[Bibliography:_Specific_energy_consumption_models#SPI.2FMinnovex_Method|Bennett, Dobby & Kosick, SAG 2001]]</sup>. Published data from Cadia confirms a 5% difference in throughput (and therefore, specific energy) when the pebble crusher is fully utilized.<sup>Hart, Valery, Clements, Reed, Song & Dunne, SAG 2001</sup>.
   
 
=== Discussion, Rowland EF factors ===
 
=== Discussion, Rowland EF factors ===
The model is only interested in the EF factors that are properties of the rock, namely the coarse feed EF4 factor and the fine product EF5 factor. The EF factors that correspond to machine characteristics are not included (EF<sub>1</sub>, EF<sub>2</sub>, EF<sub>3</sub>, EF<sub>6</sub>, EF<sub>7</sub> and EF<sub>8</sub>) as they are intended to capture the inefficiency of certain types of milling circuits. The assumption in this model is that the E<sub>ssbm</sub> is the "maximum efficiency" possible for a grinding circuit, and therefore one can neglect the machine characteristics.
+
The model is only interested in the EF factors that are properties of the rock, namely the coarse feed EF4 factor and the fine product EF5 factor. The EF factors that correspond to machine characteristics are not included (EF<sub>1</sub>, EF<sub>2</sub>, EF<sub>3</sub>, EF<sub>6</sub>, EF<sub>7</sub> and EF<sub>8</sub>) as they are intended to capture the inefficiency of certain types of milling circuits. The assumption in this model is that the E<sub>ssbm</sub> is the "maximum efficiency" possible for a grinding circuit, to which we apply an empirical adjustment for SAG milling (''Cf'') which encapsulates any other EF factors in a form suitable for a SAG process.
   
 
=== Model defaults ===
 
=== Model defaults ===

Latest revision as of 15:05, 23 September 2024

Bond/Barratt Specific Energy Consumption SABC Model

This is a SAG or AG mill plus ball mill model that estimates the overall circuit specific energy consumption using the classical Bond work index equation for multi-stage crushing and single-stage ball milling (Essbm) including Rowland efficiency factors. The circuit Etotal is equal to the Essbm plus an inefficiency factor (CF) related to the difference in grinding efficiency of the two types of circuits. The SAG mill specific energy consumption (ESAG) is calculated using the 1979 Barratt equation and the ball mill specific energy consumption (Ebm)

This model includes a phantom cyclone effect in the equations due to the ball mill being calculated by difference and not being calculated by the normal Bond equation. The resultant operating work index of the ball mill will vary according to the ratio of the ball mill, rod mill and crushing work index values and is in the range of 80% of the measured ball mill work index value.

Testwork Required

Required parameters

  • F80, µm is the 80% passing size of the fresh feed to the circuit (expected to be a Bond-compatible size distribution).
  • P80, µm is the 80% passing size of the circuit product (expected to be a Bond-compatible size distribution).
  • Availability, expressed as a decimal (0.90 = 90% availability) is used to convert t/h to t/d.
  • Ball mill Operating strategy when circuit is SAG-limited (note, fixed speed can only be vary P80)

Optional parameters

  • Description and Comment are optional text fields
  • Maximum t/h limit is a t/h throughput limit above which the mills will turn down by, for example, grinding out.
  • T80 min and T80 min override the transfer size restrictions built into the model
  • Ball mill work index adjustment used to adjust WiBM for different P80 sizes.
    • exponent (a), the fitted Hukki exponent to the adjustment equation (enter as a positive number)
  • Essbm calibration factor (CF), model tuning factor for Etotal
  • Easag crushing circuit product, model tuning factor that dictates where the Barratt Easag equation switches from using the crushing Wi to the rod mill Wi
  • Easag RM circuit product, model tuning factor that dictates where the Barratt Easag equation switches from using the rod mill Wi to the ball mill Wi
  • Essbm crushing circuit product, model tuning factor that dictates where the Bond Essbm equation (the basis for the Barratt Etotal equation) switches from using the crushing Wi to the rod mill Wi

Formulae


E_{ssbm} = Wi_{C} \times \left ( \tfrac {10}{\sqrt{ Pc_{ssbm} }} - \tfrac {10}{\sqrt{ F_{80} }} \right ) + Wi_{RM} \times \left ( \tfrac {10}{\sqrt{ Prm }} - \tfrac {10}{\sqrt{ Pc_{ssbm} }} \right )\times EF_4^{RM} + Wi_{BM} \times \left ( \tfrac {10}{\sqrt{ P_{80} }} - \tfrac {10}{\sqrt{ Prm }} \right ) \times EF_4^{BM} \times EF_5 Rowland, 2006

where

  • Pcssbm is the calibration value entered for "Essbm crushing circuit product" (default: 9400)
  • Prm is the calibration value entered for "Easag RM circuit product" (default: 2100)

 E_{total}  = E_{ssbm} \times (CF)


\begin{align}
E_{SAG} = \Big[ & Wi_{C}   \times \left ( \tfrac {10}{\sqrt{ P_{C}  }} - \tfrac {10}{\sqrt{ F_{80}  }} \right ) +  Wi_{RM} \times \left ( \tfrac {10}{\sqrt{ P_{R}  }} - \tfrac {10}{\sqrt{ P_{C}  }} \right ) \times EF_4^{RM} \\
 & + Wi_{BM} \times \left ( \tfrac {10}{\sqrt{110} } - \tfrac {10}{\sqrt{ P_{R}  }} \right ) \times EF_4^{BM} \Big] \times 1.25 - Wi_{BM} \times \left ( \tfrac {10}{\sqrt{ {110} }} - \tfrac {10}{\sqrt{ T_{80} }} \right ) 
\end{align}
Barratt, 1989

where

  • PC is the calibration value entered for "Easag crushing circuit product" (default: 18850)

 E_{bm}= E_{total} - E_{sag}

As with most models, the power split between the SAG and ball mills is evaluated and the transfer size (T80) necessary to balance the power draw between the two stages is estimated. If that transfer size is within the specified limits, then the calculation proceeds to calculate the ESAG and Ebm using that transfer size. If the transfer size to balance the power is outside the specified limits, then the transfer size size is forced to the limit and the ESAG and Ebm values are then calculated. SAG-limited circuits (where the T80 hits the upper limit) will result in the ball mill over-grinding unless the ball mill speed is reduced to control the grind size (in the case of a variable speed mill). Ball-limited circuits (where the T80 hits the lower limit) will result in the charge level in the SAG mill dropping until its power draw matches the throughput limit dictated by the ball mill. It is important to realize that once the SAG mill charge level drops to be equal to the ball charge, the model will be operating "out of range" and the predictions will not be valid.

Circuits that operate at or near the transfer size limits should be less efficient than the default formulas will estimate. If you have a design that is consistently operating around the limits, then consider adding a contingency to inflate the specific energy consumption estimates to account for this poor grinding efficiency. Alternatively, if you are designing a new circuit then choose a different combination of SAG and ball mills that operates in a more reasonable transfer size range.

The optimum feed size for the rod mill oversize feed factor (EF4RM) in the Barratt ESAG equation is calculated using the greater of the sample's rod mill or crushing work index.

F_{O}= 16000 \times \sqrt{\frac{13}{Max(Wi_{RM}, Wi_{C}) \div 1.1023} }

The EF4RM is calculated greater of the crushing work index or the rod mill work index.

EF_{4}^{RM} = \frac{RR+ \bigl(\frac{Max(Wi_C, Wi_{RM})}{1.1023} -7 \bigr) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}

The ball mill oversize feed factor (EF4BM) is always calculated with the rod mill work index for the optimum feed size and the ball mill work index for EF4BM.

F_{O}= 4000 \times \sqrt{\frac{13}{Wi_{RM} \div 1.1023} }

EF_{4}^{BM} = \frac{RR+ \bigl( \frac{Wi_{BM}}{1.1023} -7 \bigr) \times \frac {F_{80}-F_{O}}{F_{O}}}{RR}

The model assumes that the calculated value of Essbm + the specified contingency is greater than the sum of Esag + Epeb; otherwise Ebm will be negative. The model does a check to see if Ebm calculated is negative (or close to zero) and, if so, will substitute the Raw Bond Model equation estimate for Ebm.

Recommended Usage

This model should be the lower end-member of the range of Bond/Barratt style models. It is useful when interpreting pilot plant results to determine where the pilot plant results fall relative to the spectrum of models. It can be used as a realistic model for design work in softer ore types like copper porphyries. The Optimized Bond/Barratt model is best used for low competency ores, such as those found in the Western Cordillera of North America or the Andes, where the impact crushing work index is similar to, or less than the rod mill and ball mill work index values. These ores tend to demonstrate a phantom cyclone effect.

Diagram of work index by size


Discussion, pebble crushing

The pebble crusher is assumed to be part of the SAG mill for the purposes of this calculation. Epeb is calculated as the power drawn by the pebble crusher divided by the circuit throughput. This Epeb is then deducted from Esag.

The total circuit specific energy consumption is generally higher when no pebble crusher is present in a circuit. The rule of thumb is apply the CF value of 1.15 for SAB circuits (5% higher overall specific energy consumption) Bennett, Dobby & Kosick, SAG 2001. Published data from Cadia confirms a 5% difference in throughput (and therefore, specific energy) when the pebble crusher is fully utilized.Hart, Valery, Clements, Reed, Song & Dunne, SAG 2001.

Discussion, Rowland EF factors

The model is only interested in the EF factors that are properties of the rock, namely the coarse feed EF4 factor and the fine product EF5 factor. The EF factors that correspond to machine characteristics are not included (EF1, EF2, EF3, EF6, EF7 and EF8) as they are intended to capture the inefficiency of certain types of milling circuits. The assumption in this model is that the Essbm is the "maximum efficiency" possible for a grinding circuit, to which we apply an empirical adjustment for SAG milling (Cf) which encapsulates any other EF factors in a form suitable for a SAG process.

Model defaults

  • Maximum T80, µm = 6000
  • Minimum T80, µm = 400
  • PC for Essbm= 9400
  • PC for Easag= 18850
  • PR= 2100
  • CF = 0.10