• Refine Query
  • Source
  • Publication year
  • to
  • Language
  • 21
  • 11
  • 3
  • Tagged with
  • 32
  • 10
  • 9
  • 7
  • 7
  • 6
  • 5
  • 5
  • 5
  • 5
  • 4
  • 4
  • 4
  • 4
  • 4
  • About
  • The Global ETD Search service is a free service for researchers to find electronic theses and dissertations. This service is provided by the Networked Digital Library of Theses and Dissertations.
    Our metadata is collected from universities around the world. If you manage a university/consortium/country archive and want to be added, details can be found on the NDLTD website.
11

Indécomposabilité des graphes et des tournois / Indecomposability of graphs and tournaments

Belkhechine, Houmem 15 July 2009 (has links)
Cette thèse porte sur l’indécomposabilité dans les graphes et les tournois. Elle comporte cinq chapitres dont le premier est introductif. Le deuxième chapitre consiste en une étude des tournois indécomposables suivant les tournois indécomposables à 5 ou à 7 sommets qu’ils abritent [3, 2]. Le troisième chapitre est une caractérisation des tournois (-1)-critiques avec une description morphologique de ces tournois [4,5]. Le quatrième chapitre contient une caractérisation des graphes (-1)-critiques [6], répondant ainsi, dans le cas général, à un problème posé par Y. Boudabbous et P. Ille[10]. Le cinquième chapitre est consacré à une opération d’inversion dans les tournois et un invariant, l’indice d’inversion d’un tournoi, dont l’étude a été proposée par M. Pouzet. Le fait que les tournois (-1)-critiques sont d’indice entre 2 et 4 est le lien avec l’étude de la criticalité. Plusieurs propriétés de la classe des tournois d’indice au plus n sont données. / This thesis focuses on indecomposability in graphs and tournaments. It contains five chapters, The first is introductory. The second chapter is a study of indecomposable tournaments according to the indecomposable subtournaments on 5 or 7 vertices embedding in. The third chapter is a characterization of the (-1)-critical tournaments. The fourth chapter contains a characterization of the (-1)-critical graphs answering in the general case, to a problem asked by Y. Boudabbous and P. Ille. The fifth chapter is devoted to an inversion operation in tournaments and an invariant, the inversion index of a tournament, which the study was proposed by M. Pouzet. The fact that the inversion index of a (-1)-critical tournament is between 2 and 4 is the link with the study of criticality. Many properties of the class of tournaments with inversion index at most n are given
12

Etude de l'atmosphère profonde de Jupiter / Study of the deep atmosphere of Jupiter

Blain, Doriann 10 September 2018 (has links)
L'objet de cette thèse est l'étude des répartitions spatiales des abondances chimiques de la troposphère de Jupiter. Le cas du NH3 est particulièrement étudié. Deux types de données ont été utilisées : des observations au sol via l'instrument TEXES monté sur le télescope IRTF, et des observations spatiales via l'instrument JIRAM embarqué sur la sonde JUNO. Les données sont analysées par un code parallélisé développé au cours de la thèse, combinant un transfert radiatif raie-par-raie et une méthode d'inversion non-linéaire optimale. Les données TEXES ont permis d'établir une carte du NH3 troposphérique (entre 1 et 4 bar) à une précision de l'ordre de 20% sur environ trois quarts de la surface de Jupiter à une résolution spatiale d'environ 1200 km (1,4"). Cette carte met en évidence d'importantes variations latitudinales de l'abondance de NH3 troposphérique, qui peut varier de 60 ppmv dans la NEB à proximité des "hotspots", à plus de 500 ppmv dans les zones. Les résultats dans les zones souffrent toutefois d'un faible rapport signal-sur-bruit, ce qui augmente leur incertitude. La carte met également en évidence des variations longitudinale, en particulier dans la SEB. Ces résultats sont en bon accord avec ceux de JIRAM et de MWR (un autre instrument de JUNO), avec toutefois pour ce dernier une différence importante au niveau la NEZ. / The goal of this thesis is the study of the Jupiter's troposphic chemical species spatial repartition, with an emphasis on NH3. Two types of data were used: ground-based observations with the TEXES instrument mounted on the IRTF, and space-based observations with the JIRAM instrument mounted on the JUNO spacecraft. The data were analysed by a parallelised code developped during this thesis, combining a line-by-line radiative transfer and an optimal non-linear retrieval method.The TEXES data were used to draw a map of the tropospheric NH3 (between 1 and 4 bar) with an uncertainty of 20% on three quarters of Jupiter's surface with a spatial reslution of approximately 1200 km (1.4"). This map highlight important latitudinal variations of the tropospheric NH3, which can vary from 60 ppmv near the hotspots of the NEB, to more than 500 ppmv into the zones. However, the results in the zone have a low signal-to-noise ratio, increasing their uncertainties. The map also highlight longitudinal variations, particularly in the SEB. These results are in good agreement with those of JIRAM and MWR (another instrument of JUNO), albeit an important discrepancy with MWR's results in the NEZ.
13

