Search Documents

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

By R. J. Durrheim, Y. Yabe, M. S. D. Manzi, DSeis team, N. Z. Nkosi, S. B. Mngadi, Jr. H. Ogasawara

On 5 August 2014, the Orkney M5.5 earthquake, the largest in South African gold mining districts, occurred with an unusual strike-slip mechanism at great depth of 4.78 km below the surface. In a rare case, the ML 5.5 earthquake occurred within the area covered by the legacy 2D and 3D reflection seismic data acquired in the 1990s and 2000s for gold exploration and mine development. In addition, the earthquake ML 5.5 rupture plane was recorded and accurately delineated by the underground in-mine seismic network near the source region. The integration of the legacy 2D reflection seismic data and mine seismicity data allowed us to identify a near-vertical structure, striking northnorthwest-south southeast (NNW-SSE). The ICDP-DSeis team drilled three holes (Hole A, Hole B, and Hole C), and two holes (Hole B and Hole C), intersected the upper edges of the ML 5.5 rupture. These holes recovered metasediments, metabasalts, intrusive rocks including dolerite sills and lamprophyre dykes adjacent to the fault zone, and the fault gouge. Late Prof. Tullis Onstott and geomicrobiologists installed a packer in Hole A and successfully recovered saline water and detected gas (~10 MPa). Slip weakening and rupture propagation are significantly influenced by the existence of fault gouges and the production of wear material between two sliding rock surfaces. Using the fault gouge material recovered from Hole C, we conducted friction experiments under dry and room temperature conditions at high slip velocity of ~100 mm/s and normal stress of 2 MPa. The resulting steady-state frictional strength was ~0.66 over a slip weakening distance of ~9.1 m. The steady-state frictional strength was high, which may be caused by large gouge thickness, and the rate of wear generation. As a result, it is proposed that the Moab Khotsong ML 5.5 seismogenic zone is complex and could be controlled by three main processes: a) the complex structural architecture of the seismogenic zone (e.g., intersection of fault, lamprophyre dykes, and dolerite sills); b) the mechanical process induced by tectonic and/or mining related stresses; and c) the mechanical and chemical processes caused by the water and rock interaction.

Jan 12, 2025

By I. James, R. J. Durrheim, M. S. D. Manzi, O. Valeshin, J. Lepine, L. V. Socco, C. Colombero, A. Malehmir, M. K. Rapetsoa, M. Sihoyiya

The FUTURE project aims to enhance seismic surveying in South African mining by integrating advanced acquisition technologies. In this paper, we detail surface and underground seismic surveys at the South Deep gold mine as part of the project. Surface seismic acquisition used co-located broadband micro-electromechanical sensors, one and three component recorders connected to geophones, distributed acoustic sensing using a straight fibre optic cable, along with a 6 ton broadband seismic vibrator operating with 4–180 Hz linear sweeps. In addition, underground surveys employed a combination of straight, and helically wound fibre optic cables, cabled-based and nodal sensors, along with a 500 kg drophammer mounted on a skid steer bobcat. The surface and underground seismic arrays were time synchronised using a GPS-time transmitter to allow for accurate seismic data recording in a GPS-denied environment. Results show the efficacy of these techniques in acquiring high-resolution seismic data, revealing detailed subsurface geological structures and mineralisation down to 3500 m below ground surface. Surface data revealed clear seismic reflections corresponding to the Black Reef and Ventersdorp Contact Reef. An underground P-wave velocity model exhibits low velocity anomalies that are associated with fracturing, while a P-wave reflection stack section shows significant reflections and a degree of faulting below the tunnel floor that was not observed on the surface reflection seismic data. We highlight the potential of these technologies to enhance safety, efficiency, and sustainability in mining. This study sets a benchmark for future reflection seismic surveys, paving the way for safer, more efficient, deep mineral exploration and extraction practices.

Jan 12, 2025

By R. J. Durrheim, B. Manzunzu, V. Midzi

