CF Villegas Ls

UNIVERSIDAD DE CHILE
FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS

DEPARTAMENTO DE GEOFSICA
ESTRUCTURA SISMICA CORTICAL EN LOS ANDES CENTRALES

(33-34.5S): CONCENTRACIONES DE SISMICIDAD BAJO MINAS
EL TENIENTE Y DISPUTADA.
TESIS PARA OPTAR AL TTULO DE

MAGISTER EN CIENCIAS MENCIN GEOFSICA
LUIS PEDRO VILLEGAS SEGURA
PROFESOR GUA:
EMILIO VERA SOMMER
MIEMBROS DE LA COMISIN:
MARIO PARDO PEDEMONTE
EDUARDO CONTRERAS REYES
ANDRS TASSARA ODDO
SANTIAGO DE CHILE
AGOSTO 2012
RESUMEN
Los procesos sismotectnicos de la Zona Central de Chile estn relacionados con

la subduccin de la Placa de Nazca bajo la Placa Sudamericana. Sin embargo, la
sismicidad cortical superficial en la zona central ha sido poco estudiada, tanto en su
origen como en su tasa de actividad, a pesar de la existencia de 2 eventos ssmicos de
gran magnitud localizados a profundidades menores de 15 km, uno en Las Melosas en
1958 (M=6.9) y otro al oeste de Rancagua en 1987 (M=5.9). Cabe mencionar que en
este sector se concentra el 40% de la poblacin del pas, adems se encuentran 3
grandes minas cuprferas, algunas minas de yeso, plantas hidroelctricas y el
gasoducto proveniente de Argentina.
En este escenario y enmarcado en el proyecto FONDECYT 1050758, este trabajo
tiene el objetivo de crear un modelo de estructura de velocidades a nivel cortical de la
zona central cordillerana (33-34.5S), adems de estudiar las concentraciones de
sismicidad bajo las minas Los Bronces (Ex Disputada) y El Teniente. Esto ayudar a
una mejor caracterizacin y comprensin del comportamiento ssmico en esta zona.
Con 28 estaciones temporales ms las estaciones permanentes de la Universidad
de Chile se registraron, entre Noviembre del 2005 y Marzo del 2006, una serie de
eventos, siendo de gran importancia los registros provenientes de tronaduras asociadas
a las faenas de extraccin de la mina Los Bronces. Estos eventos artificiales permiten
la obtencin de datos distancia-tiempo para ondas ssmicas, donde se tiene un
conocimiento prcticamente exacto de las ubicaciones tanto de la fuente como de las
estaciones registradoras. En el caso de El Teniente los eventos son netamente de
carcter ssmico.
El modelo 1D obtenido en Los Bronces se resume en 3 capas con gradiente lineal
de velocidad. Aumentando en profundidad, los espesores de estas capas son de 4.1,
22.1 y 3.9 km y sus velocidades de 4.8 - 5.1 km/s, 5.9 - 7.0 km/s y 7.1- 7.2 km/s para
ondas P y de 2.8 km/s, 3.5 - 3.9 km/s y 3.9 - 4 km/s para ondas S. En el caso de El
Teniente el modelo es de 2 capas con espesores de 2.4 y 23.9 km y velocidades de 4.4
- 4.5 km/s y 5.8 - 7 km/s para ondas P y 2.5 - 2.6 km/s y 3.4 - 3.9 km/s para ondas S.
Ambos modelos se asemejan a modelos realizados en estudios anteriores estudios y
las diferencias estn relacionadas principalmente a la correccin de altura y a la
distribucin de las estaciones en el rea de estudio del nuevo modelo. Se debe
mencionar que en la modelacin de El Teniente, la dispersin de datos distancia-tiempo
(X-T) es alta producto de la decisin de tomar 1 punto de origen cuando en realidad
existen distintos hipocentros debido al carcter ssmico de los eventos.
Con los modelos unidimensionales se realiz un modelo bidimensional entre estas
dos zonas mineras, creando un perfil Norte-Sur de 132.8 Km de longitud ubicado en el
sector precordillerano. Para ajustar el modelo se utilizaron los datos de estaciones
cercanas al perfil y se consideraron slo los datos distancia-tiempo de los eventos
provenientes de Los Bronces. El modelo 2D obtenido alcanza profundidades cercanas a
los 25 Km y al observar el grfico Vp/Vs se puede concluir un mayor aumento en
profundidad de los valores de Vp/Vs hacia el Norte del perfil, lo que podra indicar la
ii
existencia de fases fludas (roca hidratada o magma) y roca ms seca y rgida hacia el
Sur del perfil.
Se realiz adems la relocalizacin de los eventos concentrados en las zonas
mineras de Los Bronces y El Teniente con los nuevos modelos 1D y 2D. Las
variaciones de la relocalizacin son poco significativas y no se distingue caracterstica
geomtrica que pueda ayudar a una mejor interpretacin de la sismicidad en estos
sectores. En cuanto a los esfuerzos, en Los Bronces se lograron obtener mecanismos
focales de 3 eventos, todos con fallamiento inverso de rumbo Norte-Sur. En El Teniente
se obtuvieron 17 mecanismos focales, la mayora compresionales y tambin con
rumbos aproximados Norte-Sur.
iii
A Salom Muoz
y Berta Soto
iv
AGRADECIMIENTOS
En primer lugar agradezco a mis padres, quienes priorizaron durante mi infancia

satisfacer mis curiosidades por sobre caprichos fugaces. Los logros de mi educacin
son en gran medida vuestro mrito.
A mi amada Pachi, quien comparte mis esperanzas, rabias y formas de moverme
por la vida. Gracias por seguir confiando en el arte y la ciencia como herramientas para
la construccin de una sociedad ms justa.
A mis hermanos Juan y Sebastin por su amistad.
A la gran cantidad de personas increbles que conoc durante mis aos de
Universidad, a mis amigos de Radio Integral, a mis amigos del Grupo de Cine de
Ingeniera y a mis eternos amigos de la revista El Pjaro Verde quienes me levantan
el nimo con challa, pebre y poesa.
A mis compaeros de sala en Licenciatura y de oficina en el Magister, con
quienes forj una gran amistad y un sincero respeto: Felipe, Andrei, Marianela,
Eduardo, Cindy, Hernn, Rodrigo y Francisco.
A Ernesto Castillo por su paciencia y buena disposicin en la implementacin y
mejora del programa computacional que utilic en este trabajo.
A los funcionarios del DGF por su tiempo, buena voluntad y preocupacin, en
especial a Pato, Rodrigo y Germn.
Y finalmente a muchos de los profesores del DGF, en particular a Emilio Vera y
Mario Pardo por su particular ayuda en este trabajo.
Tabla de contenidos
Capitulo 1
Introduccin .............................................................................................1
1.1
Antecedentes .....................................................................................................1
1.2
Mina Los Bronces .............................................................................................2
1.3
Mina El Teniente ................................................................................................3
1.4
Objetivos ............................................................................................................4
Capitulo 2
2.1
Marco terico ...........................................................................................5
Teora de rayos ..................................................................................................5
Capitulo 3
Datos ........................................................................................................8
3.1
Red Sismolgica ...............................................................................................8
3.2
Los Bronces.....................................................................................................10
3.3
El Teniente .......................................................................................................11
Capitulo 4
4.1
Modelos 1D ......................................................................................................12
4.1.1
4.1.2
4.1.3
4.1.4
4.1.5
4.2
Procesamiento de datos .......................................................................12
Modelo de velocidades en Los Bronces .....................................................12

Modelo de velocidades en El Teniente .......................................................17
Correccin por altura ..................................................................................22
Correccin de altura en Los Bronces..........................................................23
Correccin de altura en El Teniente ...........................................................25
Modelacin 2D .................................................................................................27
4.2.1
4.2.2
4.2.3
4.2.4
Eleccin del Perfil .......................................................................................28

Seleccin de estaciones y correccin de altura ..........................................30
Modelos 1D.................................................................................................32
Modelacin 2D Los Bronces - El Teniente ..................................................45
4.3
Relocalizacin .................................................................................................55
4.4
Mecanismos focales .......................................................................................58
Capitulo 5
5.1
Resultados .............................................................................................59
Modelo de Velocidades 1-D Los Bronces .....................................................59

vi
5.2
Modelo de Velocidades 1-D El Teniente ........................................................60
5.3
Modelo de Velocidades 2-D (Perfiles) ............................................................61
5.4
Relocalizacin .................................................................................................65
5.5
Mecanismos focales .......................................................................................72
Capitulo 6
Interpretacin y Discusin ...................................................................74
6.1
Modelos 1D ......................................................................................................74
6.2
Modelo 2D ........................................................................................................76
6.3
Sismicidad en Los Bronces y El Teniente.....................................................80
Capitulo 7
Conclusiones .........................................................................................81
Bibliografa 83
Anexo A
86
Relacin grfica entre Vm y X y clculo del valor Vm ..............................................86

Anexo B
89
Rutina en Matlab para la creacin de la figuras de la modelacin 2D. ....................89
vii
ndice de Figuras
Fig. 1.1 Zona de Falla El Teniente (ZFT). ..................................................................................... 3
Fig. 2.1 Geometra del rayo atravesando capa horizontal. ............................................................ 5
Fig. 2.2 Gradiente de velocidad en una capa. ............................................................................... 6
Fig. 3.1 Red sismolgica utilizada. Perfil W-E. .............................................................................. 8
Fig. 3.2 Zona de estudio (33 - 34.5S). ........................................................................................ 9
Fig. 3.3 Eventos seleccionados del sector de Los Bronces (33.1S - 33.2S y 70.3W - 70.2W).10
Fig. 3.4 Eventos seleccionados del sector de El Teniente (34.04S - 34.14S y 70.4W 70.3W). ...................................................................................................................................... 11
Fig. 4.1 Vista en planta de los eventos registrados en Los Bronces. ........................................... 13
Fig. 4.2 Eventos registrados en Los Bronces. ............................................................................. 13
Fig. 4.3 Clculo del tiempo mediante regresin lineal. ................................................................ 14
Fig. 4.4 Modelo de velocidades en Los Bronces sin correcciones (Onda P). .............................. 16
Fig. 4.5 Modelo de velocidades en Los Bronces sin correcciones (Onda S). .............................. 16
Fig. 4.6 Vista en planta de los eventos registrados en El Teniente. ............................................ 18
Fig. 4.7 Eventos registrados en El Teniente. ............................................................................... 18
Fig. 4.8 Clculo del tiempo mediante regresin lineal. ................................................................ 19
Fig. 4.9 Modelo de velocidades en El Teniente sin correcciones (Onda P). ................................ 21
Fig. 4.10 Modelo de velocidades en El Teniente sin correcciones (Onda S). .............................. 21
Fig. 4.11 Diagrama sobre la correccin de altura. ....................................................................... 23
Fig. 4.12 Modelo de velocidades en Los Bronces con correccin de altura (Onda P). ............... 24
Fig. 4.13 Modelo de velocidades en Los Bronces con correccin de altura (Onda S). ................ 24
Fig. 4.14 Modelo de velocidades de El Teniente con correccin de altura (Onda P). .................. 26
Fig. 4.15 Modelo de velocidades de El Teniente con correccin de altura (Onda S). .................. 26
Fig. 4.16 Celda Trapezoidal. ....................................................................................................... 28
Fig. 4.17 Perfil entre Los Bronces y El Teniente. ........................................................................ 29
Fig. 4.18 Perfil Norte-Sur. ........................................................................................................... 29
Fig. 4.19 Seleccin de estaciones cercanas al perfil. .................................................................. 31
Fig. 4.20 Perfil Norte-Sur y estaciones seleccionadas. ............................................................... 31
Fig. 4.21 Datos Distancia-Tiempo () y simulacin de disparo en Los Bronces () (onda P). ....... 33
Fig. 4.22 Datos Distancia-Tiempo () y simulacin de disparo en Los Bronces () (onda S). ....... 33
Fig. 4.25 Datos Distancia-Tiempo con correccin de altura () .................................................... 35
viii

Fig. 4.27 Datos Distancia-Tiempo () y simulacin de disparo en El Teniente () (onda P). ......... 36
Fig. 4.28 Datos Distancia-Tiempo () y simulacin de disparo en El Teniente () (onda S). ......... 36
Fig. 4.45 Modelo inicial de velocidades 2D (onda P). .................................................................. 48
Fig. 4.46 Modelo inicial de velocidades 2D (onda S). .................................................................. 48
Fig. 4.47 Estructura de capas del modelo 2D inicial (onda P). .................................................... 49
Fig. 4.48 Estructura de capas del modelo 2D inicial (onda S). .................................................... 50
Fig. 4.49 Estructura de capas del modelo 2D (onda P y S). ........................................................ 50
Fig. 4.58 Estructura de capas del modelo 2D final (onda P y S). ................................................ 55
Fig. 5.1 Modelo final 1D y relacin Vp/Vs. Los Bronces. ............................................................. 59
Fig. 5.2 Modelo final 1D y relacin Vp/Vs. El Teniente. .............................................................. 60
ix
Fig. 5.3 Modelo final 2D (onda P). ............................................................................................... 62

Fig. 5.4 Modelo final 2D (onda S). ............................................................................................... 62
Fig. 5.5 Relocalizacin con vista en planta. Modelo 1D Los Bronces. ......................................... 65
Fig. 5.6 Relocalizacin con vista del perfil A-A. Modelo 1D Los Bronces. .................................. 66
Fig. 5.7 Relocalizacin con vista del perfil B-B. Modelo 1D Los Bronces. .................................. 66
Fig. 5.8 Relocalizacin con vista en planta. Modelo 2D Los Bronces. ......................................... 67
Fig. 5.9 Relocalizacin con vista del perfil A-A. Modelo 2D Los Bronces. .................................. 67
Fig. 5.10 Relocalizacin con vista del perfil B-B. Modelo 2D Los Bronces. ................................ 68
Fig. 5.11 Relocalizacin con vista en planta. Modelo 1D El Teniente. ......................................... 69
Fig. 5.12 Relocalizacin con vista del perfil A-A. Modelo 1D El Teniente. .................................. 69
Fig. 5.13 Relocalizacin con vista del perfil B-B. Modelo 1D El Teniente. .................................. 70
Fig. 5.14 Relocalizacin con vista en planta. Modelo 2D El Teniente. ......................................... 70
Fig. 5.15 Relocalizacin con vista del perfil A-A. Modelo 2D El Teniente. .................................. 71
Fig. 5.16 Relocalizacin con vista del perfil B-B. Modelo 2D El Teniente. .................................. 71
Fig. 5.17 Mecanismos focales en Los Bronces. .......................................................................... 72
Fig. 5.18 Mecanismos focales en El Teniente. ............................................................................ 73
Fig. 6.1 Distintos modelos 1D de Los Bronces (onda P, onda S y Vp/Vs). .................................. 75
Fig. 6.2 Comparacin de modelos 1D de El Teniente Godoy (1999) y el actual (2012). ............. 76
Fig. 6.3 Vp/Vs de modelo final 2D. .............................................................................................. 77
Fig. 6.4 Anomala de velocidades de Onda P con respecto al modelo 2D inicial. ........................ 78
Fig. 6.5 Anomala de velocidades de Onda S con respecto al modelo 2D inicial. ........................ 78
Fig. 6.6 Tomografa. Vp/Vs. Pardo (2009). ................................................................................. 79
Fig. 0.1 Trazado de rayos en Los Bronces en base a modelo sin correcciones (Onda P) ........... 86
Fig. 0.2 Trazado de rayos en Los Bronces en base a modelo sin correcciones (Onda S) ........... 86
Fig. 0.3 Trazado de rayos en El Teniente en base a modelo sin correcciones (Onda P). ............ 87
Fig. 0.4 Trazado de rayos en El Teniente en base a modelo sin correcciones (Onda S). ............ 87
ndice de Tablas
Tabla 1. Modelo de velocidades de Barrientos y otros (2004) ..................................................... 15
Tabla 2. Modelo de velocidades de Los Bronces, sin correcciones. ............................................ 15
Tabla 3. Modelo de velocidades de Godoy y otros (1999)........................................................... 20
Tabla 4 Modelo de velocidades de El Teniente, sin correcciones. .............................................. 20
Tabla 5. Modelo final 1D. Los Bronces........................................................................................ 25
Tabla 6. Modelo final 1D. El Teniente. ....................................................................................... 27
Tabla 7 Modelo 1D de El Teniente con una tercera capa. ........................................................... 45
Tabla 8. Modelo inicial 2D. Onda P. ............................................................................................ 46
Tabla 9. Modelo inicial 2D. Onda S. ............................................................................................ 47
Tabla 10 Modelo ingresado en HypoDD en base a modelo 1D de Los Bronces. ........................ 56
Tabla 11 Modelo ingresado en HypoDD en base a modelo 1D de El Teniente. .......................... 57
Tabla 12 Modelo ingresado en HypoDD en base a modelo 2D de Los Bronces. ........................ 57
Tabla 13 Modelo ingresado en HypoDD en base a modelo 2D de El Teniente. .......................... 58
Tabla 14: Modelo final 1D. Los Bronces. ..................................................................................... 60
Tabla 15. Modelo final 1D. El Teniente. ..................................................................................... 61
Tabla 16. Modelo final 2D (onda P). ............................................................................................ 63
Tabla 17. Modelo final 2D (onda S). ............................................................................................ 64
xi
Capitulo 1
Introduccin
1.1 Antecedentes
Los procesos sismotectnicos de la Zona Central de Chile estn relacionados con

la subduccin de la Placa de Nazca bajo la Placa Sudamericana. La velocidad
promedio y direccin de subduccin de la Placa de Nazca varan, dependiendo el
modelo y el perodo, entre 67 y 80 mm/ao con una direccin aproximada de 78W
[Kreemer y otros, 2003; DeMets y otros, 1994]. Se observa que al sur de los 33S la
placa de Nazca subducta con un ngulo aproximado de 30E [Cahill y otros, 1992;
Fuenzalida y otros, 1992; Pardo y otros, 2002] y que la sismicidad cortical es mayor en
el antearco que en el trasarco.
La sismicidad cortical superficial (0-50 km) en la zona central ha sido poco
estudiada tanto en su origen como en su tasa de actividad, aunque sabemos que sta
ocurre principalmente en la regin cordillerana de los Andes, y est relacionada
principalmente con la deformacin andina y su consecuente alzamiento.
Es importante hacer notar que en todo este sector y sus cercanas se concentra el
40% de la poblacin del pas, adems se encuentran 3 grandes minas cuprferas (El
Teniente, Los Bronces y Andina), algunas minas de yeso, plantas hidroelctricas y el
gasoducto proveniente de Argentina, por lo que es de suma importancia conocer y
entender los procesos geodinmicos de la regin.
En esta zona se han registrado 2 eventos ssmicos de gran magnitud localizados a
profundidades menores de 15 km y que ayudan a cuantificar el riesgo: el 4 de
Septiembre de 1958 (M=6.9) [Lomnitz, 1961; Piderit,1961; P. Acevedo,1985], (M=6.3)
[Alvarado,2009] y el 13 de Septiembre de 1987 (M=5.9) [Barrientos y Eisenberg,1988].
El terremoto de 1958 produjo daos estructurales en Santiago y en Las Melosas fue
reportada una intensidad mxima de 10 en la escala de Mercali. Slo en los ltimos
aos gracias al monitoreo y estudio en este sector se ha determinado su importante
actividad ssmica y el potencial peligro que esto conlleva [Barrientos y otros, 2003].
Este trabajo tiene el objetivo de crear un modelo de estructura de velocidades a
nivel cortical en la zona central cordillerana (33-34.5S) y conseguir as una mejor
caracterizacin del comportamiento ssmico de la zona, adems de analizar las
concentraciones de sismicidad bajo las minas Los Bronces (Ex Disputada) y El
Teniente.
Todo esto se obtiene gracias al registro de eventos en la zona de estudio

capturados por estaciones temporales y permanentes, siendo de gran importancia
aquellos registros provenientes de tronaduras asociadas a las faenas de extraccin de
la mina Los Bronces, pues estos eventos artificiales permiten la obtencin de datos
distancia-tiempo para ondas ssmicas, donde se tiene un conocimiento prcticamente
exacto de las ubicaciones tanto de la fuente como de las estaciones registradoras.
Se han realizado anterior a este trabajo otros modelos de estructura de
velocidades correspondiente a esta zona: [Acevedo, 1985], [Alvarado, 1998] y [Vera y
otros, 2006]. Estos modelos sern comparados con el modelo aqu obtenido para
lograr, de las diferencias y similitudes, una mejor comprensin de la corteza ms
superficial.
El presente trabajo se enmarca dentro del proyecto FONDECYT 1050758 Crustal
seismicity and velocity structure in the principal cordillera of central Chile (33-34.5S):
Implications on andean geodynamic and seismic hazard que pretende lograr un modelo
geodinmico general de la zona (33-34.5S).
1.2 Mina Los Bronces
La compaa minera Anglo American Chile posee yacimientos en la zona

cordillerana al noreste de la Regin Metropolitana. Los Bronces (que en el perodo de la
toma de datos y de gran parte de la investigacin perteneca a la compaa Disputada
de Las Condes) es uno de los principales yacimientos cuprferos de explotacin y se
ubica a un promedio de 3500 m.s.n.m. en el sector de las nacientes de los ros San
Francisco y Blanco en la comuna de Las Condes, a 65 km de Santiago. Sus
coordenadas geogrficas son 33 09 S y 70 17 W.
La roca husped de la mineralizacin en el yacimiento es un cuerpo de brecha
hidrotermal conocido como Brecha Los Bronces. Esta brecha aflora entre las cotas
3840 m.s.n.m. y 3600 m.s.n.m. [F. Hintze, 1978].
La explotacin del mineral (Cu y Mo) se produce a rajo abierto. Y son las
tronaduras, utilizadas para la extraccin del metal en el yacimiento Los Bronces, las
principales fuentes detectadas por las estaciones sismolgicas y que sirvieron de ayuda
para la obtencin del modelo de estructura de velocidades en esta zona.
1.3 Mina El Teniente
El Teniente es un yacimiento de tipo prfido cuprfero localizado en la alta

cordillera de los Andes, en la provincia de Cachapoal, VI regin. Se encuentra 67 km al
Noreste de la ciudad de Rancagua y cerca de 70 km al Sureste de Santiago. Sus
coordenadas geogrficas son 70 21 W y 34 05 S.
El nivel ms alto de explotacin de la mina se encuentra aproximadamente a 2400
m.s.n.m. y el ms bajo a unos 2100 m.s.n.m. La profundidad, desde la superficie, a la
que existe minera activa vara, debido a la irregularidad topogrfica, entre los 450 y 900
m.
La gnesis del yacimiento est ligada al emplazamiento de varios cuerpos
intrusivos y brechas hidrotermales en rocas volcnicas andesticas de edad Miocena
[Cuadra y Flores, 1987]. Posee un amplio sistema de fracturas, fallas y vetas, la
mayora con rumbo NE en un rea de 3 a 5 km al NE y SW de la mina; siendo su edad
mnima post Formacin Farellones [Duarte, 2000].
La falla principal de El Teniente (TFZ) es sub-vertical con rumbo NEE y una
longitud de 10 km caracterizada por un movimiento de rumbo dextral (Fig. 1.1), aqu la
roca est altamente fallada, alterada y mineralizada [Garrido y otros, 1994]. El ltimo
evento tectnico de la regin es un acortamiento en direccin NNE-SSW [Lavenu y
Cembrano, 1999; Garrido y otros, 2002].
Fig. 1.1 Zona de Falla El Teniente (ZFT).
1.4 Objetivos
El objetivo general de este trabajo es obtener un modelo de estructura de

velocidades en la zona central cordillerana (33 - 34.5S) a nivel cortical y analizar las
concentraciones de sismicidad bajo las minas Los Bronces y El Teniente, con estos
resultados ser posible mejorar el entendimiento sobre este sector, afectado
anteriormente por terremotos, pudindose realizar una mejor localizacin y
caracterizacin de futuros eventos.
Para poder lograr el objetivo general se deben cumplir los siguientes objetivos
especficos:
Determinar un modelo cortical 1D de la estructura de velocidades en la zona de

estudio.
Determinar un modelo cortical 2D de la estructura de velocidades en la zona de

estudio.
Refinar la localizacin de hipocentros con la nueva estructura de velocidades,

particularmente de los eventos concentrados por debajo de las minas El Teniente
y Los Bronces; determinar mecanismos focales de los eventos corticales en la
zona.
Capitulo 2
Marco terico
2.1 Teora de rayos
La teora de rayos es muy til para entender y obtener resultados de manera

simple sobre el comportamiento de las ondas al viajar de un punto a otro a travs de un
medio, en este caso bajo la superficie terrestre.
Gracias a la Ley de Snell sabemos que sin = p V , donde V es la velocidad de
propagacin del rayo que depende de las caractersticas del medio, es el ngulo de
incidencia, ngulo entre la vertical y el rayo, y p es el parmetro del rayo que representa
la lentitud aparente del frente de ondas en direccin horizontal.
Si consideramos que el trayecto realizado por la onda es entre dos puntos
ubicados a la misma altura y en un modelo de capas horizontales, entonces definiremos
zm como la profundidad mxima que alcanza el rayo y Vm la velocidad horizontal al
alcanzar la profundidad mxima (Vm = V(z=zm)). Como = 90 en zm, entonces en
ese punto Vm ser igual a la velocidad total.
Fig. 2.1 Geometra del rayo atravesando capa horizontal.
Por geometra simple (Fig. 2.1) podemos extraer las relaciones:
Obtenemos entonces 2 ecuaciones fundamentales para llevar a cabo el clculo

de distancia recorrida horizontalmente (X) y tiempo de viaje (T) integrando en
profundidad z. Considerando que el rayo comienza en el origen tenemos expresiones
donde tanto X como T dependen del parmetro p que define el rayo:
Es de gran utilidad, para la modelacin e interpretacin de estructuras geolgicas

bajo la superficie terrestre, poder comparar los tiempos de llegada de ondas ssmicas a
estaciones sismolgicas con resultados tericos en base a modelos de estructura de
velocidades. Por ello es necesario calcular las curvas camino-tiempo que nos entregan
el tiempo de recorrido de la onda en funcin de la distancia horizontal. Un mtodo
fidedigno para hallar esta relacin es considerar capas donde la velocidad vara
linealmente con la profundidad.
V = V0 + kz
k = (V1-V0)/H
Fig. 2.2 Gradiente de velocidad en una capa.
dV = kdz
As obtenemos
V representa la velocidad del medio en el punto ms profundo del rayo, el factor

2 es para considerar simtricamente la parte en que el rayo desciende y en la que luego
asciende. Resolviendo ambas integrales, dependientes de las propiedades del medio y
por ende del comportamiento del rayo, se pueden obtener X y T que definirn las
curvas camino-tiempo.
Para modelos de ms de una capa y con gradientes de velocidad, el clculo del
tiempo recorrido y de la distancia horizontal alcanzada por un rayo, se deben ir
sumando las distancias y los tiempos de viaje del rayo en cada capa.
Capitulo 3
Datos
3.1 Red Sismolgica
Los datos fueron proporcionados por 28 estaciones sismolgicas temporales de

tres componentes: 16 banda ancha y 12 de perodo corto, ubicadas en la zona central
del territorio chileno (33- 34.5S), registrando eventos de manera continua entre el 17
de Noviembre del 2005 y el 31 de Marzo del 2006 enmarcado dentro del proyecto
FONDECYT 1050758. Adems se trabaj con los registros diarios de las estaciones
permanentes de la red del Servicio Sismolgico de la Universidad de Chile (Fig. 3.2).
Las magnitudes de los eventos registrados fueron menores a 4.0 en la escala de
Richter.
Fig. 3.1 Red sismolgica utilizada. Perfil W-E.

o estaciones sismolgicas (latitud 33 - 34.5 S)
* origen de los eventos en Los Bronces (tronaduras).
* extremo superior de la mina El Teniente considerado como origen de los eventos (sismicidad)
Fig. 3.2 Zona de estudio (33 - 34.5S).

Las estrellas indican las dos fuentes utilizadas: Los Bronces (D) y El Teniente (T).
Las estaciones son representadas por crculos (perodo corto) y rombos (banda ancha);
los azules y verdes corresponden a estaciones temporales y los rojos a la red estable del SSN. Aparecen
adems caracterizados los eventos histricos ms importantes (1958 y 1987).
3.2 Los Bronces
Los datos provenientes de la mina Los Bronces estn asociados principalmente a

tronaduras con fines de explotacin. Para poder identificar adecuadamente estas
tronaduras, son seleccionados los eventos registrados en la zona de extraccin minera
correspondiente al sector de Los Bronces, esto es entre las latitudes 33.1 y 33.2 S y
entre las longitudes 70.3 y 70.2 W, y con una profundidad menor a 20 km obteniendo
un total de 40 eventos (Fig. 4.1 y Fig. 4.2).
Tomando en cuenta la estimacin del error que poseen los datos de profundidad
en la localizacin de hipocentros (10 km), se revisaron las horas en las que ocurrieron
los eventos. Se puede observar que todos los eventos seleccionados ocurren dentro de
una franja de tiempo que puede correlacionarse a un horario laboral (10 - 19 PM). El
evento 11 es descartado por estar dudosamente fuera del lapso horario, quedan
entonces 39 eventos seleccionados en la zona de Los Bronces (Fig. 3.3).
Fig. 3.3 Eventos seleccionados del sector de Los Bronces (33.1S - 33.2S y 70.3W - 70.2W).
El grfico superior muestra la hora chilena (GMT - 4hr) en que ocurrieron los eventos.
El grfico inferior muestra la profundidad detectada de los eventos
El evento nro. 11 marcado en rojo fue eliminado por estar fuera del lapso horario de trabajo.
10
3.3 El Teniente
Los datos provenientes de la zona de El Teniente fueron obtenidos seleccionando

