23.12.2012 Views

Gotthard base tunnel rock burst phenomena in a fault zone ...

Gotthard base tunnel rock burst phenomena in a fault zone ...

Gotthard base tunnel rock burst phenomena in a fault zone ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

World Tunnel Congress 2008 - Underground Facilities for Better Environment and Safety - India<br />

<strong>Gotthard</strong> <strong>base</strong> <strong>tunnel</strong> <strong>rock</strong> <strong>burst</strong> <strong>phenomena</strong> <strong>in</strong> a <strong>fault</strong> <strong>zone</strong>, measur<strong>in</strong>g and<br />

modell<strong>in</strong>g results<br />

He<strong>in</strong>z Hagedorn & Rolf Stadelmann<br />

Amberg Eng<strong>in</strong>eer<strong>in</strong>g AG, Switzerland<br />

Stephan Husen<br />

Swiss Seismological Service (SED) ETH Zurich, Switzerland<br />

SYNOPSIS: Dur<strong>in</strong>g construction of the multifunction section Faido consist<strong>in</strong>g of up to 5 parallel <strong>tunnel</strong>s an<br />

unknown <strong>fault</strong> <strong>zone</strong> has been encountered. Large deformations occurred dur<strong>in</strong>g excavation especially <strong>in</strong> a<br />

<strong>tunnel</strong> cross<strong>in</strong>g the <strong>fault</strong> kernel. The excavation of the <strong>tunnel</strong> system has been given rise to a large area of<br />

stress redistributions, caus<strong>in</strong>g stress concentrations <strong>in</strong> the hard <strong>rock</strong> adjacent to the <strong>fault</strong>. Dur<strong>in</strong>g 2005 <strong>rock</strong><br />

<strong>burst</strong>s started to occur at the excavation faces. The <strong>in</strong>tensity of the <strong>rock</strong> <strong>burst</strong> <strong>in</strong>creased with ongo<strong>in</strong>g<br />

construction. The <strong>tunnel</strong>’s owner ATG (AlpTransit <strong>Gotthard</strong> Base Tunnel AG) decided to monitor the<br />

seismic activity by a dense network of seismic stations. The monitor<strong>in</strong>g campaign <strong>in</strong>clud<strong>in</strong>g result<br />

<strong>in</strong>terpretation has been realized by the SED (Swiss Seismological Service, Swiss Federal Institute of<br />

Technology). The seismic measur<strong>in</strong>g results enabled to determ<strong>in</strong>e the micro tremor’s sources and to derive<br />

the dynamic load to be applied on numerical models. The numerical <strong>in</strong>vestigations served to assess the<br />

additional load<strong>in</strong>g of the <strong>tunnel</strong>’s f<strong>in</strong>al concrete l<strong>in</strong><strong>in</strong>g and the possible risk aris<strong>in</strong>g due to a micro seismic<br />

event dur<strong>in</strong>g operation of the <strong>tunnel</strong>.<br />

Amberg Eng<strong>in</strong>eer<strong>in</strong>g AG is member of the Eng<strong>in</strong>eer<strong>in</strong>g JV <strong>Gotthard</strong> Base Tunnel South, GBTS, consist<strong>in</strong>g<br />

of: PÖYRY, Lombardi AG and Amberg Eng<strong>in</strong>eer<strong>in</strong>g AG.<br />

1. INTRODUCTION<br />

1.1 Overview<br />

The longest <strong>tunnel</strong> of the NEAT (New Alp<strong>in</strong>e<br />

Transverse) through the Swiss Alps is the <strong>Gotthard</strong><br />

Base Tunnel (GBT) with a length of 57 km. The<br />

<strong>tunnel</strong> consists of two parallel s<strong>in</strong>gle track tubes<br />

each with a diameter of 9.2 m. The maximum<br />

overburden is 2350 meters. Three <strong>in</strong>termediate<br />

po<strong>in</strong>ts of attack divide the <strong>tunnel</strong> <strong>in</strong>to five sections<br />

of approximately equal length. These additional<br />

attacks are necessary to atta<strong>in</strong> a reasonable<br />

construction time and for ventilation purposes.<br />

Detailed <strong>in</strong>formation of the <strong>tunnel</strong> can be found<br />

under http://www.alptransit.ch/pages/e/.<br />

The <strong>in</strong>termediate po<strong>in</strong>ts of attack at Sedrun and<br />

Faido shall serve as Multi Function Sections (MFS)<br />

dur<strong>in</strong>g operation of the <strong>tunnel</strong>. These MFS enables<br />

the tra<strong>in</strong>s to change the <strong>tunnel</strong> tube <strong>in</strong> case of<br />

ma<strong>in</strong>tenance works and the specially ventilated<br />

419<br />

emergency sections serve the rescue of people <strong>in</strong><br />

emergency cases (Figure 1). The MFS Faido had to<br />

be accessed by a 2.7 km long access <strong>tunnel</strong><br />

decl<strong>in</strong><strong>in</strong>g at 12.7% from the portal. The overburden<br />

<strong>in</strong> the MFS Faido is approximately 1500 meters.<br />

The excavation was done by Drill and Blast<br />

1.2 Geological Prediction for the MFS FAIDO<br />

The outcrops from quarries <strong>in</strong> the area of the MFS<br />

Faido, the experiences made dur<strong>in</strong>g construction of<br />

the <strong>in</strong>vestigation system for the Triassic Piora bas<strong>in</strong><br />

as well as various vertical exploration drill<strong>in</strong>gs<br />

confirmed a favorable geological section for the<br />

designed location the MFS consist<strong>in</strong>g of Levent<strong>in</strong>a<br />

Gneisses with very good quality (Figure 2).


Cross cavern<br />

East Tube EON<br />

West Tube EWN<br />

Figure 1 Schematic orig<strong>in</strong>al layout of the multifunction section at Faido.<br />

0 km<br />

57 km<br />

420<br />

MFS Faido<br />

Figure 2. Geology along the <strong>Gotthard</strong> Base Tunnel consist<strong>in</strong>g ma<strong>in</strong>ly of gneiss and granite.<br />

MFS Faido is located <strong>in</strong> the Levent<strong>in</strong>a Gneiss of very good quality<br />

1.3 Actual Geology <strong>in</strong> the area of the MFS Faido<br />

Dur<strong>in</strong>g construction of the cross cavern (Figure 1) a<br />

