Enhanced image reconstruction of electrical impedance tomography using simultaneous algebraic reconstruction technique and K-means clustering

Electrical impedance tomography (EIT), as a non-ionizing tomography method, has been widely used in various fields of application, such as engineering and medical fields. This study applies an iterative process to reconstruct EIT images using the simultaneous algebraic reconstruction technique (SART) algorithm combined with K-means clustering. The reconstruction started with defining the finite element method (FEM) model and filtering the measurement data with a Butterworth low-pass filter. The next step is solving the inverse problem in the EIT case with the SART algorithm. The results of the SART algorithm approach were classified using the K-means clustering and thresholding. The reconstruction results were evaluated with the peak signal noise ratio (PSNR), structural similarity indices (SSIM), and normalized root mean square error (NRMSE). They were compared with the one-step gauss-newton (GN) and total variation regularization based on iteratively reweighted least-squares (TV-IRLS) methods. The evaluation shows that the average PSNR and SSIM of the proposed reconstruction method are the highest of the other methods, each being 24.24 and 0.94; meanwhile, the average NRMSE value is the lowest, which is 0.04. The performance evaluation also shows that the proposed method is faster than the other methods

Applying Compactness Constraints to Differential Traveltime Tomography

Tomographic imaging problems are typically ill-posed and often require the use of regularization techniques to guarantee a stable solution. Minimization of a weighted norm of model length is one commonly used secondary constraint. Tikhonov methods exploit low-order differential operators to select for solutions that are small, flat, or smooth in one or more dimensions. This class of regularizing functionals may not always be appropriate, particularly in cases where the anomaly being imaged is generated by a non-smooth spatial process. Timelapse imaging of flow-induced velocity anomalies is one such case; flow features are often characterized by spatial compactness or connectivity. By performing inversions on differenced arrival time data, the properties of the timelapse feature can be directly constrained. We develop a differential traveltime tomography algorithm which selects for compact solutions i.e. models with a minimum area of support, through application of model-space iteratively reweighted least squares. Our technique is an adaptation of minimum support regularization methods previously explored within the potential theory community. We compare our inversion algorithm to the results obtained by traditional Tikhonov regularization for two simple synthetic models; one including several sharp localized anomalies and a second with smoother features. We use a more complicated synthetic test case based on multiphase flow results to illustrate the efficacy of compactness constraints for contaminant infiltration imaging. We conclude by applying the algorithm to a CO[subscript 2] sequestration monitoring dataset acquired at the Frio pilot site. We observe that in cases where the assumption of a localized anomaly is correct, the addition of compactness constraints improves image quality by reducing tomographic artifacts and spatial smearing of target features.Massachusetts Institute of Technology. Earth Resources Laborator

Non-Linear Constraints with Application to Self-Potential Source Inversion

We investigate the use of non-linear constraints for geophysical inverse problems, with specific examples applied to source inversion of self-potential data. Typical regularization methods often produce smooth solutions by introducing a quadratic term in the objective function that minimizes the L2 norm of a low-order differential operator applied to the model. In some cases, however, the properties of interest may not vary smoothly. Two alternative constraints are examined that provide inversion stability while allowing for solutions with non-smooth properties. One method, often referred to as ‘compactness’ or ‘minimum support’, seeks to minimize the area (in 2D) or volume (in 3D) occupied by non-zero model parameters. The second method, ‘total variation’, minimizes an approximation of the L1 norm of the gradient of the model. Both approaches involve a non-linear regularization functional, and must therefore be solved iteratively. We discuss the practical aspects of implementing these regularization methods and compare several examples using self-potential source inversion on a synthetic model. We also apply the compactness constraint for self-potential source inversion using a field data example.Kuwait-MIT Center for Natural Resources and the EnvironmentMassachusetts Institute of Technology. Earth Resources Laborator

Nonsmooth Nonconvex Variational Reconstruction for Electrical Impedance Tomography