los eventos registrados entre las latitudes 34.04 y 34.14 S y entre las longitudes 70.3
y 70.4 W, a una profundidad menor a 15 km. obteniendo un total de 89 eventos. A
diferencia que en Los Bronces, estos eventos no estn necesariamente ligados a
explosiones o tronaduras, sino que tendran, en su mayora, origen ssmico asociados a
las fallas existentes en este sector minero.
Se puede apreciar en la figura la inexistencia de un patrn horario. (Fig. 3.4)
Fig. 3.4 Eventos seleccionados del sector de El Teniente (34.04S - 34.14S y 70.4W - 70.3W).
El grfico superior muestra la hora chilena (GMT - 4hr) en que ocurrieron los eventos.
El grfico inferior muestra la profundidad detectada de los eventos.
11
Capitulo 4
Procesamiento de datos
4.1 Modelos 1D
Es posible obtener un modelo de estructura ssmica unidimensional (velocidad en

funcin de la profundidad) a partir del conocimiento exacto de la ubicacin geogrfica
de las estaciones sismolgicas y del lugar donde ocurre el (los) evento(s), adems de
poseer los tiempos de llegada de las ondas P y S a las estaciones sismolgicas.
Para llegar a este modelo es necesario hacer ciertas suposiciones y arreglos
hasta obtener un conjunto de datos que incluyan los tiempos de viaje de las ondas
ssmicas y las distancias entre evento y estacin. Con el programa mmodpg, creado por
Emilio Vera, escrito en FORTRAN para LINUX, se puede ir ajustando interactivamente
(prueba y error) el modelo de velocidades, superponiendo los datos reales sobre las
curvas tericas camino-tiempo obtenidas por teora de rayos (Cap. 2.1)
4.1.1 Modelo de velocidades en Los Bronces
Se tienen los registros de los tiempos de llegada de las ondas P y S en cada

estacin y un clculo estimado de las distancias entre las estaciones y los hipocentros,
pero es necesario establecer una distancia constante, que denominaremos distancia
real, entre origen y llegada (eventos y estaciones) adems de definir y calcular un
tiempo nico del trayecto entre estos puntos; en este caso se establece que los eventos
registrados en Los Bronces son tronaduras, por lo que se considera slo 1 punto
geogrfico como origen: (33.1515 S, 70.2802 W). Tomando en cuenta que la
superficie terrestre es esfrica (como aproximacin) y conociendo las coordenadas
geogrficas de tronaduras y estaciones podemos calcular la distancia real.
Para definir un tiempo nico, se fija un tiempo de llegada en una de las 43
estaciones con registros de los eventos ocurridos en Los Bronces. La estacin NEVA se
encuentra cercana al origen y tiene el registro de 28 de los 39 eventos registrados, por
lo que es la estacin escogida para fijar el tiempo de llegada de la onda P. Es eliminado
1 evento por tener un mal registro de datos, as que finalmente de los 40 eventos
iniciales, son 27 los seleccionados y utilizados para los clculos posteriores (Fig. 4.1 y
Fig. 4.2).
12
Fig. 4.1 Vista en planta de los eventos registrados en Los Bronces.

El rombo amarillo es el lugar escogido como origen de las tronaduras.
Los crculos azules son los eventos seleccionados y los verdes los eventos descartados.
Fig. 4.2 Eventos registrados en Los Bronces.

Los puntos azules son los eventos seleccionados
y los verdes los descartados. El crculo rojo es la mediana.
13
Calculando la regresin lineal de la distancia versus los tiempos de llegada

registrados en la estacin NEVA, se obtiene un tiempo fijo tf = 4.8024 s que
corresponde a los 23,2327 km de distancia real que existe entre esta estacin y el
origen de las tronaduras (Fig. 4.3).
Fig. 4.3 Clculo del tiempo mediante regresin lineal.

o datos registrados por la estacin NEVA
+ regresin lineal
+ tiempo calculado con una distancia de 23.27 Km
Al fijar un tiempo de referencia en la estacin NEVA, se pueden ajustar y corregir

todos los tiempos de llegada de la onda P a las distintas estaciones, luego se calcula la
mediana de los tiempos en cada estacin para establecer un tiempo nico de viaje
para cada estacin. Finalmente se obtiene un set de datos (43 puntos = 43 estaciones)
con la distancia real y el tiempo nico de viaje. En base a estos datos se crea, con el
programa mmodpg un modelo de velocidades (Fig. 4.4). Corrigiendo con la misma
metodologa y utilizando el mismo tiempo fijo de la estacin NEVA empleado para las
ondas P, se crea un set de datos para los tiempos de llegada de la onda S (Fig. 4.5). El
proceso para crear el modelo de velocidades es idntico al recin mencionado para las
ondas P (Tabla 2).
Se observa en las figuras (Fig. 4.4 y Fig. 4.5) que no existen datos para definir la
estructura de los kilmetros ms superficiales de la corteza. Por esta razn se utiliza un
modelo anterior, obtenido a partir de tronaduras en la mina Los Bronces, que cubri con
6 estaciones porttiles los 20 km ms cercanos a las detonaciones [Barrientos y otros,
2004]. (Tabla 1)
14
Tabla 1. Modelo de velocidades de Barrientos y otros (2004)
Capa
Profundidad [km]
(techo - base)
Vp [km/s]
(techo base)
0 2.2
4.75 - 5
2.2 8.9
5.8 6.2
8.9 15.0
6.6
1
2
3
Tabla 2. Modelo de velocidades de Los Bronces, sin correcciones.
Capa
1
2
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 2.2
4.75 5.09
2.678 2.872
2.2 33.41
5.83 7.36
3.358 4.132
15
Fig. 4.4 Modelo de velocidades en Los Bronces sin correcciones (Onda P).
Velocidad de Reduccin = 7 km/s.
Fig. 4.5 Modelo de velocidades en Los Bronces sin correcciones (Onda S).
Velocidad de Reduccin = 3.7 km/s.
16
4.1.2 Modelo de velocidades en El Teniente
Al igual que en Los Bronces (Captulo 4.1.1) se poseen los datos de tiempos de
llegada de las ondas P y S en cada estacin y las distancias desde las estaciones a los
hipocentros localizados. Sin embargo, el proceso para obtener la distancia real entre
origen y llegada (eventos y estaciones) es levemente diferente, principalmente por la
eleccin del punto de origen para los eventos de El Teniente, pues en Los Bronces los
eventos se consideraron como tronaduras realizadas en superficie, mientras que en El
Teniente los eventos se suponen principalmente de carcter ssmico y estaran
asociados a las zonas de esfuerzo a las que se ve expuesto este sector minero.
Para poder realizar un clculo estimado de la distancia real, es necesario
primero definir un criterio para la eleccin de un punto de origen para todos los eventos.
Al considerar que la mayora de los eventos estn asociado a una misma zona de
esfuerzos y fallamiento y teniendo adems en cuenta el error en la coordenada vertical
al localizar los hipocentros (10 km), se seleccionan los eventos confinados dentro de
un cilindro con radio basal de 5 km y de 12 km de altura. El centro del cilindro es
elegido calculando la mediana de todos los eventos seleccionados en un principio.
Finalmente los eventos que se ubican dentro del cilindro son 64. Para efectos
prcticos en el clculo de la distancia real, entre eventos y estaciones (Fig. 4.6 y Fig.
4.7), las coordenadas geogrficas del centro del cilindro (34.088 S, 70.341 W) son
consideradas como el origen fsico de todos los eventos.
17
Fig. 4.6 Vista en planta de los eventos registrados en El Teniente.

El rombo amarillo es el punto elegido como origen de los eventos.
Los crculos azules son los eventos seleccionados y los eventos verdes los descartados.
Fig. 4.7 Eventos registrados en El Teniente.

Los puntos azules son los eventos seleccionados y los verdes los descartados.
El crculo rojo es la mediana y centro del cilindro.
18
Para establecer un tiempo fijo, se utiliza la llegada de la onda P en la estacin

SATA, cercana al origen y con 59 eventos registrados. Al calcular la regresin lineal de
la distancia versus tiempo que tenemos de los registros de llegada en la estacin SATA,
se obtiene un tiempo fijo tf = 4.2379 s que corresponde a 20.7656 km, distancia real,
distancia entre las coordenadas geogrficas de esta estacin y del centro del cilindro
suponiendo una superficie esfrica (Fig. 4.8).
Fig. 4.8 Clculo del tiempo mediante regresin lineal.

o datos registrados por la estacin NEVA
+ regresin lineal
o dato eliminado
+ tiempo calculado con una distancia de 20.7656 Km
Luego de corregir todos los tiempos de llegada de la onda P, con el tiempo fijo de
la estacin SATA como referencia, se calcula la mediana de los tiempos en cada
estacin para establecer un tiempo nico de viaje para cada estacin.
Finalmente se obtuvo un set de datos (40 puntos = 40 estaciones) con la
distancia real y el tiempo nico de viaje. En base a estos datos se crea un modelo de
velocidades con el programa mmodpg (Fig. 4.9).
Utilizando el mismo tiempo fijo de la estacin SATA empleado para las ondas P
se crea un set de datos para los tiempos de llegada de la onda S (Fig. 4.10). El proceso
para crear el modelo de velocidades es idntico al recin mencionado para las ondas P
(Tabla 4).
Para definir la estructura del segmento ms superficial de la corteza se utiliza un
modelo realizado anteriormente obtenido por medio de refraccin ssmica usando
fuentes artificiales y gracias a una red temporal de 15 estaciones ubicadas hasta 40 km
al oeste de la mina [Godoy y otros, 1999] (Tabla 3).
19
Tabla 3. Modelo de velocidades de Godoy y otros (1999)
Capa
Profundidad [km]
(techo - base)
Vp [km/s]
(techo base)
0 2.3
4.5
2.3 6.0
6.2
1
2
Tabla 4 Modelo de velocidades de El Teniente, sin correcciones.
Capa
1
2
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 2.4
4.4 4.6
2.4 3.0
2.4 24.91
5.88 6.93
3.33 3.96
20
Fig. 4.9 Modelo de velocidades en El Teniente sin correcciones (Onda P).

Velocidad de Reduccin = 6.6 km/s
Fig. 4.10 Modelo de velocidades en El Teniente sin correcciones (Onda S).

21
4.1.3 Correccin por altura
Debido a que el tiempo que demora la onda entre origen y estacin fue calculado
para puntos ubicados a una misma cota, se realiza una correccin que considera la
diferencia vertical. Para esto se supone la existencia, para cada uno de los puntos, de
una capa de relleno de espesor igual a la diferencia de altura entre el origen de los
eventos y cada estacin (Fig. 4.11). En resumen, se establece un DATUM igual a la
altura del punto que se estableci como origen de los eventos. En el caso de Los
Bronces, la altura donde ocurren las tronaduras es mayor a la altura de todas las
estaciones, mientras que en el caso de El Teniente existen algunas estaciones que se
ubican sobre el punto de origen definido arbitrariamente en el captulo 4.1.2 (Fig. 3.1).
Utilizando los conocimientos vistos en el captulo 2.1 sobre la Ley de Snell ms
geometra bsica, algunas aproximaciones y tomando en cuenta que al ser un modelo
de capas planas el ngulo de incidencia con que parte el vector velocidad es igual al de
llegada, se pueden obtener para cada evento los valores de correccin de la distancia
horizontal y el tiempo: Xs y Ts.
= sin-1 (pV)
h = | H hestacin |
Xs = h tg
p =1/Vm
Ts = x/ V sin
El valor de V, la velocidad de la capa de relleno, se obtiene del promedio de

velocidad de la primera capa. Los valores del parmetro del rayo (p) se estiman a partir
de la relacin de Vm (Captulo 2.1) con la distancia horizontal que avanza el rayo al
llegar a la superficie (Anexo A).Finalmente, al agregar a la distancia (X) y tiempo (T)
original las correcciones Xs y Ts respectivamente, se obtiene como resultado la
Distancia corregida (Xc) y el Tiempo corregido (Tc) para cada evento.
Xc = X + Xs
Tc = T + Ts
22
Fig. 4.11 Diagrama sobre la correccin de altura.

Se muestra la capa de relleno y los distintos parmetros utilizados para realizar la correccin.
4.1.4 Correccin de altura en Los Bronces
El clculo de Xs y Ts para la onda P en Los Bronces es conseguido con los

siguientes valores: altura del origen de las tronaduras: H = 3.7 km, velocidad promedio
de la onda en la capa de relleno: V = 4.9 km/s, Vm = 6 km/s, distancia promedio de las
estaciones al origen: x = 87 km. Para las dos estaciones ms lejanas (SFDO y CHNG),
Xs y Ts fueron calculados considerando un Vm = 7 km/s que corresponde a un x = 180
km.
Para los datos de la onda S, el clculo de Xs y Ts es obtenido con los siguientes
valores: H = 3.7 km, V = 2.8 km/s, Vm = 3.5 km/s, x = 87 km. Para las dos estaciones
ms lejanas (SFDO y CHNG), Xs y Ts fueron calculados considerando un Vm = 3.9
km/s que corresponde a un x = 180 km. (Anexo A)
La modelacin con el programa mmodpg de los datos corregidos por altura (Xc,
Tc) pueden observarse en Fig. 4.12 y Fig. 4.13 y el modelo obtenido se puede apreciar
en la Tabla 5.
23
Fig. 4.12 Modelo de velocidades en Los Bronces con correccin de altura (Onda P).
Velocidad de Reduccin = 7 km/s
Fig. 4.13 Modelo de velocidades en Los Bronces con correccin de altura (Onda S).
24
Tabla 5. Modelo final 1D. Los Bronces.
Capa
1
2
3
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 4.1
4.76 5.07
2.8 2.83
4.1 26.21
5.94 6.95
3.45 3.86
26.21 30.09
7.09 7.17
3.91 3.99
4.1.5 Correccin de altura en El Teniente
El clculo de Xs y Ts para la onda P en El Teniente, es conseguido con los