Selecting ground motion prediction models (GMPM) is crucial in assessing seismic hazard levels in a specific region. This study investigates the influence of ground motion model (GMM) selection on seismic hazard assessment in Johannesburg, South Africa. Most seismic events affecting Johannesburg are related to mining activities and water ingress into the underground mine voids. The most seismically active regions in South Africa include areas around Klerksdorp, Welkom, and Carletonville. These regions experience a significant number of earthquakes, with mining activities being a major factor contributing to seismic events. Ten ground motion prediction equations were selected from regions of similar tectonics setting and were assessed to see if they correspond to the available local data. Six of these were also compared to the ground motions produced by the earthquake of 14 August 2014, which occurred in the Orkney area, one of the mining areas in South Africa. Sensitivity analyses were also carried out to assess the influence of ground motion prediction equations (GMPE) on the seismic hazard of Johannesburg. The results indicate that the choice of GMPM significantly affects the estimation of seismic hazard levels in Johannesburg. Different GMM can lead to varying ground shaking intensity predictions, impacting risk assessments and engineering design considerations. The findings highlight the importance of carefully selecting appropriate GMM tailored to Johannesburg’s specific geological and tectonic characteristics. Understanding how different models perform in this region is essential for accurate seismic risk mitigation strategies.

Jan 12, 2025

By S. Plaatjie, L. Linzer, M. S. D. Manzi, M. Sihoyiya

Fiber-Optic Sensing and UAV-Platform Techniques for Innovative Mineral Exploration (FUTURE) is a joint European Union and South African consortium project funded under ERA-NET Cofund ERA-MIN3. The project aims to develop cost-effective and novel seismic methods to image geological structures ahead of the mine face, ultimately leading to improved decision-making and safer mining practices. In this study we present results from numerical simulations conducted at South Deep Gold Mine for active in-mine seismic surveys. The study uses a finite differencing code known as WAVE3D, which models models the propagation of a seismic wave generated by stress pertubations that are applied to a defined grid. The numerical simulation model included two major ore bodies (Black Reef and Ventersdorp Contact Reef) and mine tunnels. The model was constrained using laboratory physical property measurements from underground borehole samples such as seismic velocities (e.g., compressional and shear waves) and bulk densities. Two case studies were created, one modelled the active shots and receivers on surface, and the other modelled the shots at depth (in a mine tunnel at ~4000 m depth) with receivers on surface. In both cases, the objective was to investigate the wavefield propagation through the rockmass between the mine tunnel and surface. The simulated wavefield produced by the source located along the mine tunnel was successfully recorded on surface. The surface-recorded wavefield exhibits clear one-way P-wave arrival times at approximately 450 ms to 500 ms, which can be used to calculate the surface-tunnel P-wave refraction tomogram. The results from the study were used to constrain the design, and for the acquisition and processing of the real surface and in-mine seismic surveys acquired under the FUTURE project.

Jan 12, 2025

By P. Lenegan, R. Masethe, R. I. L. Ferreira

Underground mining in high rock-stress environments can induce unanticipated dynamic rockmass deformation, posing a serious risk to the workforce. Data collected by seismic networks are imperative for the design-as-you-mine process. Often an inappropriate assessment of the seismic hazard is done when unnecessary data are included in analytical techniques, particularly short-term time frames. Depending on the methodology or application of data analysis and interpretation, some of the information may not be relevant and even detract from the objective. By using the spatial event clusters associated with active workplaces to define the shape and size of the polygons for data selection for subsequent analysis, a problematic subjective component can be eliminated. The mining of tabular reefs in South Africa, be they gold- or platinum group metals, is very extensive, and underground workplaces are often widespread. Production stoppages are very costly, especially for marginal mines, so there is just cause not to evacuate an entire operation’s workforce following the occurrence of a large mining-induced seismic event. A rationale based on the damaging peak particle velocities (PPV) of the past seismicity can exclude the unaffected underground personnel and allow a quick orderly withdrawal where required. Micro-deformation may precede rockbursting, which can occur on-shift close to the workface. A mine-wide array of geophones suitable for the accurate location and analysis of seismicity may not be able to detect these higher-frequency emissions. Inexpensive but fit-for-purpose accelerometers, providing real-time, continuous measurement of acoustic rockmass deformation near the workplace, integrated through underground and surface communications networks, are shown to record ‘rock-talk’, clearly distinguishable from the instrument background noise, which can be in the order of hours leading up to an imminent damaging seismic event.

Jan 12, 2025

By T. Matlou, S. Myendeki, V. Midzi, T. Zulu, H. Mabotja, T. Mulabisana, B. S. Zulu, L. Hlatshwayo, T. Makhateng

The Council for Geoscience (CGS) performed a macroseismic investigation following the occurrence of an earthquake in the Boksburg area on 11 June 2023. The event had a magnitude of ML 4.4 and a focal mechanism that showed the event as having strike-slip motion along a northwest-southeast oriented fault. The event resulted in aftershocks, some of which occurred away from the causative fault in the east in an area without a mapped fault. However, the cluster of past events in that area suggests that there is an active fault. Due to the size of the main earthquake, a macroseismic survey was conducted to determine its impact on the community in the region. Analysis of observed data obtained from questionnaires resulted in the creation of 153 intensity data points. Their distribution showed high intensity values close to the earthquake source area. However, the investigation highlighted a number of anomalies, such as the higher than expected intensity values observed far from the event epicentre. It could be that these anomalies are caused by site conditions, for example, amplification of ground motion due to thick soil layers.

Jan 12, 2025

By S. Durapraj, R. T. Masethe, U. Maqina, A. C. Adoko

This study investigates seismicity and rockburst control within the Merensky Reef and UG2 orebody of an intermediate-depth platinum mine located in South Africa's Bushveld Complex. Utilising a mine-wide seismic monitoring system, over 1,900 seismic events were recorded between April 2023 and May 2024, revealing complex seismic patterns driven by mining-induced stress changes. The largest event, with a magnitude of ML 2.0, underscores the seismic risks associated with mining at increasing depths. Key findings include the spatial and temporal variability of seismic events, with high seismic activity observed near pillars and geological structures. A significant portion of events correlated with pillar and foundation failures. Moment tensor analysis identified dominant failure modes, such as stress-driven fracture closures and dynamic pillar failures. To mitigate risks, the study emphasises advanced support system redesign and rockburst risk assessment using the rockburst damage potential (RDP) method. The research highlights the importance of monitoring, geotechnical compliance, and adaptive designs in enhancing safety and reducing seismic hazards in underground mining operations. The study applied the RDP framework, classifying damage across R1 to R5 scales, where most instances fell into low-risk categories (R1 and R3), indicating minor to moderate damage. However, data limitations and the tabular stoping method's complexity led to inconclusive results, highlighting the need for more robust data and customised risk assessment tools. Despite these challenges, the RDP approach provided valuable insights into rockburst mechanisms and the effectiveness of existing support systems.

Jan 12, 2025

By F. C. Harris, S. J. Gaab, J. Sattarvand, S. Duzgun

Underground mines are inherently dangerous due to potential hazards such as fires, explosions, and caving of rock. Although these events are rare, once they occur a speedy and safe egress of mineworkers is critical to survival. Static evacuation systems use exit signs or directional lines but can get obstructed or be hardly visible in smoky or dusty conditions. This paper presents the investigation of the effectiveness of a smart evacuation technology in underground mining, where a real-time guidance system leads individuals on an optimised path to the egress while avoiding danger zones, compared to conventional practices. Participants simulated an evacuation from a virtual reality mine, using the conventional system and the smart evacuation system with real-time guidance. The highest reduction in evacuation time using the smart method amounted to nearly 40% on average. Moreover, 83% of participants preferred the smart guidance system, and agreed that it could enhance safety.

Sep 1, 2021

By E. Moosavi, K. Tolouei

In mine production scheduling problems are well-known, where most are categorised into deterministic and non-deterministic NP-hard problems because of their complexity. Many researchers have attempted to solve the problem by applying various optimisation techniques. While using traditional methods they observed extreme difficulty in solving highly complex problems. In the 1990s, many researchers addressed long-term production scheduling (LTPS) by using simulated annealing, the genetic algorithm, the tabu search algorithm, and the ant colony algorithm. These algorithms are known as the meta-heuristic algorithms, and have proved to be the most efficient algorithms for solving LTPS thus far. The firefly algorithm is one of the best methods for visualising problems related to LTPS, and provides the best possible optimisation. It is one of the simplest methods, and easy to apply to any NP-hard problem. The experimental, and comparison with other recently published algorithms shows that the proposed algorithm is feasible, and is an effective approach for long-term production scheduling problems. The results indicated that a better solution is obtained by the proposed method, compared to other methods in terms of cumulative net present value, and average ore grade. Additionally, the CPU time via the proposed model is roughly 2.52% higher than the other method.

