Search Documents

By J Maxlow, M J. Woods, M P. Sandy

The Newcrest Mining Limited Telfer operation comprises both open pit and underground operations producing gold and copper. The underground operation is made up of three production areas, the sublevel cave (SLC), M Reefs narrow vein open stoping and Western longhole open stoping areas. With stoping situated on the abutments of a mature and active sublevel cave, the interaction between mining areas is often complex. This is particularly evident in the M Reefs orebodies with high levels of extraction in all four reefs, and reliance on pillars to maintain local and global stability. A key concern was the potential for a cascading pillar failure to develop and rapidly progress through the M Reefs stoping area. The monitoring plan needed to capture any change in overall stability and identify if deformation was caused by pillar failure within the M Reefs and hanging wall closure, or shear on the reefs caused by SLC influence.

Nov 29, 2022

By T Warren, A Dempster, N Barnett, J Oh, S Saydam

The influence of the insulating properties of lunar regolith was shown in the sub-surface temperature tests performed by the Apollo missions, with no significant temperature fluctuations occurring below 80 cm from the lunar surface even if surface temperatures fluctuated >200 K between lunar days and nights (Heiken et al, 1991). Key to the thermal insulation of the regolith is its low density near the surface having the greatest insulative effect with the lunar regolith becomes denser with depth. As the regolith becomes denser, the thermal conductivity of the regolith increases (Hayne et al, 2017). Thermal modelling software, Oxford 3D Thermal Model, which mapped the lunar surface temperatures effect on the thermal influence at the lunar south pole sub-surface over a lunar year found areas with brief illumination with surface temperatures >120 K could have sub-surface temperatures <112 K that can sustain long-term stable water ice at depths of <30 cm (Paige et al, 1992; King et al, 2020). Further to this, the observation of water molecules or other hydroxyls by the NASA/DLR Stratospheric Observatory for Infrared Astronomy (SOFIA) in non-permanently shadowed regions (PSRs) on the Moon, indicated that water could be present beyond PSRs (Honniball et al, 2021). Using the Oxford 3D Thermal Model in selected regions in the lunar south pole, can provide insights into the location and depth of cold traps that could sustain water ice.

Nov 29, 2022

By P B. Hills

A simple algebraic method of aligning large scale underground excavations such as a crusher station or workshop or indeed any large cavern or mine layout to the optimal stress condition is outlined. The method uses the stress transformation law to determine how the component stresses vary as the stress tensor is rotated around its axis. The aim is to define the orientations at which the maximum and minimum horizontal stresses occur such that the excavation can be aligned with its long axis paralleling the trend of the maximum horizontal stress leading to the minimum potential for overstressing of the side walls. The physical environment in which the stress field is acting is unchanged. At the design location, the trend and plunge of the major principal stress is unchanged. However, the coordinate system by which the location of the stress field is described is arbitrary and the stress transformation law allows that coordinate system to be rotated to a more convenient orientation and for the component stresses to be calculated within the alterative coordinate system. The stress field is typically described in terms of normal and shear stresses aligned to a north–south, east–west, vertical, coordinate system. In this scenario, there are two orientations in the 360° horizontal plane defined by the coordinate system where the horizontal normal stresses are equal, and at 45° to those orientations, two orientations where they diverge to the maximum extent. Those orientations of maximum divergence correspond to the trend of the maximum and minimum horizontal stresses, usually designated σH and σh, and that is what the algebraic manipulation seeks to determine. This technique is a first pass assessment which relies on a knowledge of the full stress tensor. It is not a replacement for numerical modelling or detailed design, but a tool to guide planning.

Nov 29, 2022

By I Canbulat, A Mcquillan