break down of f<strong>in</strong>e gra<strong>in</strong>ed quartz occurred <strong>in</strong> the<br />

cavern’s roof form<strong>in</strong>g a cavity of 8 meters <strong>in</strong> height.<br />

The results of an extensive drill holes campaign<br />

together with seismic reflection measurements<br />

effectuated dur<strong>in</strong>g construction revealed an<br />

unknown large <strong>fault</strong> system <strong>in</strong> the area of the MFS<br />

Faido. The ma<strong>in</strong> kernel of this <strong>fault</strong> strikes at an<br />

average angle of about 15° to the <strong>tunnel</strong> axis and<br />

dips at about 80°. The <strong>fault</strong> system encountered<br />

dur<strong>in</strong>g construction is shown <strong>in</strong> Figure 3. In the<br />

<strong>fault</strong>’s kernel, layers of partially completely<br />

decomposed <strong>rock</strong> (kakitrite) are embedded.<br />

Adjacent to the east of the <strong>fault</strong> hard and brittle<br />

Levent<strong>in</strong>a gneiss is located. To the west of the <strong>fault</strong><br />

the <strong>rock</strong> mass consists of hard but less brittle<br />

Lucomagno gneiss.<br />

2. DIFFICULTIES DURING CONSTRUCTION<br />

Dur<strong>in</strong>g the first rebuild<strong>in</strong>g of the EWN tube to the<br />

north of the cross cavern (Figure 1) a support<br />

consist<strong>in</strong>g of steel arches HEM 200 backfilled with<br />

40 cm concrete was <strong>in</strong>stalled immediately after each<br />

excavation step of 1 m. This rigid support was<br />

<strong>in</strong>tended to cater for the heavy pressure and<br />

especially to protect workforce from break outs at<br />

the face. In the rebuilt section, on a length of 250 m<br />

(Figure 3) the load<strong>in</strong>g of the support however gave<br />

rise to radial displacements of up to 1 m. Stra<strong>in</strong>


North<br />

EON tube<br />

EWN tube<br />

Area of <strong>in</strong>creased<br />

<strong>rock</strong> <strong>burst</strong>s<br />

800 m<br />

Cross section for<br />

2D model<strong>in</strong>g<br />

421<br />

Fault’s kernel<br />

Cross Cavern<br />

Access <strong>tunnel</strong><br />

East<br />

West<br />

Rescue Station<br />

South<br />

Figure 3. Actual geology recognized by January 2003. Fault’s kernel with <strong>in</strong>tensively jo<strong>in</strong>ted <strong>rock</strong> to both sides of the<br />

kernel. As a consequence the crossover has been shifted to the south (Figure 1).<br />

measurements revealed yield<strong>in</strong>g of the steel arches<br />

already Figure 4 days after backfill<strong>in</strong>g[1]. Severe<br />

damage of the support developed (Figure 4) and the<br />

critical section had aga<strong>in</strong> to be rebuilt with an<br />

enlarged excavation radius of 1.5 m to allow for<br />

additional displacements. A flexible support has<br />

been <strong>in</strong>stalled us<strong>in</strong>g TH profiles with slid<strong>in</strong>g<br />

connections. Shotcrete slots were left open for<br />

unh<strong>in</strong>dered slid<strong>in</strong>g of the arches and for avoid<strong>in</strong>g<br />

damage of the shorcrete l<strong>in</strong><strong>in</strong>g due to further<br />

movements (Figure 5).<br />

Figure 4. Sheared off steel arch <strong>in</strong> the <strong>in</strong>vert<br />

3. NUMERICAL SIMULATION OF THE<br />

STATIC LOAD CASE<br />

In order to better understand the failure mechanisms<br />

<strong>in</strong> the area around the <strong>tunnel</strong>s of the MFS numerical<br />

model<strong>in</strong>g was carried out. The model section is<br />

Figure 5. Rebuild<strong>in</strong>g of the critical section with enlarged<br />

profile and slots <strong>in</strong> the concrete l<strong>in</strong><strong>in</strong>g<br />

<strong>in</strong>dicated <strong>in</strong> Figure 3. From the known geology it<br />

was obvious that a discont<strong>in</strong>uum model with<br />

deformable blocks had to be used <strong>in</strong> order to<br />

recognize the <strong>in</strong>fluence of the discont<strong>in</strong>uities and<br />

the kakirite <strong>zone</strong>s on the stability of the <strong>tunnel</strong><br />

section. For this purpose the program UDEC,<br />

developed by ITASCA was used. The <strong>in</strong>vestigation<br />

was <strong>base</strong>d on a parametric study compris<strong>in</strong>g the<br />

variation of <strong>rock</strong> and jo<strong>in</strong>t properties. The model is<br />

shown <strong>in</strong> Figure 6. The different <strong>fault</strong> regions <strong>in</strong> the<br />

model have been selected accord<strong>in</strong>g to the<br />

geologist’s f<strong>in</strong>d<strong>in</strong>gs <strong>base</strong>d on <strong>in</strong>terpretation of<br />

numerous borehole results. From overburden, <strong>rock</strong><br />

mass density and gravity a primary pr<strong>in</strong>cipal stress<br />

level of σyy=42 MPa and σxx = 21 MPa was<br />

def<strong>in</strong>ed at the top of the model. In Figure 6 the<br />

vertical stresses yy are represented. The results


clearly show a considerable extension of the stress<br />

redistribution due to the excavation of the <strong>tunnel</strong>s.<br />

In the area of the middle <strong>tunnel</strong> EWN both the σxx<br />

and σyy stresses decreased to approximately 10<br />

MPa. A general stress reduction occurred around all<br />

of the three <strong>tunnel</strong>s. To the east of the <strong>fault</strong> <strong>in</strong> the<br />

hard brittle gneiss a stress concentration is<br />

extend<strong>in</strong>g beneath the <strong>tunnel</strong> along the <strong>fault</strong>’s<br />

eastern boundary The stress concentration at <strong>tunnel</strong><br />

level amounts to yy = 60 to 80 MPa (Figure 6). The<br />

area of this stress is <strong>in</strong> accordance with the<br />

hypocenters of the micro tremors (Figure 8).<br />

Kernel Intensively<br />

jo<strong>in</strong>ted<br />

Hard <strong>rock</strong>, brittle<br />

