The region of West Bohemia/Vogtland in the Czech–German border area is well known for the repeated occurrence of earthquake swarms, CO2 emanations and mofette fields. To deepen the understanding of these phenomena local earthquake tomography of the Vp and Vp/Vs structure and attenuation tomography are carried out in this study. In comparison with previous investigations the travel time tomography revealed more details of the near-surface geology, potential fluid pathways and features around and below the swarm focal zone. In the uppermost crust, for the first time the Cheb basin and the Bublák/Hartoušov mofette fields were imaged as distinct anomalies of Vp and Vp/Vs. The well-pronounced low-Vp anomaly of the Cheb basin is not continuing into the Eger rift indicating a particular role of the basin within the rift system. A steep channel of increased Vp/Vs is interpreted as the pathway for fluids ascending from the earthquake swarm focal zone up to the Bublák/Hartoušov mofette fields. As a new feature, a mid-crustal body of high Vp and increased Vp/Vs is revealed just below and north of the earthquake swarm focal zone. It may represent a solidified intrusive body which emplaced prior or during the formation of the rift system. The enhanced fluid flow into the focal zone and triggering of earthquakes could be driven by the presence of the intrusive body if cooling is not fully completed. The assumed intrusive structure is considered as a heterogeneity leading to higher stress particularly at the junction of the rift system with the basin and prominent fault structures. This may additionally contribute to the triggering of earthquakes. The three-dimensional (3-D) P-wave attenuation (Qp) model for West Bohemia is the first of its kind. Path-averaged attenuation t * is calculated from amplitude spectra of time windows around the P-wave arrivals of local earthquakes. Average value or Qp for stations close to Nový Kostel are very low (< 150) compared to that of stations located further away from the focal zone (increases up to 500 within 80 km distance). The SIMUL2000 tomography scheme is used to invert the t * for P-wave attenuation perturbation. Analysis of resolution shows that the model is wellresolved in the vicinity of earthquake swarm hypocenters. The prominent features of the model are located around Nový Kostel focal zone and its northern vicinity. Beneath Nový Kostel a vertically stretched (down to depth of 11 km) and a highly attenuating body is observed. This might be due to fracturing and high density of cracks inside the weak earthquake swarm zone in conjunction with presence of free gas/fluid. Further north of Nový Kostel two high attenuating body are located at depths between 2 to 8 km which can represents trapped laterally distributed fluids. The eastern anomaly shows a good correlation with the fluid accumulation area which was suggested in 9HR seismic profile.

In this dissertation an approach is presented for the three-dimensional electrical resistivity tomography (ERT) using unstructured discretizations. The geoelectrical forward problem is solved by the finite element method using tetrahedral meshes with linear and quadratic shape functions. Unstructured meshes are suitable for modelling domains of arbitrary geometry (e.g., complicated topography). Furthermore, the best trade-off between accuracy and numerical effort can be achieved due to the capability of problem-adapted mesh refinement. Unstructured discretizations also allow the consideration of spatial extended finite electrodes. Due to a corresponding extension of the forward operator using the complete electrode model, known from medical impedance tomography, a study about the influence of such electrodes to geoelectrical measurements is given. Based on the forward operator, the so-called triple-grid-technique is developed to solve the geoelectrical inverse problem. Due to unstructured discretization, the ERT can be applied by using a resolution dependent parametrization on arbitrarily shaped two-dimensional and three-dimensional domains. A~Gauss-Newton method is used with inexact line search to fit the data within error bounds. A global regularization scheme is applied using special smoothness constraints. Furthermore, an advanced regularization scheme for the ERT is presented based on unstructured meshes, which is able to include a-priori information into the inversion and significantly improves the resulting ERT images. Structural information such as material interfaces known from other geophysical techniques are incorporated as allowed sharp resistivity contrasts. Model weighting functions can define individually the allowed deviation of the final resistivity model from given start or reference values. As a consequent further development the region concept is presented where the parameter domain is subdivided into lithological or geological regions with individual inversion and regularization parameters. All used techniques and concepts are part of the open source C++ library GIMLi, which has been developed during this thesis as an advanced tool for the method-independent solution of the inverse problem.