Electrical Impedance Tomography is an imaging technique that aims to reconstruct the inner conductivity distribution of a medium starting from a set of measured voltages registered by a series of electrodes that are positioned on the surface of the medium. Such technique was used for the first time in geological studies in 1930 and then applied to industrial procedures. The first clinical use of EIT dates back to 1987. In 2018 EIT was validated in tissue engineering as a noninvasive and label-free imaging and monitoring tool for cell distribution (cell growth, differentiation and tissue formation) in 3D scaffolds. EIT problem can be split into a Forward and an Inverse problem. The aim of Forward EIT is to define the set of measured voltages starting from a known conductivity distribution. If the forward problem is characterized by a nonlinear mapping, called Forward Operator, from the conductivity distribution to the measured voltages, inverse EIT consists of inverting the Forward Operator. This leads to an ill-posed problem which requires regularization, either in the model or in the numerical method that is applied to define the solution. The inverse problem is modelled as a Nonlinear Least Squares problem, where one seeks to minimize the mismatch beetween the measured voltages and the ones generated by the reconstructed conductivity. Reconstruction techniques require the introduction of a regularization term which forces the reconstructed data to stick to certain prior information. In this dissertation, some state-of-the-art regularization methods are analyzed and compared via EIDORS, a specific software for EIT problems. The aim is to reconstruct the variation in conductivity within a 2D section of a 3D scaffold. Furthermore a variational formulation on a 2D mesh for a space-variant regularization is proposed, based on a combination of high order and nonconvex operators, which respectively seek to recover piecewise inhomogeneous and piecewise linear regions

Very Low Frequency Electrical Impedance Tomography Image Reconstruction System Using FPGA Software-Hardware Co-design

Electrical Impedance Tomography (EIT) is an imaging technique which is noninvasive and uses the internal conductivity distribution of the object of interest to form a tomographic image. It is performed by applying electrodes to the surface of the object. An alternating current up to frequency 10kHz is applied through a pair of electrodes, and the induced voltage is measured on other electrodes. These current and voltage values are used to reconstruct the internal conductivity distribution. The EIT imaging is increasingly getting used in clinical applications, as it is safer, portable, and low cost if compared with available imaging technologies used in clinical settings. The goal of this project is to develop a low frequency Zynq SoC-based EIT system. A Zynq 7020 device-based development board, Zedboard, interfaced with a customized hardware circuit, is used to develop a complete EIT system. A graphical user interface is developed using C# Graphic User Interface (GUI) application to control the hardware and visualize the results. It is a twelve-electrode system, and current injection and voltage measurement is performed through Zynq SoC. There are two image reconstruction algorithms developed, Gauss Newton One Step and Total variation. The algorithms are implemented in Zynq SoC using software-hardware co-design. The algorithms are also implemented in C#. The image reconstruction performance between the two algorithms is compared. The computation performance between Zynq SoC implementation and C# implementation is also compared to understand the feasibility of FPGA implementation of EIT image reconstruction algorithms

Time-Lapse Seismic and Electrical Monitoring of the Vadose Zone during a Controlled Infiltration Experiment at the Ploemeur Hydrological Observatory, France

The vadose zone is the main host of surface and subsurface water exchange and has important implications for ecosystems functioning, climate sciences, geotechnical engineering, and water availability issues. Geophysics provides a means for investigating the subsurface in a non-invasive way and at larger spatial scales than conventional hydrological sensors. Time-lapse hydrogeophysical applications are especially useful for monitoring flow and water content dynamics. Largely dominated by electrical and electromagnetic methods, such applications increasingly rely on seismic methods as a complementary approach to describe the structure and behavior of the vadose zone. To further explore the applicability of active seismics to retrieve quantitative information about dynamic processes in near-surface time-lapse settings, we designed a controlled water infiltration experiment at the Ploemeur Hydrological Observatory (France) during which successive periods of infiltration were followed by surface-based seismic and electrical resistivity acquisitions. Water content was monitored throughout the experiment by means of sensors at different depths to relate the derived seismic and electrical properties to water saturation changes. We observe comparable trends in the electrical and seismic responses during the experiment, highlighting the utility of the seismic method to monitor hydrological processes and unsaturated flow. Moreover, petrophysical relationships seem promising in providing quantitative results

