Este trabalho tem como objetivo mapear os aquíferos sedimentar e cristalino na bacia sedimentar do Paraná no Estado de São Paulo, por meio do emprego da Inversão Conjunta de dados obtidos em sondagens com os métodos TDEM (Eletromagnético no Domínio do Tempo) e Eletrorresistividade por meio da técnica de SEV (Sondagem Elétrica Vertical). Os estudos foram realizados na região de Termas de Ibirá, distante cerca de 450 km do Município de São Paulo. Nesta região os poços tubulares perfurados captam água principalmente do aquífero sedimentar raso (Aquífero Bauru) relacionado com os sedimentos arenosos da Formação Adamantina, Grupo Bauru. Por outro lado, o aquífero fraturado na camada de basalto da Formação Serra Geral (Aquífero Serra Geral) ainda é pouco conhecido, sendo que este último constitui-se num excelente alvo para exploração de água subterrânea. Com as SEVs, obtiveram-se as informações mais rasas dos arenitos do Grupo Bauru, o Aquífero Bauru e o topo dos basaltos da Fm. Serra Geral. Com o método TDEM o pacote basáltico pôde ser investigado com maior detalhe e o Aquífero Serra Geral caracterizado. A análise dos resultados foi feita por meio da Inversão Conjunta-1D de sondagens SEVs/TDEM, de modo que as informações de cada metodologia foram integradas complementarmente, isto é, todo o intervalo de profundidade pôde ser caracterizado: a parte mais rasa com os dados das SEVs e a mais profunda com os dados do TDEM. Os resultados foram integrados com informações litológicas de poços de exploração de água subterrânea catalogados no SIAGAS-CPRM. Esta integração diminuiu as ambiguidades inerentes ao processo de interpretação de dados geofísicos, permitindo a demarcação de três zonas de favorabilidade para o posicionamento de poços tubulares profundos visando a exploração de água subterrânea no Aquífero Serra Geral e cinco zonas de favorabilidade para o Aquífero Bauru. Também foi possível caracterizar o comportamento estrutural da litoestragrafia na área de estudo. Deste modo, esses resultados contribuíram para o desenvolvimento das pesquisas hidrogeofísicas no Estado de São Paulo. / This work presents the geophysical mapping of sedimentary and fractured aquifers in Parana Sedimentary Basin, in Sao Paulo State, Brazil, obtained by 1D - Joint Inversion of TEM (Transient Electromagnetic) soundings and VES (Vertical Electrical Soudings). The studies were conducted in Termas de Ibira region, in Ibira City, which is located 450 km northwest of the city of Sao Paulo. The groundwater exploration wells in the area are located mainly into the shallow sedimentary aquifer (Bauru Aquifer), whose sediments and rocks are related to Adamantina Formation, Bauru Group. On the other hand, the fractured aquifer inside the basalts of Serra Geral Formation (Serra Geral Aquifer) is still not well known and has a huge groundwater exploration potential. The VES results can investigate the shallow sandstones from Bauru Group, the Bauru Aquifer and the top of Serra Geral basalt layer. The TEM soundings investigate the deepest depth interval inside the basalt layer and can examine in detail the fractured aquifer (basalts). The result analysis were integrated by 1D - Joint Inversion of VES/TEM soundings pairs, where the information of each depth interval was used in a complementary way, i.e., the Joint Inversion used the VES to solve shallow information and the TEM soundings the deepest information. Moreover, the results were integrated with the lithological information from the groundwater exploration wells registered in SIAGAS-CPRM database. The integration between Joint Inversion results and the lithological well information could decrease the inherent ambiguities in the geophysical interpretation process, and allowed to map five groundwater exploration favorabitility zones in the Bauru Aquifer and three favorability zones in the Serra Geral Aquifer. The results also provided the structural behavior of the area lithostrigraphy. Furthermore, this study can contribute to the development of hydrogeophysical research in countryside of Sao Paulo State.

QC 20100713

