Search Documents

By S. Pule, Z. Ndabezimbi, A. Hefer, E. Cupido

This article outlines various projects implemented at Rand Refinery to reduce effluent generation and freshwater consumption. The refinery generates about 2800 m³ of effluent annually and has initiated projects such as granulation water recycle and silver leach optimisation to address these issues. The granulation water recycle project has successfully reduced effluent volumes, while the silver leach optimisation has further decreased effluent by implementing stoichiometric addition of HCl and NaClO3. The vibratory shear enhanced processing (VSEP) trial showed promising results for gold recovery from wastewater, though further trials are recommended. The water savings of the effluent reduction projects are monitored by measuring the amount of effluent pumped to service yard effluent storage facilities with the use of flow meters connected on SCADA. Since the effluent reduction initiative of granulation effluent recycle in FY23, the effluent reduced to 2676 m³. In FY24, the stoichiometric addition of HCl and NaClO3 was implemented in the silver leach process, further reducing the effluent to 2083 m³.

Mar 5, 2025

By G Capes, L Morizzi, N Mitkas, E Jones, D Grant

This paper describes the mining setting for the sub-level open stoping (SLOS) mining method used at Olympic Dam (OD), and investigates the range of stope performance conditions encountered over 30 years of production, using a qualitative description of the factors influencing performance. Utilising a continuous improvement tool, the Plan-Do-Check-Act (PDCA) process, and a forward-looking rock mechanics program, the intent is to identify the key factors driving productivity and performance to create a stable and achievable long term mine plan.

Mar 15, 2021

By G M. Zhomin, A S. Arkhipin, A Kondratiev

Viscosity of multicomponent slags is one of the crucial physicochemical properties for many industrial processes and technologies. While the viscosity-temperature relationship of completely molten slags can well be described by an Arrhenian-type law, the compositional dependence of viscosity is still not fully understood, especially in compositional regions with the lack of experimental data (eg at the glass-former concentration of 90 molar per cent or higher). In the present work major compositional relationships of the activation energy of viscous flow have been revised for numerous binary and ternary slag systems on the basis of experimental data collected and stored in the viscosity database OxiVis. For example, it has been found that the activation energy relevant to the polymerisation effect in binary systems of the MO-GFO type, where MO is the modifier oxide (ie monovalent K2O, Na2O, Li2O and bivalent CaO, MgO, BaO, MnO, PbO etc) and GFO is the glass-forming oxide (SiO2 or B2O3) is not always proportional to the third power of XGFO (contrary to the Urbain formalism) and can vary significantly dependent on the slag system. Also, the activation energy of the charge compensation effect has been revised for ternary systems of the AO-MO-GFO and AO-MO type, where AO is the amphoteric oxide (Al2O3, Fe2O3 or other). A simple polynomial model has been proposed on the basis of the carried-out analysis for estimation of viscosities of multicomponent slag systems. The model has been optimised against the experimental viscosities in the selected binary, ternary, and higher-order oxide systems and has been shown reasonable agreement with the experimental data.

Aug 21, 2024

By R W. Shaw

The MTM Process targets treating lower grade and/or problem concentrates that are not attractive to smelters. The process has been developed at laboratory scale for copper and is based on a related development for zinc that has been extensively developed both in the laboratory and through pilot testing. The basis of the process is novel ammonia-ammonium chloride leach chemistry adapted to give high copper extraction whilst being selective over iron, arsenic and many other elements. The preferred flow sheet depends on location, especially infrastructure and utilities. These impact upon the leach and crystallisation stages more than copper reduction.

Aug 24, 2022

By A J. Morrison, K Demir, J Kohmuench, K Runge, C Evans

The Eriez HydroFloat® is a flotation device capable of recovering significantly coarser and lessliberated particles than is possible using a conventional mechanical flotation cell. Originally developed for industrial and bulk mineral beneficiation, the HydroFloat® is attracting increasing interest as a means of increasing the required grind size for metalliferous ore flotation. Doing so would reduce the costs of processing complex ores, either by decreasing the energy required to grind the flotation feed or by increasing the throughput achievable with a given comminution circuit. Processing coarser particles would also result in coarser waste streams, potentially improving the mechanical and hydraulic performance of the waste dams that they are sent to. Recent work at pilot – and industrial-scale has proven the potential of the HydroFloat® technology to achieve these gains in metalliferous applications. However, in order to realise the full potential of the HydroFloat® technology, it is necessary to understand not just its performance but also the drivers for that performance. Doing so will allow the design and operation of HydroFloat® devices to be optimised for different ores and applications. For this work, a HydroFloat® pilot rig was deployed to a coppergold plant in Australia and was used to investigate the effect of stream and operating conditions for treating the tailings stream of one of the plant’s rougher flotation cells. Under each set of operating conditions, samples were taken of the HydroFloat® feed, tailings and concentrate streams, and processed to yield the grade and recovery achieved by the device. At the same time, the bubble size and bed behaviour inside the device were measured. This paper presents preliminary relationships observed between the operating conditions and the performance of the HydroFloat® cell, and thus suggests a range of possible drivers for the measured coarse particle flotation performance.

