6-26

6.4.1.3 Adquisición Gravimétrica

La adquisición de datos gravimétricos se realizó entre el 13 y el 25 de Marzo de 2009,

utilizando un gravímetro SCINTREX CG3M (Fotografía 6.2). Se registró un total de
350 estaciones, a lo largo de las cuatro transectas, en una longitud total (medidas en línea
recta) de 133.5 km. Se había planteado realizar una transecta de gravimetría a través del
golfo de Urabá, aprovechando un viaje de un barco maderero, pero no fue posible debido
a que no se contaba con el equipo adecuado para tal fin.

Fotografía 6.1 Estación gravimétrica La Fotografía 6.2 Gravímetro SCINTREX CG3M

Chinita

6.4.1.4 Procesamiento de Datos Gravimétricos

Para que los valores de gravedad medidos puedan ser interpretados, es necesario
corregir todas las variaciones del campo gravitacional de la Tierra que no dependen de las
diferencias de densidad de las rocas subyacentes. Este proceso es conocido como
reducción de la gravedad o reducción de la gravedad al geoide, donde el nivel del mar es
usualmente el datum más conveniente. Para ello se aplican las siguientes correcciones
por deriva y por factor de escala instrumental.

Corrección por Factor de Escala Instrumental

El factor de escala instrumental corrige una lectura instrumental en una lectura en

miligales, utilizando un factor de conversión.
rc = rS (r ) (Ec. 6.2)

Donde: rc = valor corregido, en miligales

r = lectura instrumental en unidades del dial del equipo
S(r) = factor de escala
 6-27

El gravímetro CG3 SCINTREX registra las medidas directamente en unidades de

miligales. El factor de escala instrumental es igual a uno (1), por lo tanto, utilizando este
gravímetro no se requiere aplicar esta corrección.

Corrección por Mareas

Las medidas de gravedad en una localidad fija varía con el tiempo debido a la variación
periódica en los efectos gravitacionales del sol y la luna, asociados con sus movimientos
orbitales, razón por la cual correcciones deben ser hechas para estas variaciones en
exploraciones de alta precisión. A pesar de su menor masa, la atracción gravitacional de
la Luna es más grande que la del Sol, debido a su proximidad. También estos efectos
gravitacionales causan que la forma de Tierra sólida varíe en la misma forma que la
atracción celestial causa las mareas en el mar. Las mareas terrestres son
considerablemente más pequeñas que las mareas oceánicas y están lejos de seguir el
movimiento lunar. Ellas causan la elevación de un punto de observación al ser alterado
unos pocos centímetros y así varían la distancia desde el centro de masa de la Tierra Las
variaciones periódicas de la gravedad causadas por los efectos combinados del Sol y la
Luna son conocidos como variaciones de marea. Ellas tienen una máxima amplitud de
algunos 3 gu y un periodo mínimo de 12 horas.

r1 = rc + g tide (Ec. 6.3)

Donde: r1 = lecturas corregidas por mareas

rc = lecturas corregidas por factor de escala, a partir de la ecuación 6.2
gtide = corrección por mareas

Si un gravímetro con una alta tasa de deriva es usado, el amarre de bases se hace en un
intervalo mucho más pequeño que el periodo de la marea y las variaciones de marea son
automáticamente removidas durante las correcciones de deriva. Si un medidor con una
baja tasa de deriva es usado, el amarre de bases es normalmente hecho únicamente al
inicio y al final del día, así que las variaciones de marea experimentan un ciclo total. En
un caso tal, una corrección separada de marea puede necesitarse. Los efectos de marea
son predecibles y publicados todos los años en la prensa geofísica. Esta corrección es
realizada automáticamente por el Software incorporado en el gravímetro SCINTREX CG3.

Altura Instrumental

Cada lectura es corregida por efecto de la altura instrumental por encima de la estación o
base a la cual la lectura gravimétrica es obtenida. Esta es una corrección positiva de
anomalía de aire libre.

rh = rt + 0.308596 hl (Ec. 6.4)

Donde rh = medida corregida por altura instrumental

rt = lectura corregida por efecto de mareas
hI = altura instrumental en metros
 6-28

Todas las mediciones de gravedad se realizaron colocando el gravímetro sobre el trípode,

de manera que la altura instrumental permaneció constante durante todo el programa de
adquisición, haciendo innecesario la realización de la corrección de la altura instrumental.

Corrección de Deriva

Las correcciones por deriva instrumental están basadas en la repetición de lecturas en

una estación base a lo largo del día. Las lecturas son graficadas y se asume que la deriva
es lineal entre lecturas de dos bases consecutivas. La corrección en el tiempo t es d, la
cual es sustraída del valor observado.

 ( r − r ) − ( gb 2 − g b1 ) 
d =  b 2 b1  (Ec. 6.5)
 tb2 − tb1 
Donde: d = deriva instrumental en miligales/hora
rb1 = lectura en la base 1
tb1 = tiempo en la base 1
gb1 = valor absoluto de gravedad en miligales en la base 1
rb2 = lectura en la base 2
tb2 = tiempo en la base 2
gb2 = valor absoluto de la gravedad en miligales en la base 2
Después de la corrección por deriva, la diferencia de gravedad entre un punto de
observación y la base es encontrada multiplicando la diferencia de la medición por el
factor de calibración del gravímetro. Conociendo las diferencias en gravedad, la gravedad
absoluta en el punto de observación gobs puede ser calculada a partir del valor de un punto
conocido en la base. Alternativamente, las lecturas pueden ser relacionadas a un datum
arbitrario, pero esta práctica no es conveniente ya que los resultados a partir de diferentes
servicios no pueden entonces ser relacionados. La deriva instrumental es corregida
automáticamente por el gravímetro SCINTREX CG3.

Gravedad Observada

La gravedad observada es la atracción gravitacional terrestre en el punto de medida,

después de aplicadas las correcciones de efecto de mareas, altura y deriva instrumental a
partir de la estación base.
g obs = g b1 + (rh − rb1 ) − (t − tb1 )d (Ec. 6.6)

Donde gobs = valor de gravedad observada, en miligales

gb1 = gravedad absoluta en miligales en la base 1
rh = lectura de estación corregida por altura instrumental
rb1 = lectura en la base 1
t = tiempo de lectura
 6-29

tb1 = tiempo en la base 1

d = deriva instrumental en miligales/hora
Una vez realizadas las correcciones anteriores se ha ajustado el valor de cada estación
de medida a valores de gravedad absoluta, a partir de las diferencias establecidas en el
ciclo de amarre a la estación gravimétrica situada en el estribo derecho del puente sobre
la quebrada La Chinita en Apartadó, Antioquia.

Corrección de Latitud

La gravedad varía con la latitud debido a que la Tierra no es exactamente de forma

esférica y debido a que la velocidad angular de un punto sobre la superficie de la Tierra
decrece desde un máximo en el ecuador a cero en los polos. La aceleración centrípeta
generada por la rotación de la Tierra tiene una componente radial negativa que
consecuentemente causa que la gravedad decrezca a partir del polo hasta el ecuador. La
verdadera forma de la Tierra es un esferoide aplanado o elipsoide aplanado en los polos,
con una diferencia en longitud de los radios ecuatorial y polar de unos 21 km.
Consecuentemente, los puntos cercanos al ecuador están más lejos del centro de masa
de la Tierra que aquellos cercanos a los polos. La amplitud de este efecto es reducido por
la diferencia en el subsuelo de la distribución de masas resultante del abultamiento
ecuatorial, ya que la masa subyacente en las regiones ecuatoriales es más grande que la
subyacente en las regiones polares. El efecto neto de estos factores es que la gravedad
en los polos excede la gravedad en el ecuador por unos 5186 miligales. La formula de
Clairaut relaciona la gravedad a la latitud sobre el esferoide de referencia de acuerdo con
una ecuación de la forma:

gφ = g 0 (1 + k1sen 2φ − k 2 sen 2 2φ ) (Ec. 6.7)

Donde gφ es el valor predicho del valor de gravedad a la latitud φ, g0 es el valor de la

gravedad en el ecuador y k1 y k2 son constantes que dependen de la forma y rapidez de la
rotación de la Tierra. La ecuación es de hecho una aproximación de una serie infinita.
Antes de 1967 fueron empleadas constantes menos precisas en la Fórmula Internacional
de 1930:

gφ = 978049 [1 + 0,0052884sen 2 (φ ) − 0,0000059sen 2 (2φ )] (Ec. 6.8)

Los valores de g0, k1 y k2 usados en la Fórmula Internacional de Gravedad de 1967 son:

gφ = 978031,846 [1 + 0,005278895sen 2 (φ ) + 0,000023462sen 4 (φ )] (Ec. 6.9)

Una forma más precisa de presentación de la Fórmula de Gravedad 1967 (Mittermayer

1969), es aquella en la cual las constantes son ajustadas para minimizar errores
resultantes de truncar la serie:

gφ = 978032,7 [1 + 0,0053024sen 2 (φ ) + 0,0000058sen 4 (φ )] (Ec. 6.10)

La corrección de latitud requiere la gravedad teórica en el sitio de la estación sobre el

esferoide terrestre. El valor gφ da el valor pronosticado de gravedad al nivel del mar en
 6-30

cualquier punto de la superficie de la Tierra y es sustraído de la gravedad observada a

corregir por variaciones de latitud.

Anomalías de Aire Libre

La corrección de aire libre se aplica debido a que las estaciones de gravimetría tienen
diferentes elevaciones y corrige el decrecimiento de la gravedad con la altura, resultante
de un incremento de la distancia desde el centro de la Tierra a la estación de medida, de
acuerdo con la Ley de Newton. Esta corrección es positiva por encima del datum utilizado,
en este caso, la elevación sobre el nivel del mar. La anomalía de Aire Libre remueve la
corrección por latitud y calcula la elevación de la estación por encima del esferoide.

g faa = g obs − g t + 0,308596hs (Ec. 6.11)

Donde gfaa = anomalía de aire libre, en miligales

gobs = gravedad observada, en miligales
gt = corrección por latitud
hS = elevación de la estación, en metros

La corrección de aire libre cuenta solo para variaciones en la distancia del punto de
observación al centro de la Tierra; no toma en cuenta los efectos de la roca presente entre
el punto de observación y el datum. La corrección de aire libre es la diferencia entre la
gravedad medida a nivel del mar y la medida a una elevación h sobre el nivel del mar, sin
tener en cuenta la masa de roca que se encuentra en medio.

La anomalía de aire libre obtenida a partir de los datos gravimétricos terrestres de la zona
de estudio se presenta en la figura 6.12. Las mayores diferencias positivas de gravedad
(anomalía de aire libre positiva) se encuentran en la franja de terreno amplia, de forma
arqueada, que cubre desde el sur hasta el noroccidente de la zona. Las anomalías de aire
libre positivas (diferencias positivas de gravedad) van disminuyendo hacia el suroccidente,
nororiente, oriente y suroriente de la zona, hasta pasar a anomalías de aire libre muy
bajas o negativas, en las esquinas nororiental y suroriental de la misma. Los sectores con
anomalías de aire libre positivas bajas o negativas, están separados hacia el suroriente
por un sector relativamente pequeño donde se presentan anomalías de aire libre
moderadas. Los valores estadísticos de la anomalía de aire libre corresponden a un
mínimo de -32.38 miligales, un máximo de 116.95 miligales, una media de 40.16 miligales
y desviación estándar de 25.97 miligales.
 6-31

Figura 6.15 Anomalía de Aire Libre. Valores estadísticos: Mínimo = -32.38 miligales,
Máximo = 116.95 miligales, Media = 40.16 miligales y Desviación Estándar = 25.97 miligales

Anomalía de Bouguer Simple

La anomalía de Bouguer es la base de la interpretación de los datos de gravedad en

tierra, que además provee una amplia evaluación del grado de compensación isostática
de un área (e.g. Bott 1982). La corrección de Bouguer tiene en cuenta el efecto
gravitacional de las rocas presentes entre el punto de observación y el datum. Esta
corrección considera las rocas por debajo del punto de medida como una capa horizontal,
con un espesor igual a la elevación del punto de observación sobre el datum y asignando
un valor de densidad, que en este caso ha sido 2.67 gr/cm3, considerado como el valor
promedio para la corteza terrestre. Esta anomalía es sustraída, ya que la roca presente
entre el punto de observación y el datum presenta una atracción adicional. La corrección
de Bouguer asume que el terreno alrededor del punto de medida es plano. La anomalía
de Bouguer (gab) corrige la anomalía de aire libre por efecto de la masa de roca existente
entre la elevación de la estación y el esferoide.

g ba = g faa − 0,0419088ρ B hs (Ec. 6.12)

Donde gba = Anomalía de Bouguer en miligales

 6-32

gfaa = Anomalía de Aire Libre en miligales

ρB = densidad de Bouguer de la roca en g/cm3
hS = elevación de la estación

Las anomalías de gravedad son convencionalmente desplegadas en perfiles o en mapas

de “isogales”. Los colores y tonalidades de las imágenes pueden revelar rasgos
estructurales que generalmente no son discernibles en mapas no procesados. La
anomalía de Bouguer obtenida a partir de los datos gravimétricos terrestres de la zona de
estudio se presenta en la Figura 6.13. Se aprecia un alto gravimétrico en forma de arco,
que se extiende desde el sur al noroccidente de la zona, alcanzando valores por encima
de 110 miligales, el cual es generado por la cercanía a la superficie de las rocas del arco
de Dabeiba (o arco de Sautatá). A parir del alto gravimétrico, la anomalía de Bouguer va
disminuyendo hasta llegar a valores negativos en las esquinas nororiental y suroriental de
la zona, alcanzando valores por debajo de -10 miligales, generados por la profundización
del basamento de la cuenca de Urabá y la presencia de rocas sedimentarias. Los valores
estadísticos de la anomalía de Bouguer corresponden a un mínimo de -46.29 miligales, un
máximo de 115.03 miligales, media de 36.55 miligales y desviación estándar de
27.28 miligales.

Figura 6.16 Anomalía de Bouguer simple. Valores estadísticos de la anomalía de Bouguer:

Mínimo = -46.29 miligales, Máximo = 115.03 miligales, Media = 36.55 miligales y Desviación
Estándar = 27.28 miligales.
 6-33

Anomalía de Bouguer Total

La corrección por efecto del terreno se realiza con el fin de considerar las irregularidades
topográficas presentes en las vecindades del punto de medida. Esta corrección es
siempre positiva. Un relieve positivo en las vecindades de un punto de medición ejerce
una atracción adicional hacia arriba, causando una disminución de la gravedad. Esta
atracción es eliminada mediante una corrección positiva de terreno. Para el caso de
relieves negativos, en vecindades del sitio de medida, se asume que estas depresiones
no existen en la corrección de Bouguer, considerando una capa plana de extensión
infinita; por lo tanto se ha sobre corregido el efecto de masas de roca, que no existe
completamente, de manera que se debe corregir esta sobre estimación mediante una
corrección positiva de terreno.
g cba = g ba + g tc (Ec. 6.13)

Donde gcba = Anomalía de Bouguer total, en miligales

gba = Anomalía de Bouguer, en miligales
gtc = Correcciones por efectos del terreno, en miligales

La corrección de Bouguer asume que la topografía alrededor de la estación de gravedad

es plana. Esto raramente ocurre y por lo tanto debe realizarse la corrección de terreno
para corregir los efectos del relieve topográfico vecino a la estación de gravedad. Esta
corrección es siempre positiva.

Los efectos de terreno son bajos en áreas de topografía baja y raramente exceden
1 miligal (10 gu) en áreas planas. En áreas de topografía montañosa los efectos de
terreno son considerablemente grandes, teniendo máximos en los flancos empinados de
los valles y en la base y cimas de las montañas. Donde los efectos de terreno son
considerablemente menores que la precisión deseada en la exploración, las correcciones
de terreno pueden ser ignoradas. Sin embargo, la usual necesidad de aplicar esta
corrección demanda un tiempo importante en la reducción de los datos de gravedad y
constituye una importante contribución al costo de una exploración gravimétrica.

La anomalía de Bouguer total obtenida a partir de los datos gravimétricos terrestres de la

zona de estudio se presenta en la Figura 6.14. Como en el caso de la anomalía de
Bouguer simple, se aprecia un alto gravimétrico que sobrepasa los 110 miligales, el cual
es generado por las rocas graníticas del arco de Dabeiba. De igual forma se aprecian dos
bajos gravimétricos situados en las esquinas nororiental y suroriental de la figura, con
valores por debajo de -10 miligales, generados por el descenso del basamento de la
cuenca de Urabá y la presencia rocas sedimentarias. Las anomalías de Bouguer simple y
de Bouguer total son similares, ya que los datos gravimétricos fueron adquiridos en una
zona prácticamente plana, retirados de las zonas montañosas adyacentes
correspondientes a las serranías del Darién, Baudó y Abibé. Los valores estadísticos de la
anomalía de Bouguer corresponden a un mínimo de -44.02 miligales, un máximo de
115.15 miligales, media de 37.15 miligales y desviación estándar de 27.04 miligales.
 6-34

Figura 6.17 Anomalía de Bouguer total. Valores estadísticos de la anomalía de Bouguer:

Mínimo = -44.02 miligales, Máximo = 115.15 miligales, Media = 37.15 miligales y
Desviación Estándar = 27.04 miligales

6.4.2 Magnetometría

El alma de los estudios magnéticos es investigar la geología del subsuelo sobre la base
de las anomalías del campo magnético de la Tierra resultante de las propiedades
magnéticas de las rocas subyacentes. Aunque la mayoría de los minerales constituyentes
de las rocas no son magnéticos, ciertos tipos de roca contienen suficiente minerales
magnéticos para producir anomalías magnéticas significativas. Similarmente, objetos
ferrosos hechos por el hombre también producen anomalías magnéticas. La exploración
magnética tiene un amplio rango de aplicaciones, desde servicios a pequeña escala
ingenieriles o arqueológicos para detectar objetos metálicos enterrados, hasta servicios a
gran escala adelantados para investigar estructuras geológicas regionales. Las
exploraciones magnéticas pueden ser realizadas en tierra, en el mar y en el aire.
Consecuentemente la técnica es ampliamente empleada y la rapidez de operación de los
servicios aéreos, hace que el método sea muy atractivo en la búsqueda de depósitos
minerales que contienen minerales magnéticos.
 6-35

El método magnético mide las variaciones del campo magnético terrestre debido al
contraste del grado de magnetización o susceptibilidad magnética de las rocas que
conforman la corteza terrestre. La prospección magnética, considerado el método más
antiguo en la exploración geofísica, se emplea en la búsqueda de petróleo y de minerales.
En la prospección de petróleo es empleado regularmente para determinar el espesor de la
secuencia sedimentaria o para cartografiar rasgos estructurales de la superficie del
basamento que podrían influenciar la estructuras suprayacentes. Las rocas sedimentarias
ejercen un efecto magnético tan pequeño, en comparación con las rocas ígneas y
metamórficas infrayacentes, que virtualmente todas las variaciones de la intensidad
magnética medibles en la superficie están asociadas a la topografía o a los cambios
litológicos del basamento.

6.4.2.1 Susceptibilidad Magnética

La susceptibilidad magnética es una propiedad intrínseca de la roca y expresa la cantidad

de magnetismo que una roca puede tener ante la presencia de un campo magnético
determinado. Las rocas están compuestas de minerales, los cuales varían en su
susceptibilidad magnética y en sus cantidades de magnetita. Las medidas de
susceptibilidad pueden ser usadas antes de los levantamientos magnetométricos, para
tener mayores elementos de juicio con respecto al conocimiento de cuales rocas pueden
ser detectables magnéticamente y en que grado. Esto puede realizarse en afloramiento, o
en el laboratorio. El grado de magnetización (“POLARIZACIÓN”), por polarización
inducida es el producto de su susceptibilidad k y el campo magnético actuante.
I = kH (Ec. 6.14)

H es el campo magnético terrestre (0.3 a 0.6 oersteds)

En general, la susceptibilidad de las rocas es una medida de su contendido de magnetita.

Para bajas concentraciones de magnetita existe una aproximada relación lineal entre el
porcentaje de ésta y k, la cual puede expresarse como:

k = 0,30 ρ (Ec. 6.15)

Donde ρ es % por volumen de magnetita

En el Cuadro 6.6 se relacionan valores de susceptibilidad magnética de algunas clases de
rocas, en unidades c.g.s y SI.
Es importante conocer las susceptibilidades magnéticas de las rocas que constituyen las
unidades litoestratigráficas de la zona de trabajo, para el modelamiento inverso de los
datos de magnetometría. Con un susceptibilimetro SM-30 del Departamento de
Geociencias de la Universidad Nacional de Colombia, se midió la susceptibilidad
magnética a aproximadamente 77 muestras de roca tomadas en afloramientos de la
Serranía de Abibé y en las costas del golfo de Urabá (Fotografía 6.3). En el Cuadro 6.7 se
relaciona los valores de susceptibilidad magnética en unidades SI, medidos en
laboratorio.
 6-36

Cuadro 6.6 Susceptibilidades magnéticas de algunas rocas comunes (Tomado de Reynolds,

1998)
Litología Rango – cgs Rango – si
Peridotita 7599.63 – 15597.15 95500 – 196000
Basalto oceánico 23.87 – 2864.78 300 – 36000
Basalto 39.79 – 14483.07 500 – 182000
Gabro 63.66 – 6047.87 800 – 76000
Pegmatita 238.73 – 5968.30 3000 – 75000
Riolita 19.89 – 3000.06 250 – 37700
Granito con magnetita 1.59 – 3978.86 20 – 50000
Granito 0.80 – 5.17 10 – 65
Gneiss 9.95 – 1989.43 125 – 25000
Pizarra 0 – 3023.94 0 – 38000
Esquisto 25.07 – 238.73 315 – 3000
Shale 4.78 – 1480.14 60 – 18600
Arenisca 0 – 1671.12 0 – 21000
Caliza 0.80 – 1989.43 10 – 25000
Dolomita impura 1591.55 20000
Dolomita pura -0.99 – 3.50 -12.5 a +44

Cuadro 6.7 Valores de susceptibilidad magnética de muestras de campo medidas en

laboratorio
Susceptibilidad Magnética
Periodo Unidad Litoestratigráfica Convención
Intervalo Promedio
Cuaternario Depósitos Volcanes de Lodo Qvl 0.00001 0.00086 0.00020
Cuaternario Depósitos Aluviales Qal 0.00002 0.00048 0.00018
Neógeno Formación Corpa Ngco 0.00001 0.00258 0.00025
Neógeno Formación Arenas Monas Ngam 0.00000 0.00362 0.00038
Neógeno Formación Morrocoy Ngmp 0.00007 0.00014 0.00010
Neógeno Formación Paujil Superior Ngpas 0.00000 0.00117 0.00026
Neógeno Formación Paujil Inferior Ngpai 0.00002 0.00017 0.00011
Neógeno Formación Pavo Superior Ngps 0.00000 0.01920 0.00155
Neógeno Formación Pavo Inferior Ngpi 0.00001 0.02570 0.00213
Unidades Sedimentarias 0.00002 0.00751 0.00068
Paleógeno Batolito de Mandé Etm 0.00789 0.09500 0.02885
Complejo Santa Cecilia-
Cretácico Ksvx 0.00065 0.05600 0.01726
La Equis

En general se observan datos de susceptibilidad altos para las muestras de roca del
Batolito de Mande, incluso por encima de las muestras de roca del Complejo Santa
Cecilia-La Equis. Teniendo en cuenta que esta medición es un indicador directo del
 6-37

volumen de magnetita en la roca, la información muestra que hay altos niveles en el

volumen de este mineral y en general de hierro en las rocas del Batolito de Mandé.

Fotografía 6.3 Toma de susceptibilidad magnética a muestra de roca

6.4.2.2 Diseño de Adquisición Magnetométrica

La adquisición magnetométrica se realizó a lo largo de tres transectas de orientación
general EW, aprovechando tres vías carreteables localizadas sobre la llanura aluvial del
valle bajo del río Atrato y piedemonte occidental de la Serranía de Abibé. La primera
transecta con una longitud de 70 km, se realizo a lo largo de la vía Caucheras-Belén de
Bajirá-Riosucio; la segunda transecta con una longitud de 40 km se realizó a lo largo de la
vía El Tigre-El Cuarenta y la tercera transecta con una longitud de 25 km se realizó a lo
largo de la vía Vereda El Silencio-Carepa-Piedras Blancas (Figura 6.18). Las estaciones
de adquisición tenían una separación de 500 m. Las coordenadas de los puntos extremos
de las transectas magnetométricas se relacionan en el Cuadro 6.8:

Cuadro 6.8 Coordenadas de los puntos extremos de las transectas de magnetometría

Coordenadas Punto Inicial Coordenadas Punto Final
Transecta Estac Long
Este Norte Elev Este Norte Elev
Caucheras-Riosucio 734767 1303483 19 671426 1312633 24 166 64
El Tigre-El Cuarenta 718053 1331186 52 680931 1338973 23 91 38
Piedras Blancas-
734597 1351074 70 700461 1351075 25 102 34
Carepa-El Silencio
Totales 136

Se establecieron dos bases de control magnetométrico ubicadas respectivamente en la

Finca La Chinita, a las afueras de Apartadó y en un predio ubicado similarmente a las
afueras de Belén de Bajirá. Por fallas del magnetómetro base, fue necesario recurrir a una
 6-38

tercera base magnetométrica en cercanías de la población de Turbo. Las coordenadas de

las bases magnetométricas se relacionan en el Cuadro 6.9:

Cuadro 6.9 Coordenadas de las bases magnetométricas

Base Magnetométrica Este Norte Elevación (m)
La Chinita 718441.79 1365829.48 34
Bajirá 707170.89 1307661.08 29

6.4.2.3 Adquisición Magnetométrica

La adquisición de datos de magnetometría se realizó entre el 28 de Enero y el 13 de
Febrero de 2009. Se realizaron tres perfiles mediante circuitos cerrados, es decir,
refiriendo las lecturas de un grupo de estaciones móviles a una sola estación base y
cerrando el circuito en la misma estación base donde se hizo la primera lectura al
comenzar el ciclo. Las mediciones de magnetometría se realizaron utilizando dos
magnetómetros de precesión protónica GEOMETRICS G-856, de sensibilidad 0.1 nT
(Fotografías 6.4 y 6.5).

Fotografía 6.4 Magnetométros GEOMETRICS Fotografía 6.5 Magnetómetro y aditamentos

G-856

Para las transectas de El Tigre-El Cuarenta y Piedras Blancas-Carepa-El Silencio se

utilizó la base magnetométrica La Chinita en Apartadó. La transecta Caucheras-Bajirá-
Riosucio usó la base magnetométrica de Bajirá (Fotografía 6.6). Debido a que el
magnetómetro base presentó problemas en su funcionamiento, fue necesario apoyarse en
la estación base ubicada en Turbo, del trabajo de adquisición adelantado por Geokinetics.
La adquisición comprendió el registró del campo magnético en un total de 359 estaciones
(Fotografía 6.7).
 6-39

