UNIVERSIDAD NACIONAL DE TUCUMAN

FACULTAD DE CIENCIAS EXACTAS Y TECNOLOGÍA

DEPARTAMENTO DE GEODESIA Y TOPOGRAFÍA

CATEDRA DE GEOFÍSICA

APUNTES DE

PROSPECCIÓN
MAGNÉTICA

PARA ALUMNOS DE INGENIERÍA GEODÉSICA Y GEOFÍSICA

DE LA FACULTAD DE CIENCIAS EXACTAS Y TECNOLOGÍA
DE LA UNIVERSIDAD NACIONAL DE TUCUMAN

Prof. Ing. Luis A. Estrada

Año 2015
Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Introducción

Ventajas y desventajas del método.

La primera desventaja es el hecho que se desprecia la imantación remanente, la que a veces
puede ser muy importante. Esto da un cierto grado de incertidumbre a la interpretación.
Otra es el amplio rango de variación de la susceptibilidad magnética para los distintos tipos de
rocas, y no hay garantías de que la magnetización esté uniformemente distribuida. Recordemos
que en un flujo basáltico la susceptibilidad es bastante alta debido a la abundante magnetita.
Sin embargo, esta se acumula a veces en bolsones que dan anomalías discontinuas cuando el
relevamiento es a gran escala o de gran detalle.
En síntesis, las desventajas del método provienen de la fuerte dependencia de la anomalía
respecto de las características propias de la anomalía y de la dirección de magnetización.
Pero la gran ventaja del método está en el relativamente bajo costo de exploración por área de
investigación. Muy especialmente cuando se trata de relevamientos aéreos. Por ello es
generalmente el primero de los métodos que se utiliza para delimitar zonas de interés, y
fundamentalmente ahorrar recursos en el uso de otros métodos más costosos.
El desconocimiento de la verdadera orientación y las propiedades de la anomalía es una
desventaja del método, pero no es grave porque es constante en toda el área de exploración.
Puesto que las susceptibilidades son tan bajas en la mayoría de las rocas, es mejor focalizar la
investigación a grandes anomalías, ya que las altas susceptibilidades se encuentran en un
reducido tipo de rocas.
Como veremos, pueden relacionarse los potenciales magnéticos y gravimétricos de manera tal
que los datos magnéticos podrán ser transformados en pseudo-gravimétricos. Concretamente,
un mapa de anomalías magnéticas de componente vertical puede transformarse en uno de
gradiente vertical de anomalías pseudo-gravimétricas. Esto hace menos complejo el análisis y
permite que se comparen con los datos gravimétricos si están disponibles. Una buena
correlación entre ambas anomalías indica una misma fuente para las dos y permite mayor
definición en la interpretación.
Muchas de las técnicas de procesamiento de los datos usadas en gravimetría, fueron
adaptadas para magnetometría, tanto los de 2D, 3D, los de separación de tendencias
regionales, segundas derivadas y los de continuación hacia arriba y hacia abajo.

Susceptibilidad de las rocas

Antes de analizar las anomalías, debemos considerar la susceptiblidiad, que a pesar de ser
adimensional, su valor se expresa en unidades electromagnéticas (uem) en el sistema cgs.
Esto es solo para llamar la atención de que el sistema a usar es el cgs y no el SI. Los valores
pueden ser convertidos al SI con solo multiplicarlos por 4π.
Algunos minerales como la magnetita, ilmenita y pirrotita tienen alta susceptibilidad, siendo la
primera la mayor. Como todas las rocas tienen algo de magnetita, la susceptibilidad de todas
ellas se mide como un porcentaje de esta sustancia. La verdadera susceptibilidad de la
magnetita varía entre 0.1 y 1.0 uem, dependiendo del tamaño del grano, forma e impurezas,
adoptándose como media k=0.35 uem. Entonces, una roca con un 1% de magnetita tendrá una
susceptibilidad k=0.003, es decir 3.10-3 cgs uem.
Los valores medios para diferentes tipos de roca son los siguientes:

Sedimentarias 5.10-5 cgs uem 0.016% magnetita

Metamórficas 3.10-4 cgs uem 0.100% ,,
Granitos y Riolitas 5.10-4 cgs uem 0.160% ,,
Gabros y Basaltos 6.10-3 cgs uem 0.200% ,,
Rocas Ultabásicas 12.10-3 cgs uem 4.000% ,,

Como dijimos, la magnetización remanente Ir no es tenida en cuenta para la interpretación de

las mediciones. Esta, como se sabe, está presente en toda roca conservando la historia
magnética, es decir el magnetismo existente al momento de la formación de la roca. A veces
suele tener una magnitud similar a la magnetización inducida Ii y hasta estar orientada en
una dirección muy diferente.

Ing. Luis Estrada - 2015 2

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

En estos casos la magnetización total, suma de ambas, será muy distinta de la inducida que
intentamos interpretar. En rocas basálticas la relación Ir/Ii suele ser a veces mayor que 10, en
cambio para rocas graníticas es aproximadamente 1, en las metamórficas <1 y en las
sedimentarias <0.1.

Procedimientos básicos de campo

La exploración magnética no requiere de tantos cuidados y correcciones como la gravimétrica.
Sin embargo, hay que seguir un procedimiento para obtener una precisión aceptable en las
mediciones, aunque la exploración desde el aire requiere de otros cuidados adicionales.
La limpieza magnética es fundamental cuando se trabaja con magnetómetros portátiles sobre
el suelo. Esto implica desprenderse de objetos metálicos aparentemente inocentes como
cortaplumas, anteojos, lapiceras, llaves, etc. Además, debe tenerse cuidado de no acercarse
demasiado a automóviles, alambrados, líneas eléctricas, tuberías, etc. El sensor siempre debe
estar a más de un metro del suelo, y si es posible a dos o tres metros, para que no afecten las
lecturas los minerales o desechos magnéticos de la superficie.
A raíz de las variaciones diurnas o de corto período del campo magnético, las mediciones
deben corregirse por sus efectos, al igual que por marea y deriva en gravedad. Una forma
consiste en reocupar una estación base cada una hora y luego cambiar de base. Otra más
práctica pero que requiere de dos magnetómetros, consiste en dejar uno fijo como base para
así conocer las fluctuaciones del campo en un punto y descontarlas en los otros.
Antes de analizar las correcciones es conveniente recordar lo elemental de la ecuación del
dipolo vista en Geomagnetismo.

ya que F/m′′ = H = m/r

2
V = F.r ó V = m/r

H H+ El Potencial V es el trabajo negativo

H-m realizado sobre un polo unitario en un
r1 P campo magnético, entonces
r
-m r
Hθ V = ∫ (m / r2).dr = m / r
H ∞
l/2
θ r2
x
La importancia del Potencial viene del
rdθ dr
l/2 P hecho de que se puede obtener el
Hr campo en una dirección dada:
+m
-m Hr = - dV / dr
r Deducimos la ecuación del dipolo
+m para describir el campo magnético que
genera un dipolo en un punto P.
El Potencial total será: V = (m / r1) – (m / r2)

θ
Reemplazando r1 y r2 en función de r y x que es igual a (l / 2)cosθ