Sep 1, 2021

By F. M. Ramadhan, M. R. Pratama

Monitoring open-pit mine road conditions has become very critical nowadays. There are several benefits to in conducting research on this field such as opportunity to improve truck’s speed which directly influences truck’s productivity. This research aims to provide a semi real-time, and sensors-based automatic road damage detector that predicts the road quality through road damage exact location tracing, road damages auto-classification, and real-time instruction in the form of work order format to mining road maintenance, and repair authorities through mobile apps and maintenance unit’s in-cabin display. It will also provide the maintenance/repair activity performance evaluation via integrated dashboards. The artificial neural network algorithm is applied during the process of data training regarding classifying road damages and severities which resulted in a valid model. The road damage classification model has a precision level of 80.5% in predicting whether the road is damaged or not. Using this approach, we are expecting to visualise a road quality map of all active roads on mine site.

Sep 1, 2021

By R. Minnitt, C. Morgan, D. Fourie

The problems associated with outliers and their adverse effects on resource estimation in deep level South African gold mines are presented, and solutions suggested in available literature are reviewed. Traditional Witwatersrand type outlier capping uses the 99th percentile grade as the upper limit, a somewhat blunt approach. The occurrence of high and low grade areas in mining operations and the estimation of resources using geostatistical techniques require that the assumption of stationarity be upheld. Honouring the stationarity of data requires that mining areas be subdivided into appropriately sized and spaced domains. In order to overcome high-grade conditionally biased block estimates, the overly strong influence of extreme grades is mitigated by using naturally occurring grade sub-domains, by identifying outliers within the grade sub-domains, allocating relevant capping limits, and weighting them in a post-estimation process to. Rather than applying a capping value across a domain it applies a maximum cap based on the extreme value and its range of influence; the method also ensures that relevant capping occurs even in lower grade sub-domains. The proposed method for capping outliers recognises that metal contained in high grade chip samples should not be over-capped, and based on the grade of surrounding samples. The post-processing method of capping outliers does not interfere with the kriging equations. The method is simple, elegant, and can be easily applied. An advantage is that the output from the kriging process is not overwritten, allowing for comparison between original kriging results and post-processing results.

Sep 1, 2021

By A. Melkumyan, K. Silversides

Chemical assays collected from blast holes have the potential to greatly improve the local accuracy of grade models. Assays must be assigned to the correct domain or they will introduce errors. Boundary surfaces created from the widely spaced exploration holes (~50m) do not provide sufficient resolution to correctly delineate domain boundaries between closely spaced (3-5m) blast holes. This paper demonstrates a probabilistic method for classifying production blast holes near domain boundaries. Gaussian processes models are applied to identify the assays and assay ratios most suited to distinguishing particular domains. These models are combined with the exploration based surfaces to classify the blast holes that can be confidently assigned to a particular domain. A reiterative method is used, where information from the blast holes closest to the exploration holes is used to label blast holes further away. This is tested on two domain boundaries in a banded iron formation-hosted iron ore mine.

Sep 1, 2021

By N. Madani, S. Abulkhair

In multivariate geostatistics, it is common to have different types of complexities between variables of interest. In this context, an inequality constraint is an example of complex bivariate relationships. Unfortunately, traditional co-kriging and co-simulation algorithms cannot reproduce this type of bivariate complexity, leading to the overestimation of disturbing elements. This paper proposes a new algorithm based on a hierarchical sequential Gaussian co-simulation framework, integrated with inverse transform sampling, to model inequality constraints between variables. First, the proposed methodology's validity was evaluated by applying it to a real case study from an iron deposit, with an inequality constraint between iron and aluminum oxide. Then the simulated results were compared with a conventional hierarchical co-simulation algorithm to investigate the effect of inverse transform sampling on the quality of the co-simulation. The results showed that the proposed algorithm can reproduce an inequality constraint between variables.

