El método magnetotelúrico es una técnica muy importante dentro de los métodos geofísicos de inducción electromagnética. La sensibilidad de las medidas a la existencia de inhomogeneidades laterales hace que la interpretación unidimensional no sea la herramienta adecuada para la interpretación de estructuras geológicas complicadas. Por dicho motivo, una de las actuales líneas actuales de investigación corresponde a la elaboración de algoritmos de modelización numérica de estructuras conductoras bi- y tridimensionales.En este trabajo, presentamos un algoritmo de inversión bidimensional cuyo objetivo fundamental es la determinación de la forma de la frontera entre los diferentes conductores que caracterizan al medio.En este proceso, para conseguir una mayor simplicidad sin perder generalidad, hemos trabajado con modelos de dos medios y una frontera entre ambos a determinar.El problema ha sido abordado imponiendo que la frontera está descrita por una función expresable de forma analítica (f(y,z;p)= 0), de manera que los parámetros p de la curva son los parámetros del modelo bidimensional. Ello permite describir el modelo mediante un número reducido de parámetros, con lo cual se consigue una mayor estabilidad en la resolución del problema inverso.El trabajo ha estado estructurado básicamente en tres unidades.En la primera se ha revisado y descrito los fundamentos físicos del método magnetotelúrico.En la segunda se presenta el proceso de elaboración de un algoritmo de modelización basado en la resolución de las ecuaciones diferenciales mediante el método de los elementos finitos (MEF) para su utilización posterior en el problema inverso. El algoritmo ha sido comprobado utilizando el modelo de control de Weaver, LeQuang y Fischer, y ofrece unos resultados satisfactorios. De forma complementaria, se ha realizado un estudio de los efectos provocados por los accidentes topográficos asociados a formas simples y, de esta manera, estimar sus consecuencias.Los accidentes topográficos considerados han sido el de una cadena montañosa y el de un valle, para ambas polarizaciones E y H, y para perfiles de resisitividad aparente y de fase. Se ha observado que sus efectos son importantes en las proximidades de los cambios de pendiente y decrecen rápidamente al alejarse de los mismos. Igualmente se ha estudiado la proximidad de un mar cerrado y se observa una región para la que hay un cambio de comportamiento de los perfiles, que está originado por el hecho de que el mar sea cerrado.La tercera unidad del trabajo corresponde al problema inverso. En ella se describe, en primer lugar, el proceso de inversión utilizado, que está basado en los estimadores bayesianos que permiten introducir de forma natural la información a priori, dada por la geología u otros datos geofísicos. La no linealidad del problema provoca que la resolución del problema inverso se realice mediante un proceso iterativo. El comportamiento del mismo ha sido ampliamente comprobado para problemas unidimensionales.La realización del proceso iterativo para la resolución del problema inverso bidimensional implica el desarrollo de dos partes importantes para el mismo: el generador de mallas y el cálculo de la matriz de sensibilidad. El primero surge del hecho de que los parámetros del modelo son los parámetros de la curva que describe la frontera, y éstos cambiarán de una iteración a otra, por consiguiente, también se modificará la geometría del problema. Para que el proceso sea automático es preciso disponer de un generador de mallas. Su elaboración ha seguido el siguiente criterio: se construye en primer lugar una malla rectangular; a continuación, se deforma en las proximidades de la frontera, de manera que la curva que describe la frontera corte a la malla únicamente en los nudos; finalmente se unen dichos nudos, con lo que la frontera queda representada por una línea quebrada. Ello permite trabajar con elementos triangulares de conductividad constante.Por otra parte, en el algoritmo de inversión empleado es necesario calcular la matriz de sensibilidad del modelo para obtener el nuevo modelo a partir del antiguo. Dado que la resolución del problema directo se ha efectuado de forma numérica, esta matriz se ha de calcular de forma igualmente numérica. El proceso consiste en la resolución de tantos sistemas de ecuaciones lineales similares al obtenido al aplicar el MEF en la resolución del problema directo, como parámetros tenga el modelo. Para cada uno de ellos, únicamente es preciso calcular el término independiente del sistema de ecuaciones, porque la matriz del sistema coincide con la del problema directo. Además, el hecho de que la curva que describe la frontera se puede expresar de forma analítica permite calcular el término independiente de la ecuación de forma asimismo analítica. Todo ello lleva a una importante reducción en el tiempo de resolución del problema inverso.Dado que la matriz de sensibilidad es una pieza fundamental en el proceso de inversión, hemos estudiado su comportamiento frente a las características de la malla (tamaño de los elementos, proximidad de la frontera lateral) y hemos comprobado que sus efectos no son importantes.Finalmente, como comprobación del algoritmo, mostramos algunos ejemplos de inversión a partir de datos generados sintéticamente. Los modelos utilizados, y que pueden simular diferentes estructuras geológicas, han sido tres: el de una frontera descrita por una curva gaussiana, el una frontera descrita por una función salto inclinado y el de la frontera descrita por una elipse. En todos los casos se ha comprobado el adecuado funcionamiento del generador de mallas. En los mismos se observa la convergencia satisfactoria del algoritmo de inversión para los diferentes modelos y su correcto funcionamiento cuando se incorpora información a priori. / The MT method studies the conductive structure of the Earth's interior. The fact that this technique is strongly influenced by lateral inhomogeneities suggests the use of two- or three-dimensional models.In this work, we present an algorithm of 2-D inversion directed to those problems where the unknown is the shape of the boundary between different conductors. The way we use the inversion is considering that this shape can be described by analytical functions. This achieves a reduction of the number of parameters and provides a better stability of the inversion algorithm. This algorithm is based on Bayessian estimators that let us to introduce a priori information in a natural way.As an important part of this process we have developed the forward problem used by the inversion algorithm. It is performed numerically by the finite element method. As a complement of it, we have studied simple topographic effects: a range, a valley and the coast effect.The automatic inversion relies on two process, that we also describe and study: the mesh generator, that builds the mesh automatically in every iteration, and the sensibility matrix, that is necessary to calculate the new models from the old ones, which is obtained numerically.We present different examples with synthetic data that allow us to study the performance of the algorithm for simple shaped boundaries between media (Gaussian, ellipse and sloping step) that simulate different geological models. The results show that the algorithm works properly with and without a priori information.