Stress<br />

concentration<br />

West East<br />

Y Medium hard<br />

<strong>rock</strong><br />

X<br />

Figure 6. Static load case: Distribution of the<br />

vertical stresses σyy<br />

4. MICRO TREMORS<br />

4.1 Development of the seismic activity<br />

In this paper, we def<strong>in</strong>e micro tremors as a seismic<br />

event generally occurr<strong>in</strong>g at a larger distance from<br />

the <strong>tunnel</strong> (few hundred meters) whereas a <strong>rock</strong><br />

<strong>burst</strong> occurs <strong>in</strong> the direct vic<strong>in</strong>ity of the <strong>tunnel</strong>.<br />

Dur<strong>in</strong>g 2004 <strong>rock</strong> <strong>burst</strong>s took place ma<strong>in</strong>ly at the<br />

face of the EON (Figure 3) to the north of the cross<br />

cavern (Figure 1). Most of these <strong>burst</strong>s (ca. 75%)<br />

occurred approximately Figure 3 hours after<br />

blast<strong>in</strong>g. Dur<strong>in</strong>g 2005 <strong>in</strong>tensive <strong>rock</strong> <strong>burst</strong>s<br />

happened <strong>in</strong>dependently of the excavation <strong>in</strong> the<br />

back area of the eastern tube giv<strong>in</strong>g rise to partially<br />

severe damage of the support. Constructional<br />

adaptations of the support and exclud<strong>in</strong>g critical<br />

<strong>tunnel</strong> sections for access were required to provide<br />

for the safety of workforce and equipment [2].<br />

Between March 2004 and June 2005 the Swiss<br />

422<br />

Seismological Service (SED) recorded an<br />

accumulation of seismic activity <strong>in</strong> the area of the<br />

MFS Faido, a region which normally shows a very<br />

low seismicity. Dur<strong>in</strong>g the above mentioned period<br />

the permanently <strong>in</strong>stalled Swiss Digital Seismic<br />

Network (SDSNet) registered 10 seismic events<br />

with local magnitudes ML between 0.9 and 1.9.<br />

With the SDSNet located at the surface, the<br />

epicenters could be associated with the area of the<br />

MFS Faido with<strong>in</strong> an accuracy of one kilometer.<br />

Together with the M1.9 tremor of 1.6.2005 a strong<br />

<strong>rock</strong> <strong>burst</strong> could be associated. The same holds for<br />

two additional tremors of similar magnitude. On the<br />

other hand no relations to <strong>rock</strong> <strong>burst</strong>s could be<br />

identified for tremors dur<strong>in</strong>g the period March to<br />

April 2004. S<strong>in</strong>ce it could not be excluded that the<br />

accumulation of the tremors correlated with the<br />

construction of the MFS Faido the owner ATG<br />

decided <strong>in</strong> July 2005 to form a Work Group ’Micro<br />

tremors’ consist<strong>in</strong>g of representatives from the<br />

Swiss Seismological Service (SED), Geology,<br />

Eng<strong>in</strong>eer<strong>in</strong>g jo<strong>in</strong>t venture GBTS, Supervision and<br />

ATG, the owner of the <strong>tunnel</strong>. On March 25, 2006<br />

the strongest micro tremor of ML 2.4 was registered.<br />

This tremor was felt by the <strong>in</strong>habitants of the village<br />

Faido close to the jobsite. This tremor triggered a<br />

stress drop nearby the EON caus<strong>in</strong>g heav<strong>in</strong>g <strong>in</strong> the<br />

<strong>in</strong>vert.<br />

4.2 Additional seismic measur<strong>in</strong>g stations<br />

With first priority, additional seismic stations were<br />

<strong>in</strong>stalled at the surface and <strong>in</strong> the MFS Faido. For<br />

precise monitor<strong>in</strong>g and location of the seismic<br />

activity a special local seismic network consist<strong>in</strong>g<br />

of n<strong>in</strong>e stations at the surface, <strong>in</strong>clud<strong>in</strong>g one station<br />

from the SDSNet, were <strong>in</strong>stalled <strong>in</strong> a circular<br />

arrangement 10 to 15 km around the MFS Faido. In<br />

addition, two additional stations were <strong>in</strong>stalled at<br />

different locations <strong>in</strong> the <strong>tunnel</strong>s of the MFS Faido.<br />

The circular position of the seismic measur<strong>in</strong>g<br />

equipment allows a precise determ<strong>in</strong>ation of the<br />

epicenters whereas the measur<strong>in</strong>g stations directly<br />

above and <strong>in</strong>side the MFS Faido serve the<br />

evaluation of the depths of the micro tremors’<br />

sources. The read<strong>in</strong>gs of the measur<strong>in</strong>g stations<br />

were <strong>in</strong>tegrated <strong>in</strong> the SED’s data acquisition<br />

system. The real time transmission of the measur<strong>in</strong>g<br />

data guaranteed a cont<strong>in</strong>uous survey of the seismic<br />

activity allow<strong>in</strong>g for an immediate alert of the<br />

responsible organizations such as ATG, supervision<br />

and authorities <strong>in</strong> case of a strong tremor. This was


of particular importance <strong>in</strong> case of the M2.4 tremor<br />

occurr<strong>in</strong>g on 25.03.2006. Accurate seismic wave<br />

velocities required for the hypocenter’s<br />

determ<strong>in</strong>ations were derived from two calibration<br />

shots carried out <strong>in</strong> the MFS Faido. An average Pwave<br />

velocity of 5.33 km/s was calculated.<br />

4.3 Chronology of the seismic events<br />

Figure 7 shows the development of the recorded<br />

micro tremors’ number and magnitudes as functions<br />

of time. The highest seismic activity took place<br />

dur<strong>in</strong>g December 2005, March 2006 and Mai 2006.<br />

The highest magnitude of 2.4 occurred on 25 th<br />

March 2006. Dur<strong>in</strong>g the period from October 2005<br />

to February 2008 a total of 112 micro tremors were<br />

recorded. The magnitudes of most of the tremors<br />

were below 1.0. With term<strong>in</strong>ation of excavation <strong>in</strong><br />

the MFS Faido the micro tremors’ number and<br />

magnitudes decreased cont<strong>in</strong>uously. S<strong>in</strong>ce<br />

September 2007 no more micro tremors have been<br />