Aug 24, 2022

By A Jyoti, J R. Black, R R. Haese, E O. Ansah

Copper (Cu) is a critical metal in the transition to a low-carbon power supply. Presently, low-grade copper sulfide ores account for almost 80 per cent of the global reserves. Large amounts of water are required to process these low-grade copper ores. The mining industry utilises lixiviants from recycled water from different sources to conserve water, whilst enhancing Cu recovery. Some utilise intermediate to pregnant leach solutions, sea water, and processed water from electrowinning circuits with highly variable concentrations of ions. Previous studies have highlighted the formation of precipitates that passivate the surfaces of minerals in these low-grade copper ores and/or block flow paths for lixiviant during heap leaching. In this study, we utilised concepts of phase replacement to devise an innovative means of controlling precipitation and increasing copper recovery using industrial brine as opposed to freshwater. We propose the formation of Al-rich phases (such as AlSO4 + or alunite) to replace iron sulfate and iron hydroxysulfates (such as jarosite) during the chalcopyrite dissolution-precipitation mechanism. The rationale is to have available Fe3+ throughout the chalcopyrite leaching process. We tested this method by carrying out various batch experiments using CaCl2, AlCl3, and NaCl industrial brine types compared to freshwater and explored the capacity to leach Cu from chalcopyrite. There is higher Cu recovery during oxidative leaching of chalcopyrite in industrial brine with 91 per cent recovery using AlCl3-rich brine, 84 per cent with NaCl-rich brine, and 56 per cent with CaCl2-rich brine as compared to 35 per cent recovery using freshwater. In the absence of Al3+, various forms of iron sulfates, and jarosite are formed, reducing the availability of Fe3+ for chalcopyrite dissolution. The introduction of Al3+ leads to the formation of AlSO4 + species, thus lowering sulfate activity and the tendency for K-jarosite formation.

Aug 24, 2022

By D Pienaar, B McFadzean, C T. O’Connor

Collector and frother molecules in sulfide mineral flotation are thought to have separate roles in the flotation process, where the collector selectively hydrophobizes the mineral surface and the frother adsorbs onto the bubble surface, inhibiting bubble coalescence. Molecular interactions have been observed between these molecules which allow them to co-adsorb at each of the respective interfaces, improving pulp phase recovery. This study investigated the molecular interactions between sodium diethyl dithiophosphate (SDEDTP) and hexanol, and the mechanism by which these interactions improve the recovery of galena. In this study each of the respective interfaces were investigated for the adsorption of SDEDTP. The adsorption of SDEDTP on galena was measured by UV-Vis spectrophotometry and isothermal titration calorimetry. At the air-water interface foam stability and pulp bubble coalescence was used to quantitatively discern if there are interactions between SDEDTP and hexanol. Microflotation and induction time experiments were performed to determine the effect of SDEDTP, hexanol and their mixtures on mineral recovery and bubble-particle attachment. The results of all these tests were collated to determine if SDEDTP-hexanol mixtures synergistically improve sulfide mineral recovery. It was found that SDEDTP reduced bubble coalescence in the pulp demonstrating activity at the airwater interface. Together SDEDTP and hexanol synergistically enhanced foam stability. This indicated that there are attractive molecular interactions between these reagents. Even though DTP did not adsorb onto galena, DTP-frother mixtures synergistically improved the recovery and attachment probability of this mineral. Adsorption studies with aeration showed that DTP in the presence of a frother transfers from the air-water interface to the galena mineral-water interface during bubble-particle collision. It was proposed that the transfer of hexanol through the thin film enhanced the kinetics of film thinning between the bubble and the particle by destabilising the thin film thus improving pulp phase mineral recovery.