No description available.

La colisión entre dos placas continentales plantea importantes cuestiones acerca de los procesos litosféricos que tienen lugar. Dependiendo del acortamiento litosférico la colisión puede tener asociada la subducción de una de las placas. La subducción no se debe tanto a efectos gravitatorios, ya que la corteza continental no es lo suficientemente densa para subducir en conjunto dentro del manto, sino más bien al campo de esfuerzos compresionales generado por el movimiento de aproximación entre placas (Ziegler, 1990). La convergencia y colisión entre las placas continentales ibérica y europea se manifestó en superficie con la edificación del orógeno pirenaico. Estudios realizados en los Pirineos basados en la compatibilidad del acortamiento de los dominios superficiales de la corteza (datos geológicos) y de los dominios profundos (datos geofísicos), indican un acortamiento superior a 100 kilómetros. Un problema todavía sin resolver es la acomodación en profundidad de este acortamiento. Una posible solución es la subducción de la corteza inferior de la placa ibérica bajo la placa europea, sin embargo los métodos geofísicos empleados hasta ahora no presentan una subducción de gran magnitud y no permiten decidir sobre esta cuestión. La conductividad eléctrica puede ayudar a solucionar las cuestiones planteadas ya que es un parámetro independiente de los estudiados anteriormente.Los métodos electromagnéticos se emplean en geofísica para determinar las propiedades de las estructuras geológicas en términos de la conductividad eléctrica ("sigma") o su inversa: la resistividad eléctrica ("ro"). Entre los métodos electromagnéticos el método magnetotelúrico (MT) es el que presenta una mayor resolución en los estudios corticales y litosféricos. Es un método de fuente natural que utiliza las fluctuaciones del campo electromagnético producido por la actividad solar en la ionosfera. A partir de la relación existente entre los campos eléctrico y magnético registrados simultáneamente en un punto de la superficie se puede inferir la distribución de resistividades eléctricas en profundidad.En la memoria de esta tesis se presenta la adquisición, procesado e interpretación de datos magnetotelúricos obtenidos en los Pirineos con el objetivo de determinar las propiedades electromagnéticas de la litosfera. Mediante el modelo de resistividades eléctricas se pueden restringir los posibles escenarios de evolución y estado actual de la cadena que permiten explicar conjuntamente los observables geofísicos y geológicos disponibles.Este trabajo está estructurado en once capítulos: en el primer capítulo se introduce el método magnetotelúrico. En el segundo capítulo se realiza un repaso de los fundamentos matemáticos de la teoría del electromagnetismo aplicada al método magnetotelúrico y se introduce el concepto de tensor de impedancias electromagnético. El tercer capítulo se centra en la caracterización y eliminación de los efectos de distorsión galvánica causada por cuerpos superficiales y locales sobre las funciones de transferencia regionales. Se presenta una de las aportaciones metodológicas de este trabajo, que consiste en la recuperación de las funciones de transferencia regionales en medios tridimensionales afectados por distorsiones galvánicas. En el capítulo cuarto se presenta la metodología de adquisición de los datos experimentales, así como el tratamiento y procesado robusto de las series temporales de los campo electromagnéticos.La aportación al conocimiento de la estructura Iitosférica de los Pirineos se realiza en los capítulos siguientes. En el capítulo quinto se describe la evolución geodinámica y formación del orógeno Pirenaico como una serie de pliegues y cabalgamientos, así como las principales estructuras geológicas. En el capítulo sexto se presenta la adquisición y procesado de los datos de Pirineos. El análisis del tensor de impedancias muestra que la estructura geoeléctrica de los Pirineos es bidimensional con una dirección de las estructuras principales este-oeste. Se realiza un estudio detallado del tensor de impedancias que muestra la presencia de tres grandes zonas o unidades que presentan unas características geoeléctricas propias.En el capítulo séptimo se presenta la interpretación bidimensional de los datos magnetotelúricos a lo largo de un perfil norte-sur cercano al perfil de sísmica de reflexión profunda ECORS-Pirineos. El resultado final de esta interpretación es un modelo bidimensional de resistividades eléctricas para la litosfera de los Pirineos. A nivel superficial, el modelo reproduce las principales estructuras geológicas: niveles de despegue entre la cobertera y el basamento en la Cuenca del Ebro y en la Cuenca de Aquitania, zonas de fractura de la Zona Axial, Falla Norpirenaica, etc... En profundidad, la disminución de la resistividad eléctrica en la placa ibérica a medida que nos acercamos al contacto entre placas culmina con la presencia de una estructura de elevada conductividad entre 30 y 80 km de profundidad en el contacto entre placas que se interpreta como fusión parcial de la corteza inferior de Iberia que subduce bajo la placa europea.La realización de nuevos registros hacia el oeste a lo largo de los Pirineos detecta la continuidad del conductor profundo. En el capítulo octavo se presenta la modelización bidimensional del perfil norte-sur de Belagua, situado al oeste del perfil ECORS. En el capítulo noveno se presenta un modelo tridimensional de resistividades eléctricas de la estructura litosférica de los Pirineos y se discuten sus interpretaciones tectónicas.En el capítulo décimo se comparan los modelos litosféricos obtenidos a partir del método magnetotelúrico con los modelos ya existentes obtenidos a partir de otros métodos geofísicos. En particular se realiza una interpretación conjunta con los resultados de tomografía sísmica obtenidos por Souriau y Granet (1995). Se discute la compatibilidad del modelo de conductividades eléctricas con los datos de la anomalía de Bouguer y de geoide. Para ello se realiza una modelización de la distribución de densidades litosféricas, a partir de la interpretación conjunta de la anomalía de Bouguer y la del geoide residual en los Pirineos. Para finalizar, en el capítulo undécimo se presentan las conclusiones generales del trabajo. / The Pyrenees is an Alpine collsional orogen formed as a result of the convergence between the European and Iberian plates. Restoration of upper crustal thrust sheets gives a minimum shortening of 150 km. How has this surface been shortening balanced by lower crustal and mantle rocks at depth? Subduction of the Iberian lower crust under the European Plate is a possible phenomenon that could explain it. However, geophysical studies up to present does not shown an important subduction. The distribution of the electrical conductivity at the lithosphere given by magnetotelluric method (MT) can help us to answer this question.This Thesis is concerned with the acquisition, processing and interpretation of MT data in the Pyrenees with the target of determine the lithospheric electrical properties and its structure. More than 50 sites have been measured are processed in order to obtain the magnerotelluric impedance tensor. The MT data soows the presence of an anomalous conductivity structure in the contact between plates at depth. This structure is related with the partial melting of me subducted lower crust. The amount of subduction obtained by MT modelling can explain the total shortening during the build up of the orogen and is compatible with:a) Restoration of upper crustal sheetsb) Numerical models of Pyrenean orogene) Negative anomaly in P-waves velocity (Seismic tomography)d) Bouguer anomalye) Geoid anomaly

