Assessment of Fault Slip Potential from Sill Pillar Mining

Problem Statement

Note

To view this project in 3DEC, use the menu command Help ► Examples…. Choose “3DEC/ ExampleApplications/ fault_slip” and select “seis.prj” to load. The main data file used is shown at the end of this example. All data files are found in the project.

A 10-meter thick vertically dipping orebody is being mined by cut-and-fill methods at a depth of 950 to 1050 meters. Mining has taken place on two different levels, resulting in the formation of a sill pillar. The plan is to mine the sill pillar from the base up, in four cut-and-fill stages. Two vertical faults intersect the sill pillar, and it is necessary to determine the potential for slip on these features as mining of the sill pillar proceeds. Figure 1 shows the shape of the orebody in plan and long section, along with the locations of the faults, and the magnitude and the orientation of the in-situ stresses.

The properties of the ore, host rock and fill are outlined:

  ore host rock backfill
density (kg/m3) 2700 2700 2000
bulk modulus (GPa) 17.6 3.2 0.03
shear modulus (GPa) 13.2 9.9 0.014
cohesion (MPa) 10 8 0.05
friction angle (degrees) 43 43 40

The properties of the faults are:

normal stiffness (N/m): 1 × 1011

shear stiffness (N/m): 1 × 1011

cohesion (MPa): 0

friction angle (degrees): 30

../../../../_images/faultslip-problem.png

Figure 1: Diagram of the orebody, sill pillar, and faults.

Modeling Procedure

The model is generated with the file “SEIS.DAT”. The first step in building the model is to define the horizontal limits of the orebody by using the block cut joint-set, block hide, and block hide off commands to cut up the model block into smaller blocks. After defining each arm of the orebody independently, the orebody as a whole (running the full vertical height of the model) is marked as a separate region from the host rock. The vertical limits of the lower stope, upper stope and sill pillar cuts are then defined and marked as separate regions by hiding the entire model except the orebody, and creating individual horizontal cuts with the block cut joint-set command. The individual stopes and cuts are marked as separate regions by using the block hide, block hide off, and block group commands. Figure 2 is a perspective plot of the upper and lower stopes in the 3DEC model.

Once the orebody and stopes are defined, the two faults are cut and the model is divided into inner, middle and outer regions for zoning purposes. Any cuts in the model must be made before zoning occurs. The inner region, corresponding to the volume surrounding the orebody, is zoned with constant-strain elements with a 5-meter edge length. The middle and outer regions are zoned with 20- and 80-meter edge length zones, respectively. After zoning the model, all of the blocks surrounding the faults are joined since they are only an artifact of the model construction process and do not represent actual discontinuities in the rock mass.

In the model, the host rock, ore and fill are all assumed to behave as Mohr-Coulomb materials. The faults are assumed to behave according to the default joint constitutive model, which involves point area contact and elastic/plastic behavior with Coulomb slip failure. All of the model boundaries are pinned by applying a zero-velocity boundary condition in all coordinate directions using the block gridpoint apply command. In-situ stresses are applied to the model to match the gravitational gradient.

After the model has been brought to an initial equilibrium state, the upper and lower stopes are excavated using the block excavate command, and the model is run to equilibrium. These stopes are then filled, and the sill pillar cuts are excavated and filled in sequence. The block fill command is used to fill stopes. As part of this command, the constitutive models and properties for the blocks and joints making up the fill may be changed from the original material. For this model, the properties are changed to that of the fill, while the joint is glued (jmodel elastic). This prevents slip from occurring on the faults where they pass through the fill. Following excavation of each cut, the model is run to equilibrium before filling, to simulate the delay in filling that occurs in reality. Following mining of the upper and lower stopes, and after each sill pillar cut, the displacements on the faults are reset to zero so that the influence of each mining step on the faults may be determined. This does not affect the solution in any way, since velocities, rather than displacements, are used in the solution process. The displacements are calculated for the convenience of the user.

The model contains 65,383 zones and requires 13 MB of RAM.

../../../../_images/seis-stopes.png

Figure 2: Perspective plot of the upper and lower stopes in the 3DEC model.

Results

Figures 3 and 4 show the predicted shear displacements and area of slip on Fault #1 and Fault #2, following mining of the upper and lower stopes. The view of the surface of Fault #1 is from the southwest, and the view of the surface of Fault #2 is from the east. Areas on the fault that have slipped are indicated by a hollow circle, and the shear displacement is contoured. As seen, areas of slip and shear displacement are predicted on the faults where they intersect both sides of the two excavated stopes. This is expected, since the blocks of rock are much freer to move relative to one another at the excavation boundary. The maximum shear displacement of 1.5 cm occurs where Fault #2 intersects the upper stope. In the 3DEC model, the upper and lower stopes are excavated in a single step, then filled. In reality, the stopes were excavated and filled in an incremental fashion. So the shear displacement predicted for this stage would be expected to occur incrementally as mining progressed. The sill cut, however, is excavated and filled in an incremental fashion in the model as well, so that displacements can be determined as a function of mining step.

Figures 5 and 10 show the predicted shear displacements and area of slip on both faults for the four sill pillar cuts. For the first 2 cuts (Figures 5 to 8), shear displacements are very small (2-8 mm) and are limited to the areas on the faults adjacent to the current cut. On the third cut (Figures 9 and 10), slip and shear displacement also occur in the ore between the current cut and the bottom of the upper stope. Shear displacements are on the order of 1 cm, but the area of slip is less than 50 m2. When the fourth and final cut is taken (Figures 11 and 12), a very large area of slip (approximately 1850 m2 in size), with a maximum shear displacement of 2.4 cm, is predicted on Fault #2 where it intersects the southern boundary of the excavation.

The direction of slip on the faults after the fourth sill cut is shown on a plan section through the model at the mid-height of the fourth sill cut in Figure 13. The sense of slip agrees with the direction of the major principal stresses shown on the same section in Figure 14.

Discussion

The potential energy output resulting from slip on Fault #2 at the southern boundary of the excavation due to mining of the fourth sill cut may be estimated from the seismic moment M0 (Aki and Richards 1980), which is defined by

\[M_0=Gu_sA\]

where \(G\) is the rock mass shear modulus, \(_s\) is the average slip and \(A\) is the area over which slip occurs. Based on a rock-mass shear modulus of 9.9 GPa, an average slip of 0.8 cm (a third of the maximum shear displacement) and a slip area of 1850 m2, a seismic moment of 1.47 × 1011 N-m is predicted.

The Richter magnitude \(M\) of a seismic event may be related to the seismic moment by the following expression (Hanks and Kanamori 1979):

\[M={{2}\over{3}}logM_0-10.7\]

The units of \(M_0\) in this equation are dyne-cm. For a seismic moment of 1.47 × 1018 dyne-cm on Fault #2, the preceding equation gives a Richter magnitude of 1.4. Given the potential for an event of this magnitude, and the proximity of the potential rockburst to the working area, the design of a support system to resist rockburst damage in this mining scenario should be carefully considered. Other mining methods could be examined with the model to determine their potential for reducing the likelihood of large fault-slip seismic events.

References

Aki, K., and P. G. Richards. Quantitative Seismology. W. H. Freeman & Co., San Francisco (1980).

Hanks, T. C., and H. Kanamori. “A Moment Amplitude Scale,” J. Geophys. Res., 84, 2348-50 (1979).

../../../../_images/stopes-fault1.png

Figure 3: 3DEC plot of shear displacement on Fault #1 following extraction of the upper and lower stopes (view is from the southwest).

../../../../_images/stopes-fault2.png

Figure 4: 3DEC plot of shear displacement on Fault #2 following extraction of the upper and lower stopes (view is from the east).

../../../../_images/sill1-fault1.png

Figure 5: 3DEC plot of shear displacement on Fault #1 following extraction of sill pillar cut #1 (view is from the southwest).

../../../../_images/sill1-fault2.png

Figure 6: 3DEC plot of shear displacement on Fault #2 following extraction of sill pillar cut #1 (view is from the east).

../../../../_images/sill2-fault1.png

Figure 7: 3DEC plot of shear displacement on Fault #1 following extraction of sill pillar cut #2 (view is from the southwest).

../../../../_images/sill2-fault2.png

Figure 8: 3DEC plot of shear displacement on Fault #2 following extraction of sill pillar cut #2 (view is from the east).

../../../../_images/sill3-fault1.png

Figure 9: 3DEC plot of shear displacement on Fault #1 following extraction of sill pillar cut #3 (view is from the southwest).

../../../../_images/sill3-fault2.png

Figure 10: 3DEC plot of shear displacement on Fault #2 following extraction of sill pillar cut #3 (view is from the east).

../../../../_images/sill4-fault1.png

Figure 11: 3DEC plot of shear displacement on Fault #1 following extraction of sill pillar cut #4 (view is from the southwest).

../../../../_images/sill4-fault2.png

Figure 12: 3DEC plot of shear displacement on Fault #2 following extraction of sill pillar cut #4 (view is from the east).

../../../../_images/sill4-sheardisp.png

Figure 13: 3DEC vector plot of shear displacement on the faults following extraction of sill pillar cut #4 (plan section at the mid-height of sill pillar cut #4).

../../../../_images/sill4-stress.png

Figure 14: 3DEC vector plot of principal stresses following extraction of sill pillar cut #4 (plan section at the mid-height of sill pillar cut #4).

Data File

seis.dat

model new
model random 10000
model large-strain on

block create brick 880 1105 900 1130 3890 4110
block group 'rock' slot 'geology'

;create orebody
block cut joint-set dip 90 dip-direction 90 ori 950 0 0
block hide range pos-x 880 950
block cut joint-set dip 90 dip-direction 90 ori 1000 0 0
block hide range plane pos-x 1000 1105
block cut joint-set dip 90 dip-direction 0 ori 0 1000 0
block hide range pos-y 900 1000
block cut joint-set dip 90 dip-direction 0 ori 950 1010 0
block hide range pos-y 1010 1130
block group 'orebody' slot 'geology'

block hide off
block cut joint-set dip 90 dip-direction 45 ori 1035 1035 0
block hide range plane dip 90 dip-direction 45 ori 1035 1035 0 above
block cut joint-set dip 90 dip-direction 225 ori 1000 1000 0
block hide range plane dip 90 dip-direction 225 ori 1000 1000 0 above
block cut joint-set dip 90 dip-direction 135 ori 1000 1000 0
block hide range plane dip 90 dip-direction 135 ori 1000 1000 0 above
block cut joint-set dip 90 dip-direction 315 ori 1028 1042 0
block hide range plane dip 90 dip-direction 315 ori 1028 1042 0 above
block group 'orebody' slot 'geology'

block hide off range group 'orebody'

;create upper stope
block cut joint-set dip 0 dip-direction 0 ori 0 0 4000
block cut joint-set dip 0 dip-direction 0 ori 0 0 4050
block group 'upper stope' slot 'excavation' range pos-z 4000 4050

;create lower stope
block cut joint-set dip 0 dip-direction 0 ori 0 0 3950
block cut joint-set dip 0 dip-direction 0 ori 0 0 3980
block group 'lower stope' slot 'excavation' range pos-z 3950 3980

;create sill pillars
block cut joint-set dip 0 dip-direction 0 ori 0 0 3985
block group 'sill pillar 1' slot 'excavation' range pos-z 3980 3985
block cut joint-set dip 0 dip-direction 0 ori 0 0 3990
block group 'sill pillar 2' slot 'excavation' range pos-z 3985 3990
block cut joint-set dip 0 dip-direction 0 ori 0 0 3995
block group 'sill pillar 3' slot 'excavation' range pos-z 3990 3995
block group 'sill pillar 4' slot 'excavation' range pos-z 3995 4000
block hide off

;create faults
block cut joint-set dip 90 dip-direction 45 ori 975 1005 0 jointset-id 100
block cut joint-set dip 90 dip-direction 90 ori 1017.5 1017.5 0 jointset-id 200

;fine zoning around orebody
block cut joint-set dip 90 dip-direction 90 ori 940 0 0
block hide range pos-x 880 940
block cut joint-set dip 90 dip-direction 90 ori 1045 0 0
block hide range pos-x 1045 1105
block cut joint-set dip 90 dip-direction 0 ori 0 960  0
block hide range pos-y 900 960
block cut joint-set dip 90 dip-direction 0 ori 0 1070 0
block hide range pos-y 1070 1130
block cut joint-set dip 0 dip-direction 0 ori 0 0 3950
block hide range pos-z 3890 3950
block cut joint-set dip 0 dip-direction 0 ori 0 0 4050
block hide range pos-z 4050 4110

block zone generate edgelength 5

;intermediate zoning
block hide off
block cut joint-set dip 90 dip-direction 90 ori 910 0 0
block hide range pos-x 880 910
block cut joint-set dip 90 dip-direction 90 ori 1075 0 0
block hide range pos-x 1075 1105
block cut joint-set dip 90 dip-direction 0 ori 0 930  0
block hide range pos-y 900 930
block cut joint-set dip 90 dip-direction 0 ori 0 1100 0
block hide range pos-y 1100 1130
block cut joint-set dip 0 dip-direction 0 ori 0 0 3920
block hide range pos-z 3890 3920
block cut joint-set dip 0 dip-direction 0 ori 0 0 4080
block hide pos-z 4080 4110
block zone generate edgelength 20

;outer zoning
block hide off
block zone generate edgelength 80

model save 'seis-zoned'

================================================================

;join all blocks surrounding faults
block hide off
block hide range plane dip 90 dip-direction 45 ori 975 1005 0 above
block hide range pos-x 1017.5 1105 
block join on
block hide off
block hide range plane dip 90 dip-direction 45 ori 975 1005 0 below
block hide range pos-x 1017.5 1105 
block join on
block hide off
block hide range plane dip 90 dip-direction 45 ori 975 1005 0 below
block hide range pos-x 880 1017.5 
block join on
block hide off
block hide range plane dip 90 dip-direction 45 ori 975 1005 0 above
block hide range pos-x 880 1017.5 
block join on
block hide off
;
;block properties
block zone cmodel assign mohr-coulomb
block zone property bulk 13.2e9 shear 9.9e9 density 2700 ...
                    cohesion 8e6 friction 43 
block zone property bulk 17.6e9 shear 13.2e9 ...
                    density 2700 cohesion 10e6 ...
                    friction 43 range group 'orebody'

;joint properties
block contact property stiffness-normal 1e11 stiffness-shear 1e11 ...
                       friction 30 
block contact material-table default property stiffness-normal 1e11  ...
                                              stiffness-shear 1e11 friction 30 

;boundary conditions
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-x 880
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-x 1105
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-y 900
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-y 1130
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-z 3890
block gridpoint apply velocity-x 0 velocity-y 0 velocity-z 0 ...
                      range position-z 4110
;stresses
model gravity 0 0 -9.81
block insitu stress -198.65e6 -264.87e6 -132.44e6 0 0 0 ...
             gradient-z 39.73e3 52.947e3 26.487e3 0 0 0
;
block mechanical damping global
model history mechanical unbalanced-maximum
;initial equilibration
model solve
model save 'seis-ini'

========================================================

;sequential excavating and filling

block mechanical damping local

block excavate range group 'upper stope' or 'lower stope'
model solve
model save 'seis-1'
;
block fill 
block contact jmodel assign elastic range group 'upper stope' or 'lower stope'
block contact property stiffness-normal 1e11 stiffness-shear 1e11 ...
              range group 'upper stope' or 'lower stope'
block zone cmodel assign mohr-coulomb range group 'upper stope' or ...
                                                  'lower stope'
block zone property bulk 30e6 shear 14e6 density 2000 cohesion 50e3 ...
                    friction 40 range group 'upper stope' or 'lower stope'

block contact reset displacement
block excavate range group 'sill pillar 1'
model solve
model save 'seis-2'

block fill
block contact jmodel assign elastic range group 'sill pillar 1'
block contact property stiffness-normal 1e11 stiffness-shear 1e11 ...
                       range group 'sill pillar 1'
block zone cmodel assign mohr-coulomb range group 'sill pillar 1'
block zone property bulk 30e6 shear 14e6 density 2000 cohesion 50e3 ...
                    friction 40 range group 'sill pillar 1'
  
block contact reset displacement
block excavate range group 'sill pillar 2'
model solve
model save 'seis-3'

block fill 
block contact jmodel assign elastic range group 'sill pillar 2'
block contact property stiffness-normal 1e11 stiffness-shear 1e11 ...
              range group 'sill pillar 2'
block zone cmodel assign mohr-coulomb range group 'sill pillar 2'
block zone property bulk 30e6 shear 14e6 density 2000 cohesion 50e3 ...
              friction 40 range group 'sill pillar 2'
  
block contact reset displacement
block excavate range group 'sill pillar 3'
model solve
model save 'seis-4'

block fill 
block contact jmodel assign elastic range group 'sill pillar 3'
block contact property stiffness-normal 1e11 stiffness-shear 1e11 ...
              range group 'sill pillar 3'
block zone cmodel assign mohr-coulomb range group 'sill pillar 3'
block zone property bulk 30e6 shear 14e6 density 2000 cohesion 50e3 ...
              friction 40 range group 'sill pillar 3'
  
block contact reset displacement
block excavate range group 'sill pillar 4'
model solve
model save 'seis-5'
;

program return