-
3.1 Localización Geográfica
El proyecto carretero Tuxpan Tampico tiene proyectado un trazo de 105 Km, une la
ciudad de Tampico, Tamaulipas, con la autopista México Tuxpan; siendo este el
único acceso a la carretera hacia naranjos y de ahí al tramo en construcción.
El área de estudio para este informe consiste en un tramo de casi 19 kilómetros con
una orientación NW-SE y se encuentra ubicada al este de los municipios de Naranjos
de Amatlán y Cerro Azul; en el estado de Veracruz, México. Estos municipios se
localizan al norte del estado. De manera general, la zona de estudio se encuentra en
parte de la Huasteca alta (Municipio de Naranjos) y la Huasteca baja (Municipio de
Cerro Azul); limitada al norte por la sierra de Tamaulipas, al sur por el macizo de
Teziutlán, al este por el golfo de México y al oeste por la Madre Sierra Oriental.
Los 17 tendidos de refracción sísmica tienen una longitud de 128 metros y se
distribuyeron a lo largo de los 19 km abarcados por el proyecto. (Figura 3.1)
Figura 3.1 Localización de los perfiles sísmicos adquiridos. (Imagen tomada de Google Earth)
10
3.2 Geología General
El área de estudio se encuentra dentro de la carta geológica-minera Tamiahua F14-9,
ubicada en la porción norte del estado de Veracruz. Fisiográficamente se ubica en la
cuenca Tampico Misantla y en la llanura del Costera del Golfo de México. (Raisz,
1964). La cuenca sedimentaria de Tampico Misantla está delimitada por tres
unidades estructurales primarias: la parte sur de la plataforma de Tamaulipas, la
Antefosa de Chicontepec y la sierra Madre Oriental.
Regionalmente, las unidades sedimentarias son generalmente de tipo arcillo-arenoso-
calcáreo con edades del Terciario que van del Eoceno al Pleistoceno; representadas
por las formaciones Palma Real (ToAr-Lu) constituida por arenisca-lutita del oligoceno
inferior medio, Mesón (ToAr-Lm) conformada por arenisca - limolita y la Formación
Tuxpan (TmAr-Cgp) constituida por arenisca- conglomerado polimíctico.
De acuerdo con la columna estratigráfica las unidades más antiguas están definidas
por lutita -arenisca del Eoceno medio definida como la Formación Guayabal (TeLu-Ar),
sobreyaciendo a esta unidad se encuentra marga-lutita denominada Formación
Chapopote -Tantoyuca (TeLu-Ar). La unidad más reciente está definida por una unidad
Cuaternaria constituida por aluvión, generado por el propio intemperismo natural de la
zona.
Geológicamente el área de estudio se encuentra asentada en un depósito Terciario
correspondiente a depósitos sedimentarios de arenisca - limolita (ToAr-Lm), al cual
sobreyacen algunos depósitos cuaternarios constituidos por roca sedimentaria aluvión.
La Figura 3.1, muestra la ubicación del área de estudio en la carta geológica así como su
columna estratigráfica correspondiente.
Tramo en estudio
Figura 3.2. Localización de la zona de estudio en la carta geológico-minera Tamiahua F14-9 (Escala 1:250,000).
12
3.3. Zonificación Geotécnica
La adquisición de los tendidos de refracción y el cálculo de sus propiedades elásticas
son la base para el cálculo del espectro de diseño por sismo. El espectro de diseño por
sismo sirve para conocer los parámetros bajo los que se debe diseñar y construir las
estructuras para resistir, como mínimo, los efectos de los movimientos sísmicos del
terrero sobre las que se desplanten. Por lo tanto es importante mencionar e incluir un
poco de información de la zonificación geotécnica de la zona que estudié, como lo
describo a continuación.
La República mexicana se encuentra divida en cuatro regiones sísmicas, esta
zonificación se elaboró para servir de apoyo en el diseño sísmico de estructuras. Dicha
clasificación se realizó en base a los registros históricos de la frecuencia de sismos y la
intensidad de los mismos, los cuales datan de aproximadamente inicios de siglo.
Estas regiones sísmicas reflejan que tan frecuentes son los temblores y cuál es la
aceleración máxima del suelo que se espera por la acción de las fuerzas sísmicas.
Como se mencionó anteriormente, nuestro país se divide en cuatro zonas o regiones
sísmicas las cuales son:
Zona A (Baja sismicidad): En esta zona no se han registrados ningún sismo de
magnitud considerable en los últimos 80 años, ni se esperan aceleraciones del
suelo mayores al 10 % de la aceleración de la gravedad.
Zona B (Media intensidad): Esta zona es de moderada intensidad, pero las
aceleraciones no alcanzan a rebasar el 70% de la aceleración de la gravedad.
Zona C (Alta intensidad): En esta zona hay más actividad sísmica que en la
zona B, aunque las aceleraciones del suelo tampoco sobrepasan el 70% de la
aceleración de la gravedad.
Zona D (Muy alta intensidad): Aquí es donde se han originado los grandes
sísmicos históricos, y la ocurrencia de sismos es muy frecuente, además de que
las aceleraciones del suelo sobrepasan el 70% de la aceleración de la gravedad.
Figura 3.2 Zonificación sísmica (Memoria del Seminario práctico de análisis de los efectos por viento y
sismo, 16 al 18 de junio de 2010, Morelia, Mich., Centro Regional de Desarrollo en Ingeniería Civil. )
De acuerdo al análisis en el software PRODISIS v2.3 (CFE), el tramo en estudio se
localiza dentro de la ZONA B (Véase Figura 2.3), de donde se pudo dividir el área de
estudió en dos partes desde el punto de vista de la aceleración del terreno.
La primera con una aceleración de 49 gals y va del kilometraje 36 al 44, mientras que
en la segunda se tiene una aceleración de 50 gals y abarca los kilometrajes 44 al 53.
14
4.1 Trabajo de campo realizado
Los trabajos de campo, procesamiento e interpretación para cada uno de los tendidos
de refracción los realice con la misma metodología la cual ejemplifico con el perfil
sísmico PS024.
El arreglo sísmico que utilicé para todos los TRS fue diseñado para alcanzar una
profundidad de investigación de 40 metros con una longitud total de 128 m, donde se
utilizaron 12 geófonos, en cuyo tendido se efectuaron tres puntos de tiro de fuente
sísmica artificial generada con un marro de acero de 20 lb, como se muestra en la
Figura 4.1.
PT : Punto de tiro
: Geófonos de 4.5 Hz
Figura. 4.1. Arreglo sísmico utilizado en el proyecto. Los geófonos se encuentran espaciados 8 y 16 m
con sus respectivas fuentes para la generación de ondas sísmicas.
Para realizar esta fase del proyecto se utilizó un sismógrafo que cuenta con un sistema
de registro múltiple marca Geometrics, modelo GEODE, de 24 canales. Este
sismógrafo se conecta directo a una PC por medio de un puerto USB, desplegando en
la PC todas las opciones del arreglo sísmico, y permitiendo introducir todos aquellos
parámetros que se deben controlar durante el levantamiento en campo.
4.2 Procesamiento de Datos Geofísicos de TRS
El procesamiento de los datos obtenidos en campo consistió en el análisis y lectura de
los sismogramas grabados digitalmente, para lo cual se utilizó el programa Tom-Time
(GeoTome, LLC) que permite visualizar el registro sísmico en pantalla y seleccionar los
(Figura 4.2)
Figura. 4.2 Picado de sismogramas para definir tiempo de arribo de la onda P para cada uno de los tres
puntos de tiro del perfil sísmico 024.
Una vez leídos los diferentes tiempos de arribo, se procedió a construir las
dromocrónicas, como se muestra en la Figura 4.3, para determinar el número de capas
sismoelásticas que conforman el subsuelo somero y sus velocidades de propagación
Vp y Vs, correspondientes. Las capas están definidas por los cambio de pendientes en
las curvas.
16
Figura 4.3 Dromocrónicas a partir del procesamiento de los datos del PS024.
4.3 Procesamiento de Datos Geofísicos de MASW1D
El procesamiento consistió en el análisis de los registros obtenidos en campo, haciendo
uso del software GEOPSY (Wathelet, 2006) para extraer las curvas de dispersión de
las ondas Rayleigh, para ello se utiliza el módulo FK-activo, el cual realiza una
transformada en el dominio del espacio para obtener una función dependiente del
número de onda k, posteriormente calcula la velocidad de fase y finalmente despliega
las curvas de dispersión; una para cada punto de excitación. Los resultados de este
proceso para el tiro 1 de la línea ejemplo, se muestran en la Figura 4. 4. En esa figura
se muestran los máximos espectrales asociados a los trenes de ondas superficiales de
PT1. En el diagrama de la misma figura he seleccionado la curva de dispersión. En esa
figura también muestro el límite de confiabilidad, como regla general, se puede derivar
los límites de longitud de onda del arreglo con la abertura seleccionada, mediante
sensor, obteniendo:
La curva de dispersión promedio y/o representativa, es invertida para obtener el perfil
unidimensional de onda de corte.
Figura 4.4 Curva de dispersión obtenida a partir del procesamiento de los datos del PS024. La
línea punteada es el picado de la curva de dispersión, mientras que la línea continua marca el
límite de confiabilidad.
El proceso de inversión tuvo lugar en el módulo Dinver de GEOPSY, el cual genera un
conjunto aleatorio de modelos de velocidad de onda S, con base al modelo obtenido en
el procesamiento del TRS, los cuales se comparan con la curva de dispersión medida
arrojando un valor de misfit (ajuste del error) y el mejor modelo.
Finalmente, determiné los módulos elásticos dinámicos de los estratos que constituyen
el subsuelo; para ello utilicé las relaciones derivadas de la propagación de ondas en
medios continuos, la cual establece que los parámetros elásticos dinámicos de los
materiales son función de las velocidades de propagación de las ondas sísmicas
longitudinal (Vp) y transversal o de cortante (Vs) que se propagan a través de él.
Las ecuaciones que me permitieron determinar estos parámetros elásticos son:
Relación de Poisson Módulo de corte
2
Vp G = Vs 2
0.5 () -1
= Vs
2
Vp Módulo de Young
( ) -1
Vs E = 2 Vs 2(1 + )
Tras el procesamiento y la obtención las velocidades VP y Vs correspondientes a cada
uno de las unidades sismoestratigraficas, determiné el grado de arabilidad del subsuelo
de acuerdo con el Manual de Rendimiento de Caterpillar, Edición 25.
18
Los resultados de este reporte consisten en el modelo 2D de velocidad de onda P y
modelo 1D o perfil de velocidad de la onda S, una tabla de valores de constantes
elásticas, y el dictamen de arabilidad basado en la tabla de arabilidad para un
desgarrador D8R del Manual de Caterpillar (Ver anexo I), que comúnmente se emplea
en la toma de decisiones en lo referente a excavaciones y movimiento de materiales.
En el perfil sísmico, el cual se muestra en la Figura 5.1, observamos que la línea de
exploración tiene una dirección SE-NW; tras su procesamiento e interpretación se
identificaron tres unidades sismoestratigráficas, las cuales se describen a continuación.
La primera unidad se asocia a materiales arcillo arenosos café y arenosos limosos del
mismo color, presenta un espesor medio de 8 metros, tiene una velocidad de onda
longitudinal de 396 m/s y una velocidad de onda de corte de 180 m/s; por lo que se
clasifica como 100% arable.
Subyaciendo de manera concordante se encuentra la segunda unidad, ubicada bajo 8
metros de profundidad, la cual se asocia a arcilla café de consistencia dura, presenta
un espesor medio de 16 metros, tiene una velocidad de onda longitudinal de 2,057 m/s
y una velocidad de onda de corte de 421 m/s; por lo que se clasifica como desgarrable
con bajos rendimientos o marginal.
Finalmente, la tercera unidad que se asocia posiblemente a roca lutita, se presenta a
partir de una profundidad media de 24 m y su espesor es indefinido por medio de esta
técnica, tiene una velocidad de onda longitudinal de 3,177 m/s y una velocidad de onda
de corte de 77 m/s; por lo que se clasifica como no desgarrable o no arable.
Con base en los resultados anteriores, se determinaron los módulos elásticos
dinámicos de los materiales presentes en el subsuelo de este punto, mismos que se
indican en la Tabla 1.
Tabla 1. Propiedades elásticas dinámicas del PS024.
Densidad Relación de Cortante Elasticidad
Capa Vp (m/s) Vs (m/s)
(ton/m3) Poisson (kg/cm2) (kg/cm2)
1 396 180 1.93 0.37 637 1,746
2 2,057 421 1.80 0.48 3,252 9,614
3 3,177 774 2.00* 0.47 12,214 35,870
*Pesos volumétricos estimados
Figura 5.1 Perfil sísmico obtenido a partir del procesamiento de los datos del PS024.
20
Con base a la experiencia que adquirí desarrollando el puesto de PASANTE DE
INGENIERA GEOFÍSICA, concluyo lo siguiente:
Es importante plantear objetivos claros, ya que la adquisición, procesamiento e
interpretación dependerán completamente de estos.
Tener antecedentes académicos robustos es parte vital para un buen desarrollo
del proyecto, en algunas ocasiones, la eficiencia con la que resuelven de
problemas de campo depende totalmente de los antecedentes teóricos que
tengamos.
La geología y la geofísica son ciencias interrelacionadas, es necesario contar
con el contexto geológico de donde queramos aplicar geofísica.
El informe es la presentación del trabajo realizado, por lo tanto es igual de
importante que la adquisición, el procesamiento e interpretación.
Como líder de brigada, el trabajo de tu equipo (bueno o malo) será visto como
tuyo. Así que, el trabajo en equipo es vital.
Finalmente, el trabajo de un Ingeniero Geofísico de campo es muy integral ya que no
solo requiere de conocimiento de geofísica sino de desarrollar habilidades tanto
técnicas como personales.