This paper presents a parametric analysis showing the impact of virgin in situ stresses on a shallow (less than 100 to 150 m) high slope in a jointed sedimentary rock environment. Slope stability is assessed at varying ratios of slope orientation and horizontal stresses. The results of this study indicate stress conditions should be included in stability assessments of excavated slopes, particularly where structure is orientated perpendicular to the direction of horizontal stress. Some consideration of stress state is applied in deep open cut pits (eg those with planned excavations of 400–1200 m) (Myrvang et al, 1993; Stacey et al, 2003; de Bruyn et al, 2014; Kozyrev et al, 2015). Dodd and Anderson (1972), Kalkani and Piteau (1976), Lee (1978), Coulthard et al (1992), Noorani et al (2011) have studied the effect or horizontal stress and its impact on slope stability, showing that tensile stresses tend to develop at the crest of the slope, with larger tensile zones for higher horizonal stresses and steeper slope angles (Stacey, 1970, 1973; Sjoberg, 2013). However, the consideration of virgin stresses on shallow slope stability is less frequent (Lynch et al, 2005; Dight, 2006; Lucas, 2006; Sjoberg, 2013). This has led many geotechnical engineers to believe the impact of virgin in situ stresses on open cut slope stability is minimal and not worth time or consideration in modelling, or their effects on slope stability are still poorly understood (Stacey et al, 2003; Sjoberg, 2013).

Nov 29, 2022

By E J. Hancock, K H. B Chu

Failures within hard rock underground excavations have historically been defined by a depth of failure, estimated in the field or through individual measurements, with a simplified failure geometry; such as 1/3 span, cubic or paraboloid. This simplified shape, volume and mass of failure is typically used to estimate in situ demand as a simplified discrete point. However, the highly variable failure geometry indicates simplified methods to estimate the mass can be erroneous in comparison to high precision LiDAR scans. Advances in technology have enabled the accurate capture of failure geometry and allow a far more accurate analysis of failures to be completed. This case study investigates a number of real-life failures to accurately define the failure geometry, the distribution of the failure parameters and compares this to traditional methodologies. This paper also outlines a suggested means of standardising failure geometry parameters whilst accounting for their variability.

Nov 29, 2022

By O Kirtel, A A. Faizan

With the rapid increase and development of high-speed trains (HSTs), the problems such as train-induced ground vibration and its dynamic effects on nearby structures have become a major environmental concern in urban areas. Ground-born vibrations produced by high-speed trains pose a great challenge for engineers to build structures in such areas that is applicable for residents. In order to study the effect of train-induced ground vibration on nearby buildings, an advanced three-dimensional finite element model for train-induced ground vibrations has been developed using PLAXIS 3D Ultimate. The results of the in situ measurements which was previously carried out by the authors were used to validate the numerical model. The field test for ground vibration due to the passage of a high-speed train was carried out on the Istanbul-Ankara high-speed railway in Turkey with a train speed of 250 km/h. An in-depth analysis was achieved to investigate the effect of different soil properties on the train-induced ground vibrations by using the verified model. According to the result of the vibration analysis, the relative vibration response curves in all directions (x, y, z) for different soil conditions were obtained comparatively. The obtained results of the vibration analysis are discussed in order to figure out some useful conclusions.

Nov 29, 2022

By M Revell, B Sainsbury, D Sainsbury

This paper extends the existing knowledge associated with the design and placement of prefabricated pastefill barricades through a consideration of asymmetric drive geometry and conditions. The load capacity of prefabricated barricades and their failure modes are considered in relation to abutment geology and barricade geometry that includes curvature and span and drive geometry that includes abutment angles, haunch angles and misalignment. Design charts are provided for typical barricade sizes to provide general guidance for design capacities and optimising location placement based on asymmetric drive geometry. A recent barricade failure is considered in relation to the drive geometry and the potential impact it had on the installed capacity.

Nov 29, 2022

By G W. Swana

