Search Documents

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 K Li, C Zhang, i, J Oh, S Saydam, H L. Ram

Rock surface roughness plays an important role in excavation engineering projects and applications related to surface support structures, such as shotcrete and thin-sprayed liner (TSL). The morphology of the substrate affects the mechanical behaviour at the interface significantly. Therefore, it is necessary to quantify the roughness to investigate the relationship between roughness and mechanical properties. The joint roughness coefficient (JRC) proposed by Barton and Choubey (1977) is widely used in practice. It initially estimated the roughness by comparing it to 10 standard profiles visually. However, this method has lots of limitations, it is not only subjective but also a two-dimensional parameter with a small lab scale (Beer, 2002). The subsequent research on JRC, from visual comparison with Barton’s standard profiles to measuring the roughness with a profilometer (Weissbach, 1978; Alameda-Hernández et al, 2014). Then, 2D profiles or 3D surfaces by different technologies with higher resolution and accuracy, such as 3D laser or optical scanners (Fardin, 2004; Jiang, 2020, 2021). These 3D scanners are expensive and only suitable for laboratory scales. LiDAR scanners (Lato, 2009; Aubertin, 2022) have been applied in industrial size, but it is difficult to achieve a millimetre-level accuracy, which is mostly used to detect the thickness of shotcrete in practice. Thus, finding a cost-effective surface roughness measurement method with high accuracy is necessary. Photogrammetry is another method to digitalise the roughness of the surface. Similar to 3D scanners, photogrammetry also generates the point cloud of the objective surface and then reconstructs the 3D model of the objects. 3D reconstruction technology is a trendy technology in recent years. The structure from motion (SfM) algorithm (Westoby, 2012) to rebuild the 3D model with a flexible workflow, aligning photos, building dense point cloud, mesh and texture, which has been proved that can rebuild an accurate model in variable fields.

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 M Khanal, J Qin, B Shen, Y Duan, X Luo

Gaining confidence in interpreting rock and strata properties of the coalmine roof is necessary for design, maintenance, modelling and optimisation of roadways supports. One of the key causes of roadway instabilities experienced in underground coalmines is an extreme variable nature (thickness, competence, discontinuities etc) of the coalmine roof. Normally, the thickness and strength of the roof are characterised with the geotechnical and geological data gathered from exploration boreholes, which are drilled often at considerable distances apart. This limited data cannot capture localised geological variations in the coalmine roof. With this background, this paper explores feasibility of using measurement while drilling concept to perform geotechnical characterisation of the coalmine roof. The parameters, for example, penetration rate, rotary speed, torque, sound, vibration, obtained during the drilling process can provide valuable insights into geotechnical properties. The paper describes the progress of the experimental work conducted at the laboratory on synthesized concrete samples representing various rock strata using a field-scale drilling machine. The samples with and without interfaces have been tested to replicate various underground roof strata. The signals returned from the drill machine, as the drill bit passes through the layers, have been analysed and interpreted to characterise whether the drill response data can differentiate various strengths associated with the rock layers. The initial laboratory results are promising as they show distinct responses of the various sample blocks in a controlled environment. The drilling data can differentiate the synthesized rock samples and detect ‘signatures’ of the roof strength variability and delineation of lithological interfaces (transitions). The weight and torque on the bit vary and show dependency on the strength of the blocks.

Nov 29, 2022

By J H. Wood, I Gray

This paper examines the occurrence of outbursts, coal bursts and rock bursts in underground mining. It then examines them from the viewpoint of the failure mechanism and velocity of particles. The presence of gas is extremely important as it can both act as the fluid pressure component of effective stress leading to failure, and as a substantial energy source in projecting material during a burst. For gas to drive a dynamic event it must be present as free gas held in voids within the rock or coal mass. This can pre-exist the burst event or be released from desorbing gas during it. The latter is typically part of a multi-stage outburst event which can reach very large proportions. Rock and coal bursts are driven by strain energy release and can vary in form. These tend to be single events associated with failure and strain energy release. It is also possible to conceptually have hybrid events where the failure of the mass and the energy release are due to both stress and expanding gas.

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 M W. Chen, I Canbulat, S Saydam, C G. Zhang