Le développement de l’énergie géothermique a conduit à l’exploitation de ressources établies dans des contextes géologiques et géodynamiques très variés. L’exploration géophysique de ces réservoirs complexes nécessite l’utilisation de plusieurs méthodes d’imagerie complémentaire. Ce travail de thèse porte sur l’exploration d’une ressource géothermique située en contexte de socle fracturé dans le Massif Central français par magnétotellurique, tomographie de bruit ambiant et gravimétrie.La magnétotellurique est une méthode d’imagerie 3D résolvante qui est sensible à la présence d’eau et aux argiles d’altération hydrothermale mais limitée par sa couverture spatiale. La tomographie de bruit sismique présente une bonne résolution verticale mais ne résout pas les variations horizontales de vitesse. Cette méthode est sensible aux variations des propriétés mécaniques des roches et donc aux milieux fracturés. Enfin la gravimétrie apporte une contrainte sur les variations lithologiques et possède une bonne résolution latérale mais une faible résolution verticale.Nous présentons une méthode d’inversion conjointe des données sismiques et gravimétriques sous contrainte d’un modèle de résistivité obtenu par inversion magnétotellurique indépendante. L’inversion conjointe nécessite de définir des couplages entre modèles. Par absence de connaissance a priori de relations pétrophysiques, nous avons couplé les modèles de densité, de résistivité et de vitesse avec une loi qui contraint les paramètres à être corrélés en moyenne. Cette stratégie vise à faire ressortir des relations caractéristiques des objets géologiques de la ressource géothermique.Cette méthodologie d’inversion conjointe a été testée sur des modèles synthétiques. L’application aux données réelles acquises dans le Massif Central a permis de définir une zone en profondeur de forte corrélation interprétée comme la transition ductile fragile. La partie intermédiaire des modèles, plus homogène, permet de distinguer différentes unités géologiques séparées par une zone de faille. Enfin la partie superficielle se distingue par une forte hétérogénéité des paramètres résultants probablement de processus d’altération de surface. / The development of geothermal energy has led to the exploitation of resources established in varied geological and geodynamic contexts. Geophysical exploration of these complex reservoirs requires the use of several complementary imaging methods. This PhD thesis focuses on the exploration of a geothermal resource located within the fractured basement in the French Massif Central using magnetotelluric, ambient noise tomography and gravimetry. Magnetotelluric is a 3D imaging method with a good resolution power that is sensitive to the presence of water and hydrothermal weathering clays but is limited by its spatial coverage. Seismic noise tomography has a good vertical resolution but does not resolve well horizontal velocity variations. This method is sensitive to variations of the mechanical properties of rocks and thus to fractured media. Finally gravimetry brings constraint on the lithological variations and has a good lateral resolution but lacks vertical resolution.We present a method of joint inversion of seismic and gravimetric data under the constraint of a resistivity model obtained by independent magnetotelluric inversion. Joint inversion requires defining model couplings. By lack of prior knowledge of petrophysical relationships, we have coupled the density, resistivity and velocity models with a law that constraints the parameters to be correlated on average.This strategy aims to bring out the characteristic relationships of the geological objects of the geothermal resource. This joint inversion methodology has been tested on synthetic models. The application to the real data acquired in the Massif Central has made it possible to define a deep zone of high correlation interpreted as the fragile ductile transition. The intermediate part of the models, more homogeneous, allows to distinguish different geological units separated by a fault zone. Finally the superficial part is distinguished by strong heterogeneity of the parameters resulting probably from surface alteration process.

