Analysis and prediction of structural instabilities in open pit mines are an important design and operational consideration for ensuring safety and productivity of the operation. Unstable wedges and blocks occurring a...Analysis and prediction of structural instabilities in open pit mines are an important design and operational consideration for ensuring safety and productivity of the operation. Unstable wedges and blocks occurring at the surface of the pit walls may be identified through three-dimensional(3D) image analysis combined with the discrete fracture network(DFN) approach. Kinematic analysis based on polyhedral modelling can be used for first pass analysis but cannot capture composite failure mechanisms involving both structurally controlled and rock mass progressive failures. A methodology is proposed in this paper to overcome such limitations by coupling DFN models with geomechanical simulations based on the discrete element method(DEM). Further, high resolution photogrammetric data are used to identify valid model scenarios. An identified wedge failure that occurred in an Australian coal mine is used to validate the methodology. In this particular case, the failure surface was induced as a result of the rock mass progressive failure that developed from the toe of the structure inside the intact rock matrix. Analysis has been undertaken to determine in what scenarios the measured and predicted failure surfaces can be used to calibrate strength parameters in the model.展开更多
Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual conne...Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.展开更多
Discontinuity waviness is one of the most important properties that influence shear strength of jointed rock masses,and it should be incorporated into numerical models for slope stability assessment.However,in most ex...Discontinuity waviness is one of the most important properties that influence shear strength of jointed rock masses,and it should be incorporated into numerical models for slope stability assessment.However,in most existing numerical modeling tools,discontinuities are often simplified into planar surfaces.Discrete fracture network modeling tools such as MoFrac allow the simulation of non-planar discontinuities which can be incorporated into lattice-spring-based geomechanical software such as Slope Model for slope stability assessment.In this study,the slope failure of the south wall at Cadia Hill open pit mine is simulated using the lattice-spring-based synthetic rock mass(LS-SRM)modeling approach.First,the slope model is calibrated using field displacement monitoring data,and then the influence of different discontinuity configurations on the stability of the slope is investigated.The modeling results show that the slope with non-planar discontinuities is comparatively more stable than the ones with planar discontinuities.In addition,the slope becomes increasingly unstable with the increases of discontinuity intensity and size.At greater pit depth with higher in situ stress,both the slope models with planar and non-planar discontinuities experience localized failures due to very high stress concentrations,and the slope model with planar discontinuities is more deformable and less stable than that with non-planar discontinuities.展开更多
Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the con...Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.展开更多
Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlight...Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.展开更多
Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, an...Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.展开更多
The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or ...The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints）. A synthetic rock mass modelling（SRM） approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method（DEM）-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the preexisting joints is generated by employing discrete fracture network（DFN） modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation（RQD）, joint spacing, areal fracture intensity（P21）, and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index（GSI）. The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.展开更多
This study at the Esmeralda Mine,part of the El Teniente Division of CODELCO,investigates optimizing hydraulic fracturing(HF)holes’spatial distribution to improve rock material production in one of the world's la...This study at the Esmeralda Mine,part of the El Teniente Division of CODELCO,investigates optimizing hydraulic fracturing(HF)holes’spatial distribution to improve rock material production in one of the world's largest copper-molybdenum deposits.Utilizing diverse data sources,including borehole,oriented borehole,and photogrammetry data,along with hang-up frequency and hydrofracturing details,we applied discrete fracture network(DFN)modeling to analyze in-situ block size distribution and fragmentation.These results are based on 12,000 realizations of discrete fracture network(DFN)models using R-Dis-Frag computer pacakge at real cave volumes of 200 m200 m200 m,with varying parameters,which significantly enhances their reliability.The incorporation of DFN modeling and geostatistical simulation allows for capturing the interaction berween several spatial variables and explaining the variations observed in the production results at the draw points.Keyfindings of spatio-statistical analysis highlight the significance of volumetric fracture intensity(P32)and extraction column height in reducing hang-up events and enhancing fragmentation efficiency.The study integrates HF-induced and natural fracture intensities,revealing that higher P32 values and higher draw columns correlate with fewer hang-ups and better fragmentation.We recommend non-regular HF patterns for high P32 zones to improve operational efficiency.This research provides insights into optimizing mining operations,acknowledging the limitations of HF propagation efficacy and paving the way for further exploration into the interplay between hydraulic fracturing and natural discontinuities.展开更多
Cleats are the dominant micro-fracture network controlling the macro-mechanical behavior of coal.Improved understanding of the spatial characteristics of cleat networks is therefore important to the coal mining indust...Cleats are the dominant micro-fracture network controlling the macro-mechanical behavior of coal.Improved understanding of the spatial characteristics of cleat networks is therefore important to the coal mining industry.Discrete fracture networks(DFNs)are increasingly used in engineering analyses to spatially model fractures at various scales.The reliability of coal DFNs largely depends on the confidence in the input cleat statistics.Estimates of these parameters can be made from image-based three-dimensional(3D)characterization of coal cleats using X-ray micro-computed tomography(m CT).One key step in this process,after cleat extraction,is the separation of individual cleats,without which the cleats are a connected network and statistics for different cleat sets cannot be measured.In this paper,a feature extraction-based image processing method is introduced to identify and separate distinct cleat groups from 3D X-ray m CT images.Kernels(filters)representing explicit cleat features of coal are built and cleat separation is successfully achieved by convolutional operations on 3D coal images.The new method is applied to a coal specimen with 80 mm in diameter and 100 mm in length acquired from an Anglo American Steelmaking Coal mine in the Bowen Basin,Queensland,Australia.It is demonstrated that the new method produces reliable cleat separation capable of defining individual cleats and preserving 3D topology after separation.Bedding-parallel fractures are also identified and separated,which has his-torically been challenging to delineate and rarely reported.A variety of cleat/fracture statistics is measured which not only can quantitatively characterize the cleat/fracture system but also can be used for DFN modeling.Finally,variability and heterogeneity with respect to the core axis are investigated.Significant heterogeneity is observed and suggests that the representative elementary volume(REV)of the cleat groups for engineering purposes may be a complex problem requiring careful consideration.展开更多
基金supported by the IMSRN French Company through a CIFRE grant No. 2012/0710CSIRO Energy Flagship+1 种基金QCAT in AustraliaThe laboratory 3SR is part of the Lab Ex Tec 21 (Investissements d’Avenir e grant agreement No. ANR-11-LABX-0030)
文摘Analysis and prediction of structural instabilities in open pit mines are an important design and operational consideration for ensuring safety and productivity of the operation. Unstable wedges and blocks occurring at the surface of the pit walls may be identified through three-dimensional(3D) image analysis combined with the discrete fracture network(DFN) approach. Kinematic analysis based on polyhedral modelling can be used for first pass analysis but cannot capture composite failure mechanisms involving both structurally controlled and rock mass progressive failures. A methodology is proposed in this paper to overcome such limitations by coupling DFN models with geomechanical simulations based on the discrete element method(DEM). Further, high resolution photogrammetric data are used to identify valid model scenarios. An identified wedge failure that occurred in an Australian coal mine is used to validate the methodology. In this particular case, the failure surface was induced as a result of the rock mass progressive failure that developed from the toe of the structure inside the intact rock matrix. Analysis has been undertaken to determine in what scenarios the measured and predicted failure surfaces can be used to calibrate strength parameters in the model.
基金sponsored by the General Program of the National Natural Science Foundation of China(Grant Nos.52079129 and 52209148)the Hubei Provincial General Fund,China(Grant No.2023AFB567)。
文摘Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.
基金Ontario Trillium Scholarship for supporting the doctorate program at Laurentian UniversityFinancial supports from the Natural Sciences and Engineering Research Council of Canada(NSERC CRD 470490-14)of Canada+1 种基金Nuclear Waste Management Organization(NWMO)Rio Tinto。
文摘Discontinuity waviness is one of the most important properties that influence shear strength of jointed rock masses,and it should be incorporated into numerical models for slope stability assessment.However,in most existing numerical modeling tools,discontinuities are often simplified into planar surfaces.Discrete fracture network modeling tools such as MoFrac allow the simulation of non-planar discontinuities which can be incorporated into lattice-spring-based geomechanical software such as Slope Model for slope stability assessment.In this study,the slope failure of the south wall at Cadia Hill open pit mine is simulated using the lattice-spring-based synthetic rock mass(LS-SRM)modeling approach.First,the slope model is calibrated using field displacement monitoring data,and then the influence of different discontinuity configurations on the stability of the slope is investigated.The modeling results show that the slope with non-planar discontinuities is comparatively more stable than the ones with planar discontinuities.In addition,the slope becomes increasingly unstable with the increases of discontinuity intensity and size.At greater pit depth with higher in situ stress,both the slope models with planar and non-planar discontinuities experience localized failures due to very high stress concentrations,and the slope model with planar discontinuities is more deformable and less stable than that with non-planar discontinuities.
基金co-supported by the National Basic Research Program of China(Grant No.2011CB201103)the National Science and Technology Major Project(GrantNo.2011ZX05004003)
文摘Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.
文摘Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.
基金the financial support from the National Institute for Occupational Safety and Health(NIOSH)(200-2014-59613)for conducting this research
文摘Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.
基金the Nuclear Waste Management Organization (NWMO) of Canadathe National Science and Engineering Research Council (NSERC)+1 种基金the Canadian Ministry of National Defence (DND)the RMC Green Team for funding this research
文摘The main objective of this paper is to examine the influence of the applied confining stress on the rock mass modulus of moderately jointed rocks（well interlocked undisturbed rock mass with blocks formed by three or less intersecting joints）. A synthetic rock mass modelling（SRM） approach is employed to determine the mechanical properties of the rock mass. In this approach, the intact body of rock is represented by the discrete element method（DEM）-Voronoi grains with the ability of simulating the initiation and propagation of microcracks within the intact part of the model. The geometry of the preexisting joints is generated by employing discrete fracture network（DFN） modelling based on field joint data collected from the Brockville Tunnel using LiDAR scanning. The geometrical characteristics of the simulated joints at a representative sample size are first validated against the field data, and then used to measure the rock quality designation（RQD）, joint spacing, areal fracture intensity（P21）, and block volumes. These geometrical quantities are used to quantitatively determine a representative range of the geological strength index（GSI）. The results show that estimating the GSI using the RQD tends to make a closer estimate of the degree of blockiness that leads to GSI values corresponding to those obtained from direct visual observations of the rock mass conditions in the field. The use of joint spacing and block volume in order to quantify the GSI value range for the studied rock mass suggests a lower range compared to that evaluated in situ. Based on numerical modelling results and laboratory data of rock testing reported in the literature, a semi-empirical equation is proposed that relates the rock mass modulus to confinement as a function of the areal fracture intensity and joint stiffness.
基金the funding of the Agencia Nacional de Investigación y 761 Desarrollo (ANID), through grant project of Fondecyt Iniciacion No. 11221093Basal Grants Center for Modeling ACE210010 and FB210005
文摘This study at the Esmeralda Mine,part of the El Teniente Division of CODELCO,investigates optimizing hydraulic fracturing(HF)holes’spatial distribution to improve rock material production in one of the world's largest copper-molybdenum deposits.Utilizing diverse data sources,including borehole,oriented borehole,and photogrammetry data,along with hang-up frequency and hydrofracturing details,we applied discrete fracture network(DFN)modeling to analyze in-situ block size distribution and fragmentation.These results are based on 12,000 realizations of discrete fracture network(DFN)models using R-Dis-Frag computer pacakge at real cave volumes of 200 m200 m200 m,with varying parameters,which significantly enhances their reliability.The incorporation of DFN modeling and geostatistical simulation allows for capturing the interaction berween several spatial variables and explaining the variations observed in the production results at the draw points.Keyfindings of spatio-statistical analysis highlight the significance of volumetric fracture intensity(P32)and extraction column height in reducing hang-up events and enhancing fragmentation efficiency.The study integrates HF-induced and natural fracture intensities,revealing that higher P32 values and higher draw columns correlate with fewer hang-ups and better fragmentation.We recommend non-regular HF patterns for high P32 zones to improve operational efficiency.This research provides insights into optimizing mining operations,acknowledging the limitations of HF propagation efficacy and paving the way for further exploration into the interplay between hydraulic fracturing and natural discontinuities.
文摘Cleats are the dominant micro-fracture network controlling the macro-mechanical behavior of coal.Improved understanding of the spatial characteristics of cleat networks is therefore important to the coal mining industry.Discrete fracture networks(DFNs)are increasingly used in engineering analyses to spatially model fractures at various scales.The reliability of coal DFNs largely depends on the confidence in the input cleat statistics.Estimates of these parameters can be made from image-based three-dimensional(3D)characterization of coal cleats using X-ray micro-computed tomography(m CT).One key step in this process,after cleat extraction,is the separation of individual cleats,without which the cleats are a connected network and statistics for different cleat sets cannot be measured.In this paper,a feature extraction-based image processing method is introduced to identify and separate distinct cleat groups from 3D X-ray m CT images.Kernels(filters)representing explicit cleat features of coal are built and cleat separation is successfully achieved by convolutional operations on 3D coal images.The new method is applied to a coal specimen with 80 mm in diameter and 100 mm in length acquired from an Anglo American Steelmaking Coal mine in the Bowen Basin,Queensland,Australia.It is demonstrated that the new method produces reliable cleat separation capable of defining individual cleats and preserving 3D topology after separation.Bedding-parallel fractures are also identified and separated,which has his-torically been challenging to delineate and rarely reported.A variety of cleat/fracture statistics is measured which not only can quantitatively characterize the cleat/fracture system but also can be used for DFN modeling.Finally,variability and heterogeneity with respect to the core axis are investigated.Significant heterogeneity is observed and suggests that the representative elementary volume(REV)of the cleat groups for engineering purposes may be a complex problem requiring careful consideration.
