Another new app has been added to the Seismic Suite for the general analysis of short-term seismic responses. There are multiple features to analyse the short-term response to different triggers (events or blasts) in time and space. In later versions, this app will replace all of the tools in the previous “Omori Analysis Tools” app but version 1 currently only replaces the old “Response to Blasting” window.
There are two main aspects of the app, the simple response viewer and trigger assessment windows. New training videos have been uploaded for the new app.
This window is for the assessment of short-term responses. Triggers must be selected (ticked) from the list and response events for the selected triggers can be assessed with a number of tools.
Responses 3D — View the trigger locations and response events within range. Adjust the spheroid controls to include/exclude events related to the trigger. “Nearby events” are events just outside the spatial or temporal range of the trigger.
Time after Trigger — Chart the events as a function of time after trigger. Can plot events as a histogram or cumulatively. The time bin used for the histogram can be adjusted. When multiple triggers have been selected, can either view responses individually or stacked together. The best fit Modified Omori Law (MOL) is calculated automatically, if the calculation fails a message will appear in the control panel. You can override the MOL parameters (p, K and c) in the control panel.
Distance from Trigger — Chart the events as a function of distance to trigger. Options for X distance, Y distance, Z distance, Horizontal distance and 3D distance. You can normalise the chart by the number of events or by the volume under consideration. Normalisation by number of events changes the Y axis to a percentage of events from 0-100%. Normalisation by volume changes the X axis so that the volume increases linearly as the distance increases (expanding sphere). If the cumulative events is linear for the chart normalised by volume, this represents a constant event density.
Density 2D — View the distribution of events in the 2D plane (XY, XZ or YZ). Events are divided into spatial bins and then ranked from highest to lowest density. The grid points are coloured by cumulative events, where the accumulation is from highest to lowest density points.
Density 3D — Similar to the Density 2D plot, uses the same grid spacing controls. Isosurfaces are used to show the cumulative events distribution in 3D.
Trigger Assessment Window
This window is for the analysis of triggers; either blasts, events or user defined points in space and time. If the response viewer is to assess the area of exclusion and time of re-entry, the trigger assessment window is to assess which blasts should have an exclusion and which events should have an evacuation.
Triggers 3D — View triggers in 3D and use the marker style colours and scales to assess the location of triggers that typically have a response.
Trigger Summary Tables — These tables summarise triggers by various characteristics. For example, you can use these tables to assess what type of blasts have the biggest response, what magnitude events typically have aftershocks, how the responses to triggers vary by depth etc. The tables currently use the following parameters to assess potential triggers for exclusion or evacuation:
The root upgrade that coincides with the software upgrade beyond v5.9 includes a new app for creating and exporting models of mine geometry. Stope, cave and development geometry is a fundamental aspect of most geotechnical analysis. Mine geometry also varies over time and capturing these changes is critical in any back analysis or numerical modelling that investigates stability or monitoring parameters over time. This is a utility app to create models of mine geometry that can be exported to facilitate a wide range of applications such as:
Live Surveys — Mine geometry can be displayed as surfaces in General Analysis and other 3D views that automatically update with the current time filters. So when you backdate your seismicity to a year ago, the surveys will only show geometry present at the time.
Volume Calculations — Mining volume can be calculated over time. This allows seismicity and other parameters to be plotted as a function of mine geometry.
Distance Calculations — It is often important to know the distance of large events or damage locations to the nearest development or stope at the relevant time. This is like the distance to survey tool except time is also taken into account. So you can compute the closest stope at the time of damage or the closest development at the time of the large event.
Numerical Modelling Inputs — The geometry models can be exported in a number of formats, including Map3D .inp format.
The export options will gradually be expanded to more modelling formats such as FLAC, RS3 and Wave. Another possible application is to create block models of lithology or geotechnical domains based on lithology contact wireframes.
Training videos have been made to walkthrough the new app. See the page below.
We have started rolling out a new version of the mXrap software and root folder. Software versions 5.9 and above include a few interface changes. The new version is available on the download site but DO NOT download it until someone from the mXrap team has upgraded your root folder. Several root folder changes need to be made simultaneously with the software upgrade. We will be in touch soon to arrange the upgrade (if we haven’t already).
From now on, when you open an app, all available windows are shown at the top of the screen as tabs. The old Windows menu, where extra windows for an app could be opened in separate windows, is now gone. The workspaces, reload data, and undo/redo/save buttons have also moved slightly.
Note there is different styling for folders. Folders can be expanded to a 2nd level with extra tabs.
Most of the apps have had a bit of a face-lift to adjust to the new window tabs. The Hazard Assessment app for example has been split into multiple windows, for easier workflow. We have added info pages to most apps to describe the app and point to relevant literature.
We also have a bit more control over panels, to show, hide and disable certain instructions and controls. We think this will further improve usability going forwards.
New and improved apps
Aside from the interface changes, the root upgrade has also got some new stuff we have been working on over the last 12 months. We will make separate blog posts on each of the new items to explain further. Training videos for the new additions are in progress.
Moment tensors in General Analysis. Moment tensors are now read in automatically with IMS event files. ESG sites can load individual data files they have received from ESG. Beach balls can be displayed in General Analysis with various decompositions and scaling. There is a separate window with stereonets and mechanism charts.
Short-term Response Analysis app. This is an upgrade to the old “Response to Blasting” window in the “Omori Analysis Tools” app. There are more spatial analysis options in this new version. There are also trigger assessment tools for large events and blasts. Previously only blasts could be chosen as a trigger.
Mine Geometry Models app. This is a new app to create 4D models of mine geometry from surveys. Models can be used in General Analysis to automatically display mine geometry and update based on the current time filters. Models can also be exported as Map3D .inp files.
Seismic Monitoring app. This app has been updated. There are now separate apps for event monitoring with and without exclusions. There is also a window to monitor activity rate.
Survey Format Converter app. New utility app to export surveys in a different format. Can also use the decimation tool to export decimated surveys. We’ve also recently made a DXF export option. DXF is an option in this app as well as the plane fitting and mine geometry models apps.
The grid-based hazard calculations in the Hazard Assessment app were discussed in a previous post. The Iso View describes the hazard at all locations within the mine but when you are considering the seismic risk for a particular work area, large events and strong ground motions may come from multiple sources. The Excavation View estimates the seismic hazard associated with working areas (minode locations) in a few different ways as described below.
P [ ML within R ]
The P [ML within R]
marker style is the probability of exceeding the design magnitude within the
design distance (R) of the minode location (per year). This is the simplest
minode hazard estimate and you will notice the other marker styles take a bit
longer to calculate because they include more complex ground motion (PPV’)
As discussed in a previous post, the grid-based hazard calculations result in a probability of exceedance assigned to each grid volume. So, for a design magnitude of ML2, the annualised probability of exceeding ML2 within each grid cell volume is computed. To compute the probability of exceeding ML2 within R of a minode, the exceedance probabilities for all the grid cells within R are integrated together.
In the 2D example below, there are seven grid cells within R
of the minode. Let’s say each of these grid cells has a probability of
exceeding ML2 of 1 %. Then, the probability of exceeding ML2
within R of the minode would be:
P [ML2 within R] = 1 – (1 – 0.01)7 = 6.8 %
The hazard for a single minode doesn’t help you much. You really want to know the hazard for the whole mine or for any possible work area. To compute the probability of exceeding ML2 within R of multiple minodes, you just integrate all the grid cell probabilities within R of any minode. This is illustrated below for a tunnel length with three minodes. Each grid cell is only counted once. In the hazard app, the footer of the Excavation View shows the P[ML within R] for all minodes. You can also use the selections tool to select any combination of minodes and the integrated hazard will be shown in the footer. The Volume Hazard tool calculates the P[ML within R] for any minode within each of your filter volumes.
As you can see from the illustrations above, the grid cells that are included “within R” are based on the distance to the grid centre point. This is why we recommend using a value for R that is larger than the grid spacing you are using. If you specify R much less than the grid spacing, you may see some odd artefacts where some minodes have much fewer grid points associated than others.
Ground Motion Hazard
The other minode marker styles in the Excavation View express the hazard based on the probabilistic strong ground motion (PPV’). As discussed in the last blog post, a Strong Ground Motion relationship is required to calculate the probability distribution of PPV’ based on the event source details and R. For each minode, the ground motion hazard, P[PPV’], is the probability of exceeding the design PPV’, anywhere within the volume associated with the minode (the tunnel section, per year).
For each minode, there is a probability of exceeding the
design PPV’ due to a large event occurring within any particular grid cell. To
compute P[PPV’] for a minode, each grid point must be considered and the probabilistic
contribution of each integrated together. The grid cells closest to the minode
(smallest R) are the most likely to contribute the most to the minode P[PPV’]
but high hazard grid cells will have a larger zone of influence because of the
potential for very large events.
The grid-minode combinations with a very small P[PPV’] are ignored in the analysis to speed up the calculation. The minode P[PPV’] is calculated from the remaining grid point combinations. For each minode, the threshold probability to ignore a grid point can be specified from the control panel. Increasing the threshold is slightly less accurate but will increase calculation speed.
The P[PPV’] for each grid-minode combination is computed by discretising the magnitude, distance and PPV’ probability distributions. From the grid-based seismic hazard calculations, the magnitude distribution is known and the probability of each magnitude bin of large events can be calculated. Then, the distance from the grid cell to the minode volume is represented by the R distribution in a similar way to the figure below.
From every combination of ML and R bin, the strong ground motion relationship is used to calculate the full PPV’ probability curve for the minode. Only the magnitude bins that can possibly result in a PPV of interest are included in the calculations. This is controlled by the PPVmin control parameter. The PPV’ probability curve is plotted above PPVmin. The P[PPV’] can be quickly calculated for any design PPV’ values above PPVmin from the probability curve.
Note the example PPV’ probability curves in the figure below. The RED line is the probability curve for a single minode but as with the P[ML within R] calcs discussed above, you will often want to evaluate larger tunnel areas, or even the whole mine (BLUE). The PURPLE line is the probability curve for a small section of tunnel with multiple minodes and indicates the yearly probability of exceedance for PPV’ anywhere along that tunnel section. The GREEN line is the probability curve for the whole mining level.
Clearly the probability of exceedance will tend to increase as longer sections of tunnel are considered. So, even though the P[PPV’] for a single minode is small, the accumulated hazard when considering the whole mine can be much higher. This is why the hazard app displays the ground motion hazard in a couple of different ways other than the individual minode P[PPV’] to make the hazard ratings more intuitive.
The equi-probability zones marker style is a simple ranking of minodes from lowest to highest P[PPV’]. The marker style value is a percentile rather than a probability. For example, let’s say you have 100 minodes, the top 5 minodes with the highest hazard will be red (0.95-1.00). The next top 5 (6-10) highest hazard minodes would be orange (0.9-0.95) etc.
Note that the equi-probability zones do not illustrate an absolute hazard rating, rather illustrates relatively low hazard areas and relatively high hazard areas. There may be cases where the top 5% of minodes are actually low hazard. Conversely, the bottom 5% minodes may still be quite a high absolute hazard. The PPV probability chart plots the curves for each equi-probability zone.
The Cumulative P[PPV] marker style also ranks minodes from
lowest to highest P[PPV’] but the marker values are probabilities accumulated
from lowest to highest. As described previously, to accumulate the hazard for
multiple minodes, the probabilities are integrated together. For example, if
your have four minodes with individual P[PPV’] of 1%, 2%, 5% and 7%, then the probability
of exceeding the design PPV’ AT ANY of the four minodes is:
The Cumulative P[PPV] marker accumulates the probabilities
in the same way, one at a time, from lowest to highest hazard. So for example:
The lowest hazard minode has the same individual P[PPV] and Cum P[PPV].
The second lowest hazard minode has a Cum P[PPV] equal to the combined hazard for the two lowest minodes.
The minode with the median individual hazard has a Cum P[PPV] equal to the combined hazard integrating all lower 50% of minodes.
The highest hazard minode has a Cum P[PPV] equal to the accumulated hazard for all minodes.
The Cum P[PPV] marker has the benefit of showing both the cold-to-hot scale of all minodes as well as an indication of the absolute hazard.
The Probability Class marker style is similar to the
Cumulative P[PPV] marker style, except you can specify your own colour scale
and classes and give them a name. This can help to communicate changing hazard
areas to operational personnel.
You can see the classes in the risk matrix. You should notice that the Cum P[PPV] corresponds to Prob Classes according to the numbers in that matrix.
For more details on the ground motion hazard calculations see Wesseloo (2018).
The Strong Ground Motion (SGM) relationship is used to calculate the Peak Particle Velocity (PPV) generated by a seismic event. You may also hear this referred to as a Ground Motion Prediction Equation (GMPE), but only the maximum velocity is estimated, i.e. the strong ground motion, rather than the full, complex wave motion.
The PPV is generally calculated for a specific location based on:
• distance to the seismic event (R)
• source magnitude (ML)
• source radius (Rs)
• static stress drop (SSD)
The source radius can be computed as a function of magnitude and the adjustment due to SSD is sometimes excluded. So PPV is often simply a function of ML and R.
You have probably seen the SGM relationship illustrated in a similar way to the figure below. What is sometimes not recognised is there is an associated uncertainty in these relationships. In the case below, the PPV values are based on a 10% chance of exceedance for a given ML and R.
The wave motion from a seismic event through the rock mass is highly complex and uncertain. So, for a given ML and R, the PPV is not a single value but a probability distribution. This is illustrated below, for a ML2 at a distance of 100 m, the PPV distribution is plotted along with the 10th, 50th and 90th percentile values for probability of exceedance.
There are many factors that contribute to the uncertainty in ground motion that results from a seismic event. The ground motion does not radiate from a seismic event uniformly. Each source mechanism has its own radiation pattern where the magnitude of the ground motion varies depending on the direction. The radiation pattern also differs for the P-wave and S-wave.
As the body waves radiate outwards from the source they attenuate but rock mass anisotropy affects the rate of attenuation. The wave attenuates faster when cutting through lamination than it does when travelling along (parallel) to the bedding or foliation plane.
Excavations, different lithologies and major contacts and discontinuities will create reflections and refractions. The wave will split into a separate P-wave and S-wave when it reflects off a boundary or refracts into a new medium. Multiple waves can superimpose to create stronger ground motions that the individual waves.
The SGM relationship does not give the PPV expected on the excavation surface. The SGM equation is based on the recorded PPV at sensors which are normally installed well into the rock mass, away from the excavation surface. So, the calculated PPV includes the uncertainty from different radiation patterns, natural variability, reflections and refractions but does not include surface effects. The strong ground motion of the body wave, without including the effect of the excavation surface, is sometimes referred to as PPV’. The true surface PPV is the PPV’ with an additional amplification factor applied.
The amplification is due to a couple of different effects of the surface. The amplitude is expected to double at the free surface due to the superimposition of the incoming wave and the reflected wave. The amplification is more than double for a corner pillar due to the closed geometry. The body waves can also interact with the free surface to form Rayleigh and Love surface waves. These waves propagate along the surface rather than through the rock mass. The low velocity fractured zone around the excavation can enhance the formation of surface waves as the seismic energy is trapped between the free surface and the fracturing boundary. Waves also increase in amplitude as they move into a lower velocity medium such as the fractured zone.
In earthquakes, surface waves cause the most damage to infrastructure. The amplitude of Rayleigh and Love waves tend to be higher than body waves. Surface waves also attenuate more slowly, i.e. travel further, since the geometric attenuation is only along the surface rather than all three dimensions. Another common observation with earthquake damage is that buildings on soft soils are more heavily damaged than buildings on solid rock. This is a similar case as the fractured zone around excavations although there is less experimental evidence of this phenomenon in underground tunnels.
If you are interested in reading more, these papers are a good place to start:
• Wesseloo (2018) – Description of the SGM relationship and hazard calculations done in mXrap.
Procedure for site specific SGM calibration
The SGM relationship is used to calculate PPV in the Hazard Assessment and Large Event Analysis apps. The default relationship is from the Canadian Rockburst Support Handbook. This relationship is mostly based on recorded ground motions at Brunswick, El Teniente and Creighton mines and may not be applicable to your site. It is fairly simple to calibrate your own site specific SGM relationship using the data recorded by your seismic system. The PPV is recorded by each sensor for each event. With this data, we have a tool to calibrate your site specific SGM relationship. If you would like to do this for your site, the procedure is as follows:
1. Export your PPV data from your seismic database. We have done this a few times for IMS data. It is called an Event-Trigger query that is a table with a row for each sensor hit per event. Contact email@example.com if you need assistance. At minimum we need the export to include the sensor location, sensor PPV, event location, event magnitude and event static stress drop.
2. Save the PPV data into the #Data folder in your root and run a default backup in mXsync. Contact firstname.lastname@example.org to let us know you would like us to calibrate your SGM relationship.
3. We will generate your site specific SGM equation and send the info back as a patch in mXsync.
4. Apply the patch in mXsync and run another default backup.