Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Particle in an electromagnetic field: The Lorentz-Dirac equation

A new method to solve the Lorentz-Dirac equation in the presence of an external electromagnetic field is presented. The validity of the approximation is discussed, and the method is applied to a particle in the presence of a constant magnetic field.

Invariant decomposition of the retarded electromagnetic field.

The integral representation of the electromagnetic two-form, defined on Minkowski space-time, is studied from a new point of view. The aim of the paper is to obtain an invariant criteria in order to define the radiative field. This criteria generalizes the well-known structureless charge case. We begin with the curvature two-form, because its field equations incorporate the motion of the sources. The gauge theory methods (connection one-forms) are not suited because their field equations do not incorporate the motion of the sources. We obtain an integral solution of the Maxwell equations in the case of a flow of charges in irrotational motion. This solution induces us to propose a new method of solving the problem of the nature of the retarded radiative field. This method is based on a projection tensor operator which, being local, is suited to being implemented on general relativity. We propose the field equations for the pair {electromagnetic field, projection tensor J. These field equations are an algebraic differential first-order system of oneforms, which verifies automatically the integrability conditions.

Interaccion among systems of finite size in predictive relativistic mechanics -II. Electromagnetic and gravitational interactions

We explicitly construct a closed system of differential equations describing the electromagnetic and gravitational interactions among bodies to first order in the coupling constants, retaining terms up to order c-2. The Breit and Barker and O'Connell Hamiltonians are recovered by means of a coordinate transformation. The method used throws light on the meaning of these coordinates.

SPATIAL DEPENDENCE OF ELECTRICAL CONDUCTIVITY AND CHEMICAL PROPERTIES OF THE SOIL BY ELECTROMAGNETIC INDUCTION

Brazilian soils have natural high chemical variability; thus, apparent electrical conductivity (ECa) can assist interpretation of crop yield variations. We aimed to select soil chemical properties with the best linear and spatial correlations to explain ECa variation in the soil using a Profiler sensor (EMP-400). The study was carried out in Sidrolândia, MS, Brazil. We analyzed the following variables: electrical conductivity - EC (2, 7, and 15 kHz), organic matter, available K, base saturation, and cation exchange capacity (CEC). Soil ECa was measured with the aid of an all-terrain vehicle, which crossed the entire area in strips spaced at 0.45 m. Soil samples were collected at the 0-20 cm depth with a total of 36 samples within about 70 ha. Classical descriptive analysis was applied to each property via SAS software, and GS+ for spatial dependence analysis. The equipment was able to simultaneously detect ECa at the different frequencies. It was also possible to establish site-specific management zones through analysis of correlation with chemical properties. We observed that CEC was the property that had the best correlation with ECa at 15 kHz.

Investigation of Electromagnetic Gauges for Determining In-Place HMA Density, Final Report, May 2007

Density is an important component of hot-mix asphalt (HMA) pavement quality and long-term performance. Insufficient density of an in-place HMA pavement is the most frequently cited construction-related performance problem. This study evaluated the use of electromagnetic gauges to nondestructively determine densities. Field and laboratory measurements were taken with two electromagnetic gauges—a PaveTracker and a Pavement Quality Indicator (PQI). Test data were collected in the field during and after paving operations and also in a laboratory on field mixes compacted in the lab. This study revealed that several mix- and project-specific factors affect electromagnetic gauge readings. Consequently, the implementation of these gauges will likely need to be done utilizing a test strip on a project- and mix-specific basis to appropriately identify an adjustment factor for the specific electromagnetic gauge being used for quality control and quality assurance (QC/QA) testing. The substantial reduction in testing time that results from employing electromagnetic gauges rather than coring makes it possible for more readings to be used in the QC/QA process with real-time information without increasing the testing costs.

Deep crustal electromagnetic structure of Bhuj earthquake region (India) and its implications

The existence of fluids and partial melt in the lower crust of the seismically active Kutch rift basin (on the western continental margin of India) owing to underplating has been proposed in previous geological and geophysical studies. This hypothesis is examined using magnetotelluric (MT) data acquired at 23 stations along two profiles across Kutch Mainland Uplift and Wagad Uplift. A detailed upper crustal structure is also presented using twodimensional inversion of MT data in the Bhuj earthquake (2001) area. The prominent boundaries of reflection in the upper crust at 5, 10 and 20 km obtained in previous seismic reflection profiles correlate with conductive structures in our models. The MT study reveals 1-2 km thick Mesozoic sediments under the Deccan trap cover. The Deccan trap thickness in this region varies from a few meters to 1.5 km. The basement is shallow on the northern side compared to the south and is in good agreement with geological models as well as drilling information. The models for these profiles indicate that the thickness of sediments would further increase southwards into the Gulf of Kutch. Significant findings of the present study indicate 1) the hypocentre region of the earthquake is devoid of fluids, 2) absence of melt (that is emplaced during rifting as suggested from the passive seismological studies) in the lower crust and 3) a low resistive zone in the depth range of 5-20 km. The present MT study rules out fluidsand melt (magma) as the causative factors that triggered the Bhuj earthquake. The estimated porosity value of 0.02% will explain 100-500 ohm·m resistivity values observed in the lower crust. Based on the seismic velocities and geochemical studies, presence of garnet is inferred. The lower crust consists of basalts - probably generated by partial melting of metasomatised garnet peridotite at deeper depths in the lithosphere - and their composition might be modified by reaction with the spinel peridotites.

Characterization and discrimination of soils by their reflected electromagnetic energy