Elastoplastic parameter is a parameter that describe how much materials behave elastically up to certain stress states and plastically thereafter (Hudson and Harrison, 1997). In the continuum analysis, two main parameters of the elastoplastic parameter are modulus elasticity and Poisson’s ratio. In the slope, elastoplastic parameter holds the main role for stability related to deformation in ductile type of rock. There’s a limitation to determine elastic or plastic deformation at the field for ductile rock type. Field manifestation commonly showed that there was no accumulation zone during and after the sliding failure mechanism, in other hand, at the research area, manifestation of deformation at the toe of slope only indicated by bulging shape. This paper will focus at the sliding failure occurred at Coal mining Pit Y owned by Arutmin Indonesia, which according to failure investigation report similar with the condition described above. This paper also will be based on two hypotheses, first hypothesis is claystone rock mass has not been reached plastic deformation yet and the second hypothesis is there is a chance that elastoplastic parameter of the rock mass has been changed during the deformation/sliding failure process. Research method of this paper are back analysis using probabilistic analysis in Finite Element Method (FEM), which conducted before and after failure. Then, the result of the displacement or deformation of the analysis will be checked with the actual deformation shape. Meanwhile, the change of the elastoplastic parameter will be analysed and interpreted. This paper hopefully able to be developed to study about elastoplastic behaviour of rock mass in actual condition and the influence of another uncertainty factor which determine the slope stability.

Nov 29, 2022

By J J. Van Wijk, A R. Penney

Successful management of geotechnical risk at a mining operation is more than just having a good ground control management plan or using advanced numerical modelling and the latest monitoring technologies. The fundamental requirements for effective geotechnical risk management include identification of potential failure modes, identifying and implementing appropriate controls, and demonstrating the effectiveness of those controls. Effective geotechnical risk management requires an integrated approach by confirming that appropriate risk management systems have been identified and implemented (management plans, procedures etc), and that these systems are executed consistently in the field across all working areas. Development of a thorough integrated auditing process executed by an experienced team of technical subject matter experts with substantial operational knowledge, provides a more complete picture of catastrophic hazard management. Following consistent evaluation criteria allows benchmarking of systems and implementation compliance between the sites being audited, and the results can then be ranked against other operations in an industry database.

Nov 29, 2022

By A Taheri, G D. Nguyen, M Karakus, S Akdag

During underground mining operations, rock masses are highly subjected to dynamic disturbance caused by blasting, mechanical drilling and earthquakes resulting in strain burst, slabbing and spalling. The dynamic fracture is a significant manifestation of rock failure in deep underground engineering, and it is of great importance to assess the dynamic fracture behaviour of the rock mass under elevated temperature conditions. The correlation between temperature – and loading rate-dependent dynamic fracture mechanism of strain burst is yet to be established. A good understanding of the dynamic fracture behaviour of brittle rock subjected to high-temperature under dynamic loading is required for the prediction of the damage extent during strain burst, and proper design as well as control of the underground rock structures. In this study, dynamic mode I fracture toughness tests are carried out on cracked chevron notched semicircular (CCNSCB) granite specimens to reveal the dynamic fracture characteristics of the brittle rock. For this purpose, a Split Hopkinson Pressure Bar (SHPB) apparatus is utilised to conduct the dynamic tests with various striker velocities on thermally degraded granite. The dynamic mechanical behaviour of granite after high-temperature treatment is examined and discussed. The dynamic fracture toughness is determined and compared for the specimens exposed to different temperatures at various dynamic impacts. The fracturing processes are recorded by a high-speed camera (HSC), and the crack propagation speeds are estimated by HSC image analysis. In addition, the dynamic fracture process and the coupled effects of temperature and loading rate on the dynamic fracture modes are identified by HSC image analysis. Energy partition characteristics in dynamic fracture of granite are quantitatively investigated.

Nov 29, 2022

By T Wang, D Wang