siguientes valores: origen de los eventos: H = 2.4 km, V = 4.4 km/s, Vm = 6 km/s, x =
51.479 km. Para la estacin ms lejana (JACH), Xs y Ts fueron calculados
considerando un Vm = 6.9 km/s que corresponde a un x = 158.03 km.
Para los datos de la onda S, el clculo de Xs y Ts se obtiene considerando los
siguientes valores: H = 2.4 Km, V = 2.5 km/s, Vm = 3.4 km/s, x = 51 km. Para la
estacin ms lejana (JACH), Xs y Ts fueron calculados considerando un Vm = 3.8 km/s
que corresponde a un x = 158 km. (Anexo A)
La eleccin arbitraria del origen de los eventos a 2.4 km de altura es debido a
que es en esa cota donde comienzan las labores de extraccin de la mina y es ms
cmodo definir un punto de origen cercano a la superficie para no sobrestimar el
espesor de la capa de relleno.
La modelacin con el programa mmodpg de los datos corregidos por altura (Xc,
Tc) pueden observarse en Fig. 4.14 y Fig. 4.15 y el modelo obtenido se puede apreciar
en la Tabla 6.
25
Fig. 4.14 Modelo de velocidades de El Teniente con correccin de altura (Onda P).
Fig. 4.15 Modelo de velocidades de El Teniente con correccin de altura (Onda S).
26
Tabla 6. Modelo final 1D. El Teniente.
Capa
1
2
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 2.4
4.4 4.5
2.53 2.57
2.4 26.31
5.81 6.97
3.35 3.94
4.2 Modelacin 2D
Para crear un perfil 2D se ha utilizado el programa RayTracing creado por Emilio

Vera escrito en FORTRAN y modificado a C++ para LINUX. Se han ocupado las
mismas ecuaciones para realizar la modelacin en 1D, pero con una complejidad mayor
en el clculo de las curvas camino-tiempo debido a la existencia de bloques verticales y
a capas con espesor variable. La incorporacin de bloques en el modelo permite
modelar fallas geolgicas o, como en este caso, la topografa del perfil.
El programa permite crear una malla compuesta por capas (horizontales) y
bloques (verticales) formando celdas, con la restriccin de que ambos lados verticales
son paralelos entre s y perpendiculares al eje x. Cada trapecio (celda) es definido por
sus cuatro nodos (vrtices), cada nodo posee informacin sobre su ubicacin espacial
(x, z) y velocidad (v). Como se puede ver en la figura (Fig. 4.16) la velocidad de un
punto P dentro del trapecio (VP) es calculada con las velocidades VP1 y VP2 de manera
proporcional a las distancias a las que se encuentra de los puntos P1 y P2
respectivamente. A la vez, VP1 es calculado con las velocidades V1A y V1B, proporcional
a las distancias a las que se encuentra el punto P1 de los puntos A y B respectivamente,
asimismo para VP2 con respecto a los puntos C y D y sus velocidades V2A y V2B.
VP1 = [(AB P1A)/AB]*V1A + [P1A/AB]*V2A

VP2 = [(CD P2C)/CD]*V1B + [P2C/CD]*V2B
VP = [(P1P2 PP1)/ P1P2]*VP1 + [PP1/ P1P2]*VP2
27
Fig. 4.16 Celda Trapezoidal.

Formada al dividir el modelo en capas (horizontales) y bloques (verticales).
La velocidad en cada punto es calculada en base
a la velocidad de los nodos que definen la celda.
El programa permite crear o modificar un modelo de velocidades en dos

dimensiones y en base a simulacin de disparos visualizar ondas ssmicas que
atraviesan las capas y bloques definidos en el modelo, calculando distancias y tiempos
de llegada de las ondas en superficie que pueden adems compararse grficamente
con algn set de datos, en este caso, los eventos en Los Bronces y El Teniente
registrados por la red de estaciones sismolgicas (Captulo 3.1). Trabajar con los
tiempos de llegada de las ondas P y S permite ir modificando el modelo (x, z, Vp y Vs)
para que las curvas tericas de camino-tiempo puedan ajustarse a los datos. Todo esto
de manera interactiva e iterativa (prueba y error).
4.2.1 Eleccin del Perfil
La extensin y orientacin del perfil utilizado para realizar la modelacin 2D se

define en base a la recta que une los puntos geogrficos correspondientes a las minas
Los Bronces (33.1515 S, 70.2802 W) y El Teniente (34.088 S, 70.341 W). (Fig.
4.17).La longitud del perfil que es ingresado al programa de modelacin posee una
longitud de 238,02 km, pues se extiende hasta las estaciones ms alejadas que
registraron eventos, tanto de Los Bronces como de El Teniente. El perfil topogrfico real
y el perfil topogrfico simplificado para definir los bloques del modelo, se muestran en la
figura (Fig. 4.18).
28
Fig. 4.17 Perfil entre Los Bronces y El Teniente.

Utilizado para modelar en 2D.
Fig. 4.18 Perfil Norte-Sur.

Topografa real (azul) y topografa utilizada en la modelacin (rojo).
29
4.2.2 Seleccin de estaciones y correccin de altura
Para comparar correctamente datos reales con resultados tericos al modelar en

2D, es necesario que las estaciones sismolgicas se encuentren cercanas al perfil para
tener as una buena representacin del sector por donde ste cruza. Se seleccionan
entonces las estaciones que se encuentran a menos de 12 km de la recta que define al
perfil obteniendo un total de 12 estaciones (Fig. 4.19). Se puede apreciar que existen 5
estaciones que se encuentran ms al sur de El Teniente, estas estaciones permiten que
se pueda obtener un perfil final de estudio con una extensin equivalente a 132,82 km,
longitud superior al tramo inicial Los Bronces El Teniente.
Para efectos de la modelacin, se considera que las estaciones se encuentran
en la recta del perfil (Fig. 4.20). Debido a esto es necesario hacer una correccin de los
tiempos de llegada en las estaciones, ya que la altura real a la que se encuentra cada
estacin es diferente a la altura que indica la topografa utilizada en la modelacin y el
tiempo terico que tarda la onda en llegar a la superficie es distinto al tiempo registrado
por la estacin.
Para corregir este desfase del tiempo, se debe sumar, o restar, dependiendo de
la ubicacin de la estacin con respecto al perfil topogrfico del modelo, el tiempo
calculado debido a la diferencia de altura y este resultado es dividido por la velocidad
de propagacin de la onda cercana a la superficie. En este caso se utiliz para la onda
P una velocidad (V) de 4.7 km/s y para la Onda S de 2.7 km/s.
dh = htopo hest
dt = dh/V
30
tf = t + dt
Fig. 4.19 Seleccin de estaciones cercanas al perfil.
Fig. 4.20 Perfil Norte-Sur y estaciones seleccionadas.

Topografa real (azul) y topografa simple utilizada para la modelacin 2D (rojo). Las letras D y T indican Los
Bronces y El Teniente. Los puntos verdes son las estaciones pertenecientes a la franja de 24 km. Las lneas
verticales marcan el corte inicial y final del perfil finalmente utilizado.
31
4.2.3 Modelos 1D
Para comenzar el modelamiento 2D se utilizan los modelos de velocidades 1D ya

revisados (Cap. 4.1). En ambos modelos se ha considerado el origen de los eventos
como la cota vertical donde se da inicio al modelo (z=0) (3700 m las tronaduras en Los
Bronces y 2400 m el origen terico de la sismicidad que coincide con la altura mxima
de extraccin en la mina El Teniente).
Es expuesta grficamente la superposicin de datos tericos con:
1) Los registros de todas las estaciones (Fig. 4.21, Fig. 4.22, Fig. 4.27, Fig. 4.28,
Fig. 4.33, Fig. 4.34, Fig. 4.39 y Fig. 4.40).
2) Los registros de las 12 estaciones seleccionadas, comprendidas en la franja de
24 Km (Fig. 4.23, Fig. 4.24, Fig. 4.29, Fig. 4.30, Fig. 4.35, Fig. 4.36, Fig. 4.41 y
Fig. 4.42).
3) Los registros de las 12 estaciones seleccionadas con su correccin de altura.
Esta correccin de los registros se mantendr finalmente durante toda la
modelacin (Fig. 4.25, Fig. 4.26, Fig. 4.31, Fig. 4.32, Fig. 4.37, Fig. 4.38, Fig.
4.43 y Fig. 4.44).
Las curvas camino-tiempo (puntos negros) indican las llegadas de las ondas a la
superficie, en base al modelo terico, gracias a la simulacin de un disparo en Los
Bronces (altura = 3700 m) y otro disparo en El Teniente (altura = 2400 m) segn sea el
caso. Los puntos rojos son los registros existentes de las llegadas de las ondas a las
estaciones en base a los eventos de Los Bronces, los puntos azules en base a los
eventos de El Teniente.
Al comparar las curvas camino-tiempo en base a los modelos unidimensionales con
los datos reales de eventos registrados, tanto en Los Bronces como en El Teniente, se
puede apreciar que el ajuste es mucho mejor con la simulacin del disparo en Los
Bronces, pues los datos provenientes de El Teniente poseen una gran dispersin, tanto
para la onda P como para la onda S.
32
Simulaciones con modelo 1D Los Bronces Disparo Los Bronces
Fig. 4.21 Datos Distancia-Tiempo () y simulacin de disparo en Los Bronces () (onda P).
Modelo 1D Los Bronces. (Velocidad de reduccin = 7 km/s).
Fig. 4.22 Datos Distancia-Tiempo () y simulacin de disparo en Los Bronces () (onda S).
Modelo 1D Los Bronces. (Velocidad de reduccin = 3.7 km/s).
33
Seleccin de estaciones cercanas al perfil (franja de 24 km).
34
Fig. 4.25 Datos Distancia-Tiempo con correccin de altura ()

y simulacin de disparo en Los Bronces () (onda P).

y simulacin de disparo en Los Bronces () (onda S).
35
Simulaciones con modelo 1D Los Bronces Disparo El Teniente
Fig. 4.27 Datos Distancia-Tiempo () y simulacin de disparo en El Teniente () (onda P).

Fig. 4.28 Datos Distancia-Tiempo () y simulacin de disparo en El Teniente () (onda S).

36


37

y simulacin de disparo en El Teniente () (onda P).

y simulacin de disparo en El Teniente () (onda S).
38
Simulaciones con modelo 1D El Teniente Disparo Los Bronces
Modelo 1D El Teniente. (Velocidad de reduccin = 7 km/s).
Modelo 1D El Teniente. (Velocidad de reduccin = 3.7 km/s).
39
40

y simulacin de disparo en Los Bronces () (onda P).

y simulacin de disparo en Los Bronces () (onda S).
41
Simulaciones con modelo 1D El Teniente Disparo El Teniente


42


43

y simulacin de disparo en El Teniente () (onda P).

y simulacin de disparo en El Teniente () (onda S).
44
4.2.4 Modelacin 2D Los Bronces - El Teniente
Para iniciar la modelacin 2D se usan como base los modelos unidimensionales

de Los Bronces y El Teniente obtenidos anteriormente. Como cada modelo se
encuentra a un extremo del sector a modelar se unen linealmente los puntos que tienen
las mismas velocidades. A veces ocurre que en las interfases de las capas definidas en
alguno de los modelos 1D existen velocidades que no se encuentran en el otro modelo,
se ha optado, en estos casos, hacer una variacin gradual de velocidad entre un
extremo y otro.
El programa de modelacin en 2D exige una misma cantidad de capas en todo el
modelo, y debido a que el modelo unidimensional de El Teniente tiene 2 capas y el de
Los Bronces 3 capas, se extiende el modelo de El Teniente a 3 capas; esta nueva capa
llega a la misma profundidad de la tercera capa de Los Bronces y aumenta
gradualmente su velocidad con la profundidad, partiendo a la velocidad con que acaba
su segunda capa hasta llegar a la misma velocidad que alcanza la de Los Bronces en
su modelo (Tabla 7).
Tabla 7 Modelo 1D de El Teniente con una tercera capa.
Capa
1
2
3
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 2.4
4.4 4.5
2.53 2.57
2.4 26.31
5.81 6.97
3.35 3.94
26.31 28.79
6.97 7.17
3.94 3.99
Dado que slo se utilizan para la modelacin los datos (distancia-tiempo) de las
12 estaciones cercanas al perfil, la extensin total del perfil y por lo tanto del modelo es
de 132.82 km, siendo el 0 la ubicacin de Los Bronces y aumentando positivamente la
distancia hacia el Sur. Se mantiene la correccin de altura realizada a los registros
(distancia-tiempo) (Captulo 4.2.2).
La grilla base, ingresada al programa, para la onda P y S del modelo inicial 2D
pueden apreciarse en Tabla 8 y Tabla 9 donde aparecen la ubicacin y propiedades de
velocidad de los nodos de los trapecios que forman el modelo.
Las figuras obtenidas del modelo inicial que se obtiene en base a los modelos 1D
de Los Bronces y El Teniente pueden apreciarse en Fig. 4.45 y Fig. 4.46.
45
Tabla 8. Modelo inicial 2D. Onda P.
z1a
v1a
z1b
v1b
v2a
z2
v2b
v3a
z3
v3b
v4a
z4
v4b
v5a
z5
v5b
0
8.16
18.16
23.66
34.16
45.16
57.76
63
78.16
104.27
111.66
122.16
126.16
132.82
-3.7
-4
-2.65
-2.3
-2.6
-1.05
-1.95
-1.1
-3.5
-2.85
-2.8
-1.4
-3
-1.13
4.76
4.76
4.76
4.76
4.76
4.76
4.76
4.76
4.51
4.42
4.42
4.42
4.42
4.42
0.4
0.37
0.33
0.31
0.27
0.23
0.18
0.16
0.1
0
0
0
0
0
5.07
5.07
4.97
4.94
4.88
4.82
4.75
4.73
4.64
4.5
4.5
4.5
4.5
4.5
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
0.4
0.58
0.8
0.92
1.15
1.39
1.66
1.78
2.11
2.68
2.68
2.68
2.68
2.68
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
22.51
22.59
22.68
22.73
22.83
22.94
23.06
23.11
23.25
23.5
23.5
23.5
23.5
23.5
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
22.51
22.62
22.75
22.83
22.97
23.12
23.29
23.36
23.56
23.91
23.91
23.91
23.91
23.91
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
6.97
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
7.09
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
7.17
46
Tabla 9. Modelo inicial 2D. Onda S.
z1a
v1a
z1b
v1b
v2a
z2
v2b
v3a
z3
v3b
v4a
z4
v4b
v5a
z5
v5b
v6a
z6
v6b
0
8.16
18.16
23.66
34.16
45.16
57.76
63
78.16
104.27
111.66
122.16
126.16
132.82
-3.7
-4
-2.65
-2.3
-2.6
-1.05
-1.95
-1.1
-3.5
-2.85
-2.8
-1.4
-3
-1.13
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.8
2.6
2.53
2.53
2.53
2.53
2.53
0.4
0.37
0.33
0.31
0.27
0.23
0.18
0.16
0.1
0
0
0
0
0
2.81
2.78
2.77
2.74
2.72
2.69
2.67
2.64
2.57
2.57
2.57
2.57
2.57
2.57
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
0.4
0.69
1.04
1.23
1.6
1.98
2.42
2.61
3.14
4.05
4.05
4.05
4.05
4.05
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
3.45
22.51
22.37
22.19
22.09
21.91
21.71
21.49
21.4
21.13
20.67
20.67
20.67
20.67
20.67
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
3.86
22.51
22.52
22.54
22.55
22.57
22.59
22.61
22.62
22.65
22.69
22.69
22.69
22.69
22.69
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
3.91
23.96
23.96
23.96
23.95
23.95
23.94
23.93
23.93
23.92
23.91
23.91
23.91
23.91
23.91
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
3.94
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
3.99
47
Modelo Inicial Los Bronces El Teniente
Fig. 4.45 Modelo inicial de velocidades 2D (onda P).
Fig. 4.46 Modelo inicial de velocidades 2D (onda S).
48
Para lograr un buen ajuste de datos y mejorar el modelo bidimensional es

