Source: http://docplayer.es/2713636-Tema-5-metodo-electromagnetico-en-el-dominio-de-las-frecuencias-fdem.html
Timestamp: 2018-06-20 19:51:38+00:00

Document:
Tema 5: MÉTODO ELECTROMAGNÉTICO EN EL DOMINIO DE LAS FRECUENCIAS (FDEM) - PDF
Download "Tema 5: MÉTODO ELECTROMAGNÉTICO EN EL DOMINIO DE LAS FRECUENCIAS (FDEM)"
Francisco José Hidalgo Cárdenas
1 Tema 5: MÉTODO ELECTROMAGNÉTICO EN EL DOMINIO DE LAS FRECUENCIAS (FDEM) 5.1 Introducción. Los métodos electromagnéticos constituyen el grupo de mayor diversidad tanto en metodologías como en instrumental, resultando ser herramientas eficaces a la hora de abordar una gran variedad de problemas: estudio de aguas subterráneas, cartografía de la intrusión marina, detección de vertidos de contaminantes, localización de cavidades, etc.. Versatilidad y gran celeridad en la toma de datos (aspecto que posibilita el reconocer extensas áreas de terreno en un corto espacio de tiempo), son las dos principales virtudes de los métodos electromagnéticos en general. Además la no necesidad de contacto físico con el terreno (no precisan electrodos), permiten incluso el realizar estudios tanto en barco como en avión. Los métodos electromagnéticos se pueden clasificar básicamente en base a dos aspectos: La naturaleza del campo electromagnético primario, distinguiendo entre los de campo natural o artificial (el hombre es quien lo genera). La distancia a la que se encuentre el receptor con respecto al origen del campo primario. En este sentido tenemos los de campo próximo o campo lejano. Si bien existen diversos métodos electromagnéticos aplicables al problema de detección de filtraciones de agua en el subsuelo, son los de campo próximo y artificial los realmente importantes, dado que en éstos es el operador quien controla el transmisor (dispositivo encargado de generar el campo electromagnético primario), y por tanto el hecho de poder definir las características de la onda primaria, nos permitirá controlar un parámetro tan importante como es la profundidad de investigación. Dentro del conjunto de técnicas electromagnéticas de campo próximo y artificial, en esta tema vamos a describir el FDEM, método electromagnético en el dominio de frecuencias en donde la profundidad de investigación es función de la frecuencia. 50
2 5. Objetivo del método. El FDEM, al igual que el resto de los métodos electromagnéticos, tiene como objetivo medir la conductividad eléctrica del terreno. Dado que la conductividad se define como la inversa de la resistividad, la detección de posibles focos de filtración de agua se reducirá simplemente a la localización de zonas, en donde se observen incrementos anómalos en su valor. 5.3 Principios teóricos básicos Teoría de la propagación de los campos electromagnéticos. Con objeto de facilitar la exposición de la teoría, trabajaremos bajo la hipótesis de medio homogéneo, lineal e isótropo. La descripción de cualquier campo electromagnético la podemos establecer a partir de las 4 ecuaciones de Maxwell: rote roth B t D J + t = (Ley de Faraday) [ 5.1] divd ρ c = (Ley de Ampere) [ 5.] = (Ley de Coulomb) [ 5.3] divb = 0 [ 5.4] En donde E (V/m) es el vector intensidad del campo eléctrico, D (C/m ) es el vector intensidad del flujo eléctrico, J (A/m ) el vector densidad de corriente, H (A/m) vector intensidad del campo magnético, ρ c (C/m 3 ) la distribución de la carga eléctrica, y B (Wb/m o Tesla) vector densidad del flujo magnético. Además de estas expresiones, para la descripción del comportamiento de los campos electromagnéticos también se precisa de las siguientes ecuaciones constitutivas. = [ 5.5] = [ 5.6] = [ 5.7] J σe D εe B µh 51
3 En donde ε (F/m) es la constante dieléctrica, µ (H/m) la permeabilidad magnética del medio, y σ (S/m) es la conductividad eléctrica (recordemos que ρ = 1/σ con unidades de Ω.m ). Sea i = 1 e ω = πf (rad/s) la frecuencia angular del campo, si ahora relacionamos todas las expresiones anteriormente descritas, y además asumimos que E y H presentan una dependencia con el tiempo de la forma H E iωt ( t) = H 0e [ 5.8] iωt ( t) = E0e [ 5.9] se obtiene que las ecuaciones básicas que describen la propagación eléctrica y magnética en un medio isótropo, lineal y homogéneo son: H = iωµσh εµω H [ 5.10] E = iωµσh εµω E [ 5.11] 5.3. Atenuación del campo electromagnético y profundidad de penetración. En el caso de estar en un medio isotrópo aislante perfecto, las ondas electromagnéticas introducidas en el terreno podrían alcanzar una profundidad indefinida. Sin embargo este estado ideal no existe, por lo que presencia de cuerpos conductores en el subsuelo conllevan la atenuación de la amplitud de las ondas electromagnéticas. Consideremos el caso de una onda plana que se propaga por un semiespacio de conductividad σ en la dirección z hacia abajo, en donde el campo eléctrico solo presenta componente x y el campo magnético solo componente y, y ambos además con una amplitud sinusoidal con respecto z. En estas condiciones si procedemos a resolver las ecuaciones [.10] 5, la solución del campo eléctrico (análogo para el campo magnético) vendrá dada por la siguiente expresión: 5 y [.11] E x ikz = E 0 e [ 5.1] 5
4 Al término K = ( β + iα) se le denomina número complejo de onda o parámetro de propagación, en donde β es el factor de fase y α el factor de atenuación. µε σ = ω 1+ 1 ωε β [ 5.13] µε σ = ω ωε α [ 5.14] En el caso de materiales conductivos el factor de atenuación α se puede simplificar de la siguiente forma: ωµσ α = [ 5.15] Partiendo de estas hipótesis, se define la profundidad nominal de penetración δ, también denominada profundidad Skin o profundidad efectiva, como aquella distancia para la cual la amplitud del campo se reduce a 1/e de su valor inicial. Esta profundidad δ es inversamente proporcional al factor de atenuación α. 1 δ = (/ µσω) α 1/ = = ( ) 1/ / σ f [ 5.16] Esquema teórico del funcionamiento del FDEM. Los FEDM son sistemas que trabajan con un transmisor y un receptor móviles de ondas electromagnéticas. A través del transmisor, que es una bobina por la que circula una corriente alterna sinusoidal de frecuencia f (Hertz), se genera un campo electromagnético primario H p, que se propaga libremente tanto por la superficie como por el interior del terreno. Este campo está formado por una componente magnética y otra eléctrica ortogonales entre si, y a su vez perpendicular a la dirección de propagación de la onda. En presencia de un cuerpo conductor, la componente magnética del campo primario provoca la aparición de unas corrientes inducidas, que a su vez generan un campo magnético secundario H. s 53
5 Este campo secundario, junto con el campo primario que circula libremente por la superficie del terreno, serán detectados por la bobina receptora dando lugar a un campo resultante, que diferirá tanto en fase como en amplitud con respecto al campo primario, aspecto que nos proporcionará información acerca del tamaño y conductividad eléctrica del elemento conductor. Fig. 1 Esquema general del campo electromagnético (Reynolds [ 19 ]) Hipótesis de trabajo para la aplicación del FDEM. Los métodos electromagnéticos en el dominio de las frecuencias (también conocidos como conductímetros), en general pueden considerarse sistemas Slingram de bobina emisora y receptora móviles si bien con una gran diferencia. Fijada una distancia entre bobinas, la frecuencia de emisión es muy pequeña de forma que la profundidad efectiva o Skin es siempre significativamente mayor que dicha distancia. Esta condición se conoce como operación a bajo número de inducción. Sea N el número de inducción y definido como el cociente entre la distancia entre la bobina transmisora y receptora S (m), y la profundidad Skin δ (m), diremos que operamos a bajo número de inducción cuando se cumpla que N << 1, en donde σ a es la conductividad aparente del medio medida en ms/ m. N = S S ( µ 0ωσ a / ) δ = [ 5.17] 54
6 Sea H y H la amplitud del campo magnético primario y secundario p s respectivamente, el operar a bajo número de inducción conlleva como principal consecuencia, el poder garantizar que existe una relación de proporcionalidad entre el cociente H / H y la conductividad eléctrica σ a según la siguiente expresión s p H H s p = iµ σ / 4 [ 5.18] 0 as Por tanto y según la expresión [ 5.19], a través de los aparatos de medida ya obtendremos directamente el valor de la conductividad eléctrica aparente (en ms/m). McNeill [ 33 ]. σ 4 H s = [ a 5.19 ] µ S 0ω H p Destacar que fijada la frecuencia de la onda primaria así como el espaciado entre bobinas, en terrenos de elevada conductividad la expresión [ 5.17] puede no ser inferior a la unidad, y por tanto no se cumple la relación de proporcionalidad entre H s / H p y σ a. En estos casos será necesario corregir las medidas de conductividad obtenidas Modos de operar y profundidad de investigación. A la hora de medir la conductividad del terreno, se distinguen dos formas de operar en función de la orientación de las bobinas con respecto a la superficie del terreno: Modo dipolo vertical Modo dipolo horizontal El modo dipolo vertical (fig.) consiste en disponer horizontalmente con respecto a la superficie del terreno, tanto la bobina transmisora como la receptora (el momento del dipolo magnético será vertical), mientras que para el dipolo horizontal (fig.3) dispondremos ambas bobinas de forma vertical, con lo que el momento del dipolo magnético será horizontal 55
7 S S Fig. Modo dipolo vertical Fig. 3 Modo dipolo horizontal La diferencia fundamental entre ambos modos de operar, reside en la sensibilidad a la hora de captar la respuesta de un cuerpo conductor frente la acción de un campo electromagnético. En este sentido, se comprueba que la profundidad efectiva de investigación del dipolo vertical es superior al del dipolo horizontal. Como valor aproximado el dipolo vertical presenta una profundidad efectiva de investigación de 1.5 veces la separación entre bobinas, mientras que para el dipolo horizontal es de 0.75 (McNeill [ 34 ]). Fig. 4 Ejemplo de las diferentes sensibilidades que presentan el dipolo vertical y el horizontal. (McNeill [ 34 ]). También se ha observado que en el caso de usar bobinas grandes (i.e. Geonics EM-34), el dipolo horizontal es más susceptible a presentar problemas por mal alineamiento de las bobinas que el dipolo vertical, por lo que se precisará de un mayor cuidado a la hora de proceder a realizar las lecturas de campo (McNeill [ 34 ]). 56
8 5.4 Equipo requerido. La diferencia más significativa del FDEM con respecto a las técnicas que hemos visto a lo largo de la tesina, es la no utilización de electrodos. Dentro de la diversidad instrumental del FDEM, el equipamiento básico consta de: Una bobina emisora, que es la encargada de emitir el campo electromagnético primario. Una bobina receptora, con la que captaremos el campo electromagnético secundario. Una o varias consolas (según modelo), en donde podremos almacenar digitalmente los datos, así como fijar los parámetros de estudio. En función del instrumento de medida, la separación entre bobinas (parámetro que condiciona la profundidad de investigación) será una variable de estudio que podrá ser modificada durante la campaña de campo, o por lo contrario será un parámetro invariable. Por lo tanto, en función del modelo escogido podremos trabajar con más o menos niveles de profundidad de investigación. A modo de ejemplo tenemos los siguientes modelos: 1. GEONICS EM-31: consta de una barra en donde en cada extremo tenemos una pequeña bobina (ver fig. 5), de forma que la distancia entre las bobinas es siempre fija (3.66 m.). Este instrumento solo precisa de un operario y fue el aparato utilizado en el estudio de Cardona.. GEONICS EM-34: consta de dos bobinas gigantes precisándose dos operarios, uno para cada bobina (ver fig. 4). Con este modelo tenemos la opción de poder trabajar con distintas separaciones entre bobinas (i.e. 10, 0 y 30m ). Fig. 5 Modelo Geonics EM-31 Fig. 6 Detalle de la consola del EM-31 57
9 5.5. Metodología de campo Planificación de la campaña de reconocimiento. Los aspectos a tener en cuenta en la planificación son análogos a los ya descritos para los métodos anteriores: Delimitar objetivos: cartografiar posibles trayectorias, monitorizar el comportamiento del flujo, etc.. Número y ubicación de los perfiles, así como el espaciado entre medidas. Profundidad de investigación requerida. Recopilar toda la información acerca del contexto geológico de la zona. Escoger el equipo de medida adecuado. Adoptar las medidas oportunas a fin de garantizar la calidad de los datos. Quizás la cuestión más importante a la hora de planificar el estudio sea escoger el equipo de medida adecuado, ya que éste condiciona tanto la profundidad de investigación, así como el número de niveles de profundidad con los que podremos trabajar. En este sentido el Geonics EM-31 presenta una distancia fija entre bobinas de 3.66 m., aspecto que conlleva el poder realizar únicamente dos lecturas sobre la vertical de cada punto de muestro (una para dipolo vertical y otra para dipolo horizontal). En cambio con el Geonics EM-34 podemos trabajar con distintas separaciones. En general es frecuente utilizar longitudes de cable de 10, 0 y 30 m. obteniéndose por tanto hasta 6 lecturas por punto de medida. Dado que a menor separación entre bobinas la profundidad de investigación disminuye, queda claro que el Geonics EM-31 es más indicado para estudios a poca profundidad, mientras que en el caso de querer alcanzar grandes niveles de profundidad, deberemos recurrir a un modelo que nos permita utilizar grandes separaciones entre bobinas (i.e. Geonics EM-34). A modo ilustrativo en la tabla 1 se muestran algunos valores de profundidad. 58
10 separación entre bobinas (m) Frecuencia (KHz) Profundidad de investigación (m) Dipolo horizontal Dipolo vertical Geonics EM Geonics EM Tabla Relación entre los diferentes factores que afectan a la profundidad de investigación. Los dos aspectos principales que caracterizan al FDEM son: 1. Una limitada resolución vertical. Gran celeridad en la toma de datos. En cuanto a la primera cuestión y a raíz de lo expuesto anteriormente, si bien el número de lecturas sobre la vertical de un punto de muestreo varía según el instrumento de medida, éste continua siendo limitado. En lo que concierne a la segunda cuestión, la gran celeridad del proceso de medición del FDEM, posibilita el poder medir en un tiempo razonable la conductividad del terreno en un gran número de puntos. Esto lo convierte en un método de gran resolución horizontal. Por tanto el responsable del proyecto podrá optar no solo por realizar un par de perfiles (de forma análoga a lo visto en Tomografía eléctrica y Potencial espontáneo), sino también y en función de los objetivos del estudio, por discretizar todo el área de estudio a modo de malla de muestreo, a fin de poder confeccionar mapas de isolineas de conductividad a distintas cotas de profundidad. Son justamente estos dos aspectos, una limitada resolución vertical y una gran resolución horizontal, lo que conlleva a que el FDEM sea más idóneo en estudios preliminares a gran escala, en donde a través del análisis de mapas de isovalores de conducitividad, detectar áreas susceptibles de presentar flujos subterráneos de agua, para posteriormente y mediante métodos con mayor resolución vertical (i.e. Tomografía eléctrica), caracterizar con más detalle el problema. 59
11 5.5. Ejecución del método. En el caso de utilizar el Geonics EM-34 (bobinas grandes), dado un punto de medida del terreno, colocaremos el dispositivo en un modo (i.e. modo dipolo vertical), fijaremos la frecuencia de emisión del campo primario así como su correspondiente longitud de cable, y procederemos a realizar la primera medida. Ahora sin cambiar de espaciado entre bobinas ni de frecuencia, cambiaremos al otro modo (dipolo horizontal) y realizaremos la otra lectura. Posteriormente se procederá a cambiar de longitud de cable y frecuencia, y se repetirá el proceso descrito anteriormente. Es frecuente trabajar con 3 longitudes de cable, de forma que el número de medidas de conductividad realizadas es de 6. Destacar que a medida que aumentemos la separación entre bobinas, la frecuencia de emisión disminuye a fin de asegurar estar operando a bajo número de inducción (ver tabla ). En el caso de utilizar el modelo Geonics EM-31, bastará simplemente con rotar 90º la barra para pasar de dipolo vertical a dipolo horizontal. Destacar que además de proceder a tomar las lecturas paso a paso, podemos optar por medir de forma continuada Procesado e interpretación de los datos. Los instrumentos de medida ya están preparados para proporcionarnos directamente los valores de conductividad aparente (en ms/m) del terreno. Sin embargo, un aspecto a tener en cuenta será el controlar que la hipótesis de operar a bajo número de inducción ( N << 1) se cumple, ya que recordemos que fijada la frecuencia del campo electromagnético primario, así como la distancia entre bobinas, la conductividad del terreno puede ser lo suficientemente elevada como para que la expresión [ 5.17] sea superior a la unidad, y en consecuencia la expresión [ 5.19] deje de tener validez. En la práctica este proceso de verificación y corrección se realiza automáticamente mediante el software correspondiente. En lo que concierne a la fase de interpretación de los datos, la localización de zonas susceptibles de presentar filtraciones de agua se reduce a la detección de áreas en donde se produzcan incrementos anómalos de conductividad. 60
12 Sin embargo ahora además del agua, deberemos tener presente que la presencia en el subsuelo de elementos antrópicos conductivos (i.e elementos metálicos, tendido eléctrico, etc.), también inducirán incrementos significativos de la conductividad eléctrica del terreno, por lo que es aconsejable el conocer previamente posibles localizaciones de estos elementos. Por último destacar que en el FDEM, estimar la profundidad a la que se encuentra el objeto responsable de la anomalía de conductividad, no será tan precisa como ocurre con Tomografía eléctrica. Ello se debe a que dada una profundidad efectiva de investigación (i.e. 6m.), el registro obtenido recoge la influencia de los diversos elementos conductivos situados por encima de este nivel, en donde lógicamente cuanto más próximos se hallen de la profundidad de 6 m., la amplitud de la anomalía medida será mayor. 5.6 Ventajas y limitaciones del método. Las ventajas que presenta el FDEM son: Es un método que no requiere de contacto con el terreno. Esto le convierte en una técnica eficaz en terrenos muy resistivios, en los que ni el uso del método del potencial espontáneo ni el de resistividad eléctrica son factibles. Es un método de gran resolución horizontal, aspecto que le convierte en una técnica ideal para estudios a gran escala. Es un método no destructivo. En cuanto a las limitaciones del método tenemos: Baja resolución vertical. el FDEM es muy sensible a los ruidos generados por la presencia de campos electromagnéticos externos. 61
TEMA 6: Detección de filtraciones de agua en Cardona mediante los métodos geofísicos de Tomografía eléctrica y FDEM. 62 6.1 Antecedentes. La zona anexa a la factoría minera de Cardona, presenta desde hace

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución