Abstract
An approach to model seismicity induced by mining was proposed by MDG Salamon and developed further by AM Linkov. We suggest a variation of the Salamon–Linkov method that does not require a priori assumptions about sizes, orientations and locations of flaws, reducing the amount of input data. The algorithm places a 3D grid over the model domain and searches the entire grid space for possible planes of failure by identifying lobes of positive Coulomb Failure Stress (CFS). The largest CFS lobes are then utilised as flaws. Another change from the original formulation is that the flaws are discretised into smaller sub flaws. This accounts for the complexity of the shapes of the flaws and provides more accurate modelling of their effect on the stress field. Some elements of the original framework of Salamon and Linkov, e.g. tensile failure and creep effects are not included into our implementation. This is done intentionally to keep the method simple at the current stage. The implemented method requires limited input data: models of mining steps, information about virgin stress, rock mass elastic properties and strength. The strength may be anisotropic and inhomogeneous, where specific Mohr–Coulomb failure criteria can be assigned to particular orientations (describing joint sets) and 3D surfaces (describing faults or dykes). The output of the method represents two catalogues - one for seismic and another for aseismic events. The catalogues contain origin times (mining steps), coordinates, seismic potencies and source mechanisms of events. The method was applied to real case studies at two mines. The Boundary Element Method (BEM) models of the mines were of the order of 200,000 Fictitious Force (FF) elements. The first model included seven mining steps and the second model consisted of 33 steps. The computations took seven and 50 hours respectively on a multi-core machine. In the first case the resulting catalogue includes more than 900 seismic events, whereas in the second case the catalogues contained approximately 3,400 seismic and 100 aseismic events. The modelled seismicity may be used as a reference for the interpretation of the observed seismicity. Analysing the discrepancies between the modelled and observed seismicity can also help to assess the validity of input parameters – characteristics of in situ stress and failure criteria. Another possible application of the suggested method is the testing of future mining scenarios. Although the prediction of seismic response for individual scenarios may be not accurate, the difference between the modelled seismicity of the scenarios may help to choose the optimal one in regards to seismic or ground motion hazard.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.