Photogrammetry-based techniques supported by traditional mapping tools are more frequently used for open pit geotechnical mapping. It fundamentally involves the collection of geological, structural and rock mass data from pit batters when safe to do so. These are essential components for refining and implementing continuous improvements to geotechnical models. The geotechnical model is used to review bench performance, optimise geotechnical design parameters, interpretation and analysis of pit slope stability. This paper is a case study of implementing different pit mapping methods at Mount Carlton gold mine, North Queensland. It critically reviews the tools required, including the Brunton compass, mobile application Clino, geotechnical data mapping software CloudCompare (2015) and Sirovision software. The review compares the image precision, structural measurements and data analysis. The review also underlines mine site management philosophy for safe geotechnical data collection to identify in-pit geotechnical hazards and the mitigation of potential slope instability. Sirovision is a remote sensing system for site geology mapping and interpreting geotechnical characteristics of exposed mining faces. The system utilises off-the-shelf digital SLR cameras to capture stereo photographs of the rock face and uses the latest image processing technology to generate accurate 3D models. These 3D models are digitally analysed using computational geometry to extract unbiased and accurate geological and geotechnical data. CloudCompare processes high accuracy 3D dense cloud rock surfaces reconstructed with photogrammetry software to extract the geological planes using the FACET or Compass plugins. The major discontinuity sets can be observed and the orientation (dip/dip direction) obtained from the algorithm is accurate and reliable for geotechnical assessment, including kinematic and rock slope stability analysis.

Nov 29, 2022

By E Shilov, R Dixon, L Tovey

BHP Western Australian Iron Ore (WAIO) operations use a suite of slope monitoring systems which are deployed to operational areas based on the ground control risk. Slope Radars are one system used within WAIO, having recently deployed and integrated a network of 22 radar systems across the operations. The integration required the development of supporting network architecture, infrastructure, virtualised hosting, security protocols and a technology support model. This has enabled data transfer between the physical radar systems located in open pits, to on-site offices, a vendor supported remote control room in Indonesia and the WAIO Integrated Remote Operations Centre for data analysis and monitoring of system and displacement alarms 24/7. The typical technical skillset of a Geotechnical Engineer does not include the knowledge and capability to understand and develop the supporting IT infrastructure required for establishing an integrated radar monitoring network. This paper will provide a discussion of the technology considerations for integrating radar systems and highlight the importance of the developed partnerships between internal geotechnical and technology teams and the external vendor team to ensure a reliable, robust and scalable slope monitoring solution.

Nov 29, 2022

By I Canbulat, J Emery, L Yao and C Zhang

This paper documents a process of data analytics and machine learning applied to coalmine roof convergence data. It follows a publication documenting the formation of the database. For each set of analyses Python scripts were utilised to extract data from an SQL database, with subsequent analyses also executed in Python, utilising open-source tools such as scikit-learn. Preliminary analyses with classic machine learning algorithms included multiple logistic regression, K-nearest neighbours, decision tree, random forest, support vector machine, kernel support vector machine, artificial neural network (ANN), naïve Bayes and multiple linear regression. Results indicated that relationships were complex and non-linear, with the artificial neural network most suitable. An ANN model was then refined through a standard process of hyperparameter optimisation and data augmentation to arrive at a final model. Due to the black box nature of the ANN further insight into how the parameters interacted was sought and found with Shapley additive explanations (SHAP values). SHAP values utilise Game Theory, where each input parameter or feature is a ‘player’ while the data set is the ‘team’. The SHAP value is the impact of each player on the target value, essentially establishing the contribution of each input parameter. The results from all analyses are delivered in the paper, together with insights, applications, strengths, and limitations of the methodology.

Nov 29, 2022

By B Banerjee, N Kumar, S Parihar

Rampura Agucha Pb-Zn Mine (RAM) is located about 200 km to the south-west of Jaipur in the state of Rajasthan, India. Open pit mining operations commenced in 1988 and continued to 2018. RAM open pit footprint is the largest in India with length of 2000 m, width 1500 m and depth of 400 m. Development of an underground mine under the operating pit commenced in 2013, through trial mining longhole open stoping with cemented backfill was established as a preferred mining method. Nine years in operations underground mine has reached the production capacity of 4.5 MT, which represent more than 70 per cent of the ore volume mined by the pit. RAM is currently the world second largest underground Pb-Zn mine.

Nov 29, 2022

By Y Tsukurimichi, D Asahina, T Endo, T Takemura