necesario tener una idntica estructura de capas (nmero de capas, espesor y
profundidad) entre los modelos de las ondas P y S. Como estas estructuras de capas
resultan ser diferentes (Fig. 4.47 y Fig. 4.48) es que deben combinarse sus
caractersticas espaciales de manera que ambas se vean representadas. Luego por
limitaciones del programa es necesario simplificar lo ms posible el modelo (Fig. 4.49).
Se entregan las comparaciones grficas (antes y despus) del ajuste de las
curvas camino-tiempo entregadas por el modelo con simulaciones de disparo en Los
Bronces y El Teniente (Fig. 4.50 a Fig. 4.57).
Fig. 4.47 Estructura de capas del modelo 2D inicial (onda P).
49
Fig. 4.48 Estructura de capas del modelo 2D inicial (onda S).
Fig. 4.49 Estructura de capas del modelo 2D (onda P y S).

Resultado de la combinacin de los modelos 2D iniciales de la onda P y S (Fig. 4.47 y Fig. 4.48).
Este es el modelo inicial utilizado para el ajuste final de las curvas tericas con los datos reales.
50
Simulaciones con modelo 2D inicial Disparo Los Bronces
Modelo 2D inicial. (Velocidad de reduccin = 7 km/s).
Modelo 2D inicial. (Velocidad de reduccin = 3.7 km/s).
51
Simulaciones con modelo 2D final Disparo Los Bronces
Modelo 2D final. (Velocidad de reduccin = 7 km/s).
Modelo 2D final. (Velocidad de reduccin = 3.7 km/s).
52
Simulaciones con modelo 2D inicial Disparo El Teniente

Modelo 2D inicial. (Velocidad de reduccin = 7 km/s).

Modelo 2D inicial. (Velocidad de reduccin = 3.7 km/s)
53
Simulaciones con modelo 2D final Disparo El Teniente

Modelo 2D final. (Velocidad de reduccin = 7 km/s).

Modelo 2D final. (Velocidad de reduccin = 3.7 km/s).
54
Al modificar el modelo, alterando los espesores de las capas y principalmente

sus velocidades, las curvas tericas generadas por el programa se modifican,
mejorando su ajuste con los datos reales, obteniendo entonces un modelo
bidimensional definitivo (Fig. 4.58).
Debida a la excesiva dispersin de los datos en El Teniente, que impide observar
una referencia fidedigna, es que la modelacin final fue realizada considerando slo los
datos en base a Los Bronces.
Se puede ver en las figuras que, tanto el modelo de velocidades de la Onda P
como el de la Onda S, es en las estaciones ms lejanas al origen de los eventos donde
existe mayor desajuste en las curvas camino-tiempo, y por lo mismo son estas las que
principalmente condicionan los cambios finales efectuados en las capas del modelo.
Fig. 4.58 Estructura de capas del modelo 2D final (onda P y S).

Resultado del ajuste de las curvas camino-tiempo con los datos reales.
4.3 Relocalizacin
Relocalizar implica calcular en base a nuevos modelos la ubicacin de hipocentros

de eventos antes ya localizados. Mientras los modelos de velocidades ms se
asemejen a la realidad, mejor ser la estimacin de los hipocentros.
Para la relocalizacin de los eventos es utilizado el programa gratuito HypoDD
(http://earthquake.usgs.gov/research/software) que vuelve a calcular la ubicacin de los
hipocentros usando el mtodo de Doble Diferencia [Walhauser y Ellsworth, 2000]. El
programa tiene la limitacin de trabajar con un modelo de capas con velocidad
constante. Debido a que en los modelos obtenidos cada capa tiene un gradiente de
55
velocidad, al programa se ingresan modelos con una gran cantidad de capas, emulando
as los cambios de velocidad.
Se realiza la relocalizacin de los eventos en base a los modelos de velocidades
1D y 2D obtenidos en captulos anteriores. En el caso 2D se ocupan las estructuras de
capas ubicadas bajo las minas Los Bronces y El Teniente que en principio sirvieron de
referencia y fueron luego modificadas al modelar el perfil. Estas pueden considerarse
como nuevos modelos de velocidades 1D, corregidos y mejorados. Los datos
proporcionados por el programa permiten adems la comparacin grfica entre la
localizacin y la relocalizacin de los eventos registrados.
Las estructuras de capas, basadas en los modelos unidimensionales y
bidimensionales, ingresadas al programa para la relocalizacin de los eventos en Los
Bronces y El Teniente pueden observarse en la Tabla 10, Tabla 11, Tabla 12 y Tabla
13.
Tabla 10 Modelo ingresado en HypoDD en base a modelo 1D de Los Bronces.
Capa
1
2
3
4
5
6
7
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vp [km/s]
promedio
0 4.1
4.76 5.07
4.915
4.1 8.522
5.94 6.142
6.041
8.522 12.944
6.142 6.344
6.243
12.944 17.366
6.344 6.546
6.455
17.366 21.788
6.546 6.748
6.647
21.788 26.21
6.748 6.95
6.849
26.21 30.09
7.09 7.17
7.13
56
Tabla 11 Modelo ingresado en HypoDD en base a modelo 1D de El Teniente.
Capa
1
2
3
4
5
6
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vp [km/s]
promedio
0 2.4
4.4 4.5
4.45
2.4 7.182
5.81 6.042
5.926
7.182 11.964
6.042 6.274
6.158
11.964 16.746
6.274 6.506
6.39
16.746 21.528
6.506 6.738
6.622
21.528 26.31
6.738 6.97
6.854
Tabla 12 Modelo ingresado en HypoDD en base a modelo 2D de Los Bronces.
Capa
1
2
3
4
5
6
7
8
9
10
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vp [km/s]
promedio
0 4.1
4.76 5
4.88
4.1 4.7
5.81 5.94
5.875
4.7 8.807
5.94 6.124
6.032
8.807 12.914
6.124 6.308
6.216
12.914 17.021
6.308 6.492
6.4
17.021 21.128
6.492 6.676
6.584
21.128 25.235
6.676 - 6.86
6.768
25.235 26.21
6.86 6.94
6.9
26.21 27.66
7.05 7.1
7.075
27.66 30.09
7.1 7.14
7.12
57
Tabla 13 Modelo ingresado en HypoDD en base a modelo 2D de El Teniente.
Capa
1
2
3
4
5
6
7
8
9
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vp [km/s]
promedio
0 2.85
4.55 4.71
4.63
2.85 4.287
5.82 5.98
5.9
4.287 8.34
5.98 6.144
6.062
8.34 12.39
6.144 6.308
6.226
12.39 16.445
6.308 6.472
6.39
16.445 20.497
6.472 6.636
6.554
20.497 24.55
6.636 - 6.8
6.718
24.55 26.76
6.8 7.09
6.945
26.76 29.24
7.09 7.14
7.115
4.4 Mecanismos focales
Para hallar los mecanismos focales que caracterizan cada evento, es utilizado el
programa FOCMEC en SEISAN que, en base a los tiempos y formas de llegada de las
ondas P y S, entrega los valores del ngulo de rumbo, manteo y deslizamiento que
definen el mecanismo. Con estos datos, las figuras se realizan de manera sencilla con
GMT.
58
Capitulo 5
Resultados
5.1 Modelo de Velocidades 1-D Los Bronces
Se puede apreciar que tanto la velocidad de la onda P como de la onda S van en

aumento a medida que la profundidad aumenta, diferencindose 3 capas. En cuanto al
Vp/Vs, se tiene un promedio de 1.76 y tambin su valor aumenta a medida que se
avanza en profundidad, se puede apreciar un cambio importante en el cambio de fase a
los 4.1 Km de profundidad. Hay que recordar que este primer cambio de fase es
heredado de modelos realizados anteriormente, pues en este estudio no existen datos
que entreguen informacin sobre los primeros kilmetros en profundidad. Los
resultados numricos son presentados en la Tabla 14 y el modelo de velocidades
unidimensional en Los Bronces se puede apreciar en la Fig. 5.1.
Fig. 5.1 Modelo final 1D y relacin Vp/Vs. Los Bronces.
59
Tabla 14: Modelo final 1D. Los Bronces.
Capa
1
2
3
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 4.1
4.76 5.07
2.8 2.83
4.1 26.21
5.94 6.95
3.45 3.86
26.21 30.09
7.09 7.17
3.91 3.99
5.2 Modelo de Velocidades 1-D El Teniente
Al igual que en Los Bronces la velocidad de la onda P como de la onda S van en

aumento a medida que la profundidad aumenta, se pueden diferenciar 2 capas. En
cuanto a los valores de Vp/Vs se tiene un promedio de 1.75 y se puede apreciar un
cambio importante en el cambio de fase a los 2.4 Km de profundidad. El cambio de fase
es heredado de un modelo realizado anteriormente, pues en este estudio no existen
datos que entreguen informacin sobre los primeros kilmetros en profundidad. Los
resultados numricos son presentados en la Tabla 15 y El modelo de velocidades
unidimensional en El Teniente se puede apreciar en la Fig. 5.2.
Fig. 5.2 Modelo final 1D y relacin Vp/Vs. El Teniente.
60
Tabla 15. Modelo final 1D. El Teniente.
Capa
1
2
Profundidad [km]
(techo-base)
Vp [km/s]
(techo base)
Vs [km/s]
(techo base)
0 2.4
4.4 4.5
2.53 2.57
2.4 26.31
5.81 6.97
3.35 3.94
5.3 Modelo de Velocidades 2-D (Perfiles)
Los modelos bidimensionales de velocidad para la Onda P y S pueden observarse

en las figuras (Fig. 5.3 y Fig. 5.4). El perfil va de Norte a Sur, Los Bronces se encuentra
en el 0 de la horizontal y El Teniente a los 104 km. Como es lgico, debido a la
construccin en base a los modelos 1D, las velocidades van aumentando a medida que
aumenta la profundidad
Con el fin de no provocar una interpretacin errada, cabe recordar que al
comenzar la modelacin 2D (Captulo 4.2.4), por requerimiento del programa, al modelo
1D de El Teniente le fue agregada una capa en su tramo ms profundo; por lo tanto
existe en el modelo final un sector calculado en base a esta capa agregada
artificialmente. Este sector correspondera a los 2.48 Km ms profundos de la zona
cercana a El Teniente.
Se entregan las tablas con la grilla base del modelo final 2D que contiene la
ubicacin y velocidad de los nodos de los trapecios que definen la estructura del
modelo (Tabla 16 y Tabla 17).
61
Modelo Final Los Bronces El Teniente
Fig. 5.3 Modelo final 2D (onda P).
Fig. 5.4 Modelo final 2D (onda S).
62
Tabla 16. Modelo final 2D (onda P).
z1a
v1a
z1b
v1b
v2a
z2
v2b
v3a
z3
v3b
v4a
z4
v4b
v5a
z5
v5b
v6a
z6
v6b
0.0
8.16
18.16
23.66
34.16
45.16
57.76
63.0
78.16
104.27
111.66
122.16
126.16
132.82
-3.7
-4
-2.65
-2.3
-2.6
-1.05
-1.95
-1.1
-3.5
-2.85
-2.8
-1.4
-3
-1.13
4.76
4.76
4.76
4.74
4.73
4.72
4.68
4.63
4.57
4.55
4.54
4.54
4.53
4.53
0.4
0.37
0.33
0.31
0.27
0.23
0.18
0.16
0.1
0
0
0
0
0
5
5
4.98
4.92
4.87
4.85
4.8
4.78
4.74
4.71
4.72
4.73
4.7
4.7
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.81
5.82
5.82
5.83
5.81
5.80
1
1.15
1.3
143.719
143.719
153.428
163.137
163.137
153.428
143.719
143.719
163.137
192.265
192.265
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.95
5.96
5.98
5.97
5.97
5.95
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.94
5.95
5.96
5.98
5.97
5.97
5.95
5.95
215.354
205.644
201.761
202.732
205.644
212.441
216.325
219.237
21.9
21.7
21.7
21.7
21.7
21.7
6.86
6.87
6.87
6.88
6.87
6.85
6.84
6.83
6.82
6.8
6.8
6.8
6.8
6.8
6.86
6.87
6.87
6.88
6.87
6.85
6.84
6.83
6.82
6.8
6.8
6.8
6.8
6.8
22.51
220.208
217.296
221.179
223.121
225.063
22.96
23.36
23.56
23.91
23.91
23.91
23.91
23.91
6.94
6.95
6.96
6.97
6.96
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
6.95
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
7.05
23.96
237.685
236.714
236.714
236.714
236.714
23.93
23.93
23.92
23.91
23.91
23.91
23.91
23.91
7.10
7.11
7.12
7.13
7.12
7.11
7.1
7.1
7.1
7.09
7.09
7.09
7.09
7.09
7.10
7.11
7.12
7.13
7.12
7.11
7.1
7.1
7.1
7.09
7.09
7.09
7.09
7.09
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
7.14
63
Tabla 17. Modelo final 2D (onda S).
z1a
v1a
z1b
v1b
v2a
z2
v2b
v3a
z3
v3b
v4a
z4
v4b
v5a
z5
v5b
v6a
z6
v6b
0.0
8.16
18.16
23.66
34.16
45.16
57.76
63.0
78.16
104.27
111.66
122.16
126.16
132.82
-3.7
-4
-2.65
-2.3
-2.6
-1.05
-1.95
-1.1
-3.5
-2.85
-2.8
-1.4
-3
-1.13
2.75
2.73
2.73
2.7
2.7
2.69
2.69
2.68
2.65
2.62
2.6
2.6
2.59
2.59
0.4
0.37
0.33
0.31
0.27
0.23
0.18
0.16
0.1
0
0
0
0
0
2.81
2.8
2.78
2.75
2.73
2.71
2.7
2.7
2.7
2.67
2.67
2.65
2.63
2.63
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
3.35
1
1.15
1.3
143.719
143.719
153.428
163.137
163.137
153.428
143.719
143.719
163.137
192.265
192.265
3.46
3.46
3.46
3.46
3.46
3.45
3.45
3.46
3.46
3.46
3.46
3.46
3.45
3.44
3.46
3.46
3.46
3.46
3.46
3.45
3.45
3.46
3.46
3.46
3.46
3.46
3.45
3.44
215.354
205.644
201.761
202.732
205.644
212.441
216.325
219.237
21.9
21.7
21.7
21.7
21.7
21.7
3.8
3.82
3.85
3.87
3.87
3.87
3.88
3.88
3.88
3.89
3.89
3.89
3.89
3.89
3.8
3.82
3.85
3.87
3.87
3.87
3.88
3.88
3.88
3.89
3.89
3.89
3.89
3.89
22.51
220.208
217.296
221.179
223.121
225.063
22.96
23.36
23.56
23.91
23.91
23.91
23.91
23.91
3.9
3.9
3.91
3.91
3.91
3.91
3.92
3.92
3.92
3.93
3.93
3.93
3.93
3.93
3.91
3.91
3.92
3.92
3.92
3.92
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
23.96
237.685
236.714
236.714
236.714
236.714
23.93
23.93
23.92
23.91
23.91
23.91
23.91
23.91
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
3.93
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
26.39
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
3.97
64
5.4 Relocalizacin
Se puede apreciar en las figuras siguientes la localizacin y relocalizacin de

los hipocentros de los eventos registrados en Los Bronces en base al modelo de
velocidades 1D entregado en Cap. 5.1. Las comparaciones son con vista en planta
(Fig. 5.5) y en cortes verticales en base a perfiles transversales entre s (Fig. 5.6 y
Fig. 5.7).
Los perfiles fueron elegidos en base a la informacin sobre el intenso
fallamiento en profundidad en Los Bronces, lugar de contacto entre la brecha Los
Bronces y el Plutn Los Bronces [F. Hintze, 1978]. El rumbo del sistema de
fallamiento en esta zona es N70E, por lo que los rumbos escogidos para los perfiles
A-A y B-B son N70E y N160E respectivamente y son los mismos utilizados en la
grfica de relocalizacin en base al modelo 2D.
Fig. 5.5 Relocalizacin con vista en planta. Modelo 1D Los Bronces.

Eventos localizados (izquierda) y relocalizados en base al modelo 1D (derecha).
65
Fig. 5.6 Relocalizacin con vista del perfil A-A. Modelo 1D Los Bronces.
Fig. 5.7 Relocalizacin con vista del perfil B-B. Modelo 1D Los Bronces.
66
En las siguientes figuras se aprecia la localizacin y relocalizacin de los

hipocentros de los eventos registrados en Los Bronces en base al modelo de
velocidades 2D entregado en Cap. 5.3. Las comparaciones son con vista en planta
Fig. 5.10), al igual que en la relocalizacin en base al modelo 1D, el perfil A-A tiene
un rumbo N70E y el perfil B-B un rumbo de N160E.
Fig. 5.8 Relocalizacin con vista en planta. Modelo 2D Los Bronces.