Aug 24, 2022

By S Natsui, Y Goto, S Kawanishi, J Takahashi, H Nogami

The demand for copper increases yearly, and further improvements in copper smelting efficiency using flash furnaces are required. This study was conducted to understand and promote the interfacial reaction associated with the collision of silica flux and iron oxide particles, ie the formation reaction of fayalite slag, in the reaction shaft of the flash furnace. High-temperature microscopy was used to observe the reaction behaviour of a 2 mm granular SiO2 reagent in contact with micrometreorder- sized FeO reagents in situ. When the amount of FeO was lower than the stoichiometric amount with SiO2, the molten slag formed at the interface between FeO and SiO2 diffused into the SiO2 particle before dissolving the remaining FeO particles. It did not contribute to the subsequent reaction, and despite having the lowest amount of FeO among the samples in this study, many unreacted FeO particles remained on the surface of the samples after testing. On the other hand, under conditions with high FeO content and close to the eutectic composition of the FeO-SiO2 system, a liquid phase formed from low temperatures, gradually increasing in volume and eventually completely dissolving the 2 mm SiO2 grain. These results indicate that the early-stage slag is essential in increasing the interfacial area between FeO and SiO2 and facilitating the supply of FeO to SiO2. In other words, the FeO/SiO2 determines the amount of liquid phase and significantly influences the subsequent interfacial reaction rate between FeO and SiO2. Therefore, the local FeO/SiO2 of the particles falling down the reaction shaft should also be optimised to ensure early slag formation. Meanwhile, the time required for 2 mm SiO2 to complete the reaction for the temperature history employed in this study was significantly longer than the particle residence time in the reaction shaft. Further detailed investigation of the effects of silica flux particle size, composition, and reaction temperature on the reaction time will help find the conditions under which the reaction of silica flux can be completed in the reaction shaft. These kinetic approaches could contribute to future innovations in the copper smelting process using flash furnaces.

Jun 19, 2024

By Q Qu

Reliable prediction of longwall gas emission sources and rate plays a critical role in mine ventilation and methane drainage capacity planning. Over the years, Australian longwall mines have been using the historic Winter method and the Flugge method in modelling specific gas emissions (SGE). The predicted results often experience high discrepancies from the operational reality. To improve gas emission prediction, post-mining gas contents were measured in several mines in Qld and NSW by safely drilling core holes into the longwall goaf. The results revealed many new insights into gas emission characteristics from different coal seams and the drawbacks of the European methods in their use in high-production gassy Australian longwall mines. Based on measurement data, a new, time-dependent gas emission prediction model has been developed. Considering several factors such as caving, desorption and production rate, the new model can predict more accurate and meaningful gas emissions to meet today’s needs in mine gas emission management. A case study is provided to demonstrate the application of the new SGE model.

Oct 12, 2022

By R G. Reddy

The complexity in treating environmentally harmful impurities such as As, Sb, Bi from the base metal matte or metal in the smelting stage are responsible for the high cost of the refining process. The impurity capacities (such as arsenic, antimony and bismuth) of slags were calculated a priori using Reddy-Blander (RB) model. The capacity predictions were for a wide range of matte and slag compositions in copper smelting conditions. The calculated impurities capacities and impurity distribution ratios results are in good agreement with the available experimental and industrial slags data. The a priori knowledge of impurities is useful for reduction of energy consumption and enhanced environmental control in the current and future non-ferrous metal processes.

Aug 21, 2024

By Q Shu, T Deng, B Yan, B Song, D Wang

In view of the lack of phase equilibria data of vanadium-bearing oxide system and the desirability of them for design or optimisation of vanadium extraction process, the phase equilibria relationships of V2O3-FeO-SiO2-CaO(15 mass per cent) system were determined experimentally and then an isothermal section phase diagram of this system was constructed. Totally, 81 samples of different composition were equilibrated at 1623 K for 48 hours in iron crucible under Ar atmosphere. The equilibrated samples were quenched quickly in cold water, and then the phase constituents of all the samples were identified by X-ray diffraction (XRD) combined with electron probe microanalysis (EPMA) measurements.

Aug 21, 2024

By K L. Avarmaa, M A. Rhamdhani, A D. P Putera, H T. B M Petrus, G A. Brooks