Scientific knowledge of Mars’ surface morphology and internal geological structure has been accumulated through satellite imagery and numerical topographic data using remote sensing technology, as well as seismic surveys by the on-site lander. Recent studies show that Marsquakes could be caused by the release of strain energy generated by the contraction of Mars. In particular, Marsquakes occur around the graben where many faults are existed and formed regularly spaced patterns. This implies that Marsquakes occur when the stress field of the Martian crust is an extension condition, and normal faults appear on the surface. We focus on the linear relationship between the fault spacing and the thickness of the mechanical layer as an approach to studying the internal structure from satellite images. Mechanical layer is a layer that composes of the same mechanical properties and is mainly bounded by lithology (Narr, 1991; Gross, 1993). The linear relationship between the thickness of the mechanical layer and the joint spacing has been studied and its slope is called the Fracture-Spacing Index (FSI) (Narr and Suppe, 1991). Soliva, Benedicto and Maerten (2006) showed that the stress drop zone (stress shadow) due to the fault activities affects the fault spacing. Figure 1 shows the approach used to estimate the thickness of the mechanical layer from the fault spacing on Mars. This paper aims to better understand the internal mechanical structure of Mars estimated by the geometrical pattern of faults on the Martian surface. We report the initial studies of this work. The lab-scale physical experiments and numerical experiments were presented to demonstrate the relationship between the thickness of the mechanical layer and the geometric pattern of fractures

Nov 29, 2022

By A Mcquillan, N Bar

Geotechnical models developed during the open pit mine planning stage are generally three-dimensional (3D) in component to capture the spatial variation in geological, structural, hydrogeological and geomechanical conditions. Yet when geotechnical analysis is completed, often the 3D geological, hydrogeological and structural models are simplified to two-dimensions (2D). This paper demonstrates this simplification, through referenced case studies, can lead to the wrong failure mechanism being analysed and/or a conservative Factor of Safety (FOS), or Strength Reduction Factor (SRF), being calculated, leading to a false sense of stability. There exists a paradigm that the FOS calculated using 2D limit equilibrium (LE), or the Strength Reduction Factor (SRF) calculated using 2D finite element (FE), modelling methods is generally conservative, ie lower, than the true FOS. However, this assumption that 3D FOS, or 3D SRF, is always higher than 2D FOS, or 2D SRF, is not always correct (Chen and Chameau, 1983; Bromhead and Martin, 2004; Herza et al, 2017; Pyke, 2017; Dana et al, 2018).

Nov 29, 2022

By C Jacobsen, N Harries

Assessing and managing instability hazards is an essential activity when working with unstable natural slopes and excavated open pit mining slopes. Slope monitoring has become the standard technique for the management of geotechnical risks associated with mine slope instability hazards. In recent decades, a proliferation of new monitoring techniques such as terrestrial radar, satellite InSAR and LIDAR have become standard operational tools, complementing conventional survey-based monitoring solutions. A reliable monitoring system should identify and record incipient anomalous slope movement. Although the most obvious purpose of a monitoring system is safety related, slope deformation and performance monitoring also enhance the understanding of slope behaviour and assists in improving system design and implementation. For a monitoring system to be considered an effective early warning tool, it needs to operate in real-time and should be linked to the mine’s emergency communication systems. Terrestrial radar systems have been developed that can provide near real-time measurements, at a high level of precision and covering broad areas in all weather conditions. This has had a positive impact on the management of geotechnical risks in operational conditions. The technology adoption of RADAR in mine monitoring has significantly progressed in the last 20 years from technology enthusiasts, to visionaries, to pragmatists to conservative users. However, the increasing adoption of RADAR cannot reduce or remove other monitoring techniques that often have important benefits; even if these technologies cannot be used as a real-time operational tool. An increasing awareness on the importance of combining data from different monitoring techniques to develop a holistic picture of all possible geotechnical hazards is a current industry trend. This is proving to be an excellent foundation for managing hazards associated with slope instability in open pit mines, and is only expected to improve with time with the adoption of new technologies based on big data and Artificial Intelligence (AI).

Nov 29, 2022

By Z R. Guo, Y Qiu, X F. Wang