L'inversion conjointe géophysique est la méthode la plus efficace pour imager l'intérieur de la Terre. En intégrant plusieurs techniques géophysiques elle permet de réduire les incertitudes inhérentes à chacune et ainsi améliorer la compréhension de la Terre. Dans cette étude, nous utilisons les techniques des fonctions récepteur (RF) en sismologie, de la magnétotellurique (MT) et de la gravimétrie qui permettent de caractériser respectivement la vitesse des ondes S, la résistivité électrique et la densité du sous-sol.Le but de ce travail de recherche se divise en deux parties: une première, méthodologique, sur le développement d'une nouvelle approche d'inversion conjointe en 3D et une deuxième avec l'application de ces techniques (en approche jointe ou séparée) sur la Divergence Nord Tanzanienne pour mieux comprendre le phénomène de rupture continentale. Pour la partie méthodologie deux approches ont été développées : une entre les données de MT et de gravimétrie avec un calcul original de l'effet gravimétrique de la topographie qui permet de réduire le nombre de mailles tout en gardant une résolution satisfaisante ; et une deuxième méthode entre les données de MT et de RF par une nouvelle approche d'extrapolation des modèles 1D de vitesse en pseudo modèle 3D de vitesse. L'application de ces techniques sur la Tanzanie a permis de mettre en évidence un certain nombre de structures lithosphériques dont deux zones majeures à faible vitesse dans la croûte inférieure et dans le manteau supérieur. Cette dernière semble refléter des interactions entre des structures héritées d'âge protérozoïque et le panache mantellique Africain. / Geophysical joint inversion attempts to reproduce as best as possible the interior of the Earth. By integrating several geophysical techniques the joint inversion reduces the uncertainties of each methods and improves our understanding of the Earth structure. In this study we use the receiver functions (RF), the magnetotelluric (MT) and the gravity methods which enable to charaterize the Swave velocity, the electrical resistivity and the density, respectively. The objective of this research work is divided in two parts; first with the development of a new 3D joint inversion approach and then with the application of these methods (on a joint or separate approach) on the North Tanzanian Divergence to better understand the continental breakup.For the methodologic part two approaches have been developed; one between the MT and gravity data with an original computation of the topographic effect which decreases the number of cells while keeping a satisfaying resolution. And a second method between the MT and RF data where pseudo 3D velocity model are created and combined with the MT models to better takes into account the physical properties of the receiver function. The application of these methods on the Tanzania highlighted several lithospheric structures and particularly two low-velocity areas in the lower crust and the upper mantle. This latter suggests interactions with Proterozoic inherited structures and the African plume material.

As sondagens elétricas verticais (SEV) e as sondagens eletromagnéticas no domínio do tempo (TDEM) são muito usadas em estudos ambientais, hidrogeológicos e em exploração mineral. A interpretação desses dados feita individualmente com modelos unidimensionais normalmente acarreta em resultados ambíguos. Isso acontece devido ao fato de que cada metodologia enxerga a subsuperfície de uma forma diferente. No caso da SEV as estruturas resistivas são bem detectadas, já no caso da sondagem TDEM as estruturas condutivas é que são detectadas com melhor precisão. Outra diferença é a capacidade da SEV conseguir identificar melhor as estruturas rasas, por outro lado, as sondagens TDEM permitem investigar as estruturas mais profundas. Nesta pesquisa foram exploradas as potencialidades das SEVs e das sondagens TDEM, visando obter uma interpretação dos dados mais consistentes. Neste sentido foi desenvolvido um programa em Matlab para a inversão conjunta 1D de dados de SEV e TDEM, na qual explora \"o melhor\" de ambos os métodos. O algoritmo desenvolvido foi testado inicialmente com dados sintéticos e posteriormente foi empregado em dados reais adquiridos na bacia sedimentar do Paraná, na região de Bebedouro no Estado de São Paulo. O resultado final da inversão conjunta de dados de SEV/TDEM permitiu obter um modelo geoelétrico que mais se assemelha às condições geológicas reais, e cujas ambigüidades, que são inerentes ao processo de interpretação, foram minimizadas. Os resultados tanto com dados sintéticos quanto com dados reais foram são bastante promissores, mostram uma melhor recuperação dos meios modelados e um grande potencial de aplicações em estudos geocientíficos, em particular em estudos hidrogeológicos. / Electrical (DC) and transient electromagnetic (TEM) soundings have been used in a great number of environmental, hydrological and mining exploration studies. The data interpretation usually is done individually by 1D models resulting in ambiguous results. This fact can be explained by how the two different methodologies observe the medium below the surface. The vertical electrical sounding (VES) is good at marking very resistive structures, while the transient electromagnetic sounding (TEM) is very sensitive to map conductive structures. Another difference is that the VES detects better shallow structures, while the TEM soundings can reach deeper ones. In this research we explore the potentials of the both VES and TEM soundings, in order to obtain a more consistent interpretation of the data. In this way, a Matlab program for the joint inversion for DC and TEM soundings was developed aiming explore the best of the both methods. Initially, the algorithm was tested with synthetic data and after were used real data from Paraná Sedimentary Basin in Bebedouro region, São Paulo State. The geoelectrical model obtained from joint inversion of DC and TEM data are more similar to the real geological condition and the ambiguities were minimized. The results with synthetic and real data shows that the joint inversion of DC/TEM is better for recovering the simulated models and shows a great potential in geological studies, particularly in hidrogeological studies.