Fig. 5.9 Relocalizacin con vista del perfil A-A. Modelo 2D Los Bronces.
67
Fig. 5.10 Relocalizacin con vista del perfil B-B. Modelo 2D Los Bronces.
En las figuras que aparecen a continuacin se aprecia la localizacin y

relocalizacin de los hipocentros de los eventos registrados en El Teniente en base
al modelo de velocidades 1D mostrado en Cap. 5.2. Las comparaciones son con
vista en planta (Fig. 5.11) y en cortes verticales en base a perfiles transversales
entre s (Fig. 5.12 y Fig. 5.13).
Los perfiles fueron elegidos en base a la informacin de la Zona de falla El
Teniente (Cap 1.3) [Garrido y otros, 2002]. Se puede ver con ms detalle la ZFT en
la figura del estudio del modelo geolgico estructural del yacimiento El Teniente (Fig.
1.1) [Garrido y otros, 1994]. Los rumbos escogidos para los perfiles A-A y B-B son
N45E y N135E respectivamente y son los mismos utilizados en la grfica de
relocalizacin en base al modelo 2D.
68
Fig. 5.11 Relocalizacin con vista en planta. Modelo 1D El Teniente.

Fig. 5.12 Relocalizacin con vista del perfil A-A. Modelo 1D El Teniente.
69
Fig. 5.13 Relocalizacin con vista del perfil B-B. Modelo 1D El Teniente.
En las siguientes figuras se aprecia la localizacin y relocalizacin de los

hipocentros de los eventos registrados en El Teniente en base al modelo de
velocidades 2D obtenido en Cap. 5.3. Las comparaciones son con vista en planta
Fig. 5.16), al igual que en la relocalizacin en base al modelo 1D, el perfil A-A tiene
un rumbo N45E y el perfil B-B un rumbo de N135E.
Fig. 5.14 Relocalizacin con vista en planta. Modelo 2D El Teniente.

70
Fig. 5.15 Relocalizacin con vista del perfil A-A. Modelo 2D El Teniente.
Fig. 5.16 Relocalizacin con vista del perfil B-B. Modelo 2D El Teniente.
En todos los casos, los cambios de ubicacin de los eventos entre la

localizacin y la relocalizacin son leves y salvo una pequea tendencia a la
alineacin vertical de eventos en El Teniente, no revelan mayor detalle.
71
5.5 Mecanismos focales
En Los Bronces se obtuvieron mecanismos focales de slo 3 eventos, y aunque

fue supuesto desde el comienzo que todos los eventos eran tronaduras, podra
considerarse a estos tres eventos de carcter ssmico. Los tres mecanismos focales
indican un movimiento de comportamiento inverso y con rumbo Norte-Sur.
Fig. 5.17 Mecanismos focales en Los Bronces.
72
En El Teniente se obtuvieron mecanismos focales de 17 eventos. Se aprecia un

comportamiento general de carcter compresional, siempre cercano a un rumbo
Norte-Sur.
Fig. 5.18 Mecanismos focales en El Teniente.
73
Capitulo 6
Interpretacin y Discusin
6.1 Modelos 1D
Para tener una referencia sobre el modelo unidimensional de velocidad obtenido

en Los Bronces, ste puede ser comparado con modelos realizados anteriormente
en base a eventos en esta misma zona (Fig. 6.1). El nuevo modelo se acerca a
grandes rasgos, a modelos pasados, aunque se pueden apreciar diferencias que
pueden asociarse a variadas causas.
En primer lugar el nuevo modelo es el nico que realiza una correccin por
altura, esta correccin puede jugar un rol no menor si se tiene en cuenta que la
diferencia de altura entre estaciones sismolgicas es en muchos casos superior a
los mil metros. Al no considerar estas diferencias de cota, los tiempos de llegada y,
por ende, la velocidad otorgada a las diferentes capas del modelo pueden entregar
distorsiones que finalicen con interpretaciones erradas. Otra razn de diferencias es
la distribucin geogrfica de estaciones sismolgicas; en el caso del modelo de Vera
(2006), muchas de las estaciones estaban en Argentina, en el trasarco, por lo tanto
el modelo nos muestra principalmente la estructura de velocidades bajo la cordillera
y alcanzando parte de Argentina; en el nuevo modelo estn prcticamente todas las
estaciones en el antearco, o sea, del lado poniente de la cordillera andina, por lo que
el sector representando es principalmente el chileno precordillerano. La siguiente
razn est adosada a la creacin del nuevo modelo, pues los ltimos 10 km
pudieron ser modelados gracias a las 2 estaciones ms lejanas al origen de los
eventos ubicadas aproximadamente a 180 km, las otras estaciones ms cercanas
estaban a 140 km, por lo que existen cerca de 40 km en superficie sin recepcin de
datos, esto disminuye la confianza del modelo en las zonas ms profundas.
El modelo 1D final para Los Bronces alcanza los 30.1 Km de profundidad a partir
de la superficie. Se pueden apreciar 3 capas. La existencia de la primera capa se
arrastra del modelo de Barrientos (2004) y delata su composicin grantica, luego en
la segunda capa se puede apreciar que a mayor profundidad hay un aumento del
valor Vp, al igual que con el valor del Vp/Vs, esto indicara una composicin ms
mfica adems de la aparicin de metamorfismo debido al aumento de presin y
temperatura.
En resumen el modelo 1D en Los Bronces es un buen modelo que representa en
buena parte la zona de estudio, y aunque la zona ms profunda que alcanza la
modelacin no es totalmente confiable, si lo es el resultado general en gran medida.
Es necesario adems destacar la importancia de la correccin de altura en zonas
topogrficas de grandes diferencias como ocurre en este caso.
74
Fig. 6.1 Distintos modelos 1D de Los Bronces (onda P, onda S y Vp/Vs).

Todos estos modelos fueron hechos en base a datos provenientes del sector de Los Bronces.
La lnea punteada de la ltima figura representa el modelo de Barrientos (2004),
el cual es tomado de base para los kilmetros ms superficiales del nuevo modelo.
En el caso de El Teniente no hay modelos anteriores basados en eventos de

este sector, salvo el de Godoy (1999) que es para los primeros kilmetros (Fig. 6.2).
Aunque las consideraciones tomadas para El Teniente son las mismas de Los
Bronces: la correccin de altura y la buena distribucin de las estaciones; la alta
dispersin de datos distancia-tiempo (Fig. 4.14 y Fig. 4.15) provoca que el modelo
obtenido, sin ser equvoco, pierda confiabilidad. Esta dispersin est ligada a la
eleccin de un punto como origen de todos los eventos, sabiendo que cada uno de
ellos ocurre a distintas profundidades en torno a la falla principal (TFZ).
Lamentablemente no es posible remediar este problema.
El modelo 1D final para El Teniente alcanza los 26.31 Km de profundidad a
partir de la superficie. Se pueden apreciar 2 capas. La primera capa es heredada del
modelo de Godoy (1999) y sus velocidades son levemente menores a las de Los
Bronces, lo que indica algn cambio en la composicin igualmente grantica. El
Vp/Vs aumenta a medida que se avanza en profundidad debido al aumento de
presin y temperatura y a una composicin cada vez ms mfica.
75
Fig. 6.2 Comparacin de modelos 1D de El Teniente Godoy (1999) y el actual (2012).

El modelo de Godoy es representado por la lnea punteada, el cual es tomado de base para los kilmetros
ms superficiales del nuevo modelo en El Teniente.
En ambos modelos 1D, para alcanzar mayores profundidades para su anlisis

es necesario colocar estaciones en lugares ms distantes al origen de lo que aqu se
presenta.
6.2 Modelo 2D
El modelo bidimensional es de una mayor complejidad y provee importante

informacin para la zona de estudio. La seleccin de estaciones cercanas al perfil
existente entre Los Bronces y El Teniente permite una mayor precisin. Aunque hay
que considerar la no utilizacin en la modelacin final de los registros pertenecientes
a El Teniente debido a su alta dispersin y poca coincidencia con la simulacin. En
las figuras con la modelacin final (Fig. 5.3 y Fig. 5.4) se puede apreciar que las
velocidades tanto de la onda P como de la onda S aumentan con la profundidad,
adems existe un salto de velocidades cerca de la cota cero, heredado de los
modelos unidimensionales. Los valores de velocidades en los primeros kilmetros
(primera capa) son influenciados por la correccin de altura lo que no hace confiable
el detalle, las tendencias y los contrastes que se observan, por lo tanto nos debe
realizarse una lectura de cuerpos o estructuras en estos primeros kilmetros.
El proceso de modelacin 2D ayuda tambin a constatar la importancia de la
correccin de altura, pues se puede notar su efecto al comparar la superposicin de
datos, con y sin esta correccin, con las simulaciones tericas que incluyen la
topografa. Se puede considerar esta correccin de gran utilidad en algn proyecto
donde se requiera mayor detalle en la modelacin bajo una zona de gran variacin
topogrfica.
76
Para comprender de mejor manera la estructura bajo el perfil se grafica

adems el Vp/Vs (Fig. 6.3) y las variaciones porcentuales (anomalas) de las
velocidades de ondas P y S con respecto al modelo 2D inicial (Fig. 6.4 y Fig. 6.5). La
anomala calculada en cada punto del modelo est descrita en la siguiente frmula:
anomala = (velocidad modelo final - velocidad modelo inicial) x 100

.
velocidad modelo inicial
Lo primero que se puede observar en los modelos es un crecimiento del valor

Vp/Vs al aumentar la profundidad, al igual que en los modelos unidimensionales.
Aunque el aumento es mucho mayor hacia el norte, cerca de Los Bronces. Esto
puede ocurrir por la existencia en esa zona de roca ms hidratada, tal vez con una
trayectoria ascendente dirigida hacia la zona del yacimiento minero.
Lamentablemente al ser una modelacin en base a capas planas es difcil detectar
claramente cuerpos con contraste de velocidad ms precisos en sus formas. Por los
valores de Vp, Vs y su cuociente. Se infiere por los valores del Vp/Vs que la roca es
principalmente grantica.
En cuanto a las figuras de anomalas, se observa en ambas figuras, un posible
artefacto bajo la mina El Teniente.
Fig. 6.3 Vp/Vs de modelo final 2D.
77
Fig. 6.4 Anomala de velocidades de Onda P con respecto al modelo 2D inicial.
Fig. 6.5 Anomala de velocidades de Onda S con respecto al modelo 2D inicial.
78
El modelo bidimensional puede compararse a la tomografa realizada por

Pardo (2009) donde se puede notar una concordancia general entre las grficas
Vp/Vs (Fig. 6.6). El modelo obtenido en este trabajo coincide en rasgos generales
con los resultados de este estudio que da cuenta de una posible migracin de
magma hacia la superficie, alcanzando a los yacimientos que es donde seguramente
la corteza posee zonas ms dbiles y fracturadas.
Fig. 6.6 Tomografa. Vp/Vs. Pardo (2009).

Corte Norte-Sur a una longitud de 70.3W.
Lamentablemente el no tener datos en El Teniente que provengan de un

nico origen, perjudica la precisin de la modelacin 2D. Debido adems a la
profundidad alcanzada por el modelo (26 km aproximadamente), no es posible
discernir sobre el origen y caractersticas geodinmicas del nivel ms profundo,
donde podra hallarse algn cuerpo en particular o alguna tendencia que se
mantenga en profundidad, pero al no tener mayor informacin slo se puede
conjeturar apoyado en otros estudios como la tomografa ya mencionada.
Es esencial mejorar el programa RayTracing para realizar modelos en 2D. Los
problemas que presenta son pequeos y corresponde a problemas con el ingreso de
datos y la interfaz grfica, stas limitaciones no debieran ser difciles de corregir,
pero deben realizarse si se quiere nuevamente trabajar de esta forma ya que son
estas trabas las que provocan que el trabajo se torne lento y engorroso en muchos
casos.
79
6.3 Sismicidad en Los Bronces y El Teniente
Las concentraciones de sismicidad ubicadas bajo las minas Los Bronces y El