Figura 6.18 Localización de transectas y bases de magnetometría

 6-40

Fotografía 6.6 Estación base Bajirá Fotografía 6.7 Toma de datos

magnetométricos

6.4.2.4 Procesamiento de Datos Magnetométricos

El procesamiento comprende la reducción de los datos magnéticos para remover todos
aquellos efectos magnéticos provenientes de fuentes externas a la corteza terrestre. Para
producir el mapa de anomalía magnética de una región, los datos magnéticos adquiridos
deben ser corregidos teniendo en cuenta las variaciones diurnas del campo magnético,
los efectos de la latitud, y en menor extensión la longitud, y la sustracción del campo
magnético terrestre teórico (Reynolds, 1997).

Corrección por variaciones diurnas

Los datos adquiridos en los estudios magnéticos requieren de varias correcciones. Un
efecto que puede ser compensado es la variación de intensidad del campo geomagnético
en la superficie de la Tierra durante el transcurso del día. Las variaciones diurnas del
campo geomagnético son debidas a que parte del campo magnético de la Tierra se
origina en la ionosfera. Las variaciones diurnas pueden ser corregidas mediante la
instalación de un magnetómetro que registre constantemente en una estación base fija
dentro del área de estudio. Alternativamente, los registros de un observatorio
geomagnético pueden ser usados, teniendo en cuenta que no se encuentre a una
distancia no mayor de 150 km del área de estudio. Los efectos de las variaciones diurnas
pueden ser corregidos o removidos de varias maneras.

Corrección geomagnética
La corrección geomagnética remueve los efectos del campo magnético terrestre de
referencia de los datos adquiridos. El método de corrección geomagnética requiere el uso
del IGRF, el cual expresa el campo magnético terrestre no perturbado en términos de un
gran número de armónicos e incluye términos temporales para las correcciones por
variaciones seculares.

Las variaciones del campo magnético con la altitud, latitud y longitud son dominadas por
las variaciones vertical y horizontal del campo del dipolo. La intensidad total del campo
 6-41

magnético Bt es obtenida por la sumatoria de sus componentes radial (Br) y tangencial

(Bθ):

 
 µ 0 m   1 + 3 cos θ 
2
2 2
Bt = Br + Bθ =   (Ec. 6.16)
 4π   r3 
 

Corrección por latitud, longitud y altitud

Como el campo magnético de la Tierra varía desde 25000 nT a en el ecuador magnético a
69000 en los polos, el incremento en magnitud con la latitud debe tenerse en cuenta.

En estudios regionales la corrección por latitud y longitud son inherentes al campo de

referencia que es sustraído. En estudios de regiones pequeñas, la corrección de latitud
es dada por el gradiente horizontal en dirección norte-sur del campo magnético, obtenido
por la diferenciación de Bt con respecto al ángulo polar θ. Esto da que Bt incrementa hacia
el norte, con el incremento de la latitud.

 2  3B senθ cos θ
 1   ∂B   µ 0 m   1   ∂ 1 + 3 cos θ = t
−   t  
= 
    r (1 + 3 cos 2 θ ) (Ec. 6.17)
 r   ∂r   4π   r 4   ∂θ
 
La corrección por latitud es cero en el polo magnético (θ = 0°) y noventa en el ecuador
magnético (θ = 90°), alcanzando un máximo valor de cerca de 5 nT por kilómetro (0.005
nT/m) en las latitudes intermedias. Esto es insignificante en estudios de pequeña escala.
La corrección de altitud es dada por el gradiente vertical del campo magnético, obtenido
por la diferenciación de la intensidad Bt respecto al radio r:

 2 
 µ m  1 + 3 cos θ  = −  3Bt 
= −3  0  
∂Bt
   (Ec. 6.18)
∂r  4π   r4  r 
 
El gradiente vertical del campo se encuentra sustituyendo r = R = 6371 km y un apropiado
valor del campo magnético total (Bt), que depende de la latitud del sitio de medida. En el
ecuador magnético, el campo magnético total alcanza un valor de 30000 nT y la
corrección de altitud es cerca de 0.015 nT/m, cerca de los polos el campo magnético total
tiene un valor de 60000 nT y la corrección de altitud es aproximadamente 0.030 nT/m. La
corrección por altitud es tan pequeña, que normalmente no se aplica. La influencia de la
topografía puede ser significativa en estudios magnéticos de grandes regiones, pero no es
completamente predecible si este depende de las propiedades magnéticas de los rasgos
topográficos. Habiendo aplicado las correcciones diurnas y geomagnéticas, todas las
variaciones del campo magnético remanente es causado solamente por las variaciones
espaciales en las propiedades magnéticas del subsuelo y estas son referidas como
anomalías magnéticas.