m m θ
m. l.cosθ
V= - =
θ
r - (l /2)cosθ θ
r + (l /2)cosθ θ )2
r2 - (l /2)2(cosθ

θ
M cosθ
Si r >> l ⇒ V ≅ (m. l l cosθ
θ) / r2, entonces, V=
r2

Como en la Prospección Magnética es conveniente calcular las componentes radial y

tangencial, serán también las negativas de las derivadas del Potencial es esas direcciones:

Ing. Luis Estrada - 2015 3

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Hr = - dV / dr θ ) / r3
= (2Mcosθ (θ = 0° ⇒ Hr = ZE = 2M/r3)

θ = (Msenθ
Hθ = - dV / rdθ θ ) / r3 (θ = 90° ⇒ Hθ = HE = M/r3)

Generalmente a las mediciones magnéticas no se les hace corrección por altitud. Calculando
el gradiente vertical desde la ecuación del dipolo, puede verse porqué no es necesaria esta
corrección:
θ / r4 = - 3Hr / r = - 3ZE / r
dZE / dr = dHr / dr = - 6Mcosθ

Por ejemplo, el cambio respecto a la altura en la componente vertical para un campo ZE de

52.950 nT, a una latitud de 42° resulta de 0,025 nT/m, que si bien es dependiente de la
posición, nunca alcanza valores que justifiquen la corrección.
Recordemos que la corrección de Aire Libre en gravimetría es de 0,3086mgal/m para un campo
total de 983.000 miligales (≈3/10.000.000). En magnetometría la variación por altura es del
mismo orden, es decir de 0,03nT/m para un campo de 70.000nT (=3/7.000.000). A pesar de
que la relación es muy parecida no se efectúa la corrección. La explicación entonces se debe a
que los 0,3mgal/m corresponden a una variación del Ecuador al Polo de unos 5.200 miligales,
en cambio los 0,03nT/m corresponden a una variación de 30.000nT. La densidad de estas
rocas es bastante parecida (2,60 a 2,80 Tn/m3), mientras la susceptibilidad magnética de una
roca sedimentaria es del orden de 0,0003 uem y la del basalto de 0,03 uem. Además, la
contribución magnética de este tipo de rocas en el campo magnético inducido es muy variable.
Todo esto hace que no valga la pena realizar la corrección.
En cuanto a la corrección por posición, corresponde realizarla pues resulta de un orden
significativo. La forma de calcularla es similar a la de altitud.

θ = (1/r)ZE / dθ
dHr / rdθ θ = (1/r) (-2Msenθ
θ) / r3 = - 2Msenθ
θ / r4 = - 2HE / r

Utilizando los mismos valores del ejemplo anterior, el cambio respecto a la posición horizontal
en la componente vertical para un campo HE de 18.200nT, resulta de 6nT/km, que es bastante
si se tiene en cuenta que medimos con una precisión de un nanoTesla.
72°45´ 72°30´
42°30´ La solución para realizar esta corrección es consultar
55.433nT 55.387nT
los cuatro valores extremos del área de trabajo en el
Campo Internacional Geomagnético de Referencia y
calcular la variación en latitud y longitud. En la figura
que sigue se muestra un área de medición de 15´ en
longitud por 30´ en latitud, con los valores del campo
7.2nT por minuto de latitud geomagnético calculados por el IGRF.
Cada punto de observación dentro del área tiene su
posición conocida en latitud y longitud.
3.1nT por minuto de longitud Designamos como valor del campo no perturbado FE
al de la esquina inferior izquierda (55.172nT) y desde
ahí el campo aumenta en la proporción indicada en la
55.218nT 55.172nT figura. De esta forma obtenemos los valores teóricos
42°00´
de todos los puntos del levantamiento.

Como ejemplo supongamos que medimos un campo total FET = 55.193nT en un punto ubicado
a 27 segundos (0.45´) al Norte y 10.8 segundos (0.18´) al Oeste del punto de base. El valor
medido y corregido será: FE = 55.172nT + 0.45x7.2nT + 0.18x3.1nT = 55.176nT
Ahora podemos restar el valor observado (55.193nT) del valor teórico (55.176nT) lo que nos da
una anomalía de 17nT. Por supuesto, si estamos midiendo la intensidad de campo vertical,
seguimos el mismo procedimiento, pero utilizando valores de componente vertical del IGRF.
Cabe señalar que como estamos midiendo valores absolutos del campo magnético, no es
necesario ligar las mediciones a una estación base como hicimos en gravimetría, donde
medimos valores relativos. Todo esto puede evitarse si contamos con GPS, pues conociendo la
latitud en cada punto, podemos obtener el campo teórico con el IGRF.

Ing. Luis Estrada - 2015 4

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Anomalías de Campo Total

Como primer paso debe definirse la terminología a usar: FE, HE y ZE son, como se vio, el campo
total y las componentes horizontal y vertical. De igual modo, FA y HA, ZA serán los elementos
correspondientes al campo magnético inducido por el campo terrestre a un determinado cuerpo
geológico, al que se le llama cuerpo anómalo y cuya geometría se trata de determinar.
En cualquier punto de la superficie terrestre donde se mide el campo magnético, el resultado
incluye tanto el campo principal como el anómalo. Como se ve en la figura, el magnetómetro
medirá 55.005 nT que es la suma de 55.000 nT del campo principal más la contribución de 5
nT del campo anómalo en la dirección del campo principal.
Obviamente, el problema que se presenta es determinar cuál es el campo anómalo total.
Primero partimos del supuesto que FE >> FA y que ambas se encuentran en el mismo plano
vertical. Esto permite encontrar ecuaciones que vinculen los elementos medidos con el campo
anómalo total.

Norte
α
HE
I
HA
FE= 55.000 nT
FE
ZE

FAT = 5 nT FET
ZA FAT
FA = 12 nT
Z
Oeste
Vertical

FA

(FE+FAT)2 = (ZE+ZA)2 + (HE+HA)2

FE2+2FATFE+FAT2 = ZE2+2ZEZA+ZA2+HE2+2HEHA+HA2

Suponiendo FE>>FA podemos despreciar FAT2, HA2 y ZA2, por lo tanto

FE2 + 2FATFE = ZE2 + 2ZEZA + HE2 + 2HEHA

Como FE2 = ZE2 + HE2 nos queda que FATFE = ZEZA + HEHA , ó

FAT = ZA(ZE/FE) + HA(HE/FE) ó

FAT = ZAsen i + HAcos i (#)

Si HA no estuviera sobre el meridiano magnético,

FAT = ZAsen i + HAcos i.cos α

Ing. Luis Estrada - 2015 5

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Efecto magnético de cuerpos simples

El cálculo e interpretación de las anomalías magnéticas es mucho más compleja que las
gravimétricas. Esto se debe fundamentalmente a que todo cuerpo anómalo tiene dos polos de
distinto signo, mientras que en gravedad podemos considerar que hay solo uno positivo o
negativo. Se debe calcular la anomalía de campo total, de modo que en el proceso de
reducción de los datos se pueda quitar la contribución del campo principal. Además la
magnetización remanente, que generalmente no es considerada, puede causar a veces un
efecto muy significativo.
Efecto magnético de un polo aislado. Es el primer +X
paso que debemos analizar para entender la Norte Magnético
interpretación magnética. Aunque el monopolo no
existe en la realidad, es un ejercicio conveniente por
su simplicidad. y
P
Partimos de las fórmulas ya conocidas, donde V es el
potencial magnético en P debido a un polo aislado de x c
intensidad m para un área A de un cuerpo con suscep-
tibilidad k dentro un campo terrestre FE.
r +Y
V = m/r z

I = kFE -m

m = IA = kFE A y r = (c2+z2)1/2 = (x2+y2+z2)1/2

k.FE .A
V=
(x2+y2+z2)1/2

Como el campo magnético en una dirección dada es la derivada negativa del Potencial en esa
dirección:

dV 2z(-1/2)(kFE A) z (kFE A)
ZA = = =
dz (x2+y2+z2)3/2 (x2+y2+z2)3/2

Siempre es conveniente orientar el sistema coordenado, de modo que la dirección x coincida

con el Norte Magnético. Las componentes x,y,z de la gráfica son positivas.

Determinamos HAx y HAy al igual que ZA

dV x (kFE A)
HAx = =
dx (x2+y2+z2)3/2

dV y (kFE A)
HAy = =
dy (x2+y2+z2)3/2

El campo total anómalo se calcula con la ecuación (#), donde HA es la componente horizontal
en la dirección del Norte Magnético, decir que HA2 = HAx2 + HAy2. Entonces,

Ing. Luis Estrada - 2015 6

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

FAT = ZAsen i + HAcos i

En la realidad puede considerarse como un polo aislado al de un imán extremadamente largo,

de modo que el otro polo no afecte por estar muy distante.

Variación de las componentes horizontal y vertical y del campo total anómalo debido
a un monopolo, para una inclinación magnética de 0° y 70°.

Efecto magnético de un dipolo. En este caso supondremos que el dipolo está magnetizado a
lo largo de su eje (paralelo a su longitud), y sin las restricciones de distancia y longitud
impuestas en las ecuaciones vistas. Si la magnetización es inducida, la orientación del dipolo
coincidirá con la del campo magnético FE. Esto también es un supuesto porque raramente se
presentará en la realidad (magnetismo remanente) pero puede haber una buena aproximación
con suerte. Lo importante de este análisis es que nos permite aprender sobre el
comportamiento de los campos magnéticos de cuerpos de interés geológico.
x=Norte Mag.
x=0 -- x -- P
φ1 φ2
θ θ zn
rn
a
-m rp zp
L b
L
FE 90-θ
+m

Ing. Luis Estrada - 2015 7

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

La intensidad del campo magnético en P debida al polo negativo y al positivo son:

RAn = -m/rn2 = (kFE A)/rn2 y RAp = +m/rp2 = (kFE A)/rp2

Las componentes horizontales y verticales de estas intensidades en P serán:

ZAn = RAn sen φ1 y ZAp = RAp sen φ2,

HAn = RAn cos φ1 y HAp = RAp cos φ2

La componente horizontal y vertical total será la suma de las parciales debidas a cada polo:

ZA = ZAn + ZAp y HA = HAn + HAp

Igual que para el monopolo, usamos la ecuación (#): FAT = ZAsen i + HAcos i, previo reemplazo
de las siguientes relaciones:

θ)
a = L cos(180-θ θ)
b = L sen(180-θ zp = zn + b

rn = (x2 + zn2)1/2 rp = ((x-a)2 + zp2)1/2

sen φ1 = zn/rn cos φ1 = x/rn sen φ2 = zp/rp cos φ2 = (x-a)/rp

Con estas relaciones se puede construir la siguiente gráfica, que nos muestra lo complejo que
se vuelve interpretar las anomalías debidas a un dipolo.

Variación de las componentes horizontal y vertical y del campo total anómalo debido a un dipolo, para una inclinación
magnética de 60° , un campo total de 55.000nT, una susceptibilidad de 0.003 cgs uem y una sección transversal de 1 m .
2

La complejidad citada puede verse claramente en la siguiente gráfica, en la cual la forma de

estas curvas cambia bastante con solo variar la inclinación magnética.

Ing. Luis Estrada - 2015 8

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Anomalías de campo total y vertical para cuatro inclinaciones magnéticas de un dipolo en la dirección del campo terrestre.

Efecto magnético de una esfera. Las ecuaciones son aún más complejas que las vistas. Por
ello iniciaremos el análisis partiendo del caso más simple que es cuando la inclinación es de
90°, es decir con la esfera imantada verticalmente. Utilizaremos la Relación de POISSON que
vincula los potenciales magnético y gravimétrico, también como Reducción al Polo.
Partiremos de un elemento de volumen dv, con masa dm y densidad δ, el que a su vez tiene
un elemento de Momento Magnético dM con imantación I y susceptibilidad uniformes:

dm = δ.dv y dM = I.dv

Utilizaremos la letra p para indicar la intensidad de polo y evitar confusión con la m de masa:

-p p r1 – r2 θ
p.d.cosθ θ
dM.cosθ
dVm = + = p = =
r1 r2 r1.r2 r2 r2

dM (z´-z)
r1-r2≅d.cosθ
θ = (z´-z)/r
cosθ ∴ dVm =
3 -p
r

Sabiendo que d Z

grad F = ∂F/∂
∂x + ∂F/∂
∂y + ∂F/∂
∂z +p
r1

r
y que
mg I r2
P
grad F = ∂F/∂
∂z Z´
X 90-θ
z X´

Ing. Luis Estrada - 2015 9

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Calculamos
∂ (1/r) ∂ (1/r) ∂r
grad(1/r) = =
z ∂z ∂r ∂z

Como r2 = (x´- x)2 + (z´- z)2, diferenciando, miembro a miembro tendremos:

2r.dr = 2(z´- z)dz´ = -2(z´- z)dz, ó dz´ = -dz.

Entonces: dr/dz´ = (z´-z)/r, dr/dz = - (z´-z)/r y ∂(1/r)/∂r = -1/r2,

Reemplazando estos valores en el gradiente,

∂ (1/r) dr (z´- z)
grad(1/r) = =
∂r dz r3

Entonces el potencial magnético es

dVm = dM.grad(1/r)
z
Para integrar toda la masa magnética con Momento dM, reemplazamos este por I.dv,

Vm = I ∫∫∫ grad(1/r) dv
z

Sabemos que el potencial gravimétrico de una masa elemental es dVg = G(dm/r)

Entonces el potencial total será Vg = G. δ ∫∫∫ (1/r) dv

Poisson derivó este potencial con respecto a z, y obtuvo lo siguiente:

∂Vg/ ∂z = G.δ δ ∫∫∫ [(∂

δ ∂ (∫∫∫∫ (1/r) dv)/∂z = G.δ ∂z] dv = G.δ
∂ (1/r))/∂ δ ∫∫∫ grad(1/r) dv
z

Resultando que ∂Vg/ ∂z = G.δ

δ.(Vm/I) ⇒ δ)(∂
Vm = (I/Gδ ∂Vg/∂
∂z) ⇒ δ)gz
Vm = (I/Gδ

Entonces el Potencial Magnético es proporcional a la derivada del Potencial Gravimétrico

en la dirección de magnetización, siempre que la densidad y la susceptibilidad sean uniformes.
Como la dirección de magnetización en nuestra derivada es la vertical, las anomalías de campo
vertical ZA y HA pueden definirse como:

∂Vm I ∂2Vg ∂Vm I ∂(∂

∂Vg)
ZA = = y HA = =
∂z δ
Gδ ∂z2 ∂x δ
Gδ ∂x ∂z

Para una esfera:

Gm πR3δ
G(4/3)π πR3δ
G(4/3)π
Vg = = =
r r (x2 + z2)½

Ing. Luis Estrada - 2015 10

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

∂Vg πR3δ
G(4/3)π ∂2Vg πR3δ
G(4/3)π
= 2z y = (2z2-x2)
∂z (x2 + z2)3/2 ∂z2 (x2 + z2)5/2

πR3I
(4/3)π πR3I
(4/3)π
ZA = (2z2-x2) y HA = 4zx
2 2 5/2
(x + z ) (x2 + z2)5/2

El caso general de una esfera uniformemente magnetizada, pero donde el campo de la Tierra
está inclinado, tiene una derivación similar pero más compleja, y resulta ser:

πR3kFE sen i
(4/3)π 3z2 3xz cotg i
ZA = - -1
(x2 + z2)3/2 (x2 + z2) (x2 + z2)

πR3kFE cos i
(4/3)π 3x2 3xz tg i
HA = -1- -1
(x2 + z2)3/2 (x2 + z2) (x2 + z2)

Igual que antes, la anomalía total será: FAT = ZA sen i + HA cos i

Si graficamos dando valores a estas ecuaciones, veremos que para altas latitudes las curvas
son similares, mientras que a bajas latitudes son bastante diferentes.

Anomalías de campo total, horizontal y vertical para 60º de inclinación magnética de una esfera uniformemente magnetizada.

Ing. Luis Estrada - 2015 11

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Anomalías de campo total y vertical para diferentes inclinaciones magnéticas de una esfera uniformemente magnetizada.

En todo el análisis precedente se simplificó bastante la situación. Se consideró que la anomalía

estaba orientada paralela al campo magnético terrestre y que este era vertical. Esto dista
mucho de la realidad, pero sirve para familiarizarse con la interpretación magnética y para
tomar conciencia que la prospección magnética tiene un alto grado de subjetividad y por lo
tanto su interpretación es fundamentalmente cualitativa.

Efecto magnético de una hoja horizontal delgada. Utilizaremos el mismo procedimiento de

Bouguer en la Prospección Gravimétrica. Al igual que para la Relación de Poisson, suponemos
que el campo terrestre es vertical. Si consideramos una cinta de polos negativos, entonces la
intensidad RA del campo magnético en un punto P debida a un área elemental dx.dy, sabiendo
que H = m/r2 e I = m/área y siguiendo la derivación de Bouguer, el campo vertical ZA debido
a la cinta expandida al infinito será:

Ing. Luis Estrada - 2015 12

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

x=+∞

(m/área) dx.dy (m/área) dx.dy

RA = y ZA = θ
cosθ
r2 r2
x=-∞

Expansión de una cinta a una hoja y a una placa infinita para calcular su efecto magnético.

π/2
cos θ = d/r ⇒ r = d/ (cos θ)
θ θ.dy
(m/área).rdθ
cos θ = r.dθ/dx ⇒ ZA = θ
cosθ
2
r -π/2 r
cos θ.dx = r.dθ d
π/2

(m/área).dy
rdθ ZA = θ.dθ
cosθ θ
θ d
-π/2
dx
2(m/área)
ZA = dy
d

El próximo paso es expandir la cinta a una hoja, recordando que se rota el punto vista en el
segundo paso, cambiando la distancia d de la figura (a) por la r de la figura (b).
y=+∞

2(m/área) dy 2(m/área) dy
ZA Hoja = = θ
cosθ
d y=-∞ r

Ing. Luis Estrada - 2015 13

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Ahora dy.cos θ = r.dθ

θ, entonces
π/2

ZA Hoja = θ
2(m/área) dθ y finalmente ZA Hoja = 2(m/área) π = 2.π
π .I
-π/2

Si la hoja tiene un espesor, la placa ya no tendrá solo polos negativos arriba, sino también los
positivos en el fondo. En ese caso la anomalía resultante será nula, es decir:

π .I - 2.π
ZA = ZA Tope - ZA Fondo = 2.π π .I = 0
Pero si truncamos la placa puede demostrarse que π, el ángulo subtendido al infinito, puede
reemplazarse por los que subtiende la placa truncada al tope y al fondo.
P
ZA = ZA Tope - ZA Fondo

θ1
θ1 - 2.I.θ
ZA = 2.I.θ θ2 = 2.I. (θ
θ 1 - θ 2)
FE
θ2
Esta relación se aplica para determinar la anomalía
---------------------- magnética de un basamento tapado por material no
magnético, típico de una gruesa secuencia de sedi-
+ + + + + + + + + + + + + + + + mentos Paleozoicos sobre Precámbricas rocas ígneas y
metamórficas.
Otro ejemplo es el de
una delgada cobertura
de sedimentos del Pleis-
toceno sobre rocas sedi-
mentarías Mesozoicas
cortadas por rocas
basálticas, la geología
del valle Conectitud de
Massachusetts.
Veamos el ejemplo
correspondiente a la
geología mencionada.
Suponemos el campo
terrestre vertical por
simplicidad. Si el
basamento fuera de
gran espesor (a), el
efecto sería como el de
una hoja con todos los
polos negativos sobre la
superficie del basa-
mento. La susceptibi-
lidad del sedimento del
Pleistoceno es tan baja
que puede despreciarse.
k’=0,0045 es la suscepti-
bilidad efectiva del
basalto.
Con la última ecuación,
y el ángulo en radianes,
podemos calcular la
anomalía vertical como:

Ing. Luis Estrada - 2015 14

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

ZA = ZA Tope - ZA Fondo ZA = 2(0.0045x60.000x2.38) - 0 = 1.286nT

De un modo similar puede calcularse un horst de solo 4 metros sobre el basamento. Este salto
es significativo en comparación al espesor del sedimento, por ello se hace necesario considerar
los polos positivos bajo dicho salto. La anomalía en P debida a este horst será:

ZA = ZA Tope - ZA Fondo = 2(0.005x60.000x2.52) - 2(0.005x60.000x2.38) = 84nT

Como un ejercicio final en obtener ecuaciones para calcular efectos magnéticos de figuras
geométricas sencillas, intentamos ahora con una losa semi-infinita utilizando dos ángulos en
razón de los polos negativos y positivos.
P θ1 = π/2 + tg-1(x/z)
x=0
θ1 θ2 = π/2 + tg-1(x/(z+t))
FE
z θ2 θ1 - θ2)
Reemplazando en la ecuación ZA = 2 I (θ
-----------------------
t tendremos la atracción de una losa semi-infinita:
+++++++++++++++++ ∞

π/2 + tg-1(x/z) - π/2 - tg-1(x/(z+t))] = 2 k FE [tg-1(x/z) - tg-1(x/(z+t))]

ZA = 2 k FE [π

Efecto de cuerpos de sección irregular. Este efecto puede calcularse igual que en
gravimetría. Para ello existe un desarrollo matemático similar al que utilizó Talwani, y que tiene
su aplicación mediante el uso de programas de computación específicos.
Además pueden adaptarse dos técnicas de interpretación directa utilizadas en gravimetría. La
primera es la técnica del medio-máximo que da una idea bastante aproximada de la
profundidad de la fuente. Partimos de la ecuación ya vista para un monopolo, y la
particularizamos para una varilla vertical delgada con el polo profundo muy lejos. Suponemos
además que nuestro perfil pasa sobre la varilla:
z z (kFE A)
ZA =
(x2+z2)3/2

El máximo se obtiene cuando x = 0, y la mitad,

dividiéndolo en dos obviamente, entonces el medio-
zmax máximo será:
(kFE A)
z(1/2)max ZA(1/2) max =
2z2

En el punto donde la mitad del máximo toma este

x valor tendremos:

r z θ
(kFE A)(z2 + x(1/2)max2)cosθ
θ ZA =
(z2 + x(1/2)max2)3/2

θ = z/r y r = (z2+x(1/2)max2)1/2, tendremos:

Como cosθ

θ
(z2 + x(1/2)max2)1/2cosθ 1
=
(z2 + x(1/2)max2)3/2 2z2

Ing. Luis Estrada - 2015 15

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Ecuación que se reduce a 2z3 = (z2 + x(1/2)max2)3/2 y finalmente queda z = 1.3 x(1/2)max

Del mismo modo puede determinarse para la esfera que resulta z = 2.0 x(1/2)max

Para la losa semi-infinita la profundidad z es igual a la mitad de la distancia entre Zmax y Zmin.

z La segunda técnica es el método de las

Máxima →
pendiente pendientes o de máximo gradiente, que
también da una idea bastante aproximada de la
½ máxima profundidad de la fuente. Consiste en crear una
pendiente línea con una pendiente igual a la mitad de la
máxima pendiente. Esta máxima pendiente está
zmax obviamente en el punto de inflexión de la curva.
La línea creada se mueve por la curva hasta que
sea tangente a la curva. Habrá dos puntos de
tangencia y la distancia horizontal d entre ellos
es proporcional a la profundidad z, es decir:
d x
z = d / 1.6

Determinación de la magnetización de una roca

Recordemos que el efecto magnético de una masa rocosa es proporcional a su magnetización

I (A/m), la cual es un vector cuyas componentes en la dirección de los tres ejes son I x, I y e I z.
La fuerza magnetizante H generará un campo magnético B, tal que B = µ H con µ la
permeabilidad magnética de la sustancia. Las unidades con que se mide cada una son las
siguientes:
H en Ampere/metro (A/m) (1 Oersted = 79,58 A/m)
B en Weber/metro2 (wb/m2) o Tesla (T)
µ en Weber/Ampere.metro (wb/Am)

En el caso del campo magnético de la Tierra Bo, µo = 4ππ x10-7 (wb/Am). Un cuerpo magnético
de volumen unitario en el campo magnético de la Tierra experimenta una fuerza magnetizante I
por unidad de volumen dentro del cuerpo, por lo que H e I son dimensionalmente equivalentes.
La fuerza magnetizante I es proporcional a H ya que I = kH con k la susceptibilidad
adimensional. El campo magnético fuera del cuerpo será por lo tanto:

B = Bo + µokH

B = µoH + µo kH = (1 + k)Bo = µo (H + I)

El Momento magnético o magnetización de la masa rocosa de un determinado volumen V será:

M = I.V (A/m2)

Nuevamente aparece una equivalencia, puesto que todos los cuerpos con el mismo momento
magnético dan el mismo efecto magnético.
La magnetización I = kH también es llamada magnetización inducida Iind, puesto que solo
existe en presencia de una fuerza magnetizante H. Sin embargo, la mayoría de las rocas
magnéticas tienen una magnetización remanente Irem aunque H = 0. La causa de esta
remanencia radica en la estructura de los dominios magnéticos de las sustancias Ferromag-
néticas (generalmente magnetita). La magnetización total de una roca será entonces:

I = Iind + Irem

Ing. Luis Estrada - 2015 16

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

La dirección de Irem en las rocas volcánicas es la del campo magnético existente al momento
que la roca se enfrió. En el caso de rocas volcánicas jóvenes, las direcciones de los campos
remanente e inducido coincidirán.
Como el campo de la Tierra es también vectorial, puede ser descompuesto en tres
componentes Box, Boy y Boz coincidentes con el Norte, el Este y Vertical respectivamente. Por
simple álgebra vectorial estas componentes pueden ser también expresadas por el módulo de
Bo, la declinación d y la inclinación i magnéticas. Las magnetizaciones total, inducida y
remanentes también pueden ser descompuestas en componentes paralelas a los ejes x,y,z.
Ejemplo de cómo evaluar I:
Calcular la magnitud de I rem y sus componentes para una roca basáltica de Auckland (New
Zealand) desde el hecho que la declinación de I es alrededor de 0º y que I ≅ 5A/m (de
mediciones aéreas). La susceptibilidad media de este basalto es de 4π πx2x10-3. Siempre es
razonable suponer que la inclinación de I es igual a la de I ind.
El Campo magnético en New Zealand es F = 55.250 nT, la inclinación i = - 62,5º y la
declinación d = 18º.
µo) entonces
Como I ind = kH = k(F/µ

π x 2 x 10-3 x 0,55250 x 10-4 (wb/m2)

4π
 I ind = = 1.105 A/m
π x 10-7 (wb/Am)
4π
Entonces las componentes de I ind serán:

I Xind  I ind cos(i) cos(d) = 1.105 x 0,462 x 0,951 0,486 A/m

I ind = I Yind =  I ind cos(i) sen(d) = 1.105 x 0,462 = 0,158 A/m
I Zind  I ind sen(i) = 1.105 x (-0,887) -0,980 A/m

Sabemos que

IX 5 cos(62,5º) cos(0º) = 5 x 0,462 = 2,309 A/m

I = IY = 5 cos(62,5º) sen(0º) = 5 x 0,0 = 0,0 A/m
IZ 5 sen(62,5º) = 5 x (-0,887) = -4,435 A/m
Entonces
NG
I X - I Xind = I Xrem = 2,309 – 0,486
I Y - I Yind = I Yrem = 0,158 - 0
I Z - I Zind = I Zrem = -4,435 + 0,980
d=18º NM
i=62,5º

I ind 1,823
I Xrem = 0,158 = 3,91 A/m
-3,455

I rem 3,91 = √ 1,8232+0,1582+3,4552

Declinación de I rem = arctg(0,158/1,823) ≅ 5º

Ing. Luis Estrada - 2015 17

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Ejemplo de intensidad de magnética debida a un dipolo inducido a 60º

Ing. Luis Estrada - 2015 18

Geofísica – FACET – UNT – Prospección Magnética para Ingenieros

Bibliografía

An Introduction to Applied and Environmental Geophysics - John M. Reynolds – Wiley - 1997

Fundamentos de Geofísica - Agustín Udias – Julio Mezcua -Alianza Universidad Textos -1997
El Magnetismo de las Rocas – Daniel Valencio – Editorial Universitaria de Buenos Aires - 1980
Geomagnetism – G.N. Petrova, V.N. Brovov and M.I. Pudovkin – U.G.G.I. – UNESCO - 1980
Tratado de Geofísica Aplicada - José Cantos Figuerola – Litoprint - 1978
Introduction to Geophysical Prospecting - Milton Dobrin - McGraw – Hill B. Company –1976
Applied Geophysics - W. M. Telford – L. P. Geldart, R. E. Sheriff, D. A. Keys - 1976
Geofísica Minera - D. S. Parasnis – Paraninfo - 1971
Introducción a la Geofísica – Benjamín F. Howell, Jr. – Ediciones Omega - 1962
Exploration Geophysics - J. J. Yakosky - Trija Publishing Company - 1957

Ing. Luis Estrada - 2015 19