La modelització numèrica és l'eina més potent per resoldre el problema directe de la prospecció elèctrica en models que presenten una distribución complicada de la resistivitat.En aquest treball es presenta el disseny d'un algorisme de modelització d'estructures bidimensionals per prospecció elèctrica en corrent continu. Les aplicacions de l'algorisme que es presenten cobreixen part del buit existent actualment en la modelització per al sondatge elèctric vertical amb dispositiu Schlumberger, que és, possiblement, el més utilitzat en la pràctica. 1) ALGORISME DE MODELITZACIÓL'algorisme de modelització que s'ha desenvolupat i el programa que s'ha construït a partir d'ell permeten l'estudi de qualsevol estructura bidimensional. L'algorisme es basa en el mètode dels elements finits i generalitza els anteriors algorismes existents basats en aquest mètode. Resol directament el problema pel potencial total i inclou la modelització d'estructures bidimensionals tant en presència d'elèctrodes de corrent lineals (problema bidimensional) com d'electròdes puntuals (problema tridimensional). A més, contempla diferents tipus possibles de condicions de contorn: condicions de Dirichlet, de Neumann i mixtes; en particular, unes mixtes basades en el comportament asimptòtic del potencial a grans distàncies de 1'elèctrode. Tanbé dóna solució a la problemàtica del càlcul pel potencial en direccions no perpendiculars a la de l 'estructura, en particular en la direcció paral.lela.L'estudi del problema tridimensional ha comportat la utilització del potencial espectral per tal d'eliminar la dependència del potencial amb la direcció de l'estructura ("strike direction"). A partir del potencial espectral es determina el potencial antritransformat. El potencial transformat s'ha estudiat en detall i s'han determinat criteris òptims per al seu càlcul numèric, que faciliten la resolució del potencial en la direcció paral.lela a l'estructura tridimensional. De l'anàlisi de les aproximacions incloses en l'algorisme i de les fonts d'errors inherents al propi mètode numèric, hem obtingut uns criteris metodològics per l'elecció de la xarxa de discretització. La validesa de l'algorisme i el correcte funcionament dels programes han estat comprovats pels models dels quals es disposa d'una solució analítica (models de control): model homogeni, models estratificats i un model dc contacte vertical entre dos medis. Pels models més complexos es realitza un estudi de la convergència i l'estabilitat de les solucions obtingudes numèricament. 2) APLICACIONSCom a aplicació de la modelització realitzada amb el nostre algorisme, es els casos més freqüents d'estructures bidimensionals, considerant a la vegada diferents configuracions del dispositiu electròdic. Aquests són:- L'estudi de l'anomalia creada per un cilindre resistitiu en un medi homogeni en 1'aproximació de potencial bidimensional i utilitzant el dispositiu gradient.- Utilitzant el dispositiu dipol-dipol s'afronta el problema de la detectabilitat d'un dic rectangular conductor en funció de les dimensions del dispositiu i de la fondària a la que es troba el dic.- S'ha dedicat especial atenció als diferents efectes i anomalies que es presenten en els sondatges elèctrics verticals amb dispositiu Schlumberger. Els models estudiats han estat:- Efectes d'inhomogeneitats a l'entorn dels elèctrodes de potencial.- Efecte d'un esglaó en la superfície.- Efectes laterals produïts per contactes verticals i desnivells.- Efectes laterals produïts per un "horst" en el basament.- Efectes topogràfics sobre medis estratificats.La resolució d'aquests models ha portat a l'elaboració d'un conjunt d 'àbacs que permeten mostrar la morfologia de les corbes de resistivitat aparent i observar la influència dels diferents paràmetres dels models. Aquests àbacs són de gran utilitat a l'hora de realitzar la interpretació i/o d'efectuar correccions d'aquests efectes.Els efectes depenen del model patricular considerat. Els resultats indiquen que aquests efectes no són en general menyspreables i s'han detenir en compte quan es fa la interpretación, tal com mostrem a l'exemple de la interpretació d'una estructura bidimensional (con és un horst) mitjancant un model unidimensional, on es tendeix a sobrevalorar la fondària del basament.Com a conclusió general es pot afirmar que s'ha donat una solució a cada problema específic de la modelització d'estructures bidimensionals per prospección elèctrica en corrent continu i que es disposa de la infrastructura per abordar cada problemàtica geològica concreta amb independencia de si és o no és aproximable al clàssic model estratificat. / In this work, an algorithm to direct current resistivity modelling for two-dimensional structures has been developed. The algorithm is based on the finite element method and it permits the use of different kinds of boundary conditions: Dirichlet conditions, Neumann conditions and mixed conditions, specially the ones using the assymptotical behaviour of the potential far to current electrode. Even though the resistivity distribution model is 2-D, due to the fact that source is a point, the problem is tridimensional and to avoid the dependence of the potential on the strike direction, the Helmholtz equation for the potential must be solved. We have studied the general behaviour of the potential's spectrum to find out some criteria to calculate the potential in any direction, in particular in the strike direction.The validity of the algorithm and the good working of the programs have been checked with models that have analytical solutions. We have studied for complex models the stability and the convergence of the numerical solution and we have made also an analysis of the different error sources and of the aproximations considered in the algorithm.Finally, we have used the algorithm to calculate apparent resistivity curves for VES over a stratified medium with several lateral and topographic effects (a hill, a valley, a step, small heterogeneities around the electrodes, a horst in the basement). We also presented other model curves for gradient array and dipol-dipol array. This systematic study can provide information about the general characteristics of the curves, the influence of the different parameters of the model and it allows have an insight to the situations where the lateral or the topographic effect must he taken into account.