recorded above the measur<strong>in</strong>g threshold of M = -1.0<br />

<strong>in</strong> the area of the MFS Faido.<br />

Figure 7. Chronological development of number and magnitudes of the micro tremors (with courtesy form SED)<br />

423<br />

4.4 Epicenters<br />

The epicenters of all registered micro tremors<br />

dur<strong>in</strong>g Oct. 2005 to Feb. 2008 are depicted <strong>in</strong><br />

Figure 8. The micro tremors are concentrated <strong>in</strong> the<br />

<strong>rock</strong> mass to the north of the MFS Faido close to<br />

the eastern part of the <strong>tunnel</strong> system. The accuracy<br />

of the epicenters’ localization is less than 100 m and<br />

less than 250 m <strong>in</strong> focal depth as determ<strong>in</strong>ed by<br />

relocation of the calibration shots. With<strong>in</strong> the error<br />

ellipsoid the tremors’ sources are at <strong>tunnel</strong> level.<br />

4.5 Impacts of <strong>rock</strong> <strong>burst</strong>s on the <strong>tunnel</strong><br />

The <strong>rock</strong> <strong>burst</strong>’s damage potential is illustrated <strong>in</strong><br />

Figure.9. The picture to the left shows the damage<br />

of the support with a shotcrete plate ejected <strong>in</strong>to the<br />

EON. This <strong>rock</strong> <strong>burst</strong> occurred together with the<br />

M1.9 micro tremor of July 2005. The EON’s <strong>in</strong>vert<br />

heave <strong>in</strong> the picture to the right was caused by the<br />

M2.4 micro tremor of 25.03.2006. Due to the<br />

seismic impact the <strong>in</strong>vert heave was smaller than it<br />

is shown <strong>in</strong> the picture taken 2 days later.


4.6 Measured particle velocities V for the design<br />

micro tremor at station MFS-A<br />

In Figure 10 the measured particle velocities V<br />

registered at the station MFS-A (Figure 3 and<br />

Figure 8) for the M2.4 micro tremor on 25.03.2006<br />

are presented. In the diagram the velocity<br />

components Vx, Vy and Vz are shown. The<br />

orientations of the X-, Y-, and Z – coord<strong>in</strong>ates<br />

correspond to those of the model coord<strong>in</strong>ates. The<br />

maximum measured amplitude amount<strong>in</strong>g to 0.024<br />

m/s was <strong>in</strong> the horizontal X – direction normal to<br />

the <strong>tunnel</strong> axis. The total measur<strong>in</strong>g <strong>in</strong>terval was 1<br />

sec. The relevant vibration starts at 0.193 sec after<br />

the measurement’s trigger<strong>in</strong>g. The dynamic load for<br />

the model computations started at that same time to<br />

reduce computation time and amount of data.<br />

4.7 F<strong>in</strong>d<strong>in</strong>gs form the seismic measurements<br />

The micro tremors <strong>in</strong> the northern part of the MFS<br />

East tube, EON<br />

West tube, EWN<br />

M2.4<br />

Measur<strong>in</strong>g station MFS-B<br />

East<br />

Measur<strong>in</strong>g station MFS-A<br />

Figure 8. Epicenters of the micro tremors from Oct. 2005 to Feb. 2008 (with courtesy form SED)<br />

424<br />

tend to form clusters i.e. the sources of several<br />

tremors are located with<strong>in</strong> the same area.<br />

Consider<strong>in</strong>g the predom<strong>in</strong>ant steep west - east<br />

dipp<strong>in</strong>g jo<strong>in</strong>t system strik<strong>in</strong>g sub parallel to the<br />

<strong>tunnel</strong> axis (Figure 3 and Figure 6) shear failures<br />

along jo<strong>in</strong>ts are most likely. The micro tremor’s<br />

source locations <strong>in</strong> the hard Levent<strong>in</strong>a gneiss to the<br />

east of the <strong>fault</strong> corresponds to the location of the<br />

vertical stress concentrations result<strong>in</strong>g from the<br />

computations of the static load case <strong>in</strong> Figure 6.<br />

There is a general tendency of the micro tremors to<br />

move together with the excavation of the <strong>tunnel</strong>s<br />

form the cross cavern’s area to the north. Very few<br />

micro tremors occurred <strong>in</strong> the southern part of the<br />

MFS Faido; these events occurred outside the<br />

region shown <strong>in</strong> Figure 8. Many of the <strong>rock</strong> <strong>burst</strong>s<br />

caus<strong>in</strong>g the support’s damages <strong>in</strong> the <strong>tunnel</strong>s were<br />

triggered by the micro tremors.<br />

Date<br />

01.06.2005 25.03.2006/ Photo: 27.03.2006<br />

Figure 9. Left: Rock <strong>burst</strong> of 01.06.05: Damage of the support’s shotcrete, EON (on this date no seismic registration<br />

was <strong>in</strong>stalled). Right: M2.4 micro tremor of 25.03.2006: <strong>in</strong>vert – heave, EON.


5. DYNAMIC MODELLING<br />

5.1 Aims of modell<strong>in</strong>g<br />

The aim of the numerical model<strong>in</strong>g was to provide<br />

basic appraisals with respect to the structural safety<br />

and usability of the <strong>tunnel</strong>s dur<strong>in</strong>g and after a<br />

seismic event. The <strong>in</strong>vestigation’s result should<br />

furthermore disclose the residual risks to be<br />

accepted and its impact potential on the <strong>tunnel</strong>s<br />

concrete l<strong>in</strong><strong>in</strong>gs. The structural design has been<br />

completed prior to the occurrence of <strong>rock</strong> <strong>burst</strong>s.<br />

Therefore, it had to be controlled whether the l<strong>in</strong><strong>in</strong>g<br />

designed for the static load case still fulfills the<br />

<strong>tunnel</strong>s’ safety and usability requirements dur<strong>in</strong>g<br />

and after a micro tremor’s impact.<br />

5.2 Computational methods and models<br />

For the <strong>in</strong>vestigation 2D and 3D models were<br />

applied us<strong>in</strong>g the discont<strong>in</strong>uum codes UDEC and<br />

3DEC (ITASCA, [3]). Discont<strong>in</strong>uum models were<br />

used to consider block movements (rigid body<br />

movements) <strong>in</strong> the jo<strong>in</strong>ted <strong>rock</strong> mass. The blocks<br />

