Arma2015 000852
Arma2015 000852
ABSTRACT: This article proposes a methodology for building a rock mass model for the Toquepala Mine. This model will be
used as an input for analyses to evaluate the stability of the current and future open pit slopes. The rock mass modeling process can
be divided into two phases: a database building phase which involves gap analysis and data collection, compilation and validation;
and an analysis phase that involves data assessment and the definition of rock mass units (RMUs). This article briefly describes
each of these phases and how they have been applied at the Toquepala Mine, with emphasis on the data assessment and definition
of the different RMUs.
1. BACKGROUND In 2011, an update to the 2000 slope design criteria was
initiated by the mine. A Geotechnical Review Team
Toquepala is a large open pit copper mine owned and
(GRT) comprised of Itasca, Piteau and Mr. Peter Stacey
operated by Southern Copper Corporation. It is located
of Stacey MGL was formed to undertake the required
in Jorge Basadre Province in the south of Peru, about
studies. Following compilation of the available
160 km northeast of Tacna. The Toquepala Mine has
information, the GRT conducted a gap analysis which
been in operation since the late 1950s. In the 1980s,
was used to guide the design of a comprehensive field
Call and Nicholas Inc. (CNI) developed an early slope
data collection program that included bench mapping,
design. Later, in 2000, updated slope design criteria
geotechnical drilling (geomechanical logging, hydraulic
were prepared by a joint collaboration of Golder
testing and piezometer installation), and laboratory
Associates Ltd. (Golder), Itasca Consulting Group Inc.
testing. The program specifically targeted areas for
(Itasca) and Piteau Associates Engineering Ltd. (Piteau).
collection of structural fabric and fault orientations to
As part of this later study, the rock mass characterization
facilitate updating of the structural database and major
was updated.
fault models, geomechanical rock mass parameters for
The rock mass model developed for the 2000 slope updating the geomechanical database and rock mass
design study was based on a series of mapping, drilling model, and permeability estimates for updating the
and laboratory testing campaigns. The laboratory testing groundwater model.
program included a limited number of triaxial
The methodology used to compile and build the rock
compressive strength tests. Following compilation of
mass model is described in the following paragraphs.
these data, rock mass units (RMUs) were defined based
on lithology and the presence or absence of calcium
sulfates (i.e., gypsum and anhydrite). 2. DEPOSIT SCALE GEOLOGY [1]
During the last 15 years, the Toquepala pit has grown. The ore deposit at Toquepala is characterized as a copper
Pit slopes are currently on the order of 700 m high, and porphyry in which copper and molybdenum
the next phase of mine development (Phase 5) will mineralization is predominantly found in breccia
increase the slope height to approximately 950 m. Given intrusions (i.e., Angular Breccia (Bx), Diorite Breccia
this significant increase in slope height, and anticipated (Dibx), Tourmaline Breccia (Bxt) and Pebble Breccia
changes to the competency of the rock mass with (Px)), hosted within two igneous rocks (i.e., Quartz
increasing depth, it was considered prudent to update the Quellaveco Porphyry, (Qq) and an intrusion of Diorite
rock mass model and slope design. (Di)). Three additional, but barren, intrusions have been
identified at the site: Dacite Porphyry (Dp), Latite 3. AVAILABLE DATA AND GAP ANALYSIS
Porphyry (Lp) and Dacite Agglomerate (Da). A
The collaborative study conducted in 2000 included
sequence of younger volcanic rocks is located on top of
compilation of the available data and definition of new
the breccias and intrusives: Toquepala Andesite (Ta),
rock mass units (RMUs). This database included
Toquepala Dolerite (Td), Toquepala Quartz Porphyry
geomechanical core logging information collected from
(Tq), Toquepala Rhyolite (Tr) and Alta Andesite (Aa).
Accurate identification of these lithologies is challenging approximately 5000 m of core, laboratory testing of
intact rock core including 358 Unconfined Compressive
due to their complex interlayering. Figure 1 shows the
Strength (UCS), 43 Triaxial Compressive Strength
stratigraphic sequence as it is currently understood at the
(TCS) and 375 Splitting Tensile Strength (STS) tests,
Toquepala deposit.
and laboratory testing of discontinuity surfaces,
The main types of hydrothermal alteration produced by including 57 Direct Shear (DS) tests. Due to the length
the intrusion of the breccias and igneous bodies are of time over which the laboratory tests were conducted,
propylitic (Prop) and quartz-sericite (QS) alteration. The much of the original data sets were lost and as a result, it
intensity of the alteration has not been investigated in was not possible to review the quality of the test results
detail. or procedures. Table 1 summarizes the results of
Below approximately 3000 masl, gypsum and anhydrite laboratory testing in terms of the number of completed
are present as discontinuity infills and disseminated in tests and mean values.
the rock matrix as the result of calcium sulfate flooding Table 1. Summary of Laboratory Testing Data Available Prior
associated with hydrothermal processes. Their presence to the Current Study
has been identified in all of the rock types except the
sequence of younger volcanic rocks. The upper limit of Lithology UCS TCS STS DS
the presence of gypsum and anhydrite is known as the 129 100 12
Di 10
Tope de Yeso (TdY). 40MPa 6MPa 28° - 59kPa
15 18 3
Da -
67MPa 6MPa 23° - 105kPa
CHRONOSTRATIGRAPHIC
DIVISIONS LITHOSTRATIGRAPHIC DIVISIONS
10 18 2
Td -
61MPa 12MPa 24° - 79kPa
TOQUEPALA GROUP
Quellaveco Formation
TOQUEPALA RHYOLITE
11 20 1
Rhyolitic flows; white, grey, brown and red in colour.
Tq -
Prominent flow banding.
TOQUEPALA ANDESITE
94MPa 9MPa 28° - 134kPa
+- 900m
UNCONFORMITY
Andesitic flows, light grey to dark grey in colour.
15 17 4
+- 75m
Generally aphanitic in texture. Middle section is
Porphyritic with phenocrysts of feldspar. Magnetic. Tr -
Di
ksTi-Ttr
KsTi-Tta
73MPa 6MPa 25° - 15kPa
+- 60m
Toquepala TOQUEPALA QUARTZ PORPHYRY 35 48 10
Lp Rhyolitic flows, grey to greyish white in colour.
Bx 10
Series
KsTi-Ttq +-100m
Porphyritic texture with phenocrysts of quartz.
Generally well banded and tuffaceous. 78MPa 5MPa 29° – 53kPa
TOQUEPALA DOLERITE 32 32 3
KsTi-Ttd Andesitic flows, grey to greenish grey in colour. Px -
Da
aphanitic to fine grained. Flow banding generally
presents Agglomeratic towards the bottom, with 56MPa 7MPa 24°- 56kPa
boulders of Qq.
UNCONFORMITY
+- 75m
Note: The first number corresponds to the number of tests conducted prior to
QUELLAVECO QUARTZ PORPHYRY
Rhyolitic flows, white to greyish white in colour; hard dense.
2000 and the second number corresponds to the mean value.
Quellaveco Porphyritic texture with small, irregularly shaped
Quartz KsTi-Tqq phenocrysts of quartz. Characteristic aphanitic,
Porphyry Dp Bxs milky white groundmass.
As indicated in Table 1 the quantity of testing data for
some lithologic units is very limited. The lack of triaxial
testing on the younger volcanic units was identified as a
Di +-150m
INTRUSIVE AND HYDROTHERMAL ROCKS
Lp
Latite Porphyry particular weakness in the database as triaxial testing
data is a key input in the evaluation of the intact and
Da
Dacite Agglomerate
Breccias
Dp
Di
Dacite Porphyry
Diorite
rock mass shear strength of the material.
Fig. 1 Stratigraphic sequence at Toquepala (modified from Some of the strength values in Table 1 appear low in
Luciano Prieto, 1977). comparison to the hardness and general rock mass
competency observed in the current pit. These reduced
strengths are attributed, in part, to sample location. Most
of the historic samples were obtained at higher
elevations and closer to the core of the deposit, and have
therefore likely been exposed to a greater degree of • Use of a secondary core orientation device (e.g.,
weathering and alteration. The current mining phase is ACT II, manufactured by Reflex™) as a
deeper, and the walls are farther away from the core of comparative check on approximately 25% of the
the deposit and are therefore less influenced by core,
weathering rock. • Geomechanical core logging for all drillholes
As one of the considerations used in the gap analysis, the using Rock Mass Rating, RMR (Bieniawski,
volume of each lithology that will be present in the pit 1976),
walls of the future mine phases (i.e., Phases 4 and 5) was • Point Load testing in all drillholes, for
estimated using Toquepala’s geological block model. correlation of Point Load Index (Is(50)) with
This allowed prioritization of field data collection laboratory UCS values,
methods to target the most predominant lithologies. The • Targeted laboratory sampling and testing
results of this assessment indicated the following units, program,
according to their priority: • Hydraulic testing at nominal 100m spacing using
double or triple packer system in each drillhole,
• Diorite, above and below TdY (38 and 28%, and
respectively), • Installation of two to three vibration wire
• All intrusive rocks, excluding Di (12%), piezometers in selected drillholes.
• All volcanic rocks (8%),
• All breccias (3%), including Pebble Breccia or
Px, which was targeted regardless its low 5. 2014 LABORATORY TESTING PROGRAM
percentage on the Phases 4 and 5 pit as it has The pre 2000 laboratory testing database and prioritized
been historically identified as an important list of lithologies expected in the Phase 4 and 5 pit walls
factor on the development of interramp level were used to prepare sampling targets for the laboratory
instabilities in the past). testing program. As part of this program, UCS, TCS,
It is important to note that Da is the most predominant STS and DS tests were identified as being required to
intrusive rock besides Di, and makes up 7% of the 12% obtain a reasonable characterization of the strength of
reported. the rock mass and joint fabric. The GRT provided the
target number of samples for each test type.
Based on review of this information, the GRT made
recommendations for a data collection program that As part of ongoing quality control efforts for the project,
targeted each rock type, the presence or absence of it was also recommended that approximately 20% of the
gypsum/anhydrite, and alteration. One of the goals of laboratory tests be completed at a different laboratory.
the program would be to determine whether alteration The primary laboratory testing program was conducted
has an important role in the mechanical behavior of the at the SGS Rock Mechanics Laboratory in Santiago,
rock mass. Chile, with the remained of the testing being conducted
at the Mecánica de Roca Ltda. laboratory in Calama,
Chile. Project engineers regularly visited the
4. 2013 DRILLHOLE PROGRAM laboratories during testing to confirm that the specified
The 2013 drillhole program consisted of 38 oriented testing procedures were being followed. [3].
boreholes distributed throughout the pit area and Following reception of laboratory test results, all results
targeted the Phase 4 and 5 pit walls, the most important were validated by reviewing the raw data, before and
lithologies, and major faults. In particular, four after photos, and comments and observations provided
drillholes were designed to help define the persistence, by the laboratories. The following list summarizes the
character and thickness of major faults. number of valid tests completed during the 2014
In addition to the proposed geotechnical drillholes, the program:
GRT also prepared recommendations to ensure the best • UCS: 154 (120 and 34 conducted at SGS and
quality of drilling. These recommendations included Mecánica de Rocas, respectively),
[2]: • TCS: 141 (133 and 8 conducted at SGS and
• Use of triple tube wireline drilling equipment to Mecánica de Rocas, respectively),
minimize core disturbance, • STS: 263 (6 and 257 conducted at SGS and
• 12hr/day drilling supervision, Mecánica de Rocas, respectively), and
• Geomechanical core logging at the drill platform • DS: 151 (121 and 30 conducted at SGS and
in core tube splits, Mecánica de Rocas, respectively).
• Geophysical surveys (i.e., acoustic and optical The current study also included an assessment of the
televiewer) to obtain discontinuity orientations impact of alteration impact on the strength of the intact
in all drillholes, rock and rock mass. As part of this assessment,
laboratory test results were initially grouped according (i.e., friction angle between 23 and 29°). As indicated
lithology, the presence or absence of gypsum/anhydrite, above, it was thought that these samples may have been
and alteration type. selected from higher pit elevations and as such, may
have been adversely affected by a higher degree of
Table 2 shows the distribution of number of valid tests
weathering, due to proximity to the surface, and
conducted for each geological combination.
alteration, due to proximity to the core of the deposit.
As part of the GRT’s recommendations, multi-stage
Table 2. Distribution of Valid Results from the 2014 direct shear testing of discontinuities at 250, 500, 750
Laboratory Testing Program and 1000 kPa normal loads, were conducted to
characterize the non-linear behavior of representative
Lith. TdY Alt. UCS TCS STS DS
discontinuities under different levels of confinement.
Prop 18 19 26 11 From a total of 160 tests, 151 were considered valid and
Above
Di
QS 2 4 5 5 were used to establish design discontinuity shear
Prop 10 9 16 10 strengths based on the following attributes:
Below
QS 11 7 18 9
Prop 7 6 7 8
• Rock Type,
Above • Presence or Absence of Gypsum/Anhydrite,
QS 1 3 3 0
Da
Below
Prop 10 6 13 9 • Alteration, and
QS 0 3 6 2 • Joint Condition (JC76).
Above QS 1 3 5 1
Dp Prop 0 2 3 1 Results of this assessment indicated that rock type and
Below
QS 7 5 14 8 the presence or absence of gypsum/anhydrite have the
Prop 7 2 6 3 biggest impact on discontinuity shear strength.
Above
QS 3 2 5 2
Lp Alteration did not appear to have a significant effect. As
Prop 0 0 2 2
Below
QS 1 2 3 0 a result of this assessment, the shear strength of joint
Prop 0 0 4 1 discontinuities was characterized according to rock type
Qq Above
QS 6 5 12 10 and location relative to the TdY.
Aa Above Prop 12 11 19 7
Ta Above Prop 5 0 7 4 Figure 2 provides an example of an analysis conducted
Td Above Prop 4 4 9 4 for Di located above the TdY. The first segment of the
Tq Above
Prop 0 0 2 0 line shown in Fig. 2 (in red) is defined by the best fit of
QS 1 2 0 0 the residual strengths from the tests conducted at the
Tr Above Prop 10 13 19 9 lowest normal load (i.e., 250kPa). Residual strengths
Prop 2 5 6 7
Above
QS 0 3 6 0
did not consider cohesion and were intended to represent
Bx the behavior of structures close to surface (i.e.,
Prop 3 3 4 8
Below equivalent to bench scale) where discontinuities tend to
QS 10 7 17 9
Dibx Below QS 5 1 6 7 be open and shear strengths are adversely affected by
Prop 8 3 6 4 blasting and excavation. The second segment of the line
Above
QS 3 6 8 1 in Fig. 2 (in blue) is defined by the best fit of the peak
Px
Prop 5 3 4 9
Below
QS 2 2 2 0
strengths from the tests conducted at all the normal
loads. For peak strength, cohesion was included and this
was considered to represent the behavior of the
As can be seen in Table 2, in comparison to the prior structures at depth, equivalent to a multi-bench scale.
testing database (Table 1) the current database is much Discontinuities at depth are less likely to be affected by
large and better distributed according to rock type. blasting and excavation.
In addition to estimating the best lineal fit for residual
6. DISCONTINUITY STRENGTH and peak strengths, lower and upper bound strengths
In rock slope design, it is common to use the Mohr were estimated by using +/- one standard deviation of
Coulomb failure criterion to express the shear strength of the data set. These estimates are intended to support the
discontinuities as a linear envelope in terms of cohesion, upcoming stability analysis stage of study.
friction angle, normal stress (σn) and shear stress (τ) [4].
The slope design studies conducted in the 1980s
included some shear strength test results for some
lithologies, based on both large and small scale shear
box testing. Results of this early testing program are
summarized in Table 1. Assessment of this original data
set indicated relatively low discontinuity shear strengths
intact rock core and rock mass characterization values
1500
DIRECT SHEAR TEST RESULTS AND BILINEAL FITS
Normal and Shear Residual Strength (from tests at 250kPa)
1400
Bench Scale Lineal Fit
Normal and Shear Peak Strength (from tests at 250, 500, 750 and 1000kPa)
(i.e., GSI) are also required to define this envelope. [5]
1300 Multi-Bench Scale Lineal Fit
1200
Range Within Friction Angle and Cohesion Could Vary
7.1. Intact Rock Strength
Characterization of the intact rock strength utilized data
1100
from the 2014 laboratory testing program. Previous data
1000
were not considered applicable to the current rock mass
900 because of the different intensity of weathering and
800 alteration between samples located close to the surface
700 and the ones located at depth. It was also not possible to
600
validate earlier test results.
500 Tables 4 to 7 summarize the mean UCS and STS results
400 (based on valid tests) by lithology, the presence or
300
absence of gypsum/anhydrite, and alteration type. For
Bench Scale:
Lower Bound: 0kPa - 39°
Best Fit: 0kPa - 43°
each of these combinations, at least three triaxial tests
200 Upper Bound: 21kPa - 48°
Multi-Bench Scale:
were conducted at three different confining stresses (i.e.,
100 Lower Bound:
Best Fit:
59kPa - 29°
135kPa - 35°
Upper Bound: 212kPa - 41°
at approximately 15, 25 and 50% of design UCS defined
0
for the 2000 study).
0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 1500
NORMAL STRESS (kPa) Table 4. Mean Values of UCS and STS from 2014 Laboratory
Testing Program for Diorite
Fig. 2 Example of bi-linear fit (bench and multi-bench
scale) for Diorite located above TdY.
Lithology TdY Alteration UCS STS
Table 3 presents the Mohr Coulomb parameters (i.e.,
friction angle and cohesion) determined from the best fit Prop 192MPa 5MPa
(bench and multi-bench scale) for the different rock Above
types above and below the TdY. QS 168MPa 5MPa
Di
Table 3. Mohr Coulomb Parameters (Best Fit) for Rock Prop 166MPa 5MPa
Type and Presence or Absence of Gypsum/Anhydrite Below
QS 132MPa 4MPa
Bench Multi-Bench
Rock Tope de Scale Scale
Type Yeso
φ Cohesion φ Cohesion
Table 5. Mean Values of UCS and STS from 2014 Laboratory
Above 42° 0kPa 36° 61kPa Testing Program for Intrusives
Diorite
Below 43° 0kPa 35° 135kPa Lithology TdY Alteration UCS STS
Prop - 7MPa
Tq Above
QS 122MPa -
Prop - 4MPa
Qq Above
QS 139MPa 5MPa
REFERENCES
1. Stacey, P.F., L.J. Lorig, and P.M. Hawley. Toquepala
Pit Slope Design Recommendations, February 2000,
Golder Associates, Itasca S.A. and Piteau Associates,
20–28.
2. Capdeville, F., and B.W. Gilmore. Recommendations
and Geotechnical Support for the 20,000m
Geotechnical Drilling Program for the Prefeasibility
Pit Slope Designs for Phases 4 and 5, 2012,
Geotechnical Review Team.
3. Geotechnical Review Team. Protocolos para el
Desarrollos del Programa de Pruebas de Laboratorio,
January 2014.
4. Wyllie, D.C and C.W. Mah. Rock Slope Engineering,
2004, ed. Spon Press, 79–80.
5. Hoek, E., C. Carranza-Torres, and B. Corkum. Hoek-
Brown Failure Criterion – 2002 Edition. In Proceedings
of the NARMS-TAC Conference, Toronto, 2002. 1, 267
– 273.