As sondagens elétricas verticais (SEV) e as sondagens eletromagnéticas no domínio do tempo (TDEM) são muito usadas em estudos ambientais, hidrogeológicos e em exploração mineral. A interpretação desses dados feita individualmente com modelos unidimensionais normalmente acarreta em resultados ambíguos. Isso acontece devido ao fato de que cada metodologia enxerga a subsuperfície de uma forma diferente. No caso da SEV as estruturas resistivas são bem detectadas, já no caso da sondagem TDEM as estruturas condutivas é que são detectadas com melhor precisão. Outra diferença é a capacidade da SEV conseguir identificar melhor as estruturas rasas, por outro lado, as sondagens TDEM permitem investigar as estruturas mais profundas. Nesta pesquisa foram exploradas as potencialidades das SEVs e das sondagens TDEM, visando obter uma interpretação dos dados mais consistentes. Neste sentido foi desenvolvido um programa em Matlab para a inversão conjunta 1D de dados de SEV e TDEM, na qual explora \"o melhor\" de ambos os métodos. O algoritmo desenvolvido foi testado inicialmente com dados sintéticos e posteriormente foi empregado em dados reais adquiridos na bacia sedimentar do Paraná, na região de Bebedouro no Estado de São Paulo. O resultado final da inversão conjunta de dados de SEV/TDEM permitiu obter um modelo geoelétrico que mais se assemelha às condições geológicas reais, e cujas ambigüidades, que são inerentes ao processo de interpretação, foram minimizadas. Os resultados tanto com dados sintéticos quanto com dados reais foram são bastante promissores, mostram uma melhor recuperação dos meios modelados e um grande potencial de aplicações em estudos geocientíficos, em particular em estudos hidrogeológicos. / Electrical (DC) and transient electromagnetic (TEM) soundings have been used in a great number of environmental, hydrological and mining exploration studies. The data interpretation usually is done individually by 1D models resulting in ambiguous results. This fact can be explained by how the two different methodologies observe the medium below the surface. The vertical electrical sounding (VES) is good at marking very resistive structures, while the transient electromagnetic sounding (TEM) is very sensitive to map conductive structures. Another difference is that the VES detects better shallow structures, while the TEM soundings can reach deeper ones. In this research we explore the potentials of the both VES and TEM soundings, in order to obtain a more consistent interpretation of the data. In this way, a Matlab program for the joint inversion for DC and TEM soundings was developed aiming explore the best of the both methods. Initially, the algorithm was tested with synthetic data and after were used real data from Paraná Sedimentary Basin in Bebedouro region, São Paulo State. The geoelectrical model obtained from joint inversion of DC and TEM data are more similar to the real geological condition and the ambiguities were minimized. The results with synthetic and real data shows that the joint inversion of DC/TEM is better for recovering the simulated models and shows a great potential in geological studies, particularly in hidrogeological studies.

No description available.