Reducciones al polo y al ecuador magnéticos

La posición y forma de las anomalías magnéticas depende no solo de las características
físicas del cuerpo que las produce, sino de su orientación y posición dentro del campo
 6-42

geomagnético (asumiendo solamente magnetización inducida). En los polos y en el

ecuador magnéticos, las anomalías se localizan encima del cuerpo causativo, pero en las
latitudes intermedias, estas poseen un máximo y un mínimo y el cuerpo se localiza
aproximadamente bajo el punto de inflexión. Para corregir la posición de las anomalías
magnéticas, de tal forma que coincida con la posición del cuerpo causativo, siempre y
cuando el magnetismo remanente sea despreciable, se aplica la reducción al polo o al
ecuador magnéticos. En el caso de las latitudes magnéticas intermedias, se aplica la
reducción al Polo (Carson Aerogravity, 2006), como es el caso de la zona de estudio,
donde el campo geomagnético presenta una inclinación magnética de 34° al norte y una
declinación magnética de -4.75° al oeste.
La reducción al polo consiste en la conversión de las anomalías observadas a su forma
equivalente en el polo norte magnético. Este proceso generalmente permite ver la
anomalía magnética como si estuviera en un campo magnético vertical. El programa de
Oasis Montaj cuenta con la herramienta MAGMAP para aplicar este filtro a una retícula de
intensidad magnética total. El filtro es aplicado en el dominio de Fourier para migrar el
campo observado a partir de la inclinación y la declinación magnética observadas, de
manera tal que el campo magnético pueda ser visto como en el polo magnético. La
reducción al polo ayuda a la interpretación de anomalías magnéticas, ya que cualquier
asimetría en el campo reducido al polo puede ser atribuida a la geometría y propiedades
magnéticas de la fuente. Una limitación de la aplicación del filtro de reducción al polo es
que este puede introducir ruido a la retícula resultante cuando esta contiene anomalías de
alta frecuencia y alta amplitud.
Por implementación de filtros tanto al espectro de amplitud y de fase de la retícula original,
las formas de las anomalías magnéticas puede ser simplificada a como ellas aparecen
como la relativamente simple anomalía positiva principalmente localizada directamente
sobre la fuente esperada para cuerpos inducidos magnéticamente en los polos
magnéticos (I = +90°). Esta fallará parcialmente donde la magnetización remanente esté
presente y las dificultades teóricas serán encontradas cuando probando la reducción al
polo de una retícula de observaciones en inclinaciones magnéticas muy bajas.
La anomalía de campo magnético total se presenta en la Figura 6.19. Se aprecia una gran
anomalía magnética positiva que cubre un amplio sector al occidente y noroccidente de la
zona, una anomalía magnética alta pero de pequeña extensión en la parte sur y una
anomalía magnética positiva moderada cubriendo los sectores norte y nororiental. De
igual manera se aprecia anomalías magnéticas negativas en los sectores central, sur,
suroccidental y suroriental de la zona. Los valores estadísticos obtenidos de la retícula de
campo magnético total corresponden a un mínimo de -712.88 nT, un máximo de 325 nT,
media de -130.54 nT y desviación estándar de 138.02 nT.
 6-43

Figura 6.19 Anomalía magnética de campo total. Valores estadísticos: Mínimo = -712.88 nT,
Máximo = 325 nT, Media = -130.54 nT y Desviación Estándar = 138.02 nT

Las anomalías de campo magnético total reducidas al polo magnético se presentan en la

figura 6.20. Se aprecian anomalías magnéticas positivas que alcanzan valores por encima
de los 200 nT al noroccidente y nororiente separadas por una anomalía magnética
negativa con valores por debajo de -400 nT situada en el sector norte. Se aprecia también
una anomalía magnética positiva que cubre una franja angosta de terreno desde el sur
hasta la parte central de la zona, la cual separa dos anomalías magnéticas negativas con
valores por debajo de -800 nT. Los valores estadísticos obtenidos de la retícula de las
anomalías de campo magnético total reducidas al polo corresponden a un mínimo
de -1439.97 nT, un máximo de 1350.38 nT, media de -133.80 nT y desviación estándar de
297.68 nT.
 6-44

Figura 6.20 Anomalía magnética reducida al polo. Valores estadísticos: Mínimo = -1439.97nT,
Máximo = 1350.38 nT, Media = -133.80 nT y Desviación Estándar = 297.68 nT