Multistage slag treatment was carried out to investigate boron removal from metallurgical grade silicon. Boron is known to be one of the most difficult elements to remove from silicon. Depending on boron initial concentration, boron can be removed from silicon into certain level; most previous studies have attained final 50–100 ppmw (parts per million weight) boron in the silicon using slag treatment. The term ‘multistage’ refers to a process where the resulting silicon from the initial slag treatment is further refined multiple times using a new slag. In this study, an initial slag composition of 45CaO-45SiO2-10Al2O3 (wt per cent) was utilised. Using this technique, boron level in the silicon was expected to be gradually lowered because the new slag will always have more capacity to reach the equilibrium distribution coefficient (LB = (ppmwB)slag/(ppmwB)Si). The current study was carried out through combined thermodynamic calculations (using FactSage™ ver 8.2, by Thermfact Ltd, thermochemical package) and high-temperature experimentation. The thermodynamic calculations were carried out for reactions at 1500, 1550, and 1600°C for three different scenarios: 1. Single stage refining with varied slag to Si ratio. 2. Multistage refining with slag to Si ratio of 1:1. 3. Multistage refining with slag to Si ratio of 2:1. The multistage experimental work was carried out at 1550°C under inert gas atmosphere each time for 2 hrs, with slag to Si ratio of 2:1. The simulation results indicate that the boron distribution coefficient stay essentially constant with varying slag to Si ratio and temperature, around 1.8 to 2, and the experimental results measured LB in the same range, between 1.5 and 2.5. By harnessing this feature, the simulation showed that boron can be removed by approximately 78 per cent in each stage. The experimental result confirmed that the boron concentration in the silicon was reduced by a total of 97.8 per cent from its initial value using four stages of slag treatment, from initial 270 ppmw down to 6 ppmw.

Aug 21, 2024

By Satyajit A. Vaidya

This paper presents a case-study in the design and construction of deep foundation element supported mats for two residential towers on an18,000-square-foot site in Brooklyn, New York. Based on site sub-surface conditions, conventional mat foundations bearing on medium dense sand, having an allowable bearing capacity of 3 tsf, were originally anticipated for the 16 and 25-story towers. However, when the tower shear walls were relocated to the edges of the mat foundations due to architectural revisions, the significantly increased mat edge contact stresses applied to the soil subgrade resulted in calculated settlements of a magnitude that would adversely impact the proposed structures and the neighboring/bordering lot-line building and rail transit subway structures. To address this issue, deep foundation element supported mats were designed for the towers. The design was performed as an iterative process involving structural SAFE and geotechnical Plaxis finite element computer analysis to develop a coordinated geo-structural model of the foundation system. Initial presumptive design values of 45 pci and 445 kips/inch were used for the modulus of subgrade reaction and minimum deep foundation element stiffness, respectively. The designed foundation system consisted of 40 and 60-inch-thick mat foundations for the 16 and 25 story towers, respectively. Initially, the deep foundation elements were designed as 50-foot-long, 14-inch-diameter elements installed to sustain a maximum individual compressive load of 250 tons. The final deep foundation element design was subsequently refined to consist of a 48-foot-long, drilled and grouted IBO Titan 103/51 element designed to sustain a maximum load of 167 tons. Subsequent field load tests on the deep foundation elements confirmed element stiffness ranging from 1,040 to 1,530 kips/inch, with measured mobilized side shear of about 29 psi. The mat foundation reinforcement was adjusted for the pile stiffness values determined in the field load tests, and the deep foundation supported mats were constructed. The two towers were subsequently constructed, and have performed successfully.

Jan 1, 2010

By Lorie M. Dilley

Alaska Village Electric Cooperative (AVEC), an electric utility cooperative providing power to 51 Alaska villages, has made a concerted effort to install wind turbines in the villages they serve in an effort to mitigate the growing costs of electrical power generation. The subsurface soils in several of these villages consist of sands and silts with ground temperatures of between 30°F and 32°F. The presence of warm, degrading permafrost in three of the villages has required innovative foundation designs, including the use of helical piers and piles to support steel and concrete wind turbine foundations. The design of the deep foundations had to account for the effects of degrading permafrost likely resulting from observed climate change patterns. Uplift forces imposed on the foundations were approximately 100 kips per pile. Uplift pile load testing was conducted on the installed piers and piles in order to determine the actual working load that each pier or pile could sustain. In addition, pile driving analyzer testing was conducted on several piles installed in the permafrost and compared to pile load tests. Results indicated that the installed piers and piles achieved the necessary capacity with sufficient factors of safety. Thermosyphons were installed near the foundations in order to limit the effects of climate change on the permafrost. The two to four wind turbines in each of the three villages are currently operational. Refinement of the design based on the pile testing from the first village to the third has decreased the cost of the tower foundations and improved our understanding of designing for warm, degrading permafrost.