Restitution des profils de dégagement de chaleur latente par radiométrie hyperfréquence. Application aux cyclones tropicaux

Burlaud, Corinne 22 December 2003 (has links) (PDF)
Les cyclones tropicaux, gigantesques machines thermiques redistribuent l'énergie entre les régions chaude et froide du globe. Cette énergie (chaleur latente) provient surtout des processus de condensation. La connaissance de la chaleur latente devrait améliorer la prévision de l'évolution de ses systèmes dévastateurs détectables par satellite. TRMM lancé en 1997, permet la restitution des contenus en hydrométéores grâce à des mesures couplées radar/radiomètre (PR/TMI). Cette thèse porte sur la restitution de la chaleur latente. On a donc développé un algorithme original de type bayésien utilisant des bases de données, construites à partir des mesures du PR. Cette méthode, appliquée au cyclone Bret (1999) a été validée par comparaison avec des résultats de mesures aéroportées et sol. On reproduit assez bien sur ce cas une structure énergétique propre à différents stade de son évolution. Une étude statistique pourrait ainsi permettre un diagnostic de l'évolution.
14

Évolution de la structure en couches dans les noyaux de masse moyenne : Recherche de l'orbitale 2d5/2 neutron dans le 69Ni

Moukaddam, Mohamad 08 March 2012 (has links) (PDF)
La fermeture de couche de l'oscillateur harmonique à N=40 dans le 68Ni est faible et perd sa rigidité après l'enlèvement (ou l'ajout) de paires de protons. Les calculs effectués dans cette région de masse prédisent un nouvel îlot d'inversion à N=40 semblable à celui à N=20 et montrant que le placement de l'orbital neutron 2d5/2 est un ingrédient essentiel pour l'interprétation de la structure nucléaire à N≈40. La différence d'énergie 1g9/2-2d5/2 a été déterminée dans le noyau 69Ni en utilisant la réaction de transfert d'un neutron d(68Ni,p) en cinématique inverse. L'expérience réalisée au GANIL utilisait un faisceau de 68Ni à 25,14 MeV/u. Les noyaux 68Ni séparés par le spectromètre LISE3 ont ensuite interagit avec une cible de CD2 d'épaisseur 2,6 mg/cm2. Le dispositif expérimental était composé principalement des détecteurs CATS/MUST2-S1/EXOGAM couplés à une chambre d'ionisation et un scintillateur plastique. Les moments angulaires et les facteurs spectroscopiques de l'état fondamental (Jπ = 9/2+) et d'un doublet d'états (Jπ = 5/2+) autour de 2,48 MeV, associés à la population des orbitales 1g9/2 et 2d5/2, ont été obtenus après la comparaison des sections efficaces différentielles et des calculs ADWA. Les spins des états observés ont été attribués après comparaison aux calculs de modèles en couches dans un grand espace de valence. La position de l'orbitale 2d5/2 dans 69Ni a été établie pour la première fois. Nos mesures confirment l'hypothèse de la faible différence d'énergie (≈2,5 MeV) entre l'orbitale neutron 2d5/2 et l'orbitale 1g9/2 et son importance pour décrire la structure des noyaux autour de N=40.
15

Quantification des paléocontraintes par l'analyse des macles de la calcite : nouvelle approche d'acquisition et d'inversion des données et mécaniques du maclage / Paleostresses quantification by calcite twin analysis : new data acquisition and inversion approach and twinning mechanics

Parlangeau, Camille 04 December 2017 (has links)
La compréhension et la simulation des mécanismes de déformation dans la croûte supérieure sont des enjeux scientifiques et techniques importants. La calcite étant un minéral commun de la croûte supérieure se déformant essentiellement par maclage sous 200 °C, c'est pour cette raison que l'on s'intéresse plus particulièrement au maclage de la calcite dans le cadre de cette thèse. L'utilisation du maclage de la calcite n'en est pas à son coup d'essai et de nombreuses méthodes d'inversions existent que ce soit pour remonter aux tenseurs des contraintes comme aux tenseurs des déformations. Cette thèse propose une nouvelle méthode d'inversion basée sur l'inversion d'Etchopar permettant de reconstruire 5 sur 6 paramètres du tenseur des contraintes avec une quantification fine des incertitudes. Cette méthode permet de détecter automatiquement l'existence d'un ou de plusieurs évènements tectoniques enregistrés par le maclage de la calcite. Un deuxième volet de la thèse consiste en l'amélioration de la méthode d'acquisition des données de macles en utilisant un EBSD (electron backscatter diffraction). En effet, l'utilisation traditionnelle de la platine universelle a des limitations techniques et amène aussi à un doute visuel sur l'appréciation du statut non-maclé de certains plans, en plus du côté long et fastidieux. Le dernier volet de la thèse consiste à déterminer le seuil de maclage de la calcite pour différentes tailles de grains à l'aide de tests mécaniques sous une presse uniaxiale. Qui plus est, le suivi en continu des déformations sur les monocristaux de calcite a permis de mettre en évidence le comportement macroscopique d'un cristal de calcite et la séquence de maclage. / The understanding and modelling of deformation mechanics in the upper crust are important scientific and technical issues. The calcite is a common mineral in the upper crust and mainly deforms by twinning under 200°C. That is why we are interested by calcite twinning as part of this thesis. It is not the first time that calcite twinning is used to reconstruct paleostress or paleostrain tensors. This thesis propose a new inversion method based on the Etchecopar’s one allowing to reconstruct 5 among 6 parameters of the stress tensor with an accurate quantification of the uncertainties. This method allows to automatically detect the realness of one or several tectonic events recorded by calcite twinning. A second part of the thesis consists in the improvement of the data acquisition by using EBSD (electron backscatter diffraction). In fact, the traditional use of the universal stage has technical limitations and brings some optical doubt about the untwinned status of few planes in addition to the long and tedious side. The last part of the thesis consists in the establishing the threshold of calcite twinning for different grain sizes by mechanical tests under a uniaxial press. Moreover, the continuous monitoring of the experiments using single crystals of calcite allowed to highlight the macroscopic behavior of a single crystal and the sequence of twinning.
16

Electromagnetic wave imaging of targets buried in a cluttered medium using an hybrid Inversion-DORT method / Imagerie d'objets enfouis en utilisant la décomposition de l'opérateur de retournement temporel.

Zhang, Ting 03 March 2014 (has links)
L'objectif de ce travail de thèse est de détecter et de caractériser des cibles tridimensionnelles dans un milieu désordonné. Ce domaine de recherche est d'intérêt pour de nombreuses applications, telles que le sondage du sous-sol, l'imagerie médicale, la détection non-destructive et l'exploration géophysique, etc. Afin de distinguer les cibles des hétérogénéités du milieu, nous proposons d'utiliser l'une des techniques de retournement temporel, à savoir la méthode DORT (Décomposition de l'Opérateur de Retournement Temporel). La méthode DORT permet de générer des ondes focalisant sélectivement sur chaque cible présente dans un environnement fortement hétérogène. Par ailleurs, la richesse de ces ondes focalisantes est combinée avec un algorithme d'inversion non-linéaire. Ceci nous permet non seulement de localiser, mais aussi de caractériser les cibles (forme et permittivité). La résolution obtenue à l'aide de cette approche est bien meilleure que celles obtenues avec la méthode DORT ou la méthode d'inversion seules, en particulier dans la direction d'illumination. Cette résolution est d'autant meilleure que les données utilisées sont vectorielles. Dans le cas spécifique d'une configuration d'objets enfouis impliquant deux semi-espaces infinis, la caractérisation s'avère problématique. Une solution est apportée en appliquant l'approche de marche récurrente en fréquences. Ces développements théoriques sont également confrontés aux données expérimentales mesurées dans le domaine optique. Une nouvelle Microscopie Tomographique par Diffraction (MTD) est mise en œuvre dans le cadre de cette thèse en tenant compte du caractère vectoriel de la lumière. Ce faisant, une résolution d'environ un quart de la longueur d'onde a été obtenue sur des échantillons en résine déposés sur un substrat de silicium. De plus, nous avons aussi appliqué avec succès la méthode DORT à la MTD afin de focaliser et caractériser de manière sélective plusieurs diffuseurs de tailles différentes.Lors de ce travail de thèse nous avons également développé des méthodes de caractérisation en régime transitoire. Les différentes méthodes d'inversion élaborées dans ce cadre ont été validées sur des données synthétiques et expérimentales dans le domaine des radio-fréquences. / The objective of this thesis work is to detect and to characterize three-dimensional targets in a disordered medium, using electromagnetic excitations. This research domain is of great interest in many applications, such as subsoil probing, medical imaging, non-destructive testing and geophysical exploration, etc. In order to extract the target information from the heterogeneities of the medium, we propose to use one of the time reversal technique, namely the DORT method (French acronym for Décomposition de l'Opérateur de Retournement Temporel). This method permits us to generate different waves that focus selectively on each target in high noisy environment. Moreover, this method is also combined with a non-linear inversion algorithm, which permits not only to localize but also to characterize the targets. The reconstruction resolution appears to be better than the ones obtained with the DORT or the inversion procedure alone, especially in the illumination direction. It is also shown that using full-polarized data is indispensable for achieving better performances rather than in scalar configuration. Moreover, in the half-space configuration, it is mandatory to use the frequency-diversity data to get an accurate reconstruction. These theoretical developments are also confronted to experimental data measured in the optical domain. A full-polarization Tomographic Diffractive Microscopy (TDM) is implemented and a resolution about one-fourth of the wavelength is thus obtained. Furthermore, the DORT method is applied in TDM to realize selective focalization and characterization. In the presence of multiple targets, selective characterization of each scatterer is achieved.This thesis work also deals with the characterization problem using transient data. Different inversion algorithms are validated using synthetic and experimental hyper-frequency data.
17

Elaboration d'hydrogels composites chitosane/charbon actif à visée cicatrisante par procédés d'inversion de phase / Elaboration of chitosan/activated carbon composites hydrogels for wound dressings applications by phase inversion processes

Venault de Bourleuf, Antoine 24 September 2010 (has links)
Des hydrogels composites de chitosane et de charbon actif ont été préparés par séparation de phase induite par (i) contact de la solution de polymère avec une solution aqueuse de non-solvant (procédé par immersion) et (ii) pénétration de vapeurs de non-solvant dans la solution de polymère (procédé VIPS). La caractérisation multi-échelle (SAXS, WAXS, MEBE, module élastique, tests de diffusion) des gels a montré que le mode de contact n'avait pas d'influence sur les propriétés fonctionnelles des gels. Les structures obtenues sont homogènes, poreuses et amorphes. Elles conviennent pour constituer la couche de pansements cicatrisants en contact avec les plaies de type escarre. Un modèle original de gélification du chitosane, couplant transferts de matière, transferts de chaleur et réactions chimiques, a été élaboré. Il a permis de déterminer les temps caractéristiques de gélification, d'approcher les profils de concentration des espèces, et de corroborer la non-influence du procédé sur la structuration des matrices finales, par l'analyse des transferts mis en oeuvre. / Chitosan/Activated carbon composite hydrogels were prepared by phase separation induced by (i) immersing the polymeric solution in a non-solvent solution (wet process) and (ii) exposing the polymeric solution to non-solvent vapors (VIPS process). The multi-scale characterization approach (SAXS, WAXS, MEBE, storage modulus, diffusion tests) did not evidence any influence of the contact mode on the functional properties of the gels. Final structures are homogeneous, porous and amorphous. They are suitable to constitute the wound dressings layer in contact with eschar-like sores. An innovative chitosan gelation model was developed coupling mass transfers, heat transfers and chemical reactions. It permitted to forecast gelation times and concentration profiles. Simulation also corroborated the non-influence of the process on the final structuring of the matrices, analyzing transfers throughout the elaboration.
18

Obtenir une représentation en continu de la lithologie et de la minéralogie. Exemples d'application du traitement statistique de données de diagraphie aux structures sédimentaires en régime de convergence de plaques (Leg ODP 134, 156 et 160)

Rabaute, Alain 26 June 1998 (has links) (PDF)
Les techniques de diagraphies sont un moyen privilégié pour étudier les formations géologiques profondes, inaccessibles à l'observation directe. Elles mesurent en continu et in situ de nombreux paramètres pétrophysiques, ainsi que les concentrations en 7 éléments majeurs (Si, Ca, Fe, Al, S, K, Ti) et 3 éléments traces (U, Th et Gd). Une connaissance préalable des différents modes de fonctionnement des outils et de leur calibration permet d'évaluer la précision de la mesure. Celle-ci est surtout fonction des conditions de mesure (diamètre et rugosité du puits, vitesse de mesure). Lorsque les données sont de bonne qualité, les méthodes de classification multivariée, comme la méthode des nuées dynamiques ou le k-means, sont employées comme premier outil de visualisation des grandeurs mesurées par diagraphie. En parallèle ou en complément, à l'aide de méthodes d'inversion linéaire ou non-linéaire, il est possible de calculer une minéralogie en continu avec une précision suffisante pour être utilisée dans des modèles pétrophysiques et géologiques. Les zones de convergence de plaques lithosphériques sont des environnements géologiques instables, dans lesquels les techniques de diagraphie par cable, classiquement utilisée dans les forages scientifiques océaniques ou continentaux, ne permettent pas d'obtenir des données de bonne qualité. La nouvelle méthode de Logging-While-Drilling, ne laisse pas à l'environnement de mesure le temps de se déteriorer et donne une mesure utilisable. L'interprétation géologique, minéralogique et faciologique des données de diagraphies apportent des indications précieuses sur la dynamique sédimentaire, ainsi que sur les éventuelles transformations minéralogiques. Le caractère in situ de la mesure rend possible la modélisation de la dynamique des fluides ou de la variation et de l'orientation des champs de contrainte, dont les environnements étudiés sont souvent le siège.
19

Évolution de la structure en couches dans les noyaux de masse moyenne : recherche de l'orbitale 2d5/2 neutron dans le noyau 69Ni

Moukaddam, Mohamad 08 March 2012 (has links) (PDF)
La fermeture de couche de l'oscillateur harmonique à N=40 dans le 68Ni est faible et perd sa rigidité après l'enlèvement (ou l'ajout) de paires de protons. Les calculs effectués dans cette région de masse prédisent un nouvel îlot d'inversion à N=40 semblable à celui à N=20 et montrant que le placement de l'orbital neutron 2d5/2 est un ingrédient essentiel pour l'interprétation de la structure nucléaire à N≈40. La différence d'énergie 1g9/2-2d5/2 a été déterminée dans le noyau 69Ni en utilisant la réaction de transfert d'un neutron d(68Ni,p) en cinématique inverse. L'expérience réalisée au GANIL utilisait un faisceau de 68Ni à 25,14 MeV/u. Les noyaux 68Ni séparés par le spectromètre LISE3 ont ensuite interagit avec une cible de CD2 d'épaisseur 2,6 mg/cm2. Le dispositif expérimental était composé principalement des détecteurs CATS/MUST2-S1/EXOGAM couplés à une chambre d'ionisation et un scintillateur plastique. Les moments angulaires et les facteurs spectroscopiques de l'état fondamental (Jπ = 9/2+) et d'un doublet d'états ( Jπ = 5/2+) autour de 2,48 MeV, associés à la population des orbitales 1g9/2 et 2d5/2, ont été obtenus après la comparaison des sections efficaces différentielles et des calculs ADWA. Les spins des états observés ont été attribués après comparaison aux calculs de modèles en couches dans un grand espace de valence. La position de l'orbitale 2d5/2 dans 69Ni a été établie pour la première fois. Nos mesures confirment l'hypothèse de la faible différence d'énergie (≈2,5 MeV) entre l'orbitale neutron 2d5/2 et l'orbitale 1g9/2 et son importance pour décrire la structure des noyaux autour de N=40.
20

Etude mathématique et numérique du problème inverse de l'électro-sismique en milieu poreux / Mathematical and numerical study of the inverse problem of electro-seismicity in porous media

Xue, Qi 20 December 2017 (has links)
Dans cette thèse, nous étudions le problème inverse du phénomène de couplage des ondes électromagnétiques (EM) et sismiques. Les équations différentielles partielles régissant le phénomène de couplage sont composées d'équations de Maxwell et de Biot. Comme le phénomène de couplage est plutôt faible, nous ne considérons que la transformation des ondes électromagnétiques en ondes sismiques. Nous utilisons le modèle électrosismique pour se référer à cette transformation. Dans le modèle, le champ électrique devient la source des équations de Biot. Un coefficient de couplage est utilisé pour désigner l'efficacité de la transformation.Chapitre 2, nous considérons l'existence et l'unicité du problème vers l'avant dans le domaine fréquentiel et dans le domaine temporel. Dans le domaine fréquentiel, nous proposons l'espace de Sobolev approprié pour considérer le problème électrocinétique. Nous prouvons que la formule faible satisfait l'inégalité de Garding en utilisant la décomposition de Helmohltz. L'alternative de Fredholm peut être appliquée, ce qui montre que l'existence est équivalente à l'unicité. Dans le domaine temporel, la solution faible est définie et l'existence et l'unicité de la solution faible est démontrée.La stabilité du problème inverse est considérée dans le chapitre 3. Nous prouvons d'abord les estimations de Carleman pour les équations de Biot et les équations électrosismiques. Basé sur les estimations de Carleman pour les équations électrosismiques, nous prouvons une stabilité de Holder pour inverser tous les paramètres dans l'équation de Maxwell et le coefficient de couplage. Pour simplifier le problème, nous utilisons des équations électrostatiques pour remplacer les équations de Maxwell. Le problème inverse est décomposé en deux étapes: le problème de source inverse pour les équations de Biot et le problème de paramètre inverse pour l'équation électrostatique. Nous pouvons prouver la stabilité du problème de source inverse pour les équations de Biot sur la base de l'estimation de Carleman pour les équations de Biot. Ensuite, la conductivité et le coefficient de couplage peuvent être reconstitués avec les informations de la première étape.Dans le chapitre 4, nous résolvons les équations électrosismiques numériquement. L'équation électrostatique est résolue par la boîte à outils Matlabe PDE. Les équations de Biot sont résolues avec une méthode de différences finies échelonnées. Pour diminuer la consommation de calcul, nous ne traitons que du problème bidimensionnel. Pour simuler des ondes se propageant dans un domaine non borné, nous utilisons le PML pour absorber les ondes atteignant la limite de coupure.Le chapitre 5 traite du problème de source inverse numérique pour les équations de Biot. La méthode que nous allons utiliser est une variante de la méthode d'inversion temporelle. La première étape de la méthode consiste à transformer le problème source en un problème de valeur initiale sans aucune source. Ensuite, l'application de la méthode d'inversion de temps récupère la valeur initiale. Des exemples numériques démontrent que cette méthode fonctionne bien même pour les équations de Biot avec un petit terme d'amortissement. Mais si le terme d'amortissement est trop grand, le processus inverse n'est pas symétrique avec le processus en avant et les résultats de la reconstruction dégénèrent. / In this thesis, we study the inverse problem of the coupling phenomenon of electromagnetic (EM) and seismic waves. Partial differential equations governing the coupling phenomenon are composed of Maxwell and Biot equations. Since the coupling phenomenon is rather weak, in low frequency we only consider the transformation from EM waves to seismic waves. We use electroseismic model to refer to this transformation. In the model, the electric field becomes the source of Biot equations. A coupling coefficient is used to denote the efficiency of the transformation.Chapter 2, we consider the existence and uniqueness of the forward problem in both frequency domain and time domain. In the frequency domain, we propose the suitable Sobolev space to consider the electrokinetic problem. We prove that the weak formula satisfies a Garding's inequality using Helmohltz decomposition. The Fredholm alternative can be applied, which shows that the existence is equivalent to the uniqueness. In the time domain, the weak solution is defined and the existence and uniqueness of the weak solution is proved.The stability of the inverse problem is considered in Chapter 3. We first prove Carleman estimates for both Biot equations and electroseismic equations. Based on the Carleman estimates for electroseismic equations, we prove a Holder stability to inverse all the parameters in Maxwell equation and the coupling coefficient. To simply the problem, we use electrostatic equations to replace Maxwell equations. The inverse problem is decomposed into two steps: the inverse source problem for Biot equations and the inverse parameter problem for the electrostatic equation. We can prove the stability of the inverse source problem for Biot equations based on the Carleman estimate for Biot equations. Then the conductivity and the coupling coefficient can be reconstructed with the information from the first step.In Chapter 4, we solve the electroseismic equations numerically. The electrostatic equation is solved by the Matlabe PDE toolbox. Biot equations are solved with a staggered finite difference method. To decrease the computation consumption, we only deal with the two dimensional problem. To simulate waves propagating in unbounded domain, we use PML to absorb waves reaching the cut-off boundary.Chapter 5 deals with the numerical inverse source problem for Biot equations. The method we are going to use is a variant of the time reversal method. The first step of the method is to transform the source problem into an initial value problem without any source. Then the application of the time reversal method recovers the initial value. Numerical examples demonstrate that this method works well even for Biot equations with a small damping term. But if the damping term is too large, the inverse process is not symmetric with the forward process and the reconstruction results degenerate.

Page generated in 0.0925 seconds