themselves are deformable. Due to the high<br />

computation time parametric studies were not<br />

possible with the 3D model. For parametric studies<br />

2D – models were used <strong>in</strong>stead. The results<br />

Particle Velocities V (m/sec)<br />

0.03<br />

0.02<br />

0.01<br />

0.00<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

-0.01<br />

-0.02<br />

-0.03<br />

Start of Computations<br />

Vy<br />

Vx<br />

Vz<br />

Time t (s)<br />

Figure 10.<br />

425<br />

presented refer solely to the dynamic load case. In a<br />

first step the static equilibrium of the supported<br />

<strong>tunnel</strong> system was computed. In a second step the<br />

f<strong>in</strong>al l<strong>in</strong><strong>in</strong>g was <strong>in</strong>serted <strong>in</strong> the model and the design<br />

micro tremor’s load superimposed to the static case.<br />

The properties for the <strong>rock</strong> and the two major jo<strong>in</strong>t<br />

systems were specified together with the geologists.<br />

A jo<strong>in</strong>t spac<strong>in</strong>g of 2 m was used for the <strong>fault</strong>’s<br />

kernel <strong>zone</strong>. The real jo<strong>in</strong>t spac<strong>in</strong>g <strong>in</strong> the kernel is<br />

much smaller which was accounted for by<br />

correspond<strong>in</strong>gly reduc<strong>in</strong>g the <strong>rock</strong>’s properties.<br />

Weak <strong>rock</strong> properties were assigned to the layers<br />

conta<strong>in</strong><strong>in</strong>g kakirites (Figure 11). For model<strong>in</strong>g<br />

support and f<strong>in</strong>al l<strong>in</strong><strong>in</strong>g <strong>in</strong> the 2D models block<br />

elements and structural bar elements were used.<br />

Special <strong>in</strong>vestigations have shown the deviation of<br />

seismic waves along weak <strong>rock</strong> formations.<br />

Therefore two different 2D models with one and<br />

with two kakirit layers (Figure 11), respectively,<br />

were <strong>in</strong>vestigated. The results from these models<br />

reveal the <strong>in</strong>fluence of a <strong>tunnel</strong>’s distance from a<br />

kakirite layer on the dynamic load<strong>in</strong>g of its l<strong>in</strong>er <strong>in</strong><br />

case of a micro tremor. With Model1 a stress drop<br />

<strong>in</strong> the vic<strong>in</strong>ity of the EON was <strong>in</strong>itiated. The results<br />

reveal the impact of such an event on the <strong>tunnel</strong><br />

l<strong>in</strong><strong>in</strong>g.<br />

Z<br />

Y<br />

X


5.3 Specification of the model design load<strong>in</strong>g wave<br />

The work group ‘Micro Tremors’ decided to<br />

consider the M2.4 micro tremor of 25 March 2006<br />

as the decisive design – tremor. For the seismic<br />

model load<strong>in</strong>g of the <strong>in</strong>put wave correspond<strong>in</strong>g to<br />

the design tremor had first to be specified. With this<br />

decision a conservative dynamic load was selected<br />

<strong>in</strong> order to control the safety and stability of the<br />

<strong>tunnel</strong>s under operation <strong>in</strong> case of a micro tremor <strong>in</strong><br />

the area of the <strong>tunnel</strong>s. The emitted wave at the<br />

tremor’s source first had to pass the <strong>tunnel</strong> system<br />

and the <strong>fault</strong> <strong>zone</strong> prior to arrive at measur<strong>in</strong>g<br />

station MFS-A. To determ<strong>in</strong>e the waves attenuation<br />

10 times amplified measured particle velocities<br />

Figure 11. 2D – Model1 and Model 2 with 2 and 1 layers of kakirite <strong>in</strong> the <strong>fault</strong>’s kernel<br />

Figure 12. Left: 3D Model. Right: Vx Horizontal cut through the model with po<strong>in</strong>t P41 simulat<strong>in</strong>g approximately<br />

the MFS-A. The curves refer to the Vx <strong>in</strong>put particle velocity and the Vx at P41. ∆t is the<br />

wave’s travel time from the model’s load boundary to Po<strong>in</strong>t P41.<br />

426<br />

were applied at the eastern boundary of the 3D<br />

model. With this load the particle velocities <strong>in</strong> a<br />

model po<strong>in</strong>t simulat<strong>in</strong>g the MFS-A were <strong>in</strong><br />

magnitude similar to those of the measured signal<br />

(Figure 12). The dynamic design load was therefore<br />

determ<strong>in</strong>ed as the 10 times amplified measured<br />

signal (Figure 10) with a peak Vx -velocity<br />

component of 0.23 m/s. This design load covers the<br />

strongest micro tremor identified s<strong>in</strong>ce start<strong>in</strong>g the<br />

extended seismic monitor<strong>in</strong>g <strong>in</strong> the MFS Faido<br />

(Figure 7) and is conservative. Figure 12 shows the<br />

3D model used. In the diagram the design load’s Vx<br />

particle velocity is compared to the signal received<br />

<strong>in</strong> model po<strong>in</strong>t P41 simulat<strong>in</strong>g the location of the<br />

MFS-A.


5.4 Results form 2D Modell<strong>in</strong>g<br />

The results presented refer to the dynamic design<br />

load only. Of particular <strong>in</strong>terest were the l<strong>in</strong>ers’<br />

concrete stresses and displacements as a function of<br />

time due to the dynamic load. For this purpose the<br />

correspond<strong>in</strong>g results <strong>in</strong> the EWN’s l<strong>in</strong>er were<br />

compared for Model_1 and Model_2. For the<br />

<strong>in</strong>terpretation a l<strong>in</strong>er segment at the eastern change<br />

from <strong>in</strong>vert to side wall was selected. The concrete<br />

l<strong>in</strong><strong>in</strong>g thickness of the EWN is 60 cm. Negative<br />

stress denotes compression. In Model 1 the EWN is<br />

embedded between two kakirite layers whereas <strong>in</strong><br />

Model 2 the location of the EWN is <strong>in</strong> front of a<br />

kakirite layer and exposed to the more or less<br />

undamped wave arriv<strong>in</strong>g from the right model side<br />

(East). The concrete stresses represented <strong>in</strong> Figure<br />

13 and Figure 14 clearly show the damp<strong>in</strong>g effect of<br />

