1.

INTRODUCCION
El electromagnetismo es la rama de la física que estudia las relaciones entre la electricidad y el
magnetismo, es decir, el campo magnético creado por la corriente eléctrica y el efecto de un
campo magnético sobre una corriente eléctrica.
Las variaciones en los campos electromagnéticos originados en la ionosfera o magnetosfera
son medidos en la superficie terrestre con el fin de inferir la impedancia característica en función
de la profundidad. Conceptos simples de teoría electromagnética, análisis de series de tiempo
y teoría de sistemas lineales se han coalesado para formar un método de exploración,
denominado “método magnetotelúrico", mediante el cual es posible determinar la estructura
geoeléctrica de la Tierra a grandes profundidades.
Los métodos electromagnéticos utilizados en geofísica de exploración tienen como objetivo
principal el cuantificar las variaciones en conductividad eléctrica de la tierra en función de la
profundidad. Avances recientes en la teoría así como en la instrumentación, han permitido
utilizar estos métodos en la exploración de zonas geotérmicas, de hidrocarburos, mantos
acuíferos y zonas de alta mineralización caracterizadas éstas por presentar cambios anómalos
en la conductividad eléctrica.
Los métodos geo eléctricos para la prospección electromagnética, calculan la resistividad
eléctrica de los materiales del subsuelo, basados en la teoría electromagnética normada por
las ecuaciones de Maxwell y están fundamentados en el fenómeno conocido como inducción.
Los métodos electromagnéticos se pueden caracterizar por el tipo de fuente que utilizan: natural
o artificial. También se pueden dividir por el tipo de operación, pueden ser en
el dominio del tiempo o en el dominio de las frecuencias. Los métodos electromagnéticos (EM)
pueden utilizar dos tipos de fuentes, las artificiales (TDEM, FDEM), donde fluye una corriente
eléctrica y esta a su vez crea un campo electromagnético entorno a los medios ilimitados en
extensión y profundidad, donde por medio de las variaciones espaciales entre la conductividad
eléctrica y la susceptibilidad magnética podemos asociar la propagación de los campos
eléctricos y magnéticos en medios homogéneos. En esta técnica electromagnética donde se
utiliza un transmisor/receptor, conectada a una espira grande de alambre (Loop), inyecta una
corriente constante en el loop produciendo un campo eléctrico primario, la corriente es
interrumpida y por Ley de Faraday se produce la inducción electromagnética de corrientes
eléctricas en el subsuelo. Estas corrientes fluyen de manera cerrada viajando en profundidad y
lateralmente, la intensidad disminuye con la profundidad, generando un campo magnético
secundario transitorio el cual induce un voltaje variable en el receptor, el decaimiento de este
voltaje contiene información de la resistividad del subsuelo.
Los métodos EM de fuente natural (MT, AMT, CSAMT), miden los campos eléctricos y
magnéticos naturales en la Tierra generados por caídas de rayos y las interacciones entre el
viento solar y la magnetosfera. Se pueden llegar a medir hasta cinco componentes, dos
eléctricas y tres magnéticas y permite explorar la Tierra hasta varias decenas de kilómetros. La
operación en campo consiste en medir campos eléctricos y magnéticos en diferentes rangos
de frecuencia con el fin de determinar la resistividad aparente del terreno.
Las técnicas electromagnéticas de prospección del subsuelo, se basan en la medida y análisis
del comportamiento de los campos electromagnéticos inducidos en el terreno mediante
impulsos de corriente de cierta intensidad, circulando por una bobina situada horizontalmente
sobre el suelo.
El equipo de prospección se compone de un transmisor de corriente con su fuente de
alimentación, una bobina transmisora, el receptor de los campos inducidos y
un registrador digital de los datos obtenidos, que permite su posterior estudio y análisis.
MÉTODOS DE PROSPECCIÓN ELECTROMAGNÉTICOS
Los dos métodos más utilizados en estudios hidrogeológicos son:
• Very Low Frequency (VLF): Medidas electromagnéticas que permiten delimitar las fracturas
o fallas de un acuífero. Particularmente útil en caso de estudio de acuíferos fracturados como
los sistemas kársticos.
• Sondeos Electromagnéticos en el dominio temporal (SEDT o TDEM en ingles)

SONDEO POR TRANSITORIO

Las técnicas electromagnéticas de exploración geofísica se utilizan cada vez más para estimar
la resistividad del subsuelo, como una alternativa a los de corriente directa. El método del
sondeo transitorio (SEMT o TDEM por sus siglas en inglés) es una técnica de exploración
geofísica usada para estimar la resistividad eléctrica del subsuelo, con aplicaciones en diversas
áreas (geohidrología, minería, geotermia, etc.).
En general, los SEMT son realizados con una unidad transmisora unida a un alambre (espira
grande, bobina o dipolo eléctrico) y un sensor (bobina o dipolo eléctrico) que recibe y manda la
señal a una unidad receptora. Al inyectar una corriente constante en el alambre transmisor se
produce un campo magnético primario.A diferencia del Método Electromagnético basado en el
dominio de frecuencias, esta metodología estudia la variación de la resistividad eléctrica con
la profundidad, en lugar de generar distribuciones horizontales. En general, los sondeos SEDT
se realizan con un transmisor constituido por una espira metálica y un receptor en el centro de
la espira. Se inyecta una corriente eléctrica primaria en la espira, la cual produce un campo
magnético primario a su alrededor. Un corte rápido de esta corriente ocasiona la interrupción
del campo magnético primario, lo que produce, para satisfacer la ley de Faraday, la inducción
electromagnética de corrientes eléctricas en el subsuelo. Estas corrientes, que fluyen en
trayectorias cerradas en el subsuelo, migran en profundidad y lateralmente, mientras que su
intensidad disminuye conforme pasa el tiempo, lo que genera a su vez un campo magnético
secundario transitorio. Este campo secundario induce un voltaje variable en tiempo en el
receptor. La forma del decaimiento de este voltaje contiene información sobre la resistividad del
subsuelo, puesto que la magnitud y distribución de las corrientes inducidas depende de la
resistividad. El carácter migratorio en profundidad de las corrientes es usado como un control de
la profundidad, es decir, los voltajes de tiempos cortos proveen información de la resistividad
somera, mientras que los de tiempos largos la dan de la resistividad a mayor profundidad. El
método de dominio del tiempo de inducción electromagnética (TDEM) que consiste en generar
un campo electromagnético que induce una serie de corrientes en la tierra en profundidades,
cada vez mayor en el tiempo. Estas corrientes, a su vez, crean campos magnéticos. Mediante
la medición de estos campos magnéticos, se obtienen las propiedades del subsuelo y sus
características. Este método permite el estudio a grandes profundidades.
La forma en que la corriente migra en profundidad ayuda a controlar la profundidad de
investigación, donde los voltajes de tiempos cortos proveen información de la resistividad
somera y los voltajes de tiempos largos otorgan información de la resistividad a mayor
profundidad.

