Tomografia Electrica - Otro Autor
Tomografia Electrica - Otro Autor
Tomografia Electrica - Otro Autor
Si bien a nivel terico no existen diferencias significativas entre ambas variantes, en esta tesina solo nos centraremos en la descripcin del Electrical imaging, ya que adems de ser la tcnica ms adecuada para el estudio de problemas de filtracin (en donde se precisan perfiles de gran longitud), es tambin la variante ms econmica, fcil y rpida de implementar en el campo.
14
La aparicin de la Tomografa elctrica ha supuesto un salto cualitativo espectacular con respecto a los mtodos de resistividad convencionales, tcnicas que aunque se han estado utilizando durante varias dcadas en estudios de filtracin de agua, su limitada resolucin 2-D les confera en general un papel secundario frente a otras tcnicas (i.e. Potencial espontneo). Este avance se debe fundamentalmente a dos razones: En primer lugar porque estos mtodos convencionales solo utilizan 4 electrodos, siendo necesario para cada medida variar manualmente sus posiciones en el terreno (proceso relativamente lento y pesado). En cambio el mtodo de Tomografa elctrica es una tcnica multielectrdica, en el que todo el proceso de adquisicin de datos est totalmente automatizado. Esto nos permitir poder realizar un gran nmero de medidas, tanto en profundidad como lateralmente, en un breve espacio de tiempo (del orden de 500 medidas en una hora y media), obteniendo por tanto modelos 2-D de gran resolucin. En segundo lugar porque el empleo de tcnicas numricas nos permitir poder procesar eficazmente todo este gran volumen de informacin.
15
3.3 Principios tericos bsicos. 3.3.1 Esquema bsico del funcionamiento del mtodo.
Independientemente de que la Tomografa elctrica sea una tcnica multielectrdica, la base terica de su funcionamiento es anloga al de los mtodos de resistividad convencionales.
Los mtodos convencionales solo precisan 4 electrodos y se basan en introducir en el terreno, un campo elctrico de corriente continua mediante dos electrodos de corriente (A,B) conectados a un miliampermetro, mientras que con los otros dos electrodos (M,N) y que estn conectados a un milivoltmetro, mediremos cual es la diferencia de potencial elctrica V entre esos dos puntos, parmetro a partir del cual podremos calcular el valor de la resistividad (ver 3.3.2) en el punto medio del dispositivo y a una profundidad determinada (ver 3.5.3). Existen diferentes configuraciones a la hora de colocar los 4 electrodos, siendo las ms utilizadas Wenner (la variante ) y Schlumberger. Tambin tenemos la denominada configuracin Wenner-Schlumberger (muy utilizada en Tomografa elctrica).
16
WENNER
SCHLUMBERGER
WENNER-SCHLUMBERGER
n =1 Wenner n 2 Schlumberger
Escogido el dispositivo y realizada la primera medida, si lo que se persigue es medir como vara la resistividad del terreno sobre la vertical del punto medio o, iremos expandiendo los electrodos que en el caso de usar el dispositivo Wenner, debern ser los 4 (para mantener la condicin de AM=MN=NB), mientras que para Schlumberger (o Wenner-Schlumberger) solo expandiremos los electrodos A y B. Este anlisis de la variacin vertical del valor de la resistividad es lo que se conoce con el nombre de sondeo elctrico vertical SEV. En cambio si lo que pretendemos es medir la variacin lateral de la resistividad del terreno, bastar con mover los 4 electrodos solidariamente (tanto para Wenner como Schlumberger) hacia uno u otro lado. Esto es lo que se conoce como calicata elctrica. Como se puede apreciar la obtencin de las medidas en los mtodos convencionales es un proceso significativamente lento y pesado.
17
Partiendo de estos conceptos bsicos, la gran innovacin del mtodo de Tomografa elctrica con respecto a los mtodos convencionales, reside en que ahora todas las medidas se realizarn de forma totalmente automatizada, es decir sin necesidad de mover manualmente ningn electrodo. Ello se debe a que por un lado trabajaremos con un gran nmero de electrodos en el terreno (dispuestos equiespaciadamente), y por otro lado a que nuestro dispositivo de medida de resistividades, se encargar de realizar automticamente toda la secuencia de medidas preestablecida, formando para ello y segn las especificaciones predefinidas, todas las posibles combinaciones de 4 electrodos. De esta forma se ir obteniendo la variacin de resistividad del subsuelo tanto en profundidad como lateralmente. En la figura 5 tenemos un ejemplo de la secuencia de medicin utilizando el dispositivo Wenner.
18
En un semiespacio homogneo e istropo de resistividad , y sometido a una intensidad I , el valor del potencial V de un punto situado a una distancia r con respecto a un electrodo puntual es
V= I 2r
[3.1]
Partiendo de la expresin [3.1] , para una disposicin electrdica como la mostrada en la figura 3, el potencial medido por los electrodos de medida M y N ser.
VM = I 1 1 AM BM ; 2 VN = I 1 1 AN BN 2
[3.2]
[3.3]
Si ahora despejamos obtendremos el valor de la resistividad para el punto medio del dispositivo. La resistividad no es ms que una medida de la dificultad que la corriente elctrica encuentra a su paso por un material determinado.
a =
V K I
[3.4]
En donde a se mide en Ohmios-metro, I en Amperios, y V en Voltios . El trmino K es el factor de penetracin o coeficiente del dispositivo y depender de las distancias entre los electrodos.
K= 2 1 1 1 1 AM BM AN + BN
[3.5]
19
Ahora bien la resistividad calculada a travs de la expresin [3.4] , en general corresponder a la resistividad aparente del terreno y no a la resistividad real. Ello se debe a que el subsuelo esta compuesto generalmente por distintos materiales, por lo que el campo elctrico introducido en el terreno afectar simultneamente a diversas capas, de forma que la resistividad medida corresponder a un valor intermedio de todas ellas. Solo para suelos homogneos o para profundidades de investigacin muy superficiales, el valor de resistividad medido coincidir con el de la resistividad real.
20
Existen algunas frmulas empricas en las que se intenta expresar la relacin entre la resistividad y algunos de estos factores, como por ejemplo la porosidad o la saturacin del terreno. Una de las ms conocida y aplicable a acuferos limpios es la de Archie [5] .
r =
a w mS n
[3.6]
En donde S es el coeficiente de saturacin, w es la resistividad del fluido de los poros, r es la resistividad del terreno, la porosidad, n 2, 1.4 m 2.20 es el coeficiente de cimentacin y 0.6 a 1 . En la siguiente figura (Orellana [6] ) se exponen algunos de los rangos tpicos de resistividad, en donde de forma implcita ya estn recogidos la influencia de algunos de los factores anteriormente citados. No obstante destacar que los rangos pueden variar significativamente en funcin de las caractersticas del problema.
21
Batera
electrodos
Conectores
22
En las figuras 8 y 9 se muestran en detalle los elementos que conforman el dispositivo de medida. Corresponden al material utilizado en Cardona, en donde el resistivmetro utilizado es el modelo compacto SYSCAL JUNIOR V.II.2a (de la casa IRIS).
Fig. 8 Detalle del resistivmetro modelo SYSCAL Jr. (derecha) y de la batera (izquierda).
23
En el caso de querer cartografiar en planta las posibles trayectorias del agua, bastar con disponer transversalmente a la traza del flujo diferentes perfiles de electrodos a lo largo de su recorrido, de forma que su correlacin nos permitir establecer una imagen cuasi 3-D del recorrido del agua. En cuanto al nmero de electrodos ste es muy variable, si bien un nmero cercano a 50 electrodos es bastante frecuente. Estos se disponen de forma equidistante oscilando de 1 a 10 m para problemas ingenieriles (Griffiths [7]), si bien para problemas de filtraciones de agua un espaciado entre 1 y 2 metros es lo ms aconsejable. Dos aspectos muy importantes a tener en cuenta, son la profundidad mxima de investigacin y el rea de definicin del mtodo, aspectos ambos estrechamente ligados. En el Electrical imaging el campo de definicin disminuye con la profundidad (ver fig. 1). En consecuencia, a la hora de planificar el trabajo de campo hemos de tener presente esta prdida de definicin en los laterales del perfil. En el caso de querer confeccionar un perfil de gran longitud (mediante la unin de varios perfiles individuales), ser necesario establecer cierto solape entre perfiles a fin de intentar reducir esta indefinicin en los extremos. En cuanto a la profundidad mxima de investigacin, sta depende esencialmente de la longitud del perfil (a mayor longitud mayor profundidad), si bien el dispositivo electrdico utilizado tambin influir en cierta medida (ver 3.5.3).
24
Ahora, a travs del porttil configuraremos todas las variables del dispositivo: nmero de electrodos utilizados, distancia entre ellos, dispositivo electrdico de medida, as como el nmero de medidas que queremos realizar (destacar que podemos suprimir aquellas medidas que creamos oportunas, aspecto que supone un ahorro de tiempo significativo). Concluida la configuracin, se transfieren estas especificaciones a la unidad central para que sta pueda realizar automticamente toda la secuencia de mediciones prefijada. Sin embargo, antes de iniciar las mediciones y tambin mediante la unidad central, se procede previamente a la verificacin de que todas las conexiones funcionan correctamente, y que la resistividad entre electrodo-suelo es suficientemente baja como garantizar unos buenos resultados (en el trabajo prctico se impuso como valor mximo de potencial aceptable 2mV). En el caso de tener valores superiores generalmente basta con clavar un poco ms los electrodos, o humedecer el terreno para mejorar la conductividad. Una vez est todo dispuesto, se inicia la secuencia de medidas que hayamos predeterminado, proceso que tardar ms o menos tiempo en funcin del nmero de medidas a realizar. En general obtener del orden de 500 medidas comportar un tiempo de clculo cercano a la hora y media. Una vez terminado todo el proceso de captura de datos, toda la informacin almacenada digitalmente se vuelca en el porttil, a fin de proceder a su procesado e interpretacin con el programa correspondiente.
25
En la prctica, y una vez definidos los parmetros de estudio, el software nos permite visualizar la distribucin lateral y en profundidad de todas las medidas a lo largo del perfil, de forma que en funcin de nuestras exigencias, podremos adecuar la profundidad mxima aparente de investigacin eliminando los puntos de medida que creamos conveniente. El dispositivo electrdico tambin influye en la profundidad mxima de investigacin. En este sentido con Wenner-Schlumberger se consigue un aumento del 10-15% con respecto al obtenido por Wenner puro.
26
No obstante y dado que en problemas de filtracin de agua en el subsuelo, las variaciones de resistividad del terreno sern graduales (en especial en los suelos), esta tcnica no es eficiente. El Smoothness-constrained inversion method (deGroot-Hedin [10] ), es en realidad el mtodo de inversin ptimo para el anlisis de problemas en donde tengamos variaciones graduales de resistividad en el terreno, como por ejemplo en filtraciones de agua o vertidos de contaminantes. Adems este mtodo (basado en el Gauss-Newton least-squares method), presenta la gran ventaja de que no requiere de un modelo inicial, a partir del cual empezar el proceso iterativo. Dado su excesivo coste computacional, junto a la poca capacidad de los ordenadores de la poca, se desarroll el mtodo Quasi-Newton least-squares (Loke and Barker [11] ), que es aproximadamente unas 10 veces ms rpido que el mtodo convencional, si bien ste precisa que la variacin de los valores de resistividad sean ms suaves, a fin de poder obtener unos buenos resultados. Para el trabajo de cardona, el programa de inversin utilizado fue el RES2DINV ver. 3.50, software que corresponde a la casa Geotomo Software. En la referencia [12] se puede consultar un manual acerca del funcionamiento del programa, no obstante y a modo de sntesis el esquema bsico de funcionamiento es el siguiente: 1. Verificar los datos de campo. El programa permite visualizar las medidas obtenidas a fin de poder eliminar posibles valores errneos (se distinguen por presentar valores muy diferenciados al de los puntos adyacentes).
2. Definir el valor de los parmetros de clculo necesarios en el proceso de inversin. Los valores con los que trabaja por defecto el programa son satisfactorios en la mayora de los problemas.
27
3. Discretizacin del terreno en bloques (manual o automtica). Por defecto el programa trabaja con diferencias finitas, con un nmero de bloques inferior al de lecturas y cuyo tamao aumenta un 10% (o 25%) con la profundidad. No obstante podemos modificar manualmente la caractersticas de la discretizacin, as como trabajar con elementos finitos (ideal en el caso de tener topografa irregular)
4. Inicio de los clculos: El programa genera un modelo homogneo de resistividades. Posteriomente calcula el modelo de resistividades aparentes que se deriva del modelo de resistividades reales. Compara el modelo aparente calculado con el medido en el campo. Comienza un proceso iterativo con el fin de encontrar un modelo de resistividades reales, tal que la diferencia entre el modelo de resistividades aparentes calculado y el medido sea mnimo. El proceso finaliza cuando se alcance un RMS inferior al predefinido. No obstante es frecuente limitar el nmero de interaciones (por defecto el programa tiene un mximo de 5 iteraciones). Visualizacin de los resultados: pseudoseccin de resistividades aparentes medida en el campo, modelo de resistividades aparentes calculado por el ordenador y el modelo de resistividades reales.
28
29
Podemos controlar la profundidad mxima de estudio. Es un mtodo no destructivo. Las principales limitaciones de este mtodo son: Disponer del material requerido. A diferencia de otros mtodos, la Tomografa elctrica precisa de un equipo bastante caro. Como pasa en todos los mtodos de resistividad elctrica, la presencia de material muy resistivo en superficie imposibilita un buen funcionamiento del mtodo, siendo necesario la utilizacin de otras tcnicas, en este caso las electromagnticas, en los que no se precisa un contacto fsico con el suelo.
30