The objective of this work was to verify if reflected energy of soils can characterize and discriminate them. A spectroradiometer (Spectral reflectance between: 400-2,500 nm) was utilized in laboratory. The soils evaluated are located in Bauru region, SP, Brazil, and are classified as Typic Argiudoll (TR), Typic Eutrorthox (LR), Typic Argiudoll (PE), Typic Haplortox (LE), Typic Paleudalf (PV) and Typic Quartzipsamment (AQ). They were characterized by their spectral reflectance as for descriptive conventional methods (Brazilian and International) according to the types of spectral curves. A method for the spectral descriptive evaluation of soils was established. It was possible to characterize and discriminate the soils by their spectral reflectance, with exception for LR and TR. The spectral differences were better identified by the general shape of spectral curves, by the intensity of band absorption and angle tendencies. These characteristics were mainly influenced by organic matter, iron, granulometry and mineralogy constituents. A reduction of iron and clay contents, which influenced higher reflectance intensity and shape variations, occurred on the soils LR/TR, PE, LE, PV and AQ, on that sequence. Soils of the same group with different clay textures could be discriminated. The conventional descriptive evaluation of spectral curves was less efficient on discriminating soils. Simulated orbital data discriminated soils mainly by bands 5 and 7.

Exposure to radio-frequency electromagnetic fields from broadcast transmitters and risk of childhood cancer: a census-based cohort study.

We investigated the association between exposure to radio-frequency electromagnetic fields (RF-EMFs) from broadcast transmitters and childhood cancer. First, we conducted a time-to-event analysis including children under age 16 years living in Switzerland on December 5, 2000. Follow-up lasted until December 31, 2008. Second, all children living in Switzerland for some time between 1985 and 2008 were included in an incidence density cohort. RF-EMF exposure from broadcast transmitters was modeled. Based on 997 cancer cases, adjusted hazard ratios in the time-to-event analysis for the highest exposure category (>0.2 V/m) as compared with the reference category (<0.05 V/m) were 1.03 (95% confidence interval (CI): 0.74, 1.43) for all cancers, 0.55 (95% CI: 0.26, 1.19) for childhood leukemia, and 1.68 (95% CI: 0.98, 2.91) for childhood central nervous system (CNS) tumors. Results of the incidence density analysis, based on 4,246 cancer cases, were similar for all types of cancer and leukemia but did not indicate a CNS tumor risk (incidence rate ratio = 1.03, 95% CI: 0.73, 1.46). This large census-based cohort study did not suggest an association between predicted RF-EMF exposure from broadcasting and childhood leukemia. Results for CNS tumors were less consistent, but the most comprehensive analysis did not suggest an association.

Investigation of Electromagnetic Gauges for Determination of In-Place Density of HMA Pavements, TR-576, 2009

This study evaluated the use of electromagnetic gauges to determine the adjusted densities of HMA pavements. Field measurements were taken with two electromagnetic gauges, the Pavement Quality Indicator (PQI) 301 and the Pavetracker Plus 2701B. Seven projects were included in the study with 3 to 5 consecutive paving days. For each day/lot 20 randomly selected locations were tested along with seven core locations. The analysis of PaveTracker and PQI density consisted of determining which factors are statistically significant, and core density residuals and a regression analysis of core as a function of PaveTracker and PQI readings. The following key conclusions can be stated: 1. Core density, traffic and binder content were all found to be significant for both electromagnetic gauges studied, 2. Core density residuals are normally distributed and centered at zero for both electromagnetic gauges, 3. For PaveTracker readings, statistically one third of the lots do not have an intercept that is zero and two thirds of the lots do not rule out a scaler correction factor of zero, 4. For PQI readings, statistically the 95% confidence interval rules out the intercept being zero for all seven projects and six of the seven projects do not rule out the scaler correction factor being zero, 5. The PQI 301 gauge should not be used for quality control or quality assurance, and 6. The Pavetracker 2701B gauge can be used for quality control but not quality assurance. This study has found that with the limited sample size, the adjusted density equations for both electromagnetic gauges were determined to be inadequate. The PaveTracker Plus 2701B was determined to be better than the PQI 301. The PaveTracker 2701B could still be applicable for quality assurance if the number of core locations per day is reduced and supplemented with additional PaveTracker 2701B readings. Further research should be done to determine the minimum number of core locations to calibrate the gauges each day/lot and the number of additional PaveTracker 2701B readings required.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

General electromagnetic simulation tool to predict the microwave nonlinear response of planar, arbitrarily-shaped HTS structures

This work describes a simulation tool being developed at UPC to predict the microwave nonlinear behavior of planar superconducting structures with very few restrictions on the geometry of the planar layout. The software is intended to be applicable to most structures used in planar HTS circuits, including line, patch, and quasi-lumped microstrip resonators. The tool combines Method of Moments (MoM) algorithms for general electromagnetic simulation with Harmonic Balance algorithms to take into account the nonlinearities in the HTS material. The Method of Moments code is based on discretization of the Electric Field Integral Equation in Rao, Wilton and Glisson Basis Functions. The multilayer dyadic Green's function is used with Sommerfeld integral formulation. The Harmonic Balance algorithm has been adapted to this application where the nonlinearity is distributed and where compatibility with the MoM algorithm is required. Tests of the algorithm in TM010 disk resonators agree with closed-form equations for both the fundamental and third-order intermodulation currents. Simulations of hairpin resonators show good qualitative agreement with previously published results, but it is found that a finer meshing would be necessary to get correct quantitative results. Possible improvements are suggested.