Los TDEM son realizados con una unidad trasmisora unida a un alambre (espira grande ó
bobina) que recibe y manda la señal a una unidad receptora. Entonces al inyectar una corriente
constante en el alambre transmisor se produce un campo magnético primario estable en el
subsuelo.
Esta técnica es principalmente usada para realizar sondeos a mayor profundidad con fines
geohidrológicos y geotérmicos, como son las correcciones estáticas (Static shift),
principalmente. Debido a su facilidad para realizar y el poco equipo que se necesita, es ideal
para zonas de difícil acceso.
Estos consisten en hacer circular
un campo eléctrico alrededor de una
bobina transmisora en periodos de
tiempo de milisegundos. Durante el
periodo de conexión se establece un
campo magnético estable en el subsuelo, y cuando este se corta de forma instantánea, el
campo electromagnético causa corrientes parásitas que se propagan en el terreno. Como
consecuencia de pérdidas de resistencia calórica estas corrientes disminuyen con el tiempo
provocando un campo magnético decreciente en la superficie, el cual puede ser medido.
Cuando en el subsuelo hay cuerpos de conductividad eléctrica elevada, la atenuación de las
corrientes parásitas es significativamente menor que en materiales dieléctricos; por tanto, la
medida de la relación de decrecimiento del campo secundario proporciona una forma de
detectar la presencia de cuerpos conductores en el subsuelo, así como estimar su
conductividad y resistividad (Casas, 2000). Es de esta forma que los TDEM muestran la
distribución del subsuelo en términos de homogeneidad, basados en la caracterización
resistiva.
El objetivo principal del método electromagnético TEM es mostrar la distribución del subsuelo;
siendo que es una técnica con gran resolución vertical, es posible obtener una imagen lo más
real posible de las condiciones del subsuelo. Los resultados de este método están en función
de la resistividad y enlazando dos o más sondeos se pueden construir secciones geoeléctricas
con las que se logran observar las heterogeneidades (zonas anómalas) debidas a estructuras
geológicas, cambios de facies y fracturamiento de la roca, donde además la presencia del agua
influye en el valor medido, provocando cambios importantes en la resistividad que es el
parámetro experimental de campo que se mide. En general, la técnica del TEM consiste en
construir una espira o bobina con un cable en forma de cuadrado, por la que se deja pasar una
corriente continua de gran intensidad que cree un campo magnético primario en la vecindad del
loop; después de cierto tiempo se interrumpe abruptamente el flujo de corriente y de acuerdo a
la Ley de Faraday, se produce la inducción electromagnética de una corriente eléctrica en el
subsuelo, que a su vez mantiene instantáneamente el campo magnético primario. Del mismo
modo, esta corriente genera un campo magnético secundario transitorio. Los cambios del
campo magnético secundario con el tiempo inducen un voltaje en la bobina receptora conocido
como transitorio. La forma del decaimiento de este voltaje contiene información sobre la
distribución de la resistividad del terreno debido a que la magnitud y la distribución de la
corriente inducida dependen de él.

En la siguiente figura se muestra una configuración típica de sondeo de resistividad TDEM,

donde se ve que el transmisor está conectado a un bucle cuadrado (generalmente de una sola
vuelta) de cable tendido en el suelo. La longitud del lado del bucle es aproximadamente igual a
la profundidad de exploración deseada, excepto que para profundidades poco profundas
(menos de 40 m), la longitud puede ser tan pequeña como 5 a 10 m en un terreno relativamente
resistente. Una bobina receptora de varias vueltas, ubicada en el centro del bucle del
transmisor, se conecta al receptor a través de un cable corto.

PRINCIPIOS FISICOS

El método TEM consiste en

medir el campo secundario inducido
en el terreno una vez que se
suprime el campo primario que lo
produjo.A grandes rasgos
,la ley que sustenta este método
de exploración es la Ley de
Faraday.Cuando una
corriente eléctrica circula
por un conductor se genera un
campo magnético ,si este campo
esta en presencia de un material
ferromagnético ,este campo se induce en el material.
El método TEM consiste en tender sobre el terreno una espira de forma cuadrada y hacer
circular corriente eléctrica por ella ,cuando la espira esta energizada se suprime la corriente y
se generara un campo magnético ,a este fenómeno se le conoce como transitorio.Este campo
magnético será inducido en el terreno e ira desapareciendo ,es decir,tendrá un decaimiento.En
este decaimiento del campo magnético secundario esta contenida información respecto a las
propiedades eléctricas del subsuelo,una de ellas es la resistividad.

APLICACIÓNES DEL SONDEO TRANSITORIO

este método en particular, podemos darle mucha utilización, como por ejemplo:

 Estratigrafía del terreno atendiendo a sus propiedades eléctricas.

 Detección de depósitos enterrados de residuos urbanos e industriales.
 Definición de zonas de intrusión marina.
 Determinación de niveles freáticos.
 Identificación de acuíferos contaminados.
 Estudio de aguas subterraneas
 Se aplica en la exploración del agua subterránea, minerales, ambiental y geotermia.

VENTAJAS DEL METODO SONDEO TRANSITORIO SOBRE ONDA CONTINUA

 Una gran sensibilidad a los terrenos conductores.
 Un poder de resolución lateral y vertical importante.
 Una buena detección delas anomalías de baja resistividad en un terreno conductor.
 Una ejecución práctica debida a la ausencia de contacto con el suelo y, por lo tanto,
rápida en un terreno despejado.
En el TEM ,el área de medición es muy enfocada ,este hecho permite la obtención de
resultados mas congruentes con las condiciones reales del terreno.Esto describe mejor
situaciones complejas tales como fallas,estructuras,cambios de facies,etc.El método TEM
tiene la capacidad de penetración en capas muy conductoras o resistivas,como podría ser el
caso de medios muy arcillosos o rocas muy compactas ,es en estos medios donde el SEV no
logra inyectar la corriente.
El equipoTEM es ligero a comparación del equipo SEV,por tanto es mas rápido y fácil la
adquisición de datos.El TEM es mas barato que el SEV. Las ventajas principales de los
métodos transitorios sobre los métodos de onda continua son que el campo primario no está
presente durante la medición del campo secundario y que las mediciones del campo
secundario en función del tiempo son equivalentes a mediciones de onda continua sobre un
amplio rango de frecuencias. La ventaja del TDEM radica en su capacidad para determinar las
profundidades y mapear características mucho más profundas. Los ensayos FDEM son
recomendables para explorar grandes superficies. Métodos transitorios se utilizan tanto para
sondeo de profundidad como para perfiles continuos.También se denomina método
electromagnético de dominio del tiempo "TDEM
LOS INCONVENIENTES SON LOS SIGUIENTES:
-Una baja resolución de los terrenos resistentes.
- Una resolución de la parte superficial bastante limitada.
- Dificultades de ejecución en áreas arboladas.
El método Transitorio Electromagnético en Dominio del Tiempo (TEM o TDEM), es una
técnica geofísica para la obtención de la resistividad eléctrica del subsuelo por medio del
fenómeno de inducción electromagnética.
Unejemplo es que en Geophysical Surveys contamos con equipo AEMR TEMFAST 48 HPC y
MONEX GEOSCOPE TERRATEM para la ejecución de sondeos TEM.

CONFIGURACIONES PARA TEM

El arreglo geométrico que se utiliza en todas las técnicas electromagnéticas inductivas en el
dominio del tiempo es una bobina con forma rectangular casi perfecta, la cual sirve como
transmisor. Existen varios tipos de arreglos para TEM, que están en función de la configuración
que presenten el transmisor y receptor, y son los siguientes:
1. Bobina simple. La configuración de este arreglo consta de una sola bobina, que funge
como receptor y transmisor, donde la bobina actúa como transmisora en el mismo
instante en que fluye la corriente a través de ella. Una vez que se apaga la corriente, la
bobina se conecta al receptor y se mide la señal transitoria durante todo el periodo en el
que el transmisor está en reposo. Las bobinas pueden presentar forma cuadrada o
rectangular donde cada lado tiene una longitud que puede variar de 5 a 300 m.
2. Bobina coincidente. Tiene la misma geometría y respuesta que el arreglo de bobina
simple, pero aquí el transmisor y receptor son bobinas independientes que coinciden en
tamaño.
3. Bobina central. Es variante de la bobina coincidente; se utiliza como un dipolo múltiple
que se coloca en el centro de la bobina transmisora.
4. Bobinas separadas. Transmisor y receptor se encuentran separados a una distancia fija
y las bobinas tienen longitudes de hasta unas cuantas decenas de metros por lado. La
variante principal en este arreglo es que el dipolo receptor permanece fijo a una distancia
en el exterior de la bobina transmisora.
5. Bobina doble. Para este arreglo se utilizan dos bobinas adyacentes y conectadas en
paralelo. Con esta configuración se reduce significantemente el ruido causado por
fuentes remotas que interfieren con la señal: sin embargo, en zonas donde la
conductividad varía lateralmente dicho ruido no se logra reducir.
6. Transmisor largo fijo con receptor móvil. Consiste en una bobina grande y fija que
funciona como transmisora, mientras que el receptor es pequeño y móvil para tomar
mediciones a lo largo de líneas perpendiculares a la bobina. La longitud de cada lado
 de la bobina va de unos cuantos metros hasta un máximo de 2 km, aunque también han
sido utilizadas bobinas de hasta 5 km.
7. TEM en pozo. Es una extensión del arreglo del transmisor fijo, en donde el dipolo
receptor se baja a lo largo de un pozo. De esta manera se puede medir a diferentes
profundidades la componente del campo magnético transitorio orientada a lo largo del
eje del pozo.