Mechanical damage leads to the increase of effective stress of sediments, and pore closure hinders the discharge of methane gas. Ignoring the influence of mechanical damage leads to an excessive prediction of hydrate gas production. Therefore, a fully coupled model considering mechanical damage-chemical decomposition-gas-water two-phase seepage-temperature was developed in this study, and the model was verified by Masuda experiment. Based on the fully coupled model, the temporal and spatial evolution of the whole process of hydrate decomposition under the influence of three initial absolute permeabilities, three pressure drops and three water saturations were considered. The results showed that: • The gas production rate reached the peak value in a short time after depressurisation, and the peak value increased with the increase of the absolute initial permeability. The larger the initial absolute permeability was, the faster the damage variable reached the peak value. The porosity increased gradually before the damage variable reached its peak, and the gas flowed more easily. • When the pressure drop was low, the hydrate remained stable in solid form. As the pressure drop increased, the hydrate gradually began to decompose. The greater the pressure drop was, the faster the initial decomposition rate of the hydrate was. So the damage variable increased rapidly and the porosity increased continuously, resulting in a large gas production rate at the initial stage of decomposition. • The larger the initial water saturation was, the smaller the effective stress wss, resulting in the slow growth of the damage variable and the small change of porosity.While the gas relative permeability was small, it resulted in a smaller gas production rate, and hydrate decomposition was slower. • The damage variable showed a process of rapid growth first, then slowly increased to the peak, and then gradually decreased to a stable state. When the damage was considered, the hydrate decomposition time to the steady state was longer, and the gas production was about 80.68 per cent of that without considering the damage.

Nov 29, 2022

By D E. Tennant

In the fields of geotechnical engineering and slope management, various software packages are used to assess the behaviour or propensity for a slope to fail. Traditional kinematic analysis calculates the probability of failure of design sectors with respect to dominant structural set orientations. A novel methodology is provided here to view potential hazardous slope behaviour in 3D, allowing engineers to make informed decisions regarding various pit slope geometries. Apparent dips are calculated for every point in space with respect to every slope orientation of the entire pit surface/design, specifically in reference to bench and inter-ramp planar failure mechanisms. Implementing a proactive approach to slope design will not only lead to early mitigation of risk but can also add direct value by enhancing or improving upon slope angles for earlier access to ore. From a safety perspective, identification of unfavourably oriented slope geometries can be dealt with before failure occurs. This methodology relies on a strong understanding of the existing rock mass anisotropy, particularly behind the face, thereby removing some of the uncertainty surrounding pit design as well as presenting engineers with an additional means to manage risk.

Nov 29, 2022

By E Ben-Dor, D Bakun-Mazor, Y Ben-Ari

Determining the mechanical properties of rocks is important for various civil engineering fields, including the mining of raw materials for aggregates used in the construction and paving industries. Traditionally, the mechanical properties of rocks are obtained through in situ and laboratory tests during geotechnical surveys. However, these time-consuming surveys involve many resources. In contrast, hyperspectral remote sensing methods make it possible to identify the mineralogical composition and crystallographic structure of the rock; properties that control the mechanical properties of the rocks. In this work, we characterise the mechanical properties of carbonate rocks by using a hyperspectral sensor in laboratory conditions. We collected about 150 cylindrical samples of carbonate rock, with a wide range of strength values from several rock outcrops in Israel. We used a point spectrometer in the range of 0.4–2.5 μm and a spectral image sensor in the range of 3.0–12.0 μm, scanning the samples to obtain the signature of their light reflections and spectral emissivity. We then measured the samples’ density, porosity, water absorption, and uniaxial compressive strength (UCS). We used sophisticated data mining to find statistical relationships between the hyperspectral signatures of the samples and their mechanical properties. We used this data to identify the most dominant wavelengths for predicting mechanical properties. We found that the density, porosity, and water absorption of carbonate rocks could be confidently predicted based on spectroscopy data, while the UCS of the rock could also be predicted, but less significantly. The results of the study pave the way for the development of measuring tools for the mechanical properties of rock, based on non-destructive tests of quarrying materials.

Nov 29, 2022