Sep 1, 2021

By S. Andika, O. A. Putri

Maintaining fleet productivity, which consists of balancing loader and truck performance, has always been a challenging task. An automatic dispatching system (ADS) is implemented to help monitoring fleet condition and assign trucks to a designated loader. Integrated data from real time Global Positioning System, Vehicle Health Monitoring System, and other production related data are acquired to support the ADS. However, the abnormal truck cycle has been the main culprit which resulted in unnecessary loader hanging time and trucks queuing time. Hanging and queuing time accounted for 2.34 and 2.07 minutes respectively per loading cycle or equal to 2,861 loose cubic metres per fleet per day production losses. Using integrated dispatching system data, the authors conducted research to study a truck’s hauling cycle that affects fleet productivity. The data is modelled to obtain the arrival rate pattern in a front loading area by adopting the queue theory algorithm. This paper aims to understand the potential pattern of a truck’s cycle abnormality and to resolve its occurrence.

Sep 1, 2021

By M. Stachowiak, M. Nowicki, K. Walas, A. Skoczylas, P. Stefaniak

In this paper the authors present a walking robot dedicated to the inspection of a belt conveyor’s transportation system in an underground mine. The developed robot can independently diagnose all of the key rotating elements of the belt conveyor (pulleys, bearings, drums, gears, engines) based on vibration and acoustic signals as well as typical RGB (red, green, blue) and thermal images. The article presents key diagnostic procedures applied to real industrial data obtained during robot testing.The work was developed as part of the ‘THING–subTerranean Haptic INvestiGator’ project. The authors presented the main assumptions of the research methodology developed for the needs of technical diagnostics of conveyor elements typical of inspection work performed by mining staff. Based on industrial data, the use of algorithms for elements such as rollers, drive stations, and the conveyor belt were shown. The main assumption of the research is to develop a robot adapted to the conditions of a deep mine, capable of carrying out an autonomous inspection of mining infrastructure. The conducted tests and the developed fault detection algorithms were aimed at preparing the robot for underground tests in the operating conditions of the mine.

Sep 1, 2021

By I. Kaabachi, F. Meunier, N. Morales

In open pit mining, ore and waste are hauled outside of the excavation area by trucks that move along a system of ramps. The design of these ramps must be carefully addressed because they have to comply with geomechanical constraints to ensure the stability of the pit walls, provide access to the different levels and phases, but also maximise the value of the extracted ore and reduce the extraction of waste. Because an optimisation approach able to deal with the design of the full ramp system seems to be out of reach at the moment, this paper addresses the problem of designing one ramp providing access to a single pushback, which has been determined beforehand. The optimisation objective is therefore to find a ramp that gives access to one given potential volume, while minimising its impact on the economic value of the volume. Two versions of the problem are formalised. One version allows the extraction of ore and waste only if this extraction is necessary for the ramp system. This version somehow minimises the intrinsic value of the ramp system. In the second version, it is possible to extract more than is needed by the ramp system. For instance, this takes into account the possibility of attaining ore on the wall opposite to the ramp as soon as the ramp reaches its level. In the setting mentioned above, the contribution of this work is three-fold. First, the computational complexity of the two different versions of the problem is presented and discussed; in particular, they are shown to be NP-hard. Second, compact integer programs modelling these two versions are proposed. Third, an efficient GRASP algorithm able to cope with the large instances of the problem is proposed. To assess the proposed approaches, the models and algorithms are applied to real instances with an off-the-shelf solver, illustrating and discussing their potential and relevance for the two versions of the problem. They are also compared with other approaches from the literature in terms of value and resulting designs.

Sep 1, 2021

By M. M. Ali, C. Musingwini, R. A. Phillips