SISTEMAS DE TDEM
El diseño modular de los Sistemas Electromagnéticos de Dominio del Tiempo (TDEM) de
PROTEM, un único receptor con múltiples transmisores intercambiables y bobinas receptoras,
proporciona la flexibilidad necesaria para adaptarse a una amplia gama de aplicaciones, desde
estudios ambientales y de ingeniería cerca de la superficie, hasta tres. Exploración de recursos
de componentes a profundidades de 1000 m o más.
 Receptor PROTEM : el componente base de cualquier sistema PROTEM, el Receptor
PROTEM proporciona una calidad de datos superior y un funcionamiento confiable
 Transmisor TEM47 : pequeño y portátil, el TEM47 a batería es apropiado para
aplicaciones cercanas a la superficie dentro de un rango de profundidad de 5-150 metros
 Transmisor TEM57-MK2 : ideal para aplicaciones de sondeo o perfilado de rango
medio, el TEM57-MK2 proporciona mediciones dentro del rango de 200-500 m de
profundidad
 Transmisor TEM67 : que proporciona la máxima profundidad de exploración, el TEM67
es adecuado para aplicaciones a profundidades de 1000 m o más
 Sistema PROTEM CM : consolidando el receptor PROTEM y los componentes del
transmisor TEM47 en una sola consola
 Sistema G-TEM : con el objetivo de reducir costos al tiempo que mantiene las
características de rendimiento principales de los sistemas PROTEM más avanzados y
agrega varias características nuevas que mejoran el rendimiento, el nuevo sistema G-
TEM representa un valor excepcional en la instrumentación electromagnética en el
dominio del tiempo (TDEM)
 SENSORES DE SUPERFICIE : están disponibles seis bobinas de inducción y un sensor
magnetómetro de flujo de tres componentes
 Transmisor EM42 : disponible solo por pedido especial, el EM42 es capaz de explorar
a profundidades del orden de 3000 metros; información adicional disponible por solicitud

Todas las configuraciones del sistema están disponibles con opciones para la recopilación de
datos secuencial o simultánea de tres componentes (xyz).
La recopilación y el procesamiento de datos asociados con cualquier sistema PROTEM están
respaldados por un sistema integrado de adquisición de datos. El software de procesamiento
incluido ( ver Software - GSPX7) facilita la revisión de datos, la edición, la presentación de
perfiles y las funciones de modelado asociadas con objetivos discretos; el software para el
modelado directo e inverso de datos de sondeo está disponible a través de terceros
comerciales.

EQUIPO

TERRATEM
Consola Receptor/transmisor TerraTEM
Capacidad de registros 400 registros completos
No. de Ventanas 53 en tres series
Rango de mediciones 50 microsegundos a 2 segundos
Promedios de señal 1-9999 muestras (seleccionable)
Tiempo de retraso 0-9999 microsegundos (seleccionable)
Filtros 50/60 Hz
Desplegado en campo Normalizado de valores de
Resistividad Aparente (en modo textual o gráfico)
Peso 10 kg

ESTUDIO GEOFÍSICO PARA DEFINIR LAS CARACTERÍSTICAS GEOHIDROLÓGICAS DEL

ACUÍFERO APAN EN UN PREDIO DELIMITADO
Esto ha conllevado a desarrollar técnicas de exploración que nos permitan obtener información
del subsuelo, que nos ayude a comprender la distribución espacial de sus estratos y las
características físicas que se puedan asociar a la presencia de agua subterránea. En la zona
de estudio se contempla realizar la construcción de un Parque Industrial que albergará grandes
empresas, que a su vez contribuirán con el desarrollo económico de la zona.
Por esas razones se encargó a GEOTEM INGENIERÍA la realización de un estudio exploratorio
con la técnica de Transitorio Electromagnético en el Dominio del Tiempo con el objetivo de
definir las características geohidrológicas del acuífero Apan y establecer zonas de interés para
la perforación de pozos exploratorios.
LOCALIZACIÓN
El primer paso, antes de comenzar con la planificación del estudio fue localizar el área de
estudio, la cual se encuentra en el estado de Hidalgo, 67 km al noreste de la Ciudad de México
y 50 Km al sureste de Pachuca. Los sondeos se encuentran distribuidos en una zona de
aproximadamente 10 km2. Se tienen un total de 15 sondeos, de los cuales 13 se encuentran
en el municipio de Apan y 2 en el municipio Emiliano Zapata.

LOCALIZACION DE LA ZONA DE ESTUDIO

 LOCALIZACION DE SONDEOS EN MAPA DIGITAL DE ELEVACIONES

MARCO FÍSICO

ACUÍFERO APAN

El Acuífero Apan cuenta con un área de 1,518 km2 aproximadamente, pertenece a la Cuenca
del Valle de México que a su vez corresponde a la Región Hidrológica XIII y a la subcuenca
Apan-Tochac. Se trata de un acuífero semiconfinado subexplotado.
MAPA DE LA ZONA DE ESTUDIO CON PROYECCION DEL ACUÌFERO DE APAN

FISIOGRAFÍA

Provincia Fisiográfica

La cuenca del Valle de México se ubica en la provincia fisiográfica denominada "Eje

Neovolcánico Transmexicano" que cruza al Territorio Nacional del Golfo al Pacifico; se divide
en tres zonas principales: Zona Meridional, Zona Septentrional y Zona Nororiental.
Litológicamente está constituida por rocas volcánicas, principalmente coladas de lava y
materiales piroclásticos que han sido depositados en ambientes lacustres hasta fechas
recientes.
El área en general del acuífero Apan corresponde a una planicie aluvial en donde se encuentran
arcillas de alta plasticidad y dispersamente arcillas cementadas. Existen también aparatos
volcánicos compuestos por lavas, brechas y cenizas cuya composición litológica va desde rocas
basálticas a riolíticas con edades que varían desde el Plioceno Tardío hasta el reciente, lo que
ha originado alteración por intemperismo en las rocas superficiales. Se tienen llanuras y
cuencas formadas por rellenos aluviales o lacustres que contienen gran variedad de rocas
mezcladas con cenizas volcánicas.