Jan 1, 2010

By P. N. Anyalebechi

A comparative study of the mechanical properties of automotive steering knuckles produced by three different shape-casting related processes has been conducted. It involved the characterization of the tensile and fatigue properties of aluminum alloy A356 knuckles produced by the vacuum/pressure riserless casting (VRC/PRC), pressure counter pressure casting (PCPC), and the hybrid Cobapress (a casting-forging) processes, in accordance with the appropriate ASTM standards. Surprisingly, the knuckles produced by the hybrid Cobapress process exhibited the lowest strength. For example, the average 0.2% yield strength of the Cobapress knuckles was 31 MPa and 40 MPa less than that of the VRC/PRC and PCPC knuckles, respectively. However, between stresses of 150-200 MPa, the fatigue lives of the VRC/PRC, PCPC, and Cobapress knuckles were comparable. But above 200 MPa, the Cobapress knuckles exhibited the lowest fatigue life. The observed differences in the mechanical properties are attributed to the inherently different microstructures of the different knuckles.

Jan 1, 2009

By Siavash Mahvelati, Alireza Kordjazi, Joseph Thomas Coe

Deep foundations are often constructed as end bearing members with a significant amount of their capacity derived from bedrock. Sites with highly variable bedrock conditions therefore present a unique challenge when attempting to estimate the design elevation of deep foundations. Typical subsurface characterization techniques such as drilling, standard penetration testing (SPT), and cone-penetration testing (CPT) may fail to adequately establish top-of-rock elevation in geologic settings prone to karst and/or significant amounts of residuum, saprolite, and deformed bedrock. Seismic surface wave geophysical methods have been increasingly used to map bedrock topography in such settings. However, the current standard of practice for these methods introduces resolution limitations in spatially variable subsurface conditions due to averaging effects from the analytical framework used to process the acquired data. This study explores the use of a tomographic approach with full waveform inversion (FWI) to map top-of-rock elevations at a numerically simulated site representative of highly variable bedrock topography. Comparisons are made with seismic refraction (SR) and a dispersion-based multichannel analysis of surface waves (MASW) approach as more commonly used in geophysical subsurface investigations. The results highlight that FWI provides a superior estimate of top-of-rock topography than SR and MASW at a cost of a significant increase in the required computational resources and time.

Sep 8, 2021

By Sarah Ramp, Daniel V. Cacciola, Justin Labrozzi, Shafiq I. Siddiqui, Faisal Ahmed

"Construction of embankments up to 34 ft in height for twelve ramps and roadways needed for the New Jersey Turnpike Interchange 14A project required ground improvement. The subsurface conditions, consisting of up to 30-ft thick layer of highly-compressible peat and organic soils, would induce up to 8 ft of primary and 4 in. of secondary consolidation settlements to the embankments over a long period of time. Innovative ground improvement solutions were needed to meet unique challenges posed by the poor subsurface conditions, space and existing utility constraints, aggressive construction schedule, and massive size of the project. A variety of ground improvement techniques were utilized in six different zones of the embankments to limit settlements and ensure stability of the embankments. These techniques included preloading with surcharge and prefabricated vertical drains (PVD), load balancing using light weight fill, and concrete column rigid inclusions (RI). Additionally, augercast piles with load transfer platform were used under an existing bridge with low vertical clearance. An innovative approach consisting of preloading with surcharge and PVDs followed by installation of RI was needed in one of the six embankment zones with the worst, i.e., very soft, soil conditions. The design of the RI was performed by a specialty ground improvement contractor (SGIC) using Plaxis 2D and 3D softwares. The performance of the ground improvement systems was verified and monitored by load tests and instrumentations. This paper presents details of the ground improvement design, construction, monitoring results, and lessons learned, with a focus on the RI.INTRODUCTIONNew Jersey Turnpike Interchange 14A is located within the cities of Bayonne and Jersey City, in Hudson County, New Jersey. The New Jersey Turnpike Authority (NJTA) identified Interchange 14A as an area of mass congestion in need of redevelopment to improve the increasing demand of traffic and roadway safety. The New Jersey Turnpike Interchange 14A project included the construction of twelve (12) new roadways, in addition to the construction of four new bridges. The limits of construction are shown in Fig. 1.The fill height required to construct the roadway embankments varied from 2 to 34 ft. The proposed embankment side slopes varied from 6H: 1V to 2H: 1V. Additionally, mechanically-stabilized-earth (MSE) walls were utilized to retain the embankment fill along some segments of the roadways. As shown in Table 1, the proposed roadways were divided into six zones, based on location, subsurface condition, and embankment height."

Jan 1, 2017

By N. Duda, P. Śliwiński, P. Stefaniak, B. Jachnik

Ventilation and air-conditioning are some of the most energy-consuming processes in deep underground mining. Additionally, the energy consumption of these processes systematically increases along with the expansion of the mine. Therefore, proper energy efficiency assessment and consumption forecasts of these systems are particularly important in modern mines for ongoing tracking and planning of the likely costs. In the article, a short-term prediction model for surface air–conditioning stations is described. This model takes into account dependencies that were identified, based on the data from the supervisory control and data acquisition (SCADA) systems. These dependencies between the response value and the external conditions are then incorporated into the model. This way, for example, the ambient temperature or the total number of specific days of the week (with sometimes significantly lower or higher than the weekly average energy consumption levels) are included. Taking into account these dependencies, and e.g. forecasted average monthly temperatures allow us to include the seasonality of energy consumption in the presented model. In the article, a sample energy efficiency analysis of the A/C system based on an example of the KGHM Polish Copper Polkowice-Sieroszowice mine is performed. The results of the solution proposed in the article are then evaluated, and compared to the currently functioning, more basic forecasting methods. However, given the general nature of the dependencies used in the model, similar solutions could be used in other case studies.

Sep 1, 2021

By M. E. Rossi, C. Badenhorst

Conditional bias occurs when the expected value of the true grade conditional on the estimated grade is not equal to the estimated grade. It is a quantification of how the average of estimates differs from the true value. Conditional bias is present in all resource estimates due to the unavoidable properties of kriging, and indeed all linear estimation (regression) methods. In a mining resource estimation context, the degree of conditional bias is dependent on the spacing of the conditioning data; more data leads to the expectation of less conditional bias This paper discusses conditional bias from two perspectives: 1. 2. Estimates seeking to eliminate, as much as possible, conditional bias; and Estimates that accept a larger amount of conditional bias in resource models, in order to achieve a specific outcome. The problem of conditional bias in estimation can only be fully appreciated if the objective or purpose of a resource model is understood. There are instances where conditional bias should be avoided because it is not reasonable to make decisions using estimates that are known to be wrong in expected value; for example, when making a final decision in grade control with respect to the destination of the material. In most instances in the mining resource estimation context, however, the estimated values will not be “final estimates”, but “interim estimates”, meaning that in the future more data will be obtained, and the ”interim estimate” will then be updated. In this context, the most useful resource model is an interim model that is able to predict the final model estimates, despite not having all the data available that will be collected in the future. This is the objective, for example, of a recoverable “long-term” model developed prior to operating a mine, for example at the feasibility stage. It attempts to predict the ore that a processing plant will receive, well into the future. In such cases, it can be unequivocally demonstrated that conditional bias is not only unavoidable, it is a necessity.

Sep 1, 2021

By Z. Saw, C. Poveda, N. L. Mai

Optimising the throughput of a mining network comprising all equipment from shovels to semi-autogenous grinding (SAG) mills is a complicated task, yet the potential benefits can be significant. Key modelling challenges stem from the sheer number of possible routes which material can travel, different equipment operating under different constraints, and the rate of SAG mills being calculated through a nonlinear function of the blended materials’ rock attributes. Nevertheless, this is an ideal onapplication area of mathematical programming, a branch of operations research. This paper describes an application of mixed integer linear programming (MILP) to optimise the mining equipment network of one of the world’s largest copper mines.The model takes as input a year-long monthly mine schedule of shovels in different pits, pushbacks, sequences and faces which produce material to be moved through a network of crushers, conveyors, silos, distributors, stockpiles and SAG mills, ultimately outputting a daily operational plan for all equipment that maximises the throughput of the network. All operational constraints and business rules are considered, including three quality attributes of copper ore. Comparing model results with actual data shows an improvement of 4.47% in network throughput tonnage. Additionally, although the model is solved using a free, open-source solver, optimal solutions are still achieved within a practical timeframe.

Sep 1, 2021