The objective of our studies is the combination of electromagnetic and direct current (DC) resistivity methods in a joint inversion approach to improve the reconstruction of a given conductivity distribution. We utilize the distinct sensitivity patterns of different methods to enhance the overall resolution power and ensure a more reliable imaging result. In order to simplify the work with more than one electromagnetic method and establish a flexible and state-of-the-art software basis, we developed new DC resistivity and electromagnetic forward modeling and inversion codes based on finite elements of second order on unstructured grids. The forward operators are verified using analytical solutions and convergence studies before we apply a regularized Gauss-Newton scheme and successfully invert synthetic data sets. Finally, we link both codes with each other in a joint inversion. In contrast to most widely used joint inversion strategies, where different data sets are combined in a single least-squares problem resulting in a large system of equations, we introduce a sequential approach that cycles through the different methods iteratively. This way, we avoid several difficulties such as the determination of the full set of regularization parameters or a weighting of the distinct data sets. The sequential approach makes use of a smoothness regularization operator which penalizes the deviation of the model parameters from a given reference model. In our sequential strategy, we use the result of the preceding individual inversion scheme as reference model for the following one. We successfully apply this approach to synthetic data sets and show that the combination of at least two methods yields a significantly improved parameter model compared to the individual inversion results. / Ziel der vorliegenden Arbeit ist die gemeinsame Inversion (\"joint inversion\") elektromagnetischer und geoelektrischer Daten zur Verbesserung des rekonstruierten Leitfähigkeitsmodells. Dabei nutzen wir die verschiedenartigen Sensitivitäten der Methoden aus, um die Auflösung zu erhöhen und ein zuverlässigeres Ergebnis zu erhalten. Um die Arbeit mit mehr als einer Methode zu vereinfachen und eine flexible Softwarebasis auf dem neuesten Stand der Forschung zu etablieren, wurden zwei Codes zur Modellierung und Inversion geoelektrischer als auch elektromagnetischer Daten neu entwickelt, die mit finiten Elementen zweiter Ordnung auf unstrukturierten Gittern arbeiten. Die Vorwärtsoperatoren werden mithilfe analytischer Lösungen und Konvergenzstudien verifiziert, bevor wir ein regularisiertes Gauß-Newton-Verfahren zur Inversion synthetischer Datensätze anwenden. Im Gegensatz zur meistgenutzten \"joint inversion\"-Strategie, bei der verschiedene Daten in einem einzigen Minimierungsproblem kombiniert werden, was in einem großen Gleichungssystem resultiert, stellen wir schließlich einen sequentiellen Ansatz vor, der zyklisch durch die einzelnen Methoden iteriert. So vermeiden wir u.a. eine komplizierte Wichtung der verschiedenen Daten und die Bestimmung aller Regularisierungsparameter in einem Schritt. Der sequentielle Ansatz wird über die Anwendung einer Glättungsregularisierung umgesetzt, bei der die Abweichung der Modellparameter zu einem gegebenen Referenzmodell bestraft wird. Wir nutzen das Ergebnis der vorangegangenen Einzelinversion als Referenzmodell für die folgende Inversion. Der Ansatz wird erfolgreich auf synthetische Datensätze angewendet und wir zeigen, dass die Kombination von mehreren Methoden eine erhebliche Verbesserung des Inversionsergebnisses im Vergleich zu den Einzelinversionen liefert.

The objective of our studies is the combination of electromagnetic and direct current (DC) resistivity methods in a joint inversion approach to improve the reconstruction of a given conductivity distribution. We utilize the distinct sensitivity patterns of different methods to enhance the overall resolution power and ensure a more reliable imaging result. In order to simplify the work with more than one electromagnetic method and establish a flexible and state-of-the-art software basis, we developed new DC resistivity and electromagnetic forward modeling and inversion codes based on finite elements of second order on unstructured grids. The forward operators are verified using analytical solutions and convergence studies before we apply a regularized Gauss-Newton scheme and successfully invert synthetic data sets. Finally, we link both codes with each other in a joint inversion. In contrast to most widely used joint inversion strategies, where different data sets are combined in a single least-squares problem resulting in a large system of equations, we introduce a sequential approach that cycles through the different methods iteratively. This way, we avoid several difficulties such as the determination of the full set of regularization parameters or a weighting of the distinct data sets. The sequential approach makes use of a smoothness regularization operator which penalizes the deviation of the model parameters from a given reference model. In our sequential strategy, we use the result of the preceding individual inversion scheme as reference model for the following one. We successfully apply this approach to synthetic data sets and show that the combination of at least two methods yields a significantly improved parameter model compared to the individual inversion results. / Ziel der vorliegenden Arbeit ist die gemeinsame Inversion (\"joint inversion\") elektromagnetischer und geoelektrischer Daten zur Verbesserung des rekonstruierten Leitfähigkeitsmodells. Dabei nutzen wir die verschiedenartigen Sensitivitäten der Methoden aus, um die Auflösung zu erhöhen und ein zuverlässigeres Ergebnis zu erhalten. Um die Arbeit mit mehr als einer Methode zu vereinfachen und eine flexible Softwarebasis auf dem neuesten Stand der Forschung zu etablieren, wurden zwei Codes zur Modellierung und Inversion geoelektrischer als auch elektromagnetischer Daten neu entwickelt, die mit finiten Elementen zweiter Ordnung auf unstrukturierten Gittern arbeiten. Die Vorwärtsoperatoren werden mithilfe analytischer Lösungen und Konvergenzstudien verifiziert, bevor wir ein regularisiertes Gauß-Newton-Verfahren zur Inversion synthetischer Datensätze anwenden. Im Gegensatz zur meistgenutzten \"joint inversion\"-Strategie, bei der verschiedene Daten in einem einzigen Minimierungsproblem kombiniert werden, was in einem großen Gleichungssystem resultiert, stellen wir schließlich einen sequentiellen Ansatz vor, der zyklisch durch die einzelnen Methoden iteriert. So vermeiden wir u.a. eine komplizierte Wichtung der verschiedenen Daten und die Bestimmung aller Regularisierungsparameter in einem Schritt. Der sequentielle Ansatz wird über die Anwendung einer Glättungsregularisierung umgesetzt, bei der die Abweichung der Modellparameter zu einem gegebenen Referenzmodell bestraft wird. Wir nutzen das Ergebnis der vorangegangenen Einzelinversion als Referenzmodell für die folgende Inversion. Der Ansatz wird erfolgreich auf synthetische Datensätze angewendet und wir zeigen, dass die Kombination von mehreren Methoden eine erhebliche Verbesserung des Inversionsergebnisses im Vergleich zu den Einzelinversionen liefert.

Nesta pesquisa foram analisados dados utilizando os métodos geofísicos Eletromagnético no Domínio do Tempo (TDEM) e Eletrorresistividade (ER) adquiridos na região central da Bacia sedimentar de Taubaté-SP com a finalidade de contribuir com estudos hidrogeológicos e realizar uma caracterização geoelétrica da área. Os dados TDEM foram obtidos por meio da técnica de caminhamento TDEM (CTDEM) usando um loop transmissor (Tx) fixo e uma bobina receptora (Rx) - 3D móvel, também conhecida como técnica fixed-loop. Os dados de ER foram coletados utilizando a técnica de Sondagem Elétrica Vertical (SEV) e perfis de Caminhamento Elétrico (CE). A interpretação das camadas geoelétricas da subsuperfície foram realizadas por meio de inversões individuais TDEM, inversões conjunta 1D de SEV/TDEM coincidentes e de pseudoseções CTDEM geradas através da interpolação dos resultados das inversões individuais TDEM, permitindo, desta forma, gerar um modelo geoelétrico para a região específica dentro do loop transmissor. A área de estudos localiza-se numa área rural da cidade de Taubaté, SP, de maneira que os resultados dos perfis de CTDEM integrados com as inversões conjuntas SEV/TDEM e informações litológicas de poços permitiram mapear o contato entre os sedimentos Quaternários e Terciários (Grupo Taubaté e Formação Pindamonhangaba). O topo do embasamento da bacia na área de estudos foi inferido com base nas informações litológicas de poços catalogados e disponibilizados pelo Sistema de Informações de Águas Subterrâneas e pelo Serviço Geológico Brasileiro (SIAGAS-CPRM). Os resultados apresentam-se de acordo com as informações geológicas e contribuem com as pesquisas hidrogeofísicas no Estado de São Paulo. / In this research the geophysics data were analyzed through the Time Domain Electromagnetic Method (TDEM) and Electrical Method (DC), which one were acquired in central region of sedimentary Taubaté basin, SP, for the purpose to contribute to hydrogeological studies and perform a geoelectrical characterization of this area. TDEM data were obtained by 3D mobile receiver coil (Rx) and transmitter fixed-loop techniques. The DC data were collected using Vertical Electric Sounding (VES) and Electrical Resistivity Tomography (ERT) techniques. The geoelectrical layers interpretation of subsurface were performed by TDEM individual invertions and 1D joint inversions VES/TDEM and the TDEM profiles were generated through interpolation of TDEM individual inversions, allowing to generate a geoelectrical model to the specific region inside the loop transmitter. The study area is located around Taubaté city, SP, in a rural area, so that TDEM profile results associated with joint inversions VES/TDEM and wells lithological informations allows mapping the contact between Quaternary and Tertiary (Taubaté Group and Pindamonhangaba Formation) sediments. The top of the basement in the studies area was inferred through of wells informations available in Groundwater Information System and Brazilian Geological Service (SIAGAS-CPRM). The results are consistent with geological informations, contributing with hydrogeophysical studies at São Paulo state.