Design the Leach Experiment: Bottle Rolls, Central Composites, and Response Surfaces
Twelve bottle-roll points scattered across dose and pH will still leave an argument about “the best condition.” A designed experiment answers a different question: which factors…
Twelve bottle-roll points scattered across dose and pH will still leave an argument about “the best condition.” A designed experiment answers a different question: which factors actually move recovery, how large the effect is, and where the model stops being trustworthy. Central-composite designs provide a useful way to test those questions; this article uses a teaching case to make the method practical.
Illustrative design points
Each marker is an illustrative bottle-roll run from a central-composite teaching design — not a drawn curve.
YX500 reaches a higher illustrative peak recovery but needs kg/t-scale dosing; the teaching economic comparison still preferred cyanide under its stated prices.
Fix the basis before you fix the design
Hold constant and document:
ore split and preparation;
grind / particle size distribution;
pulp density and solution volume;
agitation and air;
leach time and sampling schedule;
assay method and lab QC.
Vary only the factors the design assigns. If pH drifts because lime addition is inconsistent, the design is already broken.
Two illustrative designs in one picture
Design
Factors
Coded region
Centre
Illustrative max recovery
Cyanide (NaCN)
A: 147.6–232.4 ppm · B: pH 9.90–10.60
13 runs, central composite
~190 ppm · pH 10.25
92.10% at 220 ppm / pH 10.5
YX500
A: 0.14–4.66 kg/t · B: pH 9.90–10.60
13 runs, central composite
~2.4 kg/t · pH 10.25
96.23% at 4.66 kg/t
Both fit first-order models of the form Recovery = b₀ + b₁·A + b₂·B. Interactive cursor:
Loading simulation…
Read the model, then the ANOVA
Illustrative first-order fits:
Quantity
NaCN model
YX500 model
Equation
R = 55.8326 + 0.0996·A + 1.3626·B
R = −72.6799 + 3.429·A + 14.6514·B
Model F
43.50
12.07
Model p
<0.0001
0.0022
Dose term
highly significant
significant
pH term
not significant (p ≈ 0.32)
significant (p ≈ 0.02)
Lack of fit p
0.5652 (not significant)
0.1201 (not significant)
R²
0.897
0.707
Adjusted R²
0.876
0.649
Mean response
88.73%
85.73%
Interpretation discipline:
Significant model, non-significant lack of fit → the chosen form fits the design region.
pH not significant for NaCN inside 9.9–10.6 does not mean pH never matters; it means pH did not move recovery inside that narrow window once dose was modelled.
Lower R² for YX500 means more residual structure — dose and pH alone leave more unexplained variance (or centre-point scatter is larger).
Particle size and cyanide together
An illustrative tailings study varies particle size (75–212 µm) and cyanide concentration (600–1200 ppm) with bottle rolls on thickener-underflow composites. The teaching method sketch uses ~500 g solids, ~1:1 L/S, quicklime to pH 11, a 24 h leach with hourly solution samples, residue dried at 105 °C, cone-and-quartering, and XRF for elements of interest.
Illustrative economic comparison (using the example assumptions): at maximum recovery, YX500 net ~4.31 B Tzs vs cyanide ~4.76 B Tzs; near 90% recovery the gap narrows (~4.61 vs ~4.66 B). The example basis includes ~47,000 t/month, head grade ~0.87 g/t, gold price ~4 M Tzs/oz, reagent prices ~8 M Tzs/t (cyanide) and ~4 M Tzs/t (YX500).
net value = recovered gold value − reagent cost − operating cost − waste treatment − incremental capital
A 96% recovery that requires kg/t-scale reagent and different waste handling can lose to a stable 90–92% cyanide circuit.
central-composite design size in both illustrative trials
R² 0.897
NaCN first-order fit (teaching case)
4.13 pp
illustrative max recovery gap: YX500 vs NaCN
Failure modes
Design mistake
What it looks like later
No centre replicates
Cannot separate noise from curvature
Factors coupled (pH fixed by reagent)
Confounded coefficients
Unequal ore splits
Entire surface is ore variability
Optimising on residual-free extrapolation
Pilot fails outside tested dose
Recovery-only decision gate
Economics and waste surprise the board
Check yourself
Check yourselfANOVA shows a significant model and non-significant lack of fit. What may you claim?
Within the coded factor ranges, the chosen response form adequately fits the observed points relative to pure error. You may not claim global optimality, plant-scale validation, or that other factors are irrelevant.
Check yourselfWhy replicate the centre five times instead of running five new extreme points?
Centre replicates estimate pure error and detect curvature/instability without confounding the main effects; extreme points stretch the design but do not by themselves improve the error estimate.
Check yourselfYX500 shows higher max recovery but worse R². Which number drives the pilot plan?
Both: the higher peak justifies interest; the weaker fit and larger residual scatter justify more runs, not an immediate plant change.
Bottle-roll and DOE work still involve cyanide and alternate lixiviants under laboratory controls. Follow the site chemical hygiene plan; this article does not replace method statements.