GEOMORFOLOGÍA

Relieve Terrestre

En promedio la planicie de la Cuenca del Valle de México se encuentra a 2,250 msnm, sin
embargo tiene puntos naturales sobresalientes como el Cerro El Manantial al Norte, con 3,190
msnm y el Cerro San Antonio con 2,570 msnm, al este las Sierras Chichucuatlán y del Tepozán
que se derivan de la Sierra Madre Oriental en la que destacan los cerros Agua Azul con 3,020
msnm; y la Perla con 3,200 msnm; Las Tetillas 3,020 msnm y la Peñuela con 3,350 msnm.
La presencia de estas formaciones orográficas ha evitado en gran medida la erosión de las
tierras cultivadas con azolve de cauces naturales y con otras formas de deterioro de las
riquezas.
En el caso del acuífero Apan se observa una tendencia de recarga natural en las zonas de
lagunas ya que los gradientes disminuyen, razón por la cual posiblemente dichas lagunas
constituyan zonas de recarga.
ESTRATIGRAFÍA

A continuación se tiene una descripción litológica general de las unidades geológicas más
antiguas a las más jóvenes, primeramente las sedimentarias y posteriormente las ígneas.
Se tienen afloramientos de rocas del Cretácico, rocas calizas correspondientes a unidades
sedimentarias de origen marino. Del Terciario Inferior se tienen conglomerados constituidos por
clastos de caliza y pedernal mal seleccionados.
Del Terciario Superior se tienen lutitas, limolita arenisca , conglomerados con clastos de
andesita.
Del Cuaternario se tiene un depósito clástico continental, compuesto por un conglomerado
brechoide poco compactado y descansa de manera discordante sobre rocas ígneas del
Terciario presentándose en forma de abanico aluvial. Corresponden a esta época los suelos
aluviales que son depósitos clásticos no consolidados con gravas, arenas, limos y arcillas
provenientes de las rocas preexistentes constituidas en forma de abanicos y planicies.
Las rocas ígneas del terciario superior, tales como andesitas, presentan pseudo estratificación
y se intemperizan en forma de lajas presentándose morfológicamente como aparatos
volcánicos de elevación, lomas y cerros bajos.
También se encuentra andesita y brecha volcánica intermedia, unidades formadas por
derrames andesíticos compactos alternados con productos piroclásticos compuestos por
brechas, cenizas y lapilli.

MAPA GEOLOGICO DE LA ZONA DE ESTUDIO

GEOLOGÍA ESTRUCTURAL

El relieve estructural de la región está íntimamente relacionado con una intensa actividad
volcánica iniciada a principios del Terciario y desarrollada durante el Pleistoceno Inferior. El
conjunto de estructuras que caracterizan el relieve de esta provincia evolucionó sobre una
paleogeografía constituida por sedimentos mesozoicos plegados, los cuales correspondían a
la Sierra Madre Oriental.
es testimonio de la durabilidad del fenómeno durante un lapso grande de tiempo que abarcó
desde principios del Terciario hasta épocas recientes.
Existen en la región sistemas de fracturas y fallas regionales siendo tres las principales: La
primera tiene dirección Noroeste-Sureste y afectó previamente a las rocas mesozoicas y está
asociada al patrón estructural de la Sierra Madre Oriental, los otros dos sistemas están
orientados Norte Sur y Noreste-Oeste Sureste y están asociados a los esfuerzos de Tensión
que afectaron a las rocas volcánicas del Mioceno.

HIDROGRAFÍA

La red hidrográfica en la Cuenca del Valle de México está representada por un gran número de
corrientes; estas integran un sistema general bastante definido; controlado principalmente por
los sistemas de fracturamiento y fallamiento por lo que sus cauces son de dimensiones
reducidas con paredes escarpadas.
Las condiciones naturales de los escurrimientos han sido modificadas con el transcurso del
tiempo debido a los crecimientos de las áreas urbanas.
Estas corrientes no son caudalosas debido a que la precipitación pluvial en la zona es baja con
una fuerte evaporación y permeabilidad vertical alta controlada por sus características
geológicas, principalmente en las brechas volcánicas, basaltos en bloques y riolitas verdes que
han sido afectadas por el Tectonismo.
Por lo antes mencionado el Acuífero de Apan pertenece a los depósitos lacustres o de planicie
aluvial ocupada por las lagunas someras de Tecocomulco con 20 km2 de área aproximada,
Apan y Tochac con 5.5 km2 entre otras de menor cuantía.

CLIMA

El acuífero Apan, que pertenece a la Región XIII del Valle de México, se encuentra influenciado
por las características topográficas existentes y por lo tanto se tiene una gran diversidad de
climas, sin embargo dentro de esa variedad predomina el templado o mesotérmico siguiendo
en importancia los climas secos y semisecos y en menor grado los climas fríos y semifríos.

1.1 PROCESAMIENTO DE DATOS

El procesado de los registros se efectúo en las siguientes etapas, para cada sondeo
y fueron:

 Transferencia de los registros del equipo a la computadora.

 Edición de las curvas registradas con cada ganancia para obtener una
curva final de nV/AMP contra tiempo (curva de decaimiento).
 Conversión de la curva de decaimiento a Resistividad Aparente contra
 tiempo
 Cálculo del modelo estratificado
 Cálculo del modelo suavizado (inversión Occam)
 Impresión gráfica del modelo unidimensional
 Transferencia de las curvasy los modelos generados