La conservació i protecció del medi ambient estan adquirint cada cop més importància. Un dels factors que influeixen més en la seva degradació és la contaminació. Aquesta tesi es centra en la contaminació del subsòl. com a conseqüència de fuites de productes contaminants provinents de tuberies i dipòsits soterrats a poca profunditat. D'igual forma és interessant la detecció de fuites d'aigua, doncs cada vegada més s'està convertint en un bé escàs que cal preservar.Aquestes fuites poden ser detectades des de la superfície per mitjà de l'aplicació dels mètodes resistius, derivats dels de prospecció geoelèctrica. L'anàlisi de les corbes de resistivitat aparent obtingudes a la superfície del sòl i la visualització en plans verticals i horitzontals dels canvis de conductivitat elèctrica del subsòl són uns mètodes adients per poder-ne realitzar una monitorització.La tesi es desenvolupa en tres vessants: estudi del problema directe, estudi de la instrumentació per realitzar mesures i estudi del problema invers amb la conseqüent obtenció d'imatges.L'estudi del problema directe és el primer pas cap a la detecció de l'anomalia, ja que ens permet obtenir la distribució de potencials a la superfície del terreny tot i coneixent les resistivitats del terreny i de l'objecte que es vol detectar.Una forma de modelar els objectes és suposar que presenten la forma esfèrica o cilíndrica. La fuita, que es pot associar a una anomalia extensa, ja que es va estenent en el terreny, la suposarem concèntrica amb la canonada o dipòsit. Resolem l'equació de Laplace en un medi homogeni en el que hi ha un elèctrode puntual de corrent i N+1 cilindres concèntrics o N+1 esferes concèntriques de diferents conductivitats. Aquest estudi permet dissenyar una sèrie d'algorismes que generen diferents configuracions electròdiques.La particularització del problema directe a la configuració Schlumberger permet caracteritzar les anomalies per un moment bipolar normalitzat que respon a una fórmula de recurrència on es fa l'aproximació que es troba sotmesa a l'acció d'un camp elèctric uniforme. A partir de la resistivitat aparent corresponent a la configuració Schlumberger es pot discernir la grandària de la fuita i la seva profunditat.S'ha de tenir present que quasevol error a les tensions mesurades a la superficie poden afectar la resistivitat aparent, i per tant a la correcta detecció de la fuita. Alguns d'aquests errors són ocasionats pel posicionament dels elèctrodes, el soroll tel·lúric (sobretot en corrent continu) i l'acoblament electromagnètic en mesures en alterna. Per reduir aquest acoblament electromagnètic, els cables injectors i detectors haurien d'estar molt allunyats entre si, la qual cosa no és possible en totes les configuracions electròdiques d'interès i sobretot en les mesures de camp. Per aquesta raó hem considerat la possibilitat d'injectar un senyal quadrat de 1 kHz i una amplitud de pic a pic entre 20 V i 100 V.El sistema portàtil que hem contribuït a dissenyar, es basa en la injecció de senyals quadrats. La detecció es basa en un desmodulador síncron totalment diferencial. Es realitza un mostreig síncron en la part plana del senyal, de forma que es redueixen considerablement les interferències.Els algorismes dissenyats per reconstruir en forma d'imatges els canvis de conductivitat del subsòl són d'un sol pas. El mal condicionament de la matriu de sensibilitat requereix aplicar tècniques de regularització per realitzar la seva inversió. S'han estudiat els algorismes de retroprojecció i Marquardt-Levenberg, els quals s'han aplicat a dades sintètiques i a dades experimentals de laboratori i de camp. / La conservación y la protección del medio ambiente están adquiriendo cada vez más importancia. Uno de los factores que más influyen en su degradación es la contaminación. Esta tesis se centra en la detección de la contaminación del subsuelo como consecuencia de fugas de productos contaminantes provenientes de tuberías y depósitos enterrados a poca profundidad. También resulta muy importante la detección de fugas de agua, pues es un bien escaso que es necesario preservar. Estas fugas pueden ser detectadas desde la superficie por medio de la aplicación de los métodos resistivos. El análisis de las curvas de resistividad aparente y la visualización en planos verticales y horizontales de los cambios de conductividad eléctrica del subsuelo son unos métodos idóneos para poder realizar su monitorización. La tesis se desarrolla en tres vertientes: estudio del problema directo, estudio de la instrumentación para poder realizar las medidas y estudio del problema inverso para la obtención de imágenes. El modelo propuesto para la anomalía extensa ha consistido en solucionar la ecuación de Laplace en un medio homogéneo e isotrópico en el que se encuentra un electródo puntual de corriente y N+1 cilindros concéntricos o N+1 esferas concéntricas de diferentes conductividades. Este estudio permite diseñar una serie de algoritmos para generar diferentes configuraciones electródicas. La particularización del problema directo a la configuración Schlumberger permite caracterizar la anomalía por un momento dipolar normalizado, que responde a una fórmula de recurrencia, donde se ha considerado que se encuentra sometida a la influencia de un campo eléctrico uniforme. El estudio de los puntos característicos de las curvas de resistividad aparente permite obtener la profundidad y el tamaño de los objetos enterrados. Se ha contribuido a diseñar la instrumentación necesaria para realizar las medidas de campo de forma automática, la cual se basa en la inyección de una señal cuadrada de 1 kHz y una amplitud de pico a pico entre 20 V y 100 V. El muestreo síncrono de la parte plana de la señal reduce las interferencias por acoplamiento electromagnético entre el circuito inyector y el circuito detector, el cual realiza una demodulación síncrona totalmente diferencial de la señal detectada. Los algoritmos diseñados para reconstruir en forma de imágenes los cambios de conductividad del subsuelo son de un sólo paso. El mal condicionamiento de la matriz de sensibilidad requiere el uso de técnicas de regularización para realizar su inversión. Se han estudiado los algoritmos de retroproyección y de Marquardt-Levenberg, los cuales se han aplicado a datos sintéticos y a datos experimentales de laboratorio y de campo. / Electrical resistivity surveying is a versatile prospection method aplicable to groundwater location, archeological evaluation and contaminant monitoring.In this thesis we propose a methodology to detect leaks from buried pipes and underground storage tanks by means of resistive methods. The study consists on the analysis of the apparent resistivity curves and subsurface resistivity imaging. The main aspects are: resolution of the direct problem, design of the instrumentation needed to obtain experimental measurements and resolution of the inverse problem.We have performed a theoretical analysis of the electromagnetic problem where extensive anomalies are characterized by N+1 concentric cylinders or N+1 concentric spheres of different conductivities. This analysis allows us to design several algorithms to generate synthetic data.In order to obtain field measurements we have designed an instrumentation system. A portable personal computer controls current injection, electrode switching and voltage detection, which enables us to easily test different electrode arrays. The system injects 1 kHz, 20 V peak to peak square waveforms. Synchronously sampling the flat zone of square signals reduces electromagnetic coupling and electrode polarisation effects.We have developed several single step algorithms to image vertical and horizontal cross sections of the subsoil. Because the corresponding sensivity matrices are ill posed, we have applied several regularization techniques. The reconstruction algorithms that we have applied have been: retroprojection and Marquardt- Levenberg.First, the reconstruction algorithms have been applied to image synthetic data coming from an approximate analytical solution for a concentric spherical object embedded in a homogenous medium. This procedure has revealed as a powerful method to study the influence in the reconstructed image of the number of electrodes, number of pixels, depth of the anomaly, etc.Experimental measurements involving spherical and cylindrical objects immersed in a water tank validate the algorithms. Field measurements have allowed us to detect a simulated conductive leak from an insulating tube when voltages measured before the leak are used as reference.