The solution to a set of simultaneous equations is of the form A-1 B = X where A-1 is the inverse of A in the equation AX= B. The purpose of this study is to obtain an exact A-1 through the use of rational arithmetic, and to study the behavior of rational numbers when used in arithmetic calculations. This study describes a matrix inversion program written in SPS II, utilizing the concept of rational arithmetic. This program, using the Gaussian elimination matrix inversion method, is compared to the same method written in Fortran. Gaussian elimination was used by this study because of its simplicity and speed of inversion. The Adjoint method was ruled out because of its complexity and relative lack of speed when compared with Gaussian elimination. The Fortran program gives only an approximate inverse due to the rounding error while the rational arithmetic program gives an exact inverse.

A three-dimensional inversion was implemented for magnetotellurics, which is a passive electromagnetic method in geophysics. It exploits natural electromagnetic fields of the Earth, which function as sources. Their interaction with the conductive parts of the subsurface are registered when components of the electric and the magnetic field are measured and evaluated. The all-at-once approach is an inversion scheme that is relatively new to geophysics. In this approach, the objective function – the basis of each inversion – is called the Lagrangian. It consists of three parts: (i) the data residual norm, (ii) the regularisation part, and (iii) the forward problem. The latter is the significant difference to conventional inversion approaches that are built up of a forward calculation part and an inversion part. In the case of all-at-once, the forward problem is incorporated in the objective function and is therefore already taken into account in each inversion iteration. Thus, an explicit forward calculation is obsolete. As an objective function, the Lagrangian shall reach a minimum and therefore its first and second derivatives are evaluated. Hence, the gradient of the Lagrangian and its Hessian are constituent parts of the KKT system – the Newton-type system that is set up in the all-at-once inversion. Conventional inversion approaches avoid the Hessian because it is a large, dense, not positive definite matrix that is challenging to handle. However, it provides additional information to the inversion, which raises hope for a high quality inversion result. As a first step, the inversion was programmed for the more straightforward one-dimensional magnetotelluric case. This was particularly suitable to become familiar with sQMR – a Krylov subspace method which is essential for the three-dimensional case to be able to work with the Hessian and the resulting KKT system. After the implementation and validation of the one-dimensional forward operator, the Lagrangian and its derivatives were set up to complete the inversion, which successfully solved the KKT system. Accordingly, the three-dimensional forward operator also needed to be implemented and validated, which was done using published data from the 3D-2 COMMEMI model. To realise the inversion, the Lagrangian was assembled and its first and second derivatives were validated with a test that exploits the Taylor expansion. Then, the inversion was initially programmed for the Gauss-Newton approximation where second order information is neglected. Since the system matrix of the Gauss-Newton approximation is positive definite, the solution of this system of equations could be carried out by the conventional solver pcg. Based on that, the complete KKT system (Newton\\\'s method) was set up and preconditioned sQMR solved this system of equations.

Ce travail porte sur la description et la compréhension de l’inversion de phase catastrophique utilisée pour l’émulsification de produits visqueux, à travers l’analyse des effets de formulation et de procédé sur la fraction de phase dispersée à laquelle le processus se produit et sur les mécanismes mis en jeu. Les suivis rhéologique et conductimétrique simultanés in situ de l’émulsification ont permis, du point de vue procédé, de mettre en avant l’influence du débit d’addition de la phase aqueuse sur la formation d’émulsions multiples du type e/H/E lesquelles, en augmentant notablement la fraction de phase dispersée apparente, sont responsables de l’inversion dès de faibles fractions de phase dispersée ajoutée. Au niveau formulation, l’augmentation de la viscosité de l’huile induit de manière remarquable la tendance de cette phase à devenir le milieu dispersé, conduisant à une inversion pour de très faibles fractions de phase aqueuse et donc à des émulsions finales très concentrées (de 80 à 95% en volume). Le suivi au microscope du phénomène d’inversion de phase par l’intermédiaire d’un écoulement de type « squeezing flow », a permis d’établir les conditions et les mécanismes conduisant à une inversion complète ou seulement partielle. La viscosité relative des phases aqueuse et huileuse est responsable d’une inversion catastrophique suivant un mécanisme de type agglomération – coalescence plutôt que de type inclusion/fuite tel que généralement admis. L’établissement d’un modèle mathématique basé sur les bilans de population et le caractère fractal du phénomène a permis de décrire l’évolution de la taille des gouttes multiples ainsi que la fraction de phase dispersée ajoutée à laquelle l’inversion se produit / This study deals with the description and understanding of catastrophic phase inversion during high viscous oil emulsification, through the analysis of the effects of formulation and process variables on the dispersed phase fraction at which the inversion is triggered and on the involved mechanisms. The simultaneous follow – up in situ of viscosity and conductivity measurements allowed, from a process point of view, to emphasize on the effect that the aqueous phase addition rate has on the formation of multiple w/O/W emulsions. Due to the formation of w/o/W emulsions, the volume of the effective dispersed phase greatly increases while at the same time, if the aqueous phase is added by very small fractions inversion of the w/o/W system can occur. In relation with formulation, an increase in oil viscosity greatly increased the tendency of the oily phase to become the dispersed phase. At the same time, it promoted the formation of highly concentrated emulsions (about 80 to 95 % in volume) after the inversion had occurred. The microscopical follow-up of emulsion morphology by means of squeezing flow, allowed us to establish the conditions and mechanisms that lead to partial or complete inversion. Even though literature sources lead us believe that inversion will occur through the “inclusion/escape” mechanism, experimental results showed that the relative viscosity between the phases promoted inversion through the mechanism of “agglomeration – coalescence” rather than “inclusion/escape”. A mathematical model based on population balances and on the fractal nature of multiple emulsions allowed us to describe the multiples drop size and effective dispersed phase evolution until inversion phenomena

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.

In this thesis, acoustic inversion methods are employed to estimate array element locations and the geoacoustic properties of the seabed using measured acoustic data consisting of noise from a surface ship in the Gulf of Mexico. The array element localization utilizes relative travel-time information obtained by cross-correlating the recorded time series of ship noise received at spatially separated hydrophones. The relative travel-time data are used in an inversion, based on the regularized least-squares method and the acoustic ray tracing equations, to obtain improved estimates of the receiver and source positions and their uncertainties. Optimization and Bayesian matched-field inversion methods are employed to estimate seabed geoacoustic properties and their uncertainties in the vicinity of a bottom-moored vertical line array using the recorded surface ship noise. This study is used to test the feasibility of matched-field methods to detect temporal changes in the geoacoustic properties of the seabed near a known gas hydrate mound in the Gulf of Mexico. Finally, a synthetic study is performed that demonstrates how ignoring environmental range dependence of seabed sound speed and water depth in matched-field inversion can lead to biases in the estimated geoacoustic parameters. The study considers the distributions of optimal parameter estimates obtained from a large number of range-independent inversions of synthetic data generated for random range-dependent environments. Range-independent Bayesian inversions are also performed on selected data sets and the marginal parameter distributions are examined. Both hard- and soft-bottom environments are examined at a number of scales of variability in sound speed and water depth.

L'augmentation des concentrations de méthane dans l'atmosphère, directement imputable aux activités anthropiques, induit une accentuation de l'effet de serre et une dégradation de la qualité de l'air. Il existe encore à l'heure actuelle de grandes incertitudes concernant les estimations des émissions des dfférentes sources de méthane à l'échellelocale. Une meilleure caractérisation de ces sources permettrait de mettre en place des politiques d'adaptation et d'att énuation efficaces afin de réduire ces émissions. Nous avons développé une nouvelle méthode de quantificationdes émissions de méthane à l'échelle locale basée sur la combinaison de mesures atmosphériques mobiles et d'un modèle gaussien dans le cadre d'une inversion statistique. Les concentrations atmosphériques du méthane sont mesuréesainsi que celles d'un gaz traceur émis à un flux connu. Ces concentrations en gaz traceur sont utilisées pour sélectionnerla classe de stabilité représentant le mieux les conditions atmosphériques dans le modèle gaussien ainsi qu'à paramétrerl'erreur associée aux mesures et au modèle dans l'inversion statistique. Dans un premier temps, cette nouvelle méthoded'estimation des émissions de méthane a été testée grâce à des émissions contrôlées de traceur et de méthane dontles sources ont été positionnées suivant différentes configurations. J'ai ensuite appliqué cette méthode à deux sites réels connus pour leurs émissions de méthane, une exploitation agricole et une installation de distribution de gaz, afin de tester son applicabilité et sa robustesse dans des conditions plus complexes de répartition des sources de méthane. Cette méthode a permis d'obtenir des estimations des émissions totales des sites robustes prenant en compte la localisation du traceur par rapport aux sources de méthane. L'estimation séparéedes émissions des différentes sources d'un site s'est révélée fortement dépendante des conditions météorologiques durant les mesures. Je me suis ensuite focalisé sur les émissions de méthane associées au secteur des déchets en réalisant un certain nombre de campagnes de mesures au sein d'installations de stockagedes déchets non dangereux et de stations d'épuration. Les résultats obtenus pour ces différents sites montrent la grandevariabilité des émissions de méthane dans le secteur des déchets. / The increase of atmospheric methane concentrations since the beginning of the industrial era is directly linked to anthropogenic activities. This increase is partly responsible for the enhancement of the greenhouse effect leading to a rise of Earth's surface temperatures and a degradation of air quality. There are still considerable uncertainties regarding methane emissions estimates from many sources at local scale. A better characterization of these sources would help the implementation of effective adaptation and mitigation policies to reduce these emissions.To do so, we have developed a new method to quantify methane emissions from local sites based on the combination of mobile atmospheric measurements, a Gaussian model and a statistical inversion. These atmospheric measurements are carried out within the framework of the tracer method, which consists in emitting a gas co-located with the methane source at a known flow. An estimate of methane emissions can be given by measuring the tracer and methane concentrations through the emission plume coming from the site. This method presents some limitations especially when several sources and/or extended sources can be found on the studied site. In these conditions, the colocation of the tracer and methane sources is difficult. The Gaussian model enables to take into account this bad collocation. It also gives a separate estimate of each source of a site when the classical tracer release method only gives an estimate of its total emissions. The statistical inversion enables to take into account the uncertainties associated with the model and the measurements.The method is based on the use of the measured tracer gas concentrations to choose the stability class of the Gaussian model that best represents the atmospheric conditions during the measurements. These tracer data are also used to parameterize the error associated with the measurements and the model in the statistical inversion. We first tested this new method with controlled emissions of tracer and methane. The tracer and methane sources were positioned in different configurations in order to better understand the contributions of this method compared to the traditional tracer method. These tests have demonstrated that the statistical inversion parameterized by the tracer gas data gives better estimates of methane emissions when the tracer and methane sources are not perfectly collocated or when there are several sources of methane.In a second time, I applied this method to two sites known for their methane emissions, namely a farm and a gas distribution facility. These measurements enabled us to test the applicability and robustness of the method under more complex methane source distribution conditions and gave us better estimates of the total methane emissions of these sites that take into account the location of the tracer regarding methane sources. Separate estimates of every source within the site are highly dependent on the meteorological conditions during the measurements. The analysis of the correlations on the posterior uncertainties between the different sources gives a diagnostic of the separability of the sources.Finally I focused on methane emissions associated with the waste sector. To do so, I carried out several measurement campaigns in landfills and wastewater treatment plants and I also used data collected on this type of sites during other projects. I selected the most suitable method to estimate methane emissions of each site and the obtained estimates for each one of these sites show the variability of methane emissions in the waste sector.

This thesis demonstrates that the composite cascode differential stage, operating in the subthreshold region, can form the basis of a high gain (113 dB) and low-power op amp (28.1 µW). The circuit can be fabricated without adding a compensation capacitance. The advantages of this architecture include high voltage gain, low bandwidth, low harmonic distortion, low quiescent current and power, and small chip area. These advantages suggest that this design might be well-suited for biomedical applications where low power, low noise bio-signal amplifiers capable of amplifying signals in the millihertz-to-kilohertz range is required.