The height of fracturing (HoF) measures the propagation distance of continuous fracturing from the mining horizon and is also termed the height of cracking (Hebblewhite, 2020; Khanal et al, 2019) or hydrologically the height of groundwater drainage (Tammetta, 2013). In longwall practice, HoF is an important predictive indicator for assessing the dewatering effect of panel extraction on the shallow aquifer and surface flows. The HoF determination is basically dependent on field measurement works, typically using surface extensometer, borehole televiewer imaging, piezometer monitoring, packer testing, borehole breakout logging, and isotopic tracer monitoring (Brown and Walsh, 2022; Corbett, 2022; Hebblewhite, 2020; Holla and Armstrong, 1986; Mills and O’Grady, 1998; Walsh et al, 2022a, 2022b). These methods or tools have been widely applied to observe rock fracturing and groundwater responses to longwalls at Australian coalmines, for example, Dendrobium Mine (Brown and Walsh, 2022; Walsh et al, 2022b), Appin Mine (Walsh et al, 2022a), and Springvale Mine (Corbett, 2022). In Chinese mines, borehole televiewer imaging and packer testing are more frequently used due to convenient instrumentation, quick testing, and not requiring long-term maintenance of the instrumentation boreholes. The HoF determination also involves further raw data interpretations, for which specialist software for data extraction and quantitative criteria for bounding the fracture transition and termination are necessary. By combining the measured fracture distribution and longwall parameters regarding panel geometry and geological settings, several HoF models have been developed for conceptualising the understanding of strata depressurisation and providing mine engineers fast calculation tools for HoF prediction. The Mackie model considered the height of fracturing as the upper bound of fractured zone and revealed that within the fractured zone there are highly connected fractures and accordingly porosity enhancement and groundwater depressurisation, where the fracture network can be less transmissive from the mining horizon up (Hebblewhite, 2020). The fractured zone shares similar fracturing regime but different heights in varying models: HoF is 30 to 58 times of the mining height (t) based on Australian, Soviet and UK experience (Kendorski, 1993, 2006), and 21 to 33 t according to the field measurements in NSW Central Coast, Australia (Galvin, 2016). Some Australian longwall mines adopt Ditton model (Ditton and Merrick, 2014) (Equations 1.1 and 1.2 in Table 1) sourcing from NSW coalfield experience for HoF prediction, while Liu model (Liu, 1981) (Equations 2.1 to 2.4 in Table 1) developed from 27 coalmines in China is frequently used by Chinese coalmines.

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 X Z. Sun, N Jiang, J H. Zhao, C X. Wang, L Y. Bai

In backfill grouting for underground coalmines, the filling materials strengthen the caving rock and support the overlying strata to reduce surface subsidence. The broken roof strata fail and collapse during mining operations without appropriate supporting measures being taken. It is difficult to perform continuous backfill mining on the working face of such roofs using the existing mining technology. Therefore, in order to solve the above problems, fly ash and mine water are considered as filling materials, and flow characteristics of fly-ash slurry are investigated through laboratory experiments and theoretical analyses. Laws governing the diffusion of fly-ash slurry in the void of caving rock masses and in the void between a caving rock mass and a basic roof are obtained and verified. Based on the results obtained from the above analyses and field conditions at the Zhaoguan coalmine, Shandong Province, China, a cave backfill grouting system of the hauling pipeline is developed and successfully tested at the 1703 working face. The results demonstrate that a filling rate of 43.46 per cent is achieved, and the surface subsidence coefficient of the grouting process is found to be 0.475. Compared to the total caving method, the proposed system is found to achieve a reduction rate of 40.63 per cent. This effectively helps in lowering the value of the surface subsidence coefficient. Fly ash and mine water, considered as primary materials in this study, also play a significant role in improving the air quality and water environment.

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 H A. Linn, K Winn, C W. Lee

The Q-slope method for rock slope engineering was developed to enable engineering geologists and rock engineers to rapidly assess the stability of excavated rock slopes in the field and make optimal adjustments to slope angles as rock mass conditions become apparent during construction of road cuts or benches. An empirical relationship between Q-slope and long-term stable slope angle is now supported through over 500 cases studies from Asia, Australia, the Americas and Europe. This paper presents about the stability analysis by first-ever applying Q-slope method on abandoned quarry wall in Singapore. The 30 m high rock slope which is composed of monzogranite rock was severely damaged by blasting operation. As the rock slope is intended for public view the objective was to retain the natural outlook and without any artificial support application thus manual scaling was carried out to clear the loosed rock blocks. In the end stable slope angle of clean rock face is determined once again by Q-slope method.

Nov 29, 2022