al paquete WINGLINK
 En Winglink, formación de la base de datos
 Elaboración de plantas y secciones de isorresistividad e impresión de
resultados gráficos.

Todas las etapas del proceso anterior se realizaron en la oficina, y los datos de
campo fueron analizados el mismo día que fueron levantados, para conocer la
calidad de la información.

Para la ubicación de los sitios de medición, se utilizó un geo posicionamiento tipo

GPS, empleando para ello un GPS de mano marca “Magellan”, modelo “Explorist
510”, el cual nos proporciona para cada sitio de medición sus coordenadas (X,Y),
mientras que la elevación fue tomada de las cartas topográficas editadas por el
INEGI.
En el proceso de los datos se siguieron los estándares internacionales que se
aplican para este fin, utilizando software y equipos de patente, de los cuales se
tienen versiones originales, para ello podemos citar el siguiente equipo y software:

 Equipo Terratem con número de serie 1503

 Programa TEMIXS para el proceso primario de los datos
 Programa WINGLINK para utilizarlo como base de datos y formar secciones
y mapas de resistividad.

El proceso de los datos, para cada sondeo, consiste en:

Obtener un modelo unidimensional, partiendo de la curva de resistividad aparente

la cual es procesada en una primera etapa mediante el modelado por capas planas,
siguiendo el método clásico de prueba y error, en el que se propone un cierto
número de capas, dependiendo de las inflexiones de las curvas y los tiempos en los
que ocurren las inflexiones.

Una segunda etapa en el proceso de la curva de resistividad aparente es aplicar

una regresión tipo Occam con lo cual se logra obtener un modelo multicapas (17
capas) para integrar secciones y simular un comportamiento bidimensional de la
sección. A continuación se describe cada proceso.

MODELOS SUAVIZADOS

El procesado de los datos consistió en obtener los modelos suavizados mediante

una regresión tipo Occam, en la que se propone una resistividad inicial y una final,
basada en los datos de campo.
El algoritmo requiere de un número de capas inicial, un valor de ajuste al que se
quiere llegar y un número máximo de iteraciones, estos tres últimos parámetros no
intervienen en la obtención del modelo suavizado, por lo que este sólo depende de
los valores de resistividad de campo. Entonces la técnica se basa en el precepto
filosófico dictado por William de Occam en los albores del siglo XIV que dice: "Es
en vano hacer con más lo que puede ser hecho con lo menos", y en el uso de la
regresión Ridge, para encontrar un modelo que cumpla con el error marcado en el
ajuste o con el mínimo error posible.

MODELO UNIDIMENSIONAL

También los datos fueron procesados por el modelado por capas o unidimensional,
se hizo siguiendo el método clásico de prueba y error, en el que se propone un cierto
número de capas, dependiendo de las inflexiones de la curva de campo, las
pendientes de la curva y los tiempos a los que ocurren las inflexiones.

Se aprovechó la facilidad del paquete TEMIXS que permite la interacción con la

computadora en un proceso conocido como interpretación interactiva. Este proceso
posibilita que el intérprete modifique los parámetros del modelo y presenta el cálculo
de las curvas teóricas de las modificaciones realizadas, esta característica hace que
el ajuste del modelo sea rápido.

Una vez que se ha obtenido un ajuste aceptable, el modelo generado se somete a

una regresión automática con las técnicas de Ridge y de Marquardt con algoritmos
de aproximación de tipo Anderson para obtener, por inversión, el mejor ajuste bajo
el criterio del error cuadrático medio mínimo. Los modelos obtenidos, son sometidos
a un nuevo proceso para obtener una serie de modelos equivalentes que cumplan
con el mismo error de ajuste. Los modelos equivalentes posibilitan
establecer el intervalo de variación de los espesores y resistividades en cada capa
del modelo original con fines de correlación electro estratigráfica.

Para la obtención de un modelo unidimensional es necesario interpretar la curva de

resistividad aparente mediante un modelado por capas planas, siguiendo la
metodología de prueba y error. Seguido se aplica una regresión tipo Occam para
obtener un modelo multicapas (17 capas) para integrar secciones y simular un
comportamiento bidimensional de la sección.

Posteriormente se lleva a cabo un proceso llamado Inversión Unidimensional

(Modelos 1D), que consiste en obtener una curva de resistividad contra profundidad
mediante una regresión tipo Occam, en la que se propone un valor de resistividad
inicial, un número de hasta 17 capas y una profundidad estimada; esta curva de
resistividad calculada es comparada con los datos reales de campo. El algoritmo
requiere de un número de capas inicial, un valor de ajuste al que se quiere llegar y
un número máximo de iteraciones; estos parámetros no intervienen en la obtención
del modelo suavizado, por lo que la obtención de un modelo representativo de la
curva medida sólo depende de los valores de resistividad de campo.
 FIGURA 9. AJUSTE DE LA INVERSIÓN UNIDIMENSIONAL

1.2 RESULTADOS

1.2.1 CLASIFICACIÓN GEOELÉCTRICA

A partir de los contrastes de resistividad observados en los perfiles geoeléctricos y

la información geológica disponible, fue posible definir 3 unidades geoeléctricas
principales, de las cuales, una fue subdividida para lograr un mejor análisis. Estas
unidades representan variaciones litológicas o facies de los diversos cuerpos
estratigráficos que componen el subsuelo. La siguiente tabla muestra la clasificación
geofísica obtenida.
 TABLA 1. CLASIFICACIÓN GEOFÍSICA

INTERVALO
UNIDAD RESISTIVO DESCRIPCIÓN
(Ohm*m)

Capa de cobertura correspondiente a depósitos aluviales compuestos

U1 5 a 30
por arenas, gravas, limos y arcillas.