Un sitio arqueológico urbano es una realidad diversa y compleja, cuyo conocimiento, protección y puesta en valor necesariamente debe de realizarse dentro de la concepción unitaria e integral que ofrece la Gestión del Patrimonio. En este marco general, la figura del parque ecoarqueológico se ha convertido en un modelo en donde es posible interrelacionar de forma eficaz los planes de gestión y los de investigación, con objeto de alcanzar una finalidad social. En esta tesis se defiende una propuesta metodológica integral para la investigación de un sitio arqueológico que se basa en la aplicación sistemática, secuencial y sumatoria de un conjunto de estudios fundamentados en la prospección arqueológica. Se ha tomado como ejemplo la serie de trabajos que se han desarrollado y puesto en práctica en el Parque Ecoarqueológico de Xoclán, en la ciudad de Mérida, Yucatán, que alberga vestigios arqueológicos de un sitio maya prehispánico. La aplicación de estos estudios ha permitido realizar, en un primer momento, una aproximación al sitio de carácter diagnóstico; y, posteriormente, obtener un conocimiento preciso del mismo con el menor impacto posible sobre el patrimonio cultural y natural del parque. De este modo, se ha proporcionado una información cuantitativa y cualitativa con respecto al sitio arqueológico y al espacio en donde éste se halla, la cual ha podido ser utilizada por los responsables del parque tanto para la gestión inmediata del mismo como para el diseño de futuras actuaciones. / An urban archaeological site is a diverse and complex reality, whose knowledge, protection and valorization must necessarily be conducted within the unitary and integral design that offers the Heritage Management. Within this general framework, the figure of the Ecoarchaeological Park has become a model where it’s possible to interact effectively management and research goals, in order to achieve a social objective. In this thesis, we defend an integral methodological proposition for the investigation of an archaeological site, based on the systematic, sequential and summation implementation of a set of studies based on the archaeological prospecting. The series of works which have been developed and implemented in the Xoclán Ecoarchaeological Park, in the city of Mérida, Yucatán, which hosts archaeological vestiges of a pre-Columbian maya site, had taken as an example. The application of these studies has allowed to realize, at first, a diagnostic approximation to the site; and, later, to obtain a precise knowledge of the same one with the least possible impact on the cultural and natural heritage of the Park. In this way, a quantitative and qualitative information on the archaeological site and the space where it is has been provided, which has been able to be used by the Park Administrations for the immediate management actions and for the design of future activities.