the kakirite layer <strong>in</strong> front of the EWN (Figure 13).<br />

The concrete stresses due to the dynamic<br />

design load are small compared to those determ<strong>in</strong>ed<br />

for the static case. The stress peaks occur at high<br />

frequencies. The evaluation of the displacement’s<br />

corner frequencies revealed a displacement’s<br />

decrease approximately at 60 Hz. Figure 13 and<br />

Figure 14 conta<strong>in</strong> the low pass filtered stresses at 60<br />

Hz. From the curve of the filtered stresses can be<br />

deduced that the design stresses to be considered are<br />

considerably smaller than the stress peaks. Due to<br />

an irreversible shear movement along a jo<strong>in</strong>t <strong>in</strong> the<br />

<strong>tunnel</strong>’s vic<strong>in</strong>ity the l<strong>in</strong>er stresses for Model 2<br />

reveal a small permanent rema<strong>in</strong><strong>in</strong>g compression<br />

stress after the dynamic impact. The shear<br />

movement <strong>in</strong> the jo<strong>in</strong>t was caused by the<br />

unh<strong>in</strong>dered wave’s deviation towards the <strong>tunnel</strong><br />

along the weak kakirite layer. Figure 15 and Figure<br />

16 show the EWN’s displacements evaluated for the<br />

same case and at the same l<strong>in</strong>er location as used for<br />

the concrete stresses’ evaluation. In accordance<br />

with the stresses the deformations generally are<br />

small. Similar to the concrete stresses, the<br />

displacement of the EWN’s l<strong>in</strong>er embedded<br />

between two kakirite layers is smaller compared to<br />

the case of Model 2.<br />

427<br />

Concrete stress <strong>in</strong> l<strong>in</strong>er (MPa)<br />

Concrete stress <strong>in</strong> l<strong>in</strong>er (MPa)<br />

L<strong>in</strong>er displacements d (mm)<br />

2<br />

1<br />

0<br />

-1<br />

0.0 0.1 0.2 0.3 0.4 0.5<br />

-2<br />

-3<br />

-4<br />

-5<br />

Time t (s)<br />

EWN<br />

Kakirite <strong>zone</strong>s<br />

unfiltered<br />

60 Hz LP filter<br />

Figure 13 Concrete stresses (M+N) <strong>in</strong> l<strong>in</strong>er of <strong>tunnel</strong><br />

EWN for Mod_1 with 2 kakirite <strong>zone</strong>s<br />

2<br />

1<br />

0<br />

0.0<br />

-1<br />

-2<br />

-3<br />

-4<br />

-5<br />

0.1 0.2 0.3 0.4 0.5<br />

Time t (s)<br />

EWN<br />

Kakirite <strong>zone</strong><br />

unfiltered<br />

60 Hz LP filter<br />

Figure 14 Concrete stresses (M+N) <strong>in</strong> l<strong>in</strong>er of <strong>tunnel</strong><br />

EWN for Mod_2 with 1 kakirite <strong>zone</strong><br />

2<br />

1<br />

0<br />

0.0<br />

-1<br />

0.1 0.2 0.3 0.4 0.5<br />

-2<br />

-3<br />

-4<br />

Time t (s)<br />

Y<br />

dy<br />

dx<br />

X<br />

EWN<br />

Kakirite <strong>zone</strong>s<br />

Figure 15 L<strong>in</strong>er displacement of <strong>tunnel</strong> EWN for Mod_1<br />

with 2 kakirite <strong>zone</strong>s


Dur<strong>in</strong>g the dynamic computation with Model 1<br />

a ‘<strong>rock</strong> <strong>burst</strong>’ was triggered beneath the EON<br />

<strong>tunnel</strong>. This <strong>burst</strong> occurred due to a sudden stress<br />

drop <strong>in</strong> a jo<strong>in</strong>t. Figure 17 shows the stress drop’s<br />

L<strong>in</strong>er displacement d (mm)<br />

2<br />

1<br />

0<br />

0.0<br />

-1<br />

0.1 0.2 0.3 0.4 0.5<br />

-2<br />

-3<br />

-4<br />

Time t (s)<br />

Y<br />

EWN<br />

Kakirite <strong>zone</strong><br />

dy<br />

dx<br />

Figure 16 L<strong>in</strong>er displacement of <strong>tunnel</strong> EWN for Mod_2<br />

with 1 kakirite <strong>zone</strong><br />

development at different times represented by<br />

particle velocities’ Vx components. The disturbance<br />

X<br />

Figure Figure 2 Development 17 Development of a stress of a stress drop <strong>in</strong> drop the <strong>in</strong> vic<strong>in</strong>ity the vic<strong>in</strong>ity of <strong>tunnel</strong> of <strong>tunnel</strong> EON. EON. Represented Represented are the are particle the particle velocities velocities Vx. Vx –<br />

Vx. Vx Range: – Range: -0.2 ÷ -0.2 +0.2 ÷ +0.2 m/s. The m/s. represented The represented Vx - Vx <strong>in</strong>terval - <strong>in</strong>terval range range is 0.05 is 0.05 m/s. m/s.<br />

Figure Figure 2 Left: 18 Left: Development Development of Vx of around Vx around the stress the stress drop drop area. area. The The represented represented Vx - Vx <strong>in</strong>terval - <strong>in</strong>terval range range is 0.10 is 0.10 m/s. m/s. Right:<br />

Stress drop Right: of 17 Stress MPa drop <strong>in</strong> the of jo<strong>in</strong>t 17 MPa <strong>in</strong>dicated <strong>in</strong> the by jo<strong>in</strong>t the <strong>in</strong>dicated shear movement by the shear <strong>in</strong> the movement left hand picture. <strong>in</strong> the left The hand residual picture. shear stress<br />

The residual amounts shear stress to 14.5 amounts MPa. to 14.5 MPa.<br />

428<br />

trigger<strong>in</strong>g the drop was small amount<strong>in</strong>g to<br />

0 < Vx < 0.05 m/s and occurred prior to the arrival<br />

of the <strong>in</strong>put wave’s first Vx - peak. The drop’s shear<br />

direction and the magnitude of the drop’s stress<br />

release of 17 MPa are shown <strong>in</strong> Figure 18 Due to<br />

the stress drop, Vx <strong>in</strong> the <strong>tunnel</strong> area amounts to<br />

more than 0.40 m/s (Figure 18). The disturb<strong>in</strong>g<br />

<strong>in</strong>fluence of the stress drop on the <strong>in</strong>put wave is<br />

depicted <strong>in</strong> Figure 17, at time t = 0.1532 s i.e.<br />

0.0179 s after the stress drop’s trigger<strong>in</strong>g. From<br />

Figure 17 the above mentioned deviation of the<br />

<strong>in</strong>put’s Vx - peak along the eastern kakirite layer is<br />

obvious at t = 0.1532 s. The stress drop occurred<br />

only <strong>in</strong> Model 1 with two kakirite layers.<br />

The impact potential of the stress drop on the<br />

l<strong>in</strong><strong>in</strong>g was <strong>in</strong>vestigated us<strong>in</strong>g the particle velocities<br />

PV accord<strong>in</strong>g to 1) with the components aris<strong>in</strong>g at<br />

the same time t.<br />

...(1)<br />

In Figure 19 the PV’s for a po<strong>in</strong>t <strong>in</strong> the <strong>rock</strong><br />

near the l<strong>in</strong><strong>in</strong>g of the EON are plotted. Compared


are the PV’s for the cases with and without stress<br />

drop i.e. with Model 2 and Model 1. In case of the<br />

stress drop the peak PV <strong>in</strong> the <strong>rock</strong> near the l<strong>in</strong><strong>in</strong>g<br />

amounts to almost 900 mm/s. The dynamic load<strong>in</strong>g<br />

on the l<strong>in</strong><strong>in</strong>g <strong>in</strong> this case is similar to a shock impact<br />

[4], [5]. Without stress drop the correspond<strong>in</strong>g<br />

maximum PV amounts to 200 mm/s. The frequency<br />

of the stress drop <strong>in</strong>duced PV is bigger than 300 Hz<br />

whereas the frequency for the peak PV without<br />

stress drop is <strong>in</strong> the range of 50 to 60 Hz. The PV’s<br />

of the l<strong>in</strong><strong>in</strong>g segment under consideration are shown<br />

<strong>in</strong> Figure 20 for the cases with and without stress<br />

drop. In the models the seal<strong>in</strong>g membrane between<br />

<strong>rock</strong> and l<strong>in</strong><strong>in</strong>g was simulated by means of an<br />

<strong>in</strong>terface. The transmissibility of this <strong>in</strong>terface<br />

contributes to the vibration behavior of the l<strong>in</strong><strong>in</strong>g<br />

[6]. For the case with stress drop the damp<strong>in</strong>g of the<br />

l<strong>in</strong><strong>in</strong>g’s peak PV is obvious. However, due to the<br />

damp<strong>in</strong>g’s frequency reduction <strong>in</strong> both cases with<br />

and without stress drop a longer duration of high<br />

PV <strong>in</strong> the l<strong>in</strong><strong>in</strong>g can be derived from Figure 20.<br />

429<br />

Figure 21 and Figure 22 show the stresses <strong>in</strong><br />

the concrete l<strong>in</strong>er and the l<strong>in</strong>er’s dx – displacement<br />

respectively, both for the cases with and without<br />

stress drop. In the case of the stress drop the<br />

concrete stress amounts to more than 8 MPa and the<br />

l<strong>in</strong>er’s displacements to 6 mm. From the above<br />

mentioned results must be concluded, that <strong>in</strong> case of<br />

a stress drop as developed <strong>in</strong> Model 1 damage of the<br />

l<strong>in</strong><strong>in</strong>g cannot be excluded. The probability of the<br />

occurrence of the design micro tremor and the<br />

correspond<strong>in</strong>g trigger<strong>in</strong>g of a stress drop has been<br />

assessed by experts as very low. This assessment is<br />

<strong>base</strong>d on the fact that after term<strong>in</strong>ation of the<br />

excavation <strong>in</strong> the MFS no more micro tremors<br />

occurred and <strong>in</strong> addition a clear correlation between<br />

excavation activity and micro tremors was<br />

ascerta<strong>in</strong>ed. Furthermore the seismic activity <strong>in</strong> the<br />

area of the MFS Faido is very low. On 21.01.2008 a<br />

M4.0 earthquake occurred at a distance of<br />

approximately 50 km from the MFS Faido. This<br />

earthquake was recognized by the measur<strong>in</strong>g<br />

Figure 19 Particle velocities PV of a po<strong>in</strong>t <strong>in</strong> the <strong>rock</strong> Figure 20. Particle velocities PV of the l<strong>in</strong><strong>in</strong>g segment.<br />

mass near the l<strong>in</strong><strong>in</strong>g. With and without stress drop With and without stress drop<br />

Figure 21 L<strong>in</strong><strong>in</strong>g’s concrete stresses (M+N) of <strong>tunnel</strong> EON. Figure 22 L<strong>in</strong> <strong>in</strong>g displacements dx of <strong>tunnel</strong> EON.<br />

With and without stress drom With and without stress drop


stations <strong>in</strong> the MFS. No trigger<strong>in</strong>g of a micro tremor<br />

could be identified. Based on above mentioned<br />

aspects the case of a stress drop near a <strong>tunnel</strong> <strong>in</strong> the<br />

MFS Faido has been accepted as a residual risk.<br />

S<strong>in</strong>ce a residual risk cannot be completely excluded,<br />

the owner of the <strong>tunnel</strong>, ATG decided to <strong>in</strong>stall<br />

seismic measur<strong>in</strong>g equipments and vibration sensors<br />

<strong>in</strong> the l<strong>in</strong><strong>in</strong>gs of the <strong>tunnel</strong>s <strong>in</strong> the MFS Faido for a<br />

permanent seismic monitor<strong>in</strong>g dur<strong>in</strong>g operation.<br />

6. CONCLUSIONS<br />

Dur<strong>in</strong>g <strong>tunnel</strong><strong>in</strong>g at great depth <strong>in</strong> geological<br />

conditions as encountered <strong>in</strong> the MFS Faido micro<br />

tremors trigger<strong>in</strong>g <strong>rock</strong> <strong>burst</strong>s are likely to occur.<br />

The stress redistribution due to <strong>tunnel</strong><strong>in</strong>g and<br />

the stress concentration <strong>in</strong> hard <strong>rock</strong> <strong>in</strong> comb<strong>in</strong>ation<br />

with an exist<strong>in</strong>g <strong>zone</strong> of weakness (<strong>fault</strong>) favors the<br />

occurrence of micro tremors.<br />

The micro tremors did clearly correlate with the<br />

excavation activities. After term<strong>in</strong>at<strong>in</strong>g the excavation<br />

no more micro tremors have been identified.<br />

Micro tremors cannot be avoided. Dur<strong>in</strong>g<br />

construction precaution measures such as clos<strong>in</strong>g of<br />

critical sections and flexible support consist<strong>in</strong>g of<br />

flexible <strong>rock</strong> bolts and steel arches are to be applied.<br />

A seismic wave is deviated by a weak <strong>zone</strong>.<br />

With the orientation of the weak <strong>zone</strong>s <strong>in</strong> the MFS<br />

Faido a seismic wave is deviated towards the <strong>tunnel</strong>s.<br />

The dynamic impact on a l<strong>in</strong><strong>in</strong>g of a <strong>tunnel</strong> <strong>in</strong> front of<br />

a weak <strong>zone</strong> and exposed to a more or less unh<strong>in</strong>dered<br />

micro tremor’s wave is considerably higher compared to the<br />

impact on a l<strong>in</strong>er of a <strong>tunnel</strong> <strong>in</strong> the ‘shelter’ of a weak <strong>zone</strong>.<br />

The specification of the design micro tremor<br />

for the MFS Faido is conservative.<br />

Exclud<strong>in</strong>g the additional load<strong>in</strong>g due to a<br />

spontaneous ‘stress drop’ <strong>in</strong> the direct vic<strong>in</strong>ity of a<br />

<strong>tunnel</strong> there has been no need for improv<strong>in</strong>g<br />

(thickness, additional re<strong>in</strong>forc<strong>in</strong>g) the l<strong>in</strong><strong>in</strong>gs<br />

designed for the static load case <strong>in</strong> the MFS Faido.<br />

ACKNOWLEDGEMENT<br />

The authors thank the AlpTransit <strong>Gotthard</strong> for the<br />

permission to publish the paper.<br />

REFERENCES<br />

1. Hagedorn, H., Rehbock-Sander, M., Flury,<br />

S.,<strong>Gotthard</strong> Base Tunnel: State of the works and<br />

Special Aspects, Proc. Int. Symposium on<br />

430<br />

2.<br />

Construction and Operation of Long Tunnels,<br />

Taipei, Taiwan, 7. - 10.11.2005<br />

Rehbock-Sander, M., Rock Bursts Experience<br />

ga<strong>in</strong>ed <strong>in</strong> M<strong>in</strong>es and Deep Tunnels, Proc. 4 th Asian<br />

Rock Mechanics Symposium (ARMS 2006),<br />

S<strong>in</strong>gapore, 8.-10.11.2006<br />

3. ITASCA, UDEC – Manual, Theory and<br />

Background, Itasca Consult<strong>in</strong>g Group, Inc.<br />

M<strong>in</strong>neapolis, M<strong>in</strong>nesota, USA<br />

4. Kaiser, P.K., Vasak,P., Suor<strong>in</strong>emi, F.T. and<br />

Thibodeau, D. (2005). New dimensions <strong>in</strong> seismic<br />

data <strong>in</strong>terpretation with 3-D virtual reality<br />

5.<br />

visualization <strong>in</strong> <strong>burst</strong>-prone ground, RaSiM6, Perth,<br />

Australia, 33 – 47<br />

Heuze, F.E., Morris, J.P., Insights <strong>in</strong>to ground shock <strong>in</strong><br />

jo<strong>in</strong>ted <strong>rock</strong>s and the response of structures there-<strong>in</strong>, Int.<br />

Journal of Rock Mechanics & M<strong>in</strong><strong>in</strong>g Sciences 44 (2007)<br />

647-676<br />

6. Harris, C.M., Crede, E.C., Shock and Vibration<br />

Handbook, McGraw-Hill, Book Company, ISBN 0-<br />

07-026799-5<br />

BIOGRAPHICAL DETAILS OF THE AUTHORS<br />

He<strong>in</strong>z Hagedorn completed his study<br />

<strong>in</strong> civil eng<strong>in</strong>eer<strong>in</strong>g at the Swiss<br />

Federal Institute of Technology (ETH)<br />

with a M.Sc. degree <strong>in</strong> 1969. From<br />

1970 to 1979 he was work<strong>in</strong>g at ETH<br />

with the <strong>rock</strong> mechanics <strong>in</strong>stitute of<br />

Prof. K. Kovari. The activity consisted<br />

<strong>in</strong> develop<strong>in</strong>g a F<strong>in</strong>ite Element<br />

program designed for underground<br />

structures and research activity <strong>in</strong> material test<strong>in</strong>g and<br />

field measurement. S<strong>in</strong>ce 1979 he has been work<strong>in</strong>g with<br />

Amberg Eng<strong>in</strong>eer<strong>in</strong>g Ltd., Switzerland as a geotechnical<br />

expert. For numerous large underground constructions<br />

such as <strong>tunnel</strong>s, caverns etc. <strong>in</strong> Switzerland and abroad he<br />

provided geotechnical expertise and consult<strong>in</strong>g for design<br />

and dur<strong>in</strong>g construction.<br />

Rolf Stadelmann obta<strong>in</strong>ed his M. Sc <strong>in</strong><br />

Civil Eng<strong>in</strong>eer<strong>in</strong>g from the ETH<br />

Zurich, Switzerland <strong>in</strong> 1992. He<br />

specializes <strong>in</strong> general stability<br />

analyses, temporary & permanent<br />

excavation support systems, as well as<br />

<strong>in</strong> structural design of all types of<br />

<strong>tunnel</strong> l<strong>in</strong><strong>in</strong>g. His core competence is<br />

also <strong>in</strong> the quality management for<br />

geotechnical analyses as well as <strong>in</strong> the development of<br />

Contractor’s alternatives. He has <strong>in</strong>-depth experiences <strong>in</strong><br />

underground space design for all k<strong>in</strong>ds of <strong>rock</strong> conditions.<br />

S<strong>in</strong>ce 2003 he is Head of Geotechnical Department of<br />

Amberg Eng<strong>in</strong>eer<strong>in</strong>g Ltd <strong>in</strong> Switzerland.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!