Estrato formado por materiales de origen volcánico de buena

U2a 5 a 20 permeabilidad, tales como tobas con posible fracturamiento y
posibilidades de contener saturación.

Capa compuesta por materiales piroclásticos como tobas y brechas

U2b 20 a 40
volcánicas con poco fracturamiento y regular permeabilidad.

Estrato formado por materiales de origen volcánico con intercalaciones

U3a 20 a 50
de arenas. Se considera de regular permeabilidad y posible saturación.

Unidad integrada por rocas de mayor compacidad, de origen volcánico,

U3b Mayor a 50 tales como andesitas y basaltos.

Unidad 1: Capa de cobertura

La unidad U1 corresponde a la capa de cobertura, la cual está integrada por

depósitos aluviales tipo arcilla que presentan valores resistivos de entre 5 y 30
Ohm*m y un espesor promedio de 20 m sin interrupciones en su continuidad lateral.

Unidad 2: Material Volcánico

U2 ha sido subdividida para un mejor análisis, donde la subunidad U2a es un estrato

formado por materiales de origen volcánico que presumen tener buena
permeabilidad y posibilidades de saturación por agua subterránea, determinado por
sus valores resistivos que oscilan entre 5 y 20 Ohm*m, que pueden asociarse a
tobas con un grado considerable de fracturamiento. La subunidad U2b presenta
valores de resistividad de 20 a 40 Ohm*m, que se asocian a materiales
piroclásticos como tobas y brechas volcánicas con poco fracturamiento y regular
permeabilidad.

Unidad 3: Rocas Volcánicas

U3 también ha sido subdividida; la subunidad U3a muestra un intervalo resistivo de

entre 20 y 50 Ohm*m, valores que se asocian a materiales de origen volcánico con
intercalaciones de arenas, que se consideran de regular permeabilidad y posible
saturación. U3b muestra valores resistivos mayores a 50 Ohm*m, que se asocian a
rocas de mayor compacidad que de igual forma son de origen volcánico, como
andesitas y basaltos.

1.2.2 PERFILES GEOELÉCTRICOS

A partir de los 15 sondeos por Transitorio Electromagnético realizados en el área,

se trazaron 7 perfiles geoeléctricos, que se obtuvieron de los modelos suavizados.
Son de gran utilidad, pues conjuntan las variaciones de resistividad en el plano
vertical y ayudan a definir las características de los materiales, además de mostrar
la distribución geoespacial de la resistividad del subsuelo como respuesta de los
materiales rocosos que conforman los estratos en conjunto con posibles fluidos
emplazados en su interior.
 FIGURA 10. MAPA DIGITAL DE ELEVACIONES CON PROYECCIÓN DE SONDEOS Y PERFILES

Los perfiles P-1, P-2 y P-3 son paralelos entre sí, tienen una dirección general SW-
NE y longitudes de 3230 m, 2765 m y 5710 m respectivamente. El Perfil P-1 está
formado por los sondeos PI-01 a PI-04; el Perfil P-2 por los sondeos PI-05 a PI-08
y el Perfil P-3 por los sondeos PI-09 a PI-15. Los perfiles P-4, P-5, P-6 y P- 7 tienen
dirección NW-SE y una longitud promedio de 2600 m; cada uno está conformado
por 3 sondeos cubriendo la zona central del área de exploración. De esta forma se
tienen múltiples vistas de la distribución georresistiva del subsuelo y
por lo tanto es posible definir el mejor sitio para el emplazamiento de un pozo
exploratorio.

En base a la clasificación geoeléctrica, se tiene una distribución de unidades

georresistivas con buena continuidad lateral en casi todos sus estratos y perfiles.
Variando únicamente en sus espesores y profundidades. El análisis de los perfiles
se realiza en 2 grupos principales, de acuerdo a su orientación.

Perfiles P-1, P-2 y P-3

Los 3 perfiles muestran una capa de cobertura U1 correspondiente al aluvión que

predomina en la superficie de la zona, con un espesor promedio de 20 m y que se
extiende a lo largo de los perfiles sin interrupciones. Subyaciendo se encuentra la
subunidad U2b que presenta un espesor considerable, superior a los 150 m y bajo
algunos sondeos alcanza los 350 m de espesor, como es el caso de los sitios PI-
04 y PI-08. Debajo de U2b se encuentra un estrato de valores resistivos menores a
20 Ohm*m, lo cual es indicador de posible saturación por líquido subterráneo; esta
subunidad nombrada U2a corresponde a materiales de origen volcánico de buena
permeabilidad, posiblemente tobas dacíticas y andesíticas con un grado
considerable de fracturamiento. Muestra un espesor variable, en el perfil P-1 oscila
entre 250 y 400 m, para P-2 es de entre 150 y 350 m y para P-3 es de entre 120 y
250 m, siendo la zona sur la de menor espesor en la zona saturada.

A los 550 m de profundidad aproximadamente, comienzan a incrementar

nuevamente la resistividad de los materiales, por lo que se marca un contacto con
un estrato de valores resistivos que oscilan entre 20 y 50 Ohm*m, nombrado U3a,
que corresponde a rocas de origen volcánico con posible fracturamiento y regular
permeabilidad. Tiene un espesor de entre 50 y 150 m y una continuidad lateral
homogénea. Subyaciendo a U3a se observa el estrato más profundo del perfil,
que tiene valores resistivos mayores a 50 Ohm*m, los cuales se relacionan con
roca de origen volcánico de mayor compacidad y baja permeabilidad.
 FIGURA 11. PERFILES GEOELÉCTRICOS P-1, P-2 Y P-3

Perfiles P-4, P-5, P-6 y P-7

