Zhou, H, Amodio, A & Boylan, N 2019, 'Informed mine closure by multi-dimensional modelling of tailings deposition and consolidation', in AB Fourie & M Tibbett (eds), Proceedings of the 13th International Conference on Mine Closure
, Australian Centre for Geomechanics, Perth, pp. 287-302, https://doi.org/10.36487/ACG_rep/1915_23_Zhou
This paper introduces an advanced geotechnical numerical modelling approach, which can be used to simulate the gradual deposition and large-strain consolidation of tailings in a multi-dimensional space. The findings from the modelling can be used not only to inform planning and design at mine closure but also help the management of tailings impoundments. For example, the results can be used to determine the settlement of tailings with time and thus inform backfilling planning; or inform on the tailings settlement and thereby assist with the design of an effective drainage network to divert surface water.
Following a state-of-practice review of tailings consolidation modelling, the recently developed multidimensional modelling technique (the Norwegian Geotechnical Institute [NGI] model) is described in the paper, with validations against available analytical solutions and comparisons with commonly used predictions presented by Townsend & McVay (1990). An example of the application of the NGI model is presented to demonstrate its capability and performance in modelling a full-scale scenario.
The NGI modelling approach is built on commercially available specialised geotechnical modelling software, FLAC (and FLAC3D), through its embedded programming language, FISH. The NGI model extends FLAC’s existing capability of large-strain consolidation calculation to simulate the gradual deposition process of the tailings. The deposition of the tailings slurry is divided into many discontinuous layers, and these layers are activated one after another from the bottom up. Activation of each new layer (on top of the existing tailings surface) is followed by a large-strain consolidation stage, with the consolidation time being determined as a function according to the volume of the layer and discharge history. Rock backfilling can be modelled in a similar fashion or can be customised.
A user-defined constitutive model (as part of the NGI model) has been developed to reproduce the key characteristics of the tailings during consolidation, including the variation of compressibility and permeability with reducing voids ratio. The consolidation of the tailings is modelled in a large-strain mode (i.e. the coordinates of the grid are updated frequently) in order to capture its effect on the consolidation behaviour and the deformation occurred prior to addition of a new layer.
The NGI model is also capable of performing complex three-dimensional problems accounting for varying consolidation boundary conditions, non-uniformity of the tailings material, and irregular pit geometries. As illustrated by the example application, this approach can be used to predict the development of tailing consolidation settlement with time, the amount of water expressed during consolidation, the capacity of the pit for tailings storage and the required amount of rock for backfilling.
Further development is ongoing in order to expand its modelling capability, such as prediction of increase in tailings strength with consolidation, modelling drying and consolidation of tailings, simulation of tailings dam construction process to improve prediction of tailings dam stability, and so forth.
Keywords: numerical modelling, tailings, large strain, consolidation
Barron, RA 1948, ‘Consolidation of fine-grained soils by drain wells’, Transactions of the ASCE, vol. 113, iss. 1, pp. 718–742.
Gibson, RE 1958, ‘The progress of consolidation in clay layer increasing in thickness with time’, Géotechnique, vol. 8, no. 4,
Gibson, RE, England, GL & Hussey, MJL 1967, ‘The theory of one-dimensional consolidation of saturated clays’, Géotechnique, vol. 17, no. 3, pp. 261–273,
Gjerapic, G, Johnson, J, Coffin, J & Znidarcic, D 2008, ‘Determination of tailings impoundment capacity via finite-strain consolidation models’, in AN Alshawabkeh, KR Reddy & MV Khire (eds), Proceedings of GeoCongress 2008, American Society Of Civil Engineers, Reston, pp. 798–805,
Hansbo, S 1981, ‘Consolidation of fine grained soils by prefabricated drains’, Proceedings of the 10th Conference on Soil Mechanics and Foundation Engineering, vol. 3, A.A. Balkema, Rotterdam, pp. 677–682.
Imai, G 1981, ‘Experimental studies on sedimentation mechanism and sediment formation of clay materials’, Soils and Foundations, vol. 21, no. 1, pp. 7–20,
Itasca 2016, FLAC 8.0 (Fast Lagrangian Analysis of Continua) User Manual, Itasca Consulting Group, Inc., Minneapolis.
Itasca 2019, FLAC3D 7.0 (Fast Lagrangian Analysis of Continua in 3 Dimensions) User Manual, Itasca Consulting Group, Inc., Minneapolis.
Nguyen, B-P, Yun, D-H & Kim, Y-T 2018, ‘An equivalent plane strain model of PVD-improved soft deposit’, Computers and Geotechnics, vol. 103, pp. 32–42,
Tito, LAA 2015, Numerical Evaluation of One-Dimensional Large-Strain Consolidation of Mine Tailings, master’s thesis, Colorado State University, Fort Collins.
Townsend, FC & McVay, MC 1990, ‘SOA: Large strain consolidation predictions’, Journal of Geotechnical Engineering, vol. 116, no. 2, pp. 222–243,
University of Florida 1987, Symposium on Consolidation and Disposal of Phosphatic and Other Waste Clays, Department of Civil Engineering, University of Florida, Gainesville.
Znidarcic, D 1999, ‘Predicting the behavior of disposed dredging soils’, in FBJ Barends, J Lindenberg, HJ Luger, L de Quelerij & A Verruijt (eds), Geotechnical Engineering for Transportation Infrastructure: Theory and Practice, Planning and Design, Construction and Maintenance, A.A. Balkema, Rotterdam, pp. 877–886.