The long-term open-pit mine production scheduling problem (OPMPSP) is typically solved in two broad stages. The first stage entails optimisation of the ultimate pit limit, which delimits the economical final pit and maximises the undiscounted value of the deposit. The second stage entails optimisation of the long-term open-pit mine production scheduling, which maximises the discounted value of the deposit. The long-term OPMPSP is a type of an NP-hard problem, such as the travelling salesman problem (TSP), because it is intrinsically complex to solve since it incorporates both hard and soft constraints. Therefore, metaheuristic approaches have been used to find global optimum solutions to the long-term OPMPSP. This paper presents a metaheuristic approach utilising a dual interchange algorithm (DIA) to solve the long-term OPMPSP. The DIA is a parent algorithm which supervises two subservient algorithms, namely, simulated annealing (SA) and Tabu search (TS). The DIA selects the better solution generated in each run from the two algorithms, and simultaneously penalises the algorithm that produced an inferior solution, while rewarding the one that produced a better solution. The DIA was tested on the three-dimensional Arizona Copper Deposit (KD) dataset, which is publicly available for research purposes from the online MineLib repository. It was set to generate net present values (NPVs) in five-minute intervals and the solution improved until it plateaued at about 30 minutes of runtime. The NPVs compared favourably to those obtained on the same KD dataset but obtained over 24-hour runtimes, when solved using simulated annealing and Tabu search algorithms, and the NPV generated from the work of Muñoz (2012) using the modified TopoSort as found on the MineLib repository at the time of this study. Therefore, this paper demonstrates that DIAs can be developed and used to solve the three-dimensional long-term OPMPSP.

Sep 1, 2021

By J. Saunders, N. Pretorius, G. Krafft, T. P. Katakwa

Imagine being able to see all the future planned mining activities in real time. This is exactly what SIMS (Short Interval Management System) will enable you to do. SIMS is a real time short interval control system that continuously measures and sets optimised plans for mining production processes. The system utilises algorithm-based software that takes resources (people, equipment, tasks, and workplaces) and gives you the best achievable plan. The algorithm was developed using methods found in the field of metaheuristics. A metaheuristic is a strategy to designed to efficiently search for optimal or near optimal solutions to problems that are very complex computationally. (Glover, 1986) There are several algorithms that have been developed to solve optimisation problems which are difficult to solve using traditional techniques. This has given rise to a field called evolutionary computing. (Simon, 2013) Examples of these techniques are genetic algorithms, particle swarm optimisation, and ant colony optimisation. These methods are used to solve problems in a variety of fields. As an example, Kok and Lane (2012) used particle swarm optimisation for equipment allocation in a mechanised underground mining environment. The first SIMS deployment took place in a platinum underground tabular mechanised mine near Burgersfort (province?) from 2017-2019. This included the development from a paper-based system to a fully operational online system. SIMS then moved to a second site and is currently being tested in a second platinum underground tabular mechanised mine in South Africa. Full deployment of the system is currently being phased in. SIMS was built in response to issues that have been plaguing the mining industry for years. These issues are all aligned into one major problem statement: how can we optimise mining production processes without compromising safety? This product was developed by combining years of mining experience with software development. SIMS is based on a four-pillar concept. Safety, productivity, accountability, and business improvement. SIMS uses a set of pre-configured mining and production rules, and with the resources that are available, to determine valid short term production plans. A background algorithm then determines the optimum plan to be executed. The system continuously produces a set of job cards that are pushed to respective operators on their allocated tablets. The management portal component allows you to view and manage all current and planned work on any connected device available. This allows real time visibility of completed, current, and planned work to operators, service departments and management, i.e., seeing into the future.

Sep 1, 2021

By I. O. Flores, A. Anani, H. Li

Current technological advancement has increased the exploitation of deeper ore bodies using underground methods. In certain cases, the mineral deposit starts near the surface and continues deep within the earth's crust. Such deposits are best mined with both underground and surface methods. Two main techniques are used to determine the optimal transition depth from surface to underground - a disintegrated approach and an integrated approach. Both methods have limitations in their use to determine the optimal transition depth. The objective of this research is to investigate the conditions in which an integrated approach to the transition depth problem is superior to a disintegrated approach. Three binary integer linear programming models are developed for this purpose: a surface and underground model for the disintegrated instance and an integrated model. A case study is performed to validate the models, with a comparative sensitivity analysis to obtain operational insights.

Sep 1, 2021