Estos perfiles muestran una distribución de capas similar a los perfiles anteriores,
donde todos los estratos tienen buena continuidad lateral y únicamente difieren en
sus espesores. Analizando la profundidad y espesor de cada unidad se puede
observar que la subunidad U2b aumenta su espesor hacia el oriente de la zona,
mientras que la subunidad U2a disminuye en la misma dirección. Siendo el perfil P-
5 donde se observa que el estrato conductor U2a que está relacionado con la
posible zona saturada, muestra su mayor espesor, por lo tanto se considera la parte
NW de la zona de estudio como la mejor zona para realizar la perforación de un
pozo exploratorio.
FIGURA 12. PERFILES P-4, P-5, P-6 Y P-7
 1.3 CONCLUSIONES Y RECOMENDACIONES

De acuerdo al análisis de los perfiles geoeléctricos y la información geológica y

geohidrológica recabada se concluye lo siguiente:

 El acuífero Apan cuenta con un área de 1,518 km2 aproximadamente,

pertenece a la Cuenca del Valle de México, Región Hidrológica XIII,
subcuenca Apan-Tochac.

 El área del acuífero Apan corresponde en su mayoría a una planicie aluvial,

formada por arcillas de alta plasticidad y arcillas cementadas. Se tiene
presencia de aparatos volcánicos compuestos por lavas, brechas y cenizas
cuya composición litológica va desde rocas basálticas a riolíticas con edades
que varían desde el Plioceno Tardío hasta el reciente.

 Con la finalidad de conocer la distribución de resistividad en el subsuelo y

conocer la características georresistivas del acuífero Apan bajo la zona de
estudio, se realizó un estudio geofísico en la modalidad de Transitorio
Electromagnético (TEM). Se definieron 15 sitios de estudio con loops de
300x300 m con los que se alcanzó una profundidad de 700 m
aproximadamente.

 Del procesamiento y análisis de la información se definieron 3 unidades

geoeléctricas principales, 2 de ellas con subdivisiones; de las cuales las
subunidades U2a y U3a se consideran con potencial geohidrológico
importante. La subunidad U2a tiene un intervalo resistivo de entre 5 y 20
Ohm*m, valores asociados a materiales de origen volcánico de buena
permeabilidad como tobas con posible fracturamiento y posibilidades de
presentar saturación por agua subterránea. También se considera como
estrato de interés a la subunidad U3a que subyace a U2a, que presenta
valores resistivos de entre 20 y 50 Ohm*m, asociados a materiales de
 origen volcánico con intercalaciones de arenas y se considera de regular
permeabilidad con posible saturación.

 Se crearon 7 perfiles geoeléctricos con los que se cubre la totalidad del área
de estudio, 3 de ellos con dirección NW-SE y 4 con dirección SW-NE.

De acuerdo con los resultados encontrados en el presente trabajo y lo concluido

anteriormente, se observa que en general se muestran condiciones geológicas
similares en los 15 sondeos TEM realizados, donde los espesores de cada unidad
es el principal factor para la elección del mejor sitio para la perforación de un pozo,
por lo tanto se recomienda lo siguiente:

 El sondeo que muestra el mayor espesor de U2a es el sondeo PI-02, por lo

tanto se considera este sitio como el más recomendable para la realización
de un pozo exploratorio, que deberá alcanzar una profundidad del orden de
500 m para asegurar que quede emplazado dentro del acuífero profundo de
Apan. De este modo se podrá asegurar el gasto requerido para el Parque
Industrial.

 Se recomienda realizar un registro geofísico de pozo con una sonda

multiparámetros que pueda aportar información de radiación gamma natural,
resistencia del fluido, SP, curvas de resistividad y temperatura para realizar
un dictamen más integral, en cuanto a las posibles zonas permeables.
FIGURA 13. PERFIL P-1 CON UBICACIÓN DE POZO PROPUESTO A 500 M DE PROFUNDIDAD
BIBLIOGRAFÍA

CONSTABLE S.C., Parker R.L., Constable C.G., "Occam's Inversion: A

practical algorithm for generating smooth models from electromagnetic
sounding data", Geophysics, v 52, n3, Marzo 1987.
APPLIED GEOPHYSICS SECOND EDITION, W.M. Telford, L.P. Geldart,
R.E. Sheriff, Cambridge University Press, United Kingdom 1990.
Flores, L, C., La exactitud del problema directo de sondeos
electromagnéticos transitorios, GEOS, Vol 20, 2, 70-88.
Mc Neill, J.D., 1994, Principles and application of time domain
electromagnetic techniques for resistivity sounding, Technical Note TN-27, 1-
15.
Nabighian, M.N., and Macnae, J.C., 1991, Time domain electromagnetic
prospecting methods, in EMAG, v.2, Applications, Part A, Nabighian (ed.),
427- 479.
Orellana, S. E., 1982, Prospección geoeléctrica en corriente continua,
PARANINFO, 2da Edición.
Subdirección General Técnica, Gerencia de Aguas Subterráneas
Subgerencia (2008). Determinación de la Disponibilidad de agua en el
acuífero 3218 Cedros, Estado de Zacatecas. México: CONAGUA.
Hidrología de la Cuenca del Valle de México SRH-CHCUM año 1964.
Estudio de Geología en el Valle de México con fines Geohidrológicos
realizado a contrato en 1987 para CAVM por el Instituto de Ingeniería de la
UNAM.

SOFTWARE

WINGLINK, User´s Manual, Edition 4, © Geosystem srl, Milan, Italia, 1994.

©Interpex Limited, Golden Co. 1993.
OFFICE 2013
TEMIX S
SURFER v12.0
AUTOCAD 2014