Teniente estn asociadas principalmente a las mismas fallas que dan origen a los
yacimientos. Los estudios hasta ahora realizados indican que las labores de
extraccin de las minas slo gatillaran microfracturas y microsismos (del orden de 1
o 2 grados de magnitud Richter) y nicamente al interior de las zonas de trabajo.
Entender el origen y los mecanismos de la sismicidad existente en estos lugares
debieran ayudar a aclarar esta interrogante; pero no es mucho lo que se pude decir
respecto al tema con los datos y resultados aqu presentados, esto debido en mayor
medida al no poder determinar con precisin los hipocentros.
La relocalizacin realizada tanto en Los Bronces como en El Teniente gracias
a los nuevos modelos 1D y 2D no arrojan mayores pistas sobre la estructura de
estos lugares. En Los Bronces asumimos desde un principio que la mayora de los
eventos son tronaduras, as que tampoco se esperaba hallar mayor informacin con
la relocalizacin o con el clculo de mecanismos focales; aunque pudieron
obtenerse los mecanismos de 3 eventos (Fig. 5.17), estos marcaban un fallamiento
inverso con rumbo de tendencia Norte-Sur.
Con la relocalizacin en El Teniente podemos divisar apenas una tendencia a
la alineacin vertical (Fig. 5.12 , Fig. 5.13, Fig. 5.15 y Fig. 5.16), que coincidira con
la verticalidad de la falla principal datada para el sector minero, pero esta visin no
es tan evidente y slo confirmara la mejora del modelo. Se obtuvieron mecanismos
focales de 17 eventos, los cuales tambin siguen una tendencia de rumbo Norte-Sur
y responden, en su mayora, a un rgimen compresivo, lo que concuerda con las
caractersticas de esfuerzos a escala regional en la zona [McKinnon, S. 2003].
Adems de las concentraciones de sismicidad halladas bajo las minas,
existen otros eventos, a lo largo de la franja ligada al perfil, registrados por el
conjunto de estaciones que forman la red utilizada en este trabajo. Estos eventos no
pueden relocalizarse, pues para el clculo habra que utilizar nuevamente las
mismas estaciones, stas poseen una distribucin geogrfica extensa que difiere del
sector representado por el perfil, por lo que el nuevo modelo perdera validez. Y si
slo se utilizaran las estaciones confinadas a la franja, existira menos confiabilidad
de la relocalizacin al existir clculos en base a una menor cantidad de datos y sin
un registro panormico.
80
Capitulo 7
Conclusiones
Mediante una red de estaciones sismolgicas, temporales y permanentes, se

han obtenido modelos unidimensionales de estructura de velocidades en el sector de
Los Bronces y El Teniente en base a eventos registrados en estas zonas mineras.
En el caso de Los Bronces, los eventos corresponden a tronaduras asociadas al
proceso productivo, y normalmente se llevan a cabo cercanas a la superficie. El
anlisis de estos datos produce un modelo hasta 30 km de profundidad, conformado
por 3 capas de 4.1, 22.1 y 3.9 km de espesor y velocidades entre 4.8 - 5.1 km/s, 5.9
- 7.0 km/s y 7.1 - 7.2 km/s para la onda P, y 2.8 km/s, 3.5 - 3.9 km/s y 3.9 - 4.0 km/s
para la onda S. En el caso de El Teniente existe una gran dispersin de datos
distancia-tiempo (X-T) producto del carcter ssmico de los eventos que poseen
variados hipocentros. que no poseen un origen nico, hallndose sus hipocentros a
distintas profundidades. El modelo alcanza en este caso una profundidad cercana a
los 26 km, con 2 capas de 2.4 y 23.9 km de espesor y velocidades entre 4.4 - 4.5
km/s y 5.8 - 7.0 km/s para la onda P y 2.5 - 2.6 km/s y 3.4 - 3.9 km/s para la onda S.
Estos modelos se asemejan en general a modelos realizados en anteriores estudios
y las diferencias estn relacionadas principalmente a la correccin de altura y a la
distribucin de las estaciones en el rea de estudio del nuevo modelo.
Tomando como base los modelos unidimensionales de Los Bronces y El
Teniente se cre un modelo bidimensional entre estas dos zonas mineras, creando
un perfil Norte-Sur de 132.8 Km de longitud ubicado en el sector precordillerano.
Para perfeccionar el ajuste del modelo se utilizaron los datos de estaciones
instaladas dentro de una franja de 24 Km en torno al perfil y se consideraron
finalmente slo los eventos provenientes de Los Bronces debido a la gran dispersin
de los eventos de El Teniente.
Para la modelacin de los datos (X-T) se ha usado un programa que traza
rayos para un modelo bidimensional de velocidades y que permite comparar
directamente los tiempos de llegada modelados con aquellos correspondientes a los
datos reales. Variando el modelo de velocidades se mejora iterativamente el ajuste a
los datos reales y con ello, el modelo mismo. El modelo de velocidades considera un
conjunto de capas, cada una de las cuales se compone a su vez de una serie de
bloques trapezoidales donde existe una variacin bilineal de velocidad. Se puede
ajustar iterativamente las caractersticas de estos bloques permitiendo, entre otras
ventajas, considerar la topografa del perfil. El modelo bidimensional obtenido
alcanza profundidades cercanas a los 25 Km. Para interpretar el modelo no deben
considerarse los kilmetros sobre la cota cero, pues arrastran problemas
relacionados a la topografa al interpolar linealmente los modelos unidimensionales
en un comienzo, adems de estar afectado por la correccin de altura. Tambin
existe imprecisin en los 5 kilmetros ms profundos, pues estn basados en slo 2
estaciones lejanas que extienden el modelo en profundidad.
Observando las grficas de Vp/Vs se puede concluir un mayor aumento en
profundidad de los valores de Vp/Vs en el extremo Norte del perfil, bajo Los Bronces,
lo que indicara la existencia de fases fludas (roca hidratada o magma) hacia el
norte y roca ms seca y rgida hacia el Sur del perfil, bajo El Teniente. En las figuras
81
sobre anomalas con respecto al modelo 2D inicial se observa una zona bajo la mina
El Teniente que podra ser un artefacto que coincidiera con un sector ms rgido de
la corteza. Al comparar la grfica Vp/Vs de este modelo con una tomografa del 2009
donde se aprecia de manera ms detallada la existencia de cuerpos, se observa una
concordancia general entre ambas representaciones lo que confirmara tendencias y
estructuras en lneas muy generales.
Es importante mencionar la relevancia de la correccin por altura realizada en
los modelos 1D y 2D, especialmente en esta zona de estudio con significativos
cambios topogrficos, los datos utilizados provienen de estaciones ubicadas tanto en
zonas costeras como en cordilleranas, con diferencias de cota de hasta
aproximadamente 3 mil metros.
Se realiz adems la relocalizacin de los eventos concentrados en las zonas
mineras de Los Bronces y El Teniente con los nuevos modelos 1D y 2D. Las
variaciones de localizacin fueron poco significativas y no se distingue caracterstica
geomtrica alguna que pueda ayudar a una mejor interpretacin de la sismicidad en
estos sectores. En cuanto a los esfuerzos, en Los Bronces se lograron obtener
mecanismos focales de 3 eventos, todos con fallamiento inverso de rumbo NorteSur. En El Teniente se obtuvieron 17 mecanismos focales, la mayora
compresionales y tambin con rumbos aproximados Norte-Sur.
82
Bibliografa
Acevedo, P.; 1985. Estructura cortical y estudio sismotectnico de Chile
central entre las latitudes 32.08S y 34.58S. Tesis para optar al grado de Magister,
Facultad de Ciencias Fsicas y Matemticas, Universidad de Chile, 172pp.
Alvarado, P.; 1998. Sismicidad superficial de los Andes Centrales. Tesis para
optar al ttulo de Magister, Facultad de Ciencias Fsicas y Matemticas, Universidad
de Chile, 161 pp.
Alvarado, P.; Barrientos, S.; Saez, M.; Astroza, M.; Beck, S.; 2009. Source
study and tectonic implications of the historic 1958 Las Melosas cristal earthquake,
Chile, compared to earthquake damage. Physics of the Earth and Planetary Interiors.
175, 26 36.
Barrientos, S.; Eisenberg, A. 1988. Secuencia ssmica en la zona cordillerana
al interior de Rancagua. V Congreso Geolgico Chileno, Santiago, F121-F132.
Barrientos, S.; Vera, E.; Alvarado, P.; Monfret, T. 2004. Crustal seismicity in
central Chile. Journal of South American Earth Sciences. 16, 759 - 768.
Cahill, T.; B. Isacks, 1992. Seismicity and shape of the subducted Nazca
plate. Journal. Geophysical Research, 97, 17503 - 17529.
Christensen, N.; 1996. Poissons ratio and cristal seismology. Journal of
Geophysical Research, 101, 3139 3156.
Contreras, E.; 2003. Caracterizacin ssmica del monte submarino OHiggins.
Tesis para optar al ttulo de Magister, Facultad de Ciencias Fsicas y Matemticas,
Universidad de Chile, 118 pp.
DeMets, C.; Gordon, R.G.; Argus, D.F.; Stein, S., 1994. Effect of recent
revisions to the geomagnetic reversal time scale on estimates of current plate
motions., Geophys. Res. Lett., 21, 2191 2194.
Duarte, P.; 2000. Caracterizacin geolgica y geotectnica del prfido dactico
El Teniente. Memoria para optar al ttulo de Gelogo, Facultad de Ciencias Fsicas y
Matemticas, Universidad de Chile, 81pp.
Fuenzalida, A.; Pardo, M.; Cisternas, A.; Dorbath, L.; Dorbath, C.; Comte, D.;
Kausel, E. 1992. On the geometry of the Nazca plate subducted under central Chile
(32-34.5S) as inferred from microseismic data. Tectonophysics. 205, 1 - 11.
Garrido, I., Cembrano, J., Sia, A., Stedman, P., Yaez, G., 2002. High
magma oxidation state and bulk crustal shortening: key factors in the genesis of
Andean prophyry copper deposits, central Chile (31 34 S). Revista Geolgica de
Chile Vol. 29, N 1, p. 43 - 54.
83
Garrido, I., Riveros, M., Cladouhos, T., Espieira, D., Allmendinger, R., 1994.
Modelo geolgico estructural yacimiento El Teniente: 78 Congreso Geolgico
Chileno, Concepcin. v. II, pp. 1553 - 1558.
Godoy, E., Yaez, G., Vera, E., 1999. Inversion of an Oligocene volcanotectonic basin and uplifting of its superimposed Miocene magmatic arc in the Chilean
Central Andes: first seismic and gravity evidences. Tectonophysics. 306, 217 - 236.
Gutscher, M.A.; Spakman, W.; Bijwaard, H.; Engdahl, E.R. 2000.
Geodynamics of the flat subduction:Seismicity and tomographic constraint from the
Andean margin. Tectonics. 19: 814 - 833.
Hintze, F.; 1978. Geologa, zonacin, mineralizacin nivel pommerantz mina
Los Bronces de la Compaa Minera Disputada de las Condes S.A. Area
metropolitana. Memoria para optar al ttulo de Gelogo, Facultad de Ciencias Fsicas
y Matemticas, Universidad de Chile.
Kreemer, C.; Holt, W.E.; Haines, A.J., 2003. An integrated global model of
present-day plate motions and plate boundary deformation, Geophys. J. Int. 154, 8 34.
Lavenu, A., Cembrano, J., 1999. Compressional and transpressional stress
pattern for Pliocene and Quaternary brittle deformation in fore arc and intra-arc zones
(Andes of Central and Southern Chile). Journal of Structural Geology. 21, 1669 1691.
Lomnitz, C.1961. Los terremotos del 4 de septiembre en el Cajn del Maipo.
Anales de la Facultad de Ciencias Fsicas y Matemticas Publicacin 3, Depto. De
Geofsica, 18 pp. 279 - 306.
Mapa geolgico de Chile: Versin digital. Publicacin geolgica digital n4,
2003. Base geolgica escala 1:1.000.000. Servicio Nacional de Geologa y Minera.
Subdireccin nacional de Geologa. Gobierno de Chile.
McKinnon, S; Garrido de La Barra, I. 2003. Stress field analysis at the El
Teniente Mine: evidence for N S compression in the modern Andes. Journal of
Structural Geology. 25, 2125 - 2139.
Pardo, M.; Comte, D.; Monfret, T. 2002. Seismotectonic and stress distribution
in the Central Chile subduction zone. J. South Am. Earth Sciences. 15, 11 - 22.
Pardo, M.; Monfret, T.; Vera, E.; Eisenberg, A.; Gaffet, S.; Lorca, E.; Prez, A.
2002. Flat-Slab subduction zone in Central Chile-Argentina: Seismotectonic and
body-wave tomography from local data. ISAG, Fifth International Symposium on
Andean Geodynamics. Toulouse, France, September 16-18, 2002. Extended
Abstracts, ISAG Andean Geodynamics. 469 - 472.
Pardo, M.; Vera, E.; Yaez, G.; Monfret, T. 2009. Tomografa ssmica bajo los
Andes de Chile Central (33 - 34.5S): Implicaciones sismotectnicas. XII Congreso
Geolgico chileno. Santiago 22-26 Noviembre, 2009.
Piderit, E. 1961. Estudio de los sismos del Cajn del Maipo el ao 1958.
Memoria para optar al ttulo de Ingeniero Civil, Facultad de Ciencias Fsicas y
Matemticas, Universidad de Chile, 125pp.
84
Vera, E.; Pardo, M.; Monfret, T.; Eisenberg, A.; Yaez, G.; Triep, E. 2006.
Eventos ssmicos corticales en los Andes centrales de Chile y Argentina. Congreso
Geolgico Chileno 2006.
Vera, E.; Pardo, M.; Comte, D.; Yaez, G. 2009. Estructura ssmica bajo los
Andes de Chile Central: Perfil de Refraccin Los Bronces El Teniente. XII
Congreso Geolgico chileno. Santiago 22-26 Noviembre, 2009.
Waldhauser, F.; Ellsworth, W. 2000. A Double-Difference Earthquake Location
Algorithm: Method and Application to the Northern Hayward Fault, California. Bulletin
of the Seismological Society of America. 90, 6, pp. 13531368.
Zelt, C.A, and Smith, R.B. 1992. Seismic traveltime inversion for 2-D crustal
velocity structure. Geophys. J. Int, 108, pp. 16 - 34.
85
Anexo A
Relacin grfica entre Vm y X y clculo del valor Vm
Figuras donde se aprecia la relacin de Vm y X para realizar la correccin de
altura en Los Bronces y El Teniente.
Fig. 0.1 Trazado de rayos en Los Bronces en base a modelo sin correcciones (Onda P)
Fig. 0.2 Trazado de rayos en Los Bronces en base a modelo sin correcciones (Onda S)
86
Fig. 0.3 Trazado de rayos en El Teniente en base a modelo sin correcciones (Onda P).
Fig. 0.4 Trazado de rayos en El Teniente en base a modelo sin correcciones (Onda S).
87
Para estimar Vm se realizan aproximaciones y supuestos que ayudan a

simplificar el clculo. El valor de Vm sirve para encontrar los valores finales de
distancia y tiempo al realizar la correccin de altura (Captulos 4.1.4 y 4.1.5).
Primero se calcula la distancia (X) de la primera capa.
Con
de la primera capa, pero usando el H hasta la segunda capa
Luego se considera el trayecto total
Utilizando ahora el k de la segunda capa:
88
Anexo B
Rutina en Matlab para la creacin de la figuras de la modelacin 2D.
%
clear all;
%
format short g
load mod2dipB3.dat;
xzv1=mod2dipB3(:,:);
load mod2disB3.dat;
xzv2=mod2disB3(:,:);
load mod2dfpf.dat;
xzv3=mod2dfpf(:,:);
load mod2dfsf.dat;
xzv4=mod2dfsf(:,:);
load polarr2dvpvs.dat;
cp2dvpvs = polarr2dvpvs(:,:);
load seis2dp.dat;
c62dp = seis2dp(:,:);
load seis2ds.dat;
c62ds = seis2ds(:,:);
nnx = 150;
nnz = 150;
[x, z, vpini] = fngrillapiB3 (xzv1, nnx, nnz);
[x, z, vsini] = fngrillasiB3 (xzv2, nnx, nnz);
[x, z, vpfin] = fngrillapsf (xzv3, nnx, nnz);
[x, z, vsfin] = fngrillapsf (xzv4, nnx, nnz);
xp =xzv1(:,1);
zp =xzv1(:,2);
%%% paleta de colores
cp2dvpvs = cp2dvpvs/255;
c62dp = c62dp/255;
c62ds = c62ds/255;
89
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo inicial 2D. Onda P.
figure;
colormap (c62dp);
fill (x,z,vpini);
view(0,270);
colorbar;
shading interp;
title ('Modelo inicial 2D. Onda P.','FontSize',18)
xlabel('Distancia','FontSize',16);
ylabel('Profundidad','FontSize',16);
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo inicial 2D. Onda S.
figure;
colormap (c62ds);
fill (x,z,vsini);
view(0,270);
colorbar;
shading interp;
title ('Modelo inicial 2D. Onda S.','FontSize',18)
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo final 2D. Onda P.
figure;
colormap (c62dp);
fill (x,z,vpfin);
view(0,270);
colorbar;
shading interp;
title ('Modelo final 2D. Onda P.','FontSize',18)
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo final 2D. Onda S.
figure;
colormap (c62ds);
fill (x,z,vsfin);
view(0,270);
colorbar;
shading interp;
title ('Modelo final 2D. Onda S.','FontSize',18)
90
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
vpvsini = vpini./vsini;
vpvsfin = vpfin./vsfin;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo inicial 2D. vp/vs
figure;
colormap (cp2dvpvs);
fill (x,z,vpvsini);
caxis([1.65, 1.85]);
view(0,270);
colorbar;
shading interp;
title ('Modelo inicial 2D. Vp/Vs','FontSize',18)
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Modelo final 2D. vp/vs
figure;
colormap (cp2dvpvs);
fill (x,z,vpvsfin);
caxis([1.65, 1.85]);
view(0,270);
colorbar;
shading interp;
title ('Modelo final 2D. Vp/Vs','FontSize',18)
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
anomp=((vpfin-vpini)./vpini)*100;
anoms=((vsfin-vsini)./vsini)*100;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% modelo de anomalia en p
figure;
colormap (c62dp);
fill (x,z,anomp);
caxis([-6, 6]);
view(0,270);
colorbar;
91
shading interp;
title ('Anomala de velocidades. Onda P.','FontSize',16)
grid on;
hold on;
plot(xp,zp,'k-');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% modelo de anomalia en s
figure;
colormap (c62ds);
fill (x,z,anoms);
caxis([-6, 6]);
view(0,270);
colorbar;
shading interp;
title ('Anomala de velocidades. Onda S.','FontSize',16)
grid on;
hold on;
plot(xp,zp,'k-');
Funcin solicitada para crear la grilla y dar los valores a los nodos.
%%%%%%%%%%%%%%%
Funcin para crear grilla
%%%%%%%%%%%
function [xx,zz,vv] = fngrillapsf (modv, nnx, nnz)

n=length(modv(:,1));
xmax=max(modv(:,1));
xmin=min(modv(:,1));
nx = linspace (xmin, xmax, nnx); %en km
zmax=max(modv(:,19));
zmin=min(modv(:,2));
nz = linspace (zmin, zmax, nnz); %en km
nn = (nnx)*(nnz); % nro de nodos total de la grilla
nv = zeros(nnx,nnz); %matriz de velocidades inicialmente con zeros
v1izq
v1der
v2izq
v2der
=
=
=
=
zeros(nnx,nnz);
zeros(nnx,nnz);
zeros(nnx,nnz);
zeros(nnx,nnz);
z1izq
z2izq
z1der
z2der
=
=
=
=
zeros(nnx,nnz);
zeros(nnx,nnz);
zeros(nnx,nnz);
zeros(nnx,nnz);
for s = 1 : nnz
for m = 1 : nnx
92
q=1;
% q indica la columna o bloque vertical (n-1)

donde
se encuentra el nodo
while nx(m) >= modv(q,1) & q < n

q =q+1;
end
%%% profundidad de las interfaces
ztop(m,s) = modv(q-1,2) + ((modv(q,2)-modv(q1,2))/(modv(q,1)-modv(q-1,1)))*nx(m) - ((modv(q,2)-modv(q-1,2))/(modv(q,1)modv(q-1,1)))*modv(q-1,1);
zc1c2(m,s) = modv(q-1,4) + ((modv(q,4)-modv(q1,4))/(modv(q,1)-modv(q-1,1)))*nx(m) - ((modv(q,4)-modv(q-1,4))/(modv(q,1)modv(q-1,1)))*modv(q-1,1);
zc2c3(m,s) = modv(q-1,7) + ((modv(q,7)-modv(q1,7))/(modv(q,1)-modv(q-1,1)))*nx(m) - ((modv(q,7)-modv(q-1,7))/(modv(q,1)modv(q-1,1)))*modv(q-1,1);
zc3c4(m,s) = modv(q-1,10) + ((modv(q,10)-modv(q1,10))/(modv(q,1)-modv(q-1,1)))*nx(m) - ((modv(q,10)-modv(q1,10))/(modv(q,1)-modv(q-1,1)))*modv(q-1,1);
%zc7c8(m,s) = modv(q-1,22) + ((modv(q,22)-modv(q1,22))/(modv(q,1)-modv(q-1,1)))*nx(m) - ((modv(q,22)-modv(q1,22))/(modv(q,1)-modv(q-1,1)))*modv(q-1,1);
%%%%%%%%%% velocidades %%%%%%%%%%
if nz(s) >= ztop(m,s) && nz(s) < zc1c2(m,s)
elseif
v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
=
=
=
=
modv(q-1,3);
modv(q,3);
modv(q-1,5);
modv(q,5);
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
=
=
=
=
modv(q-1,2);
modv(q-1,4);
modv(q,2);
modv(q,4);
nz(s) >= zc1c2(m,s) && nz(s) < zc2c3(m,s)

v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
=
=
=
=
modv(q-1,6);
modv(q,6);
modv(q-1,8);
modv(q,8);
93
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
elseif
elseif
elseif
modv(q-1,4);
modv(q-1,7);
modv(q,4);
modv(q,7);

v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
=
=
=
=
modv(q-1,9);
modv(q,9);
modv(q-1,11);
modv(q,11);
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
=
=
=
=
modv(q-1,7);
modv(q-1,10);
modv(q,7);
modv(q,10);

v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
elseif
=
=
=
=
=
=
=
=
=
=
=
=
modv(q-1,12);
modv(q,12);
modv(q-1,14);
modv(q,14);
modv(q-1,10);
modv(q-1,13);
modv(q,10);
modv(q,13);

v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
=
=
=
=
modv(q-1,15);
modv(q,15);
modv(q-1,17);
modv(q,17);
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
=
=
=
=
modv(q-1,13);
modv(q-1,16);
modv(q,13);
modv(q,16);
nz(s) >= zc5c6(m,s) && nz(s) <= modv(q-1,19)

v1izq(m,s)
v1der(m,s)
v2izq(m,s)
v2der(m,s)
=
=
=
=
modv(q-1,18);
modv(q,18);
modv(q-1,20);
modv(q,20);
z1izq(m,s)
z2izq(m,s)
z1der(m,s)
z2der(m,s)
=
=
=
=
modv(q-1,16);
modv(q-1,19);
modv(q,16);
modv(q,19);
end
% agregar los valores de v a cada nodo, esto en fn de los trapecios
% grandes definidos en xzv
va(m,s) = v1der(m,s)*((nx(m)-modv(q-1,1))/(modv(q,1)modv(q-1,1))) + v1izq(m,s)*((modv(q,1)-nx(m))/(modv(q,1)-modv(q-1,1)));
94
vb(m,s) = v2der(m,s)*((nx(m)-modv(q-1,1))/(modv(q,1)modv(q-1,1))) + v2izq(m,s)*((modv(q,1)-nx(m))/(modv(q,1)-modv(q-1,1)));

za(m,s) = z1der(m,s)*((nx(m)-modv(q-1,1))/(modv(q,1)modv(q-1,1))) + z1izq(m,s)*((modv(q,1)-nx(m))/(modv(q,1)-modv(q-1,1)));
zb(m,s) = z2der(m,s)*((nx(m)-modv(q-1,1))/(modv(q,1)modv(q-1,1))) + z2izq(m,s)*((modv(q,1)-nx(m))/(modv(q,1)-modv(q-1,1)));
nv(m,s) = va(m,s)*((zb(m,s)-nz(s))/(zb(m,s)-za(m,s))) +
vb(m,s)*((nz(s)-za(m,s))/(zb(m,s)-za(m,s)));
if nv (m,s) == 0
nv (m,s) =nan;
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%GRILLA PARA COLOREAR
nnn = (nnx-1)*(nnz-1);
u=1;
for w = 1 : (nnz-1)
for t = 1 : (nnx-1)
if u <= nnn
xx (1,u) = nx(t);
zz (1,u) = nz(w);
vv (1,u) = nv(t,w);
xx (2,u) = nx(t+1);
zz (2,u) = nz(w);
vv (2,u) = nv(t+1,w);
xx (3,u) = nx(t+1);
zz (3,u) = nz(w+1);
vv (3,u) = nv(t+1,w+1);
xx (4,u) = nx(t);
zz (4,u) = nz(w+1);
vv (4,u) = nv(t,w+1
end
u=u+1;
end
end
95

CF Villegas Ls

Cargado por

Copyright:

Formatos disponibles

CF Villegas Ls

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

CF Villegas Ls

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD DE CHILE

FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS

ESTRUCTURA SISMICA CORTICAL EN LOS ANDES CENTRALES

TESIS PARA OPTAR AL TTULO DE

Los procesos sismotectnicos de la Zona Central de Chile estn relacionados con

En primer lugar agradezco a mis padres, quienes priorizaron durante mi infancia

Mina Los Bronces .............................................................................................2

Mina El Teniente ................................................................................................3

Marco terico ...........................................................................................5

Teora de rayos ..................................................................................................5

Red Sismolgica ...............................................................................................8

Procesamiento de datos .......................................................................12

Modelo de velocidades en Los Bronces .....................................................12

Eleccin del Perfil .......................................................................................28

Mecanismos focales .......................................................................................58

Modelo de Velocidades 1-D Los Bronces .....................................................59

Modelo de Velocidades 1-D El Teniente ........................................................60

Modelo de Velocidades 2-D (Perfiles) ............................................................61

Mecanismos focales .......................................................................................72

Interpretacin y Discusin ...................................................................74

Sismicidad en Los Bronces y El Teniente.....................................................80

Relacin grfica entre Vm y X y clculo del valor Vm ..............................................86

Rutina en Matlab para la creacin de la figuras de la modelacin 2D. ....................89

Fig. 4.26 Datos Distancia-Tiempo con correccin de altura () .................................................... 35

Fig. 5.3 Modelo final 2D (onda P). ............................................................................................... 62

Los procesos sismotectnicos de la Zona Central de Chile estn relacionados con

Todo esto se obtiene gracias al registro de eventos en la zona de estudio

1.2 Mina Los Bronces

La compaa minera Anglo American Chile posee yacimientos en la zona

1.3 Mina El Teniente

El Teniente es un yacimiento de tipo prfido cuprfero localizado en la alta

Fig. 1.1 Zona de Falla El Teniente (ZFT).

El objetivo general de este trabajo es obtener un modelo de estructura de

Determinar un modelo cortical 1D de la estructura de velocidades en la zona de

Determinar un modelo cortical 2D de la estructura de velocidades en la zona de

Refinar la localizacin de hipocentros con la nueva estructura de velocidades,

2.1 Teora de rayos

La teora de rayos es muy til para entender y obtener resultados de manera

Fig. 2.1 Geometra del rayo atravesando capa horizontal.

Por geometra simple (Fig. 2.1) podemos extraer las relaciones:

Obtenemos entonces 2 ecuaciones fundamentales para llevar a cabo el clculo

Es de gran utilidad, para la modelacin e interpretacin de estructuras geolgicas

Fig. 2.2 Gradiente de velocidad en una capa.

V representa la velocidad del medio en el punto ms profundo del rayo, el factor

3.1 Red Sismolgica

Los datos fueron proporcionados por 28 estaciones sismolgicas temporales de

Fig. 3.1 Red sismolgica utilizada. Perfil W-E.

Fig. 3.2 Zona de estudio (33 - 34.5S).

3.2 Los Bronces

Los datos provenientes de la mina Los Bronces estn asociados principalmente a

Los datos provenientes de la zona de El Teniente fueron obtenidos seleccionando

Es posible obtener un modelo de estructura ssmica unidimensional (velocidad en

4.1.1 Modelo de velocidades en Los Bronces

Se tienen los registros de los tiempos de llegada de las ondas P y S en cada

Fig. 4.1 Vista en planta de los eventos registrados en Los Bronces.

Fig. 4.2 Eventos registrados en Los Bronces.

Calculando la regresin lineal de la distancia versus los tiempos de llegada

Fig. 4.3 Clculo del tiempo mediante regresin lineal.

Al fijar un tiempo de referencia en la estacin NEVA, se pueden ajustar y corregir

VP1 = [(AB P1A)/AB]V1A + [P1A/AB]V2A