Time-lapse cross-hole electrical resistivity tomography (CHERT) for monitoring seawater intrusion dynamics in a Mediterranean aquifer

Surface electrical resistivity tomography (ERT) is a widely used tool to study seawater intrusion (SWI). It is noninvasive and offers a high spatial coverage at a low cost, but its imaging capabilities are strongly affected by decreasing resolution with depth. We conjecture that the use of CHERT (cross-hole ERT) can partly overcome these resolution limitations since the electrodes are placed at depth, which implies that the model resolution does not decrease at the depths of interest. The objective of this study is to test the CHERT for imaging the SWI and monitoring its dynamics at the Argentona site, a well-instrumented field site of a coastal alluvial aquifer located 40¿km NE of Barcelona. To do so, we installed permanent electrodes around boreholes attached to the PVC pipes to perform time-lapse monitoring of the SWI on a transect perpendicular to the coastline. After 2 years of monitoring, we observe variability of SWI at different timescales: (1) natural seasonal variations and aquifer salinization that we attribute to long-term drought and (2) short-term fluctuations due to sea storms or flooding in the nearby stream during heavy rain events. The spatial imaging of bulk electrical conductivity allows us to explain non-monotonic salinity profiles in open boreholes (step-wise profiles really reflect the presence of freshwater at depth). By comparing CHERT results with traditional in situ measurements such as electrical conductivity of water samples and bulk electrical conductivity from induction logs, we conclude that CHERT is a reliable and cost-effective imaging tool for monitoring SWI dynamics.This work was funded by the project CGL2016-77122-C2-1-R/2-R of the Spanish Government. This project also received funding from the European Commission, Horizon 2020 research and innovation programme (Marie Sklodowska-Curie (grant no. 722028)). The author Albert Folch is a Serra Húnter Fellow.Peer ReviewedPostprint (published version

Regularisation methods for imaging from electrical measurements

In Electrical Impedance Tomography the conductivity of an object is estimated from boundary measurements. An array of electrodes is attached to the surface of the object and current stimuli are applied via these electrodes. The resulting voltages are measured. The process of estimating the conductivity as a function of space inside the object from voltage measurements at the surface is called reconstruction. Mathematically the ElT reconstruction is a non linear inverse problem, the stable solution of which requires regularisation methods. Most common regularisation methods impose that the reconstructed image should be smooth. Such methods confer stability to the reconstruction process, but limit the capability of describing sharp variations in the sought parameter. In this thesis two new methods of regularisation are proposed. The first method, Gallssian anisotropic regularisation, enhances the reconstruction of sharp conductivity changes occurring at the interface between a contrasting object and the background. As such changes are step changes, reconstruction with traditional smoothing regularisation techniques is unsatisfactory. The Gaussian anisotropic filtering works by incorporating prior structural information. The approximate knowledge of the shapes of contrasts allows us to relax the smoothness in the direction normal to the expected boundary. The construction of Gaussian regularisation filters that express such directional properties on the basis of the structural information is discussed, and the results of numerical experiments are analysed. The method gives good results when the actual conductivity distribution is in accordance with the prior information. When the conductivity distribution violates the prior information the method is still capable of properly locating the regions of contrast. The second part of the thesis is concerned with regularisation via the total variation functional. This functional allows the reconstruction of discontinuous parameters. The properties of the functional are briefly introduced, and an application in inverse problems in image denoising is shown. As the functional is non-differentiable, numerical difficulties are encountered in its use. The aim is therefore to propose an efficient numerical implementation for application in ElT. Several well known optimisation methods arc analysed, as possible candidates, by theoretical considerations and by numerical experiments. Such methods are shown to be inefficient. The application of recent optimisation methods called primal- dual interior point methods is analysed be theoretical considerations and by numerical experiments, and an efficient and stable algorithm is developed. Numerical experiments demonstrate the capability of the algorithm in reconstructing sharp conductivity profiles

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