Cette thèse est consacrée à l'étude mathématique d'un problème inverse en hydrogéologie : le but est d'identifier des fractures en milieu poreux, connaissant des mesures de l'écoulement dans le sous-sol. Le nombre, la localisation et les paramètres physiques des fractures sont recherchés. Ce problème est formulé comme la minimisation au sens des moindres carrés d'une fonctionnelle évaluant l'écart entre les mesures et les résultats du modèle direct. L'écoulement est celui d'un fluide monophasique incompressible (loi de Darcy). Un modèle traitant les fractures comme des interfaces est utilisé. Le problème direct est le modèle de fracture discrétisé par la méthode des éléments finis mixtes hybrides.Pour résoudre ce problème inverse, un nouvel algorithme itératif a été développé, basé sur l’utilisation d’indicateurs de fractures mis au point pendant la thèse. Ces indicateurs donnent une information au premier ordre concernant l'effet de l'ajout d'une nouvelle fracture. Comme ces indicateurs sont peu coûteux, un grand nombre de configurations de fractures sont testées à chaque itération. L’algorithme a été programmé, validé puis testé numériquement dans des situations variées, en utilisant des mesures synthétiques. Il donne des résultats très satisfaisants, bien que ce problème soit réputé difficile.Enfin, l’étude de l’identifiabilité du problème inverse a été amorcée. Pour un modèle simplifié de fractures (failles très perméables, cas le plus courant dans le sous-sol), on a montré que le problème. / This PhD is dedicated to the mathematical study of an inverse problem in hydrogeology: the goal is to identify fractures in porous medium, knowing measurements of the underground flow. The number, the location and the physical parameters of the fracture are looked for. This problem is formulated as the least squares minimization of a function evaluating the misfit between measurements and the result of the direct model. We used a model describing the flow of a monophasic incompressible fluid (Darcy's law), in a porous medium containing some fractures represented by interfaces. The direct problem is the fracture model discretized by the mixed hybrid finite element method. To solve this inverse problem, we developed an iterative algorithm, which is based on the use of fracture indicators that have been developed durig the thesis. These indicators give a first order information concerning the effect of the addition of a new fracture. As these indicators are inexpensive, a large number of configurations of new fractures is tested at each iteration. The algorithm was programmed, validated and tested numerically in various situations, using synthetic measurements. It gives very satisfactory results, although this problem is considered difficult. Finally, an early study of identifiability of the inverse problem of fractures in porous medium has been achieved. It allowed to prove the identifiability for a simplified model (very permeable faults, which is common in the underground). The question of identifiability for the full fracture model remains open.

Afin d’étudier les classes d'isotopie rigide des courbes rationnelles nodales de degré 5 dans RPP, nous associons à chaque quintique avec un point double réel marque une courbe trigonale dans la surface de Hirzebruch Sigma3 et le dessin reel nodal correspondant dans CP/(z mapsto bar{z}). Les dessins sont des versions réelles, proposées par S. Orevkov dans cite{Orevkov}, des dessins d'enfants de Grothendieck. Un dessin est un graphe contenu dans une surface topologique, muni d'une certaine structure supplémentaire. Dans cette thèse, nous étudions les propriétés combinatoires et les recompositions des dessins correspondants aux courbes trigonales nodales C subset Sigma dans les surfaces réglées réelles Sigma . Les dessins uninodaux sur une surface a bord quelconque et les dessins nodaux sur le disque peuvent être décomposés en blocs correspondant aux dessins cubiques sur le disque D2 , ce qui conduit a une classification des ces dessins. La classification des dessins considérés mène à une classification à isotopie rigide des courbes rationnelles nodales de degré 5 dans RPP. / In order to study the rigid isotopy classes of nodal rational curves of degree $5$ in $\RPP$, we associate to every real rational quintic curve with a marked real nodal point a trigonal curve in the Hirzebruch surface $\Sigma_3$ and the corresponding nodal real dessin on~$\CP/(z\mapsto\bar{z})$. The dessins are real versions, proposed by S. Orevkov~\cite{Orevkov}, of Grothendieck's {\it dessins d'enfants}. The {\it dessins} are graphs embedded in a topological surface and endowed with a certain additional structure. We study the combinatorial properties and decompositions of dessins corresponding to real nodal trigonal curves~$C\subset \Sigma$ in real ruled surfaces~$\Sigma$. Uninodal dessins in any surface with non-empty boundary and nodal dessins in the disk can be decomposed in blocks corresponding to cubic dessins in the disk~$\mathbf{D}^2$, which produces a classification of these dessins. The classification of dessins under consideration leads to a rigid isotopy classification of real rational quintics in~$\RPP$.

Dans cette thèse nous étudions le problème de Cauchy en relativité générale. Motivés par la conjecture de censure cosmique faible formulée par Penrose, nous analysons le problème aux données initiales pour les équations d'Einstein dans le vide en faible régularité. Nous démontrons les deux résultats suivants. o Premièrement, nous nous intéressons aux équations de contrainte pour les données initiales et mettons en place une procédure de prolongement. Plus précisément, étant donné des données initiales pour les équations d'Einstein sur la boule unité dans R3, nous les prolongeons de manière continue en des données globales, asymptotiquement plates sur R3. Les équations de contrainte forment un système couplé d'équations non-lineaires sous-determinées géométriques. La preuve de notre procédure de prolongement repose sur un schéma iteratif où nous séparons ce système en deux problèmes de prolongement decouplés et solubles. Enfin, le résultat de prolongement pour les équations de contrainte est obtenu par un argument de point fixe. o Deuxièment, nous prouvons une version localisée du théorème de courbure L2 de Klainerman-Rodnianski-Szeftel. Nous montrons que, étant données des données initiales pour les équations d'Einstein sur une variété compacte avec bord, le temps d'existence de la solution des équations d'Einstein dans le domaine de dépendance de ces données initiales ne dépend que de normes de basse régularité des données initiales. En particulier, notre résultat est un critère localisé de continuité pour les équations d'Einstein. Notre preuve utilise un argument de localisation où, tout d'abord, nous généralisons la théorie de Cheeger-Gromov de convergence pour les variétés Riemanniennes à notre cas de régularité faible, et ensuite nous appliquons la procédure de prolongement pour les équations de contrainte mentionnée ci-dessus avec un argument de changement d’échelle. / In this thesis we study the Cauchy problem of general relativity. Motivated by the weak cosmic censorship conjecture formulated by Penrose, we analyse the initial value problem for the Einstein vacuum equations in low regularity. We prove the following two results. First, we consider the constraint equations of the initial data and demonstrate an extension procedure. More precisely, given small initial data for the Einstein equations on the unit ball in R3, we continuosly extend it to global, asymptotically flat initial data on R3. The constraint equations for the Einstein vacuum equations are a coupled system of non-linear under-determined geometric elliptic equations. The proof of our extension procedure is based on an iterative scheme where we split this system into two decoupled, solvable extension problems. The extension result for the constraint equations follows then by a fix point argument. Second, we prove a localised version of the bounded L2-curvature theorem by Klainerman-Rodnianski-Szeftel. We show that given low regularity initial data to the Einstein equations on a compact manifold with boundary, the time of existence of the solution to the Einstein equations in the domain of dependence of the initial data depends only on low regularity geometric data. In particular, this result is a localised continuation criterion for the Einstein vacuum equations. Our proof uses a localisation argument where we first generalise the known Cheeger-Gromov convergence theory for Riemannian manifolds to our low regularity setting, and then apply the above extension procedure for the constraint equations with a scaling argument.

Dans cette thèse, nous étudions l'équation de la chaleur couplée avec la loi de Darcy à travers de la viscosité non-linéaire qui dépend de la température pour les dimensions d=2,3 (Hooman et Gurgenci ou Rashad). Nous analysons ce problème en introduisant la formulation variationnelle équivalente et en la réduisant à une simple équation de diffusion-convection pour la température où la vitesse dépend implicitement de la température.Nous démontrons l'existence de la solution sans la restriction sur les données par la méthode de Galerkin et du point fixe de Brouwer. L'unicité globale est établie une fois la solution est légèrement régulière et les données se restreignent convenablement. Nous introduisons aussi une formulation variationnelle alternative équivalente. Toutes les deux formulations variationnelles sont discrétisées par quatre schémas d'éléments finis pour un domaine polygonal ou polyédrique. Nous dérivons l'existence, l'unicité conditionnée, la convergence et l'estimation d'erreur a priori optimale pour les solutions des trois schémas. Par la suite, ces schémas sont linéarisés par des algorithmes d'approximation successifs et convergentes. Nous présentons quelques expériences numériques pour un problème modèle qui confirme les résultats théoriques de convergence développées dans ce travail. L'estimation d'erreur a posteriori est établie avec deux types d'indicateurs d'erreur de linéarisation et de discrétisation. Enfin, nous montrons des résultats numériques de validation. / In this thesis, we study the heat equation coupled with Darcy's law by a nonlinear viscosity depending on the temperature in dimension d=2,3 (Hooman and Gurgenci or Rashad). We analyse this problem by setting it in an equivalent variational formulation and reducing it to an diffusion-convection equation for the temperature where the velocity depends implicitly on the temperature.Existence of a solution is derived without restriction on the data by Galerkin's method and Brouwer's Fixed Point. Global uniqueness is established when the solution is slightly smoother and the dataare suitably restricted. We also introduce an alternative equivalent variational formulation. Both variational formulations are discretized by four finite element schemes in a polygonal or polyhedral domain. We derive existence, conditional uniqueness, convergence, and optimal a priori error estimates for the solutions of the three schemes. Next, these schemes are linearized by suitable convergent successive approximation algorithms. We present some numerical experiments for a model problem that confirm the theoretical rates of convergence developed in this work. A posteriori error estimates are established with two types of errors indicators related to the linearisation and discretization. Finally, we show numerical results of validation.

Cette thèse est consacrée à l'analyse des singularités apparaissant dans des équations différentielles partielles elliptiques non linéaires découlant de la physique mathématique, de la biologie mathématique, et de la géométrie conforme. Les thèmes abordés sont le modèle de supraconductivité de Ginzburg-Landau, le problème de Lin-Ni-Takagi, le modèle de Keller-Segel de la chimiotaxie, et le problème de courbure scalaire prescrite. Le modèle de Ginzburg-Landau est une description phénoménologique de la supraconductivité. Une caractéristique essentielle des supraconducteurs de type II est la présence de vortex, qui apparaissent au-dessus d'une certaine valeur de la force du champ magnétique appliqué, appelée premier champ critique. Nous nous intéressons au régime de epsilon petit, où epsilon est l'inverse du paramètre de Ginzburg-Landau (une constante du matériau). Dans ce régime, les vortex sont au premier ordre des singularités topologiques de co-dimension 2. Nous fournissons une construction quantitative par approximation de vortex en dimension trois pour l'énergie de Ginzburg-Landau, ce qui donne une approximation des lignes de vortex ainsi qu'une borne inférieure pour l'énergie, qui est optimale au premier ordre et vérifiée au niveau epsilon. En utilisant ces outils, nous analysons ensuite le comportement des minimiseurs globaux en dessous et proche du premier champ critique. Nous montrons que, en dessous de cette valeur critique, les minimiseurs de l'énergie de Ginzburg-Landau sont des configurations sans vortex et que les minimiseurs, proche de cette valeur, ont une vorticité bornée. Le problème de Lin-Ni-Takagi apparait comme l'ombre (dans la littérature anglaise ``shadow'') du système de Gierer-Meinhardt d'équations de réaction-diffusion qui modélise la formation de motifs biologiques. Ce problème est celui de trouver des solutions positives d'une équation critique dans un domaine régulier et borné de dimension trois, avec une condition de Neumann homogène au bord. Dans cette thèse, nous construisons des solutions à ce problème présentant un comportement explosif en un point du domaine, lorsqu'un certain paramètre converge vers une valeur critique. La chimiotaxie est l'influence de substances chimiques dans un environnement sur le mouvement des organismes. Le modèle de Keller-Segel pour la chimiotaxie est un système de diffusion-advection composé de deux équations paraboliques couplées. Ici, nous nous intéressons aux états stationnaires radiaux de ce système. Nous sommes alors amenés à étudier une équation critique dans la boule unité de dimension 2, avec une condition de Neumann homogène au bord. Dans cette thèse, nous construisons plusieurs familles de solutions radiales qui explosent à l'origine de la boule, et se concentrent sur le bord et/ou sur une sphère intérieure, lorsqu' un certain paramètre converge vers zéro. Enfin, nous étudions le problème de la courbure scalaire prescrite. Étant donnée une variété Riemannienne compacte de dimension n, nous voulons trouver des métriques conformes dont la courbure scalaire soit une fonction prescrite, qui dépend d'un petit paramètre. Nous supposons que cette fonction a un point critique qui satisfait une hypothèse de platitude appropriée. Nous construisons plusieurs métriques, qui explosent lorsque le paramètre converge vers zéro, avec courbure scalaire prescrite. / This thesis is devoted to the analysis of singularities in nonlinear elliptic partial differential equations arising in mathematical physics, mathematical biology, and conformal geometry. The topics treated are the Ginzburg-Landau model of superconductivity, the Lin-Ni-Takagi problem, the Keller-Segel model of chemotaxis, and the prescribed scalar curvature problem. The Ginzburg-Landau model is a phenomenological description of superconductivity. An essential feature of type-II superconductors is the presence of vortices, which appear above a certain value of the strength of the applied magnetic field called the first critical field. We are interested in the regime of small epsilon, where epsilon is the inverse of the Ginzburg-Landau parameter (a material constant). In this regime, the vortices are at main order co-dimension 2 topological singularities. We provide a quantitative three-dimensional vortex approximation construction for the Ginzburg-Landau energy, which gives an approximation of vortex lines coupled to a lower bound for the energy, which is optimal to leading order and valid at the epsilon-level. By using these tools we then analyze the behavior of global minimizers below and near the first critical field. We show that below this critical value, minimizers of the Ginzburg-Landau energy are vortex-free configurations and that near this value, minimizers have bounded vorticity. The Lin-Ni-Takagi problem arises as the shadow of the Gierer-Meinhardt system of reaction-diffusion equations that models biological pattern formation. This problem is that of finding positive solutions of a critical equation in a bounded smooth three-dimensional domain, under zero Neumann boundary conditions. In this thesis, we construct solutions to this problem exhibiting single bubbling behavior at one point of the domain, as a certain parameter converges to a critical value. Chemotaxis is the influence of chemical substances in an environment on the movement of organisms. The Keller-Segel model for chemotaxis is an advection-diffusion system consisting of two coupled parabolic equations. Here, we are interested in radial steady states of this system. We are then led to study a critical equation in the two-dimensional unit ball, under zero Neumann boundary conditions. In this thesis, we construct several families of radial solutions which blow up at the origin of the ball and concentrate on the boundary and/or an interior sphere, as a certain parameter converges to zero. Finally, we study the prescribed scalar curvature problem. Given an n-dimensional compact Riemannian manifold, we are interested in finding bubbling metrics whose scalar curvature is a prescribed function, depending on a small parameter. We assume that this function has a critical point which satisfies a suitable flatness assumption. We construct several metrics, which blow-up as the parameter goes to zero, with prescribed scalar curvature.

Dans la classe de complexité QMA – la généralisation quantique de la classe NP – un état quantique est fourni comme preuve à un algorithme de vérification pour l’aider à résoudre un problème. Cette classe de complexité a un problème complet naturel, le problème des Hamiltoniens locaux. Inspiré par la Physique de la matière condensée, ce problème concerne l’énergie de l’état fondamental d’un système quantique. Dans le cadre de cette thèse, nous étudions quelques problèmes liés à la classe QMA et au problème des Hamiltoniens locaux. Premièrement, nous étudions la différence de puissance si au lieu d’une preuve quantique, l’algorithme de vérification quantique reçoit une preuve classique. Nous proposons un cadre intermédiaire à ces deux cas, où la preuve consiste en un état quantique “plus simple” et nous arrivons à démontrer que ces états plus simples sont suffisants pour résoudre tous les problèmes dans QMA. À partir de ce résultat, nous obtenons un nouveau problème QMA-complet et nous étudions aussi la version de notre nouvelle classe de complexité avec erreur unilatérale. Ensuite, nous proposons le premier schéma de délégation vérifiable relativiste de calcul quantique. Dans ce cadre, un client classique délègue son calcul quantique à deux serveurs quantiques intriqués. Ces serveurs peuvent communiquer entre eux en respectant l’hypothèse que l’information ne peut pas être propagé plus vite que la vitesse de la lumière. Ce protocole a été conçu à partir d’un jeu non-local pour le problème des Hamiltoniens locaux avec deux prouveurs et un tour de communication. Dans ce jeu, les prouveurs exécutent des calculs quantiques de temps polynomiaux sur des copies de l’état fondamental du Hamiltonien. Finalement, nous étudions la conjecture PCP quantique, où l’on demande si tous les problèmes dans la classe QMA acceptent un système de preuves où l’algorithme de vérification a accès à un nombre constant de qubits de la preuve quantique. Notre première contribution consiste à étendre le modèle QPCP avec une preuve auxiliaire classique. Pour attaquer le problème, nous avons proposé une version plus faible de la conjecture QPCP pour ce nouveau système de preuves. Nous avons alors montré que cette nouvelle conjecture peut également être exprimée dans le contexte des problèmes des Hamiltoniens locaux et ainsi que dans lecadre de la maximisation de la probabilité de acceptation des jeux quantiques. Notre résultat montre la première équivalence entre un jeu multi-prouveur et une conjecture QPCP. / In QMA, the quantum generalization of the complexity class NP, a quantum state is provided as a proof of a mathematical statement, and this quantum proof can be verified by a quantum algorithm. This complexity class has a very natural complete problem, the Local Hamiltonian problem. Inspired by Condensed Matters Physics, this problem concerns the groundstate energy of quantum systems. In this thesis, we study some problems related to QMA and to the Local Hamiltonian problem. First, we study the difference of power when classical or quantum proofs are provided to quantum verification algorithms. We propose an intermediate setting where the proof is a “simpler” quantum state, and we manage to prove that these simpler states are enough to solve all problems in QMA. From this result, we are able to present a new QMA-complete problem and we also study the one-sided error version of our new complexity class. Secondly, we propose the first relativistic verifiable delegation scheme for quantum computation. In this setting, a classical client delegates her quantumcomputation to two entangled servers who are allowed to communicate, but respecting the assumption that information cannot be propagated faster than speed of light. This protocol is achieved through a one-round two-prover game for the Local Hamiltonian problem where provers only need polynomial time quantum computation and access to copies of the groundstate of the Hamiltonian. Finally, we study the quantumPCP conjecture, which asks if all problems in QMA accept aproof systemwhere only a fewqubits of the proof are checked. Our result consists in proposing an extension of QPCP proof systems where the verifier is also provided an auxiliary classical proof. Based on this proof system, we propose a weaker version of QPCP conjecture. We then show that this new conjecture can be formulated as a Local Hamiltonian problem and also as a problem involving the maximum acceptance probability of multi-prover games. This is the first equivalence of a multi-prover game and some QPCP statement.

L’objectif de cette thèse est l’étude de la fonction zêta des hauteurs motivique associée à un problème de comptage de courbes sur les compactifications équivariantes d’espaces affines, résolvant au chapitre 6 l’analogue motivique de la conjecture de Manin pour celles-ci. La fonction zêta des hauteurs provenant du problème de comptage considéré est récrite convenablement à l’aide d'une formule de Poisson motivique démontrée au cinquième chapitre, qui généralise celle de Hrushovski-Kazhdan. Chaque terme est alors décomposé sous la forme d'un produit eulérien motivique, dont la définition et les propriétés sont établies au chapitre 3. La convergence de ces produits eulériens doit être comprise pour une topologie des poids que nous introduisons au quatrième chapitre et qui repose d'une part sur la théorie des modules de Hodge de Saito, et d'autre part sur une mesure motivique sur l’anneau de Grothendieck des variétés avec exponentielles, construite dans le chapitre 2 à l’aide de la notion de cycles évanescents motiviques. On en déduit ainsi une description de l'asymptotique d'une proportion positive des coefficients du polynôme de Hodge-Deligne des espaces de modules des courbes sur la compactification équivariante donnée, lorsque le degré tend vers l'infini. / The goal of this thesis is the study of the motivic height zeta function associated to the problem of counting curves on equivariant compactifications of vector groups, solving in chapter 6 the motivic analogue of Manin's conjecture for such varieties.The motivic height zeta function coming from this counting problem is rewritten in a convenient way using a Poisson summation formula proved in chapter 5, and which generalises Hrushovski and Kazhdan's motivic Poisson formula. Each term is then expressed as a motivic Euler product, the definition and properties of the latter being established in chapter 3. The convergence of these Euler products must be understood for a weight topology which we introduce in the fourth chapter and which relies both on Saito's theory of mixed Hodge modules and on a motivic measure on the Grothendieck ring of varieties with exponentials, constructed in chapter 2 using the notion of motivic vanishing cycles. We deduce from this a description of the asymptotic of a positive proportion of the coefficients of the Hodge-Deligne polynomial of the moduli spaces of curves on the given equivariant compactification, when the degree goes to infinity.

Le but de cette thèse est d'étudier et de développer des méthodes permettant de résoudre deux types de problèmes d'identification portant sur des équations elliptiques. Ces problèmes étant connus pour leur caractère fortement instable, les méthodes proposées s'accompagnent de procédures de régularisation, qui permettent d'assurer que la solution obtenue conserve un sens physique.Dans un premier temps, on étudie la résolution du problème de Cauchy (identification de conditions aux limites) par la méthode de Steklov-Poincaré. On commence par proposer quelques améliorations basées sur le solveur de Krylov utilisé, en introduisant notamment une méthode de régularisation consistant à tronquer la décomposition de Ritz de l'opérateur concerné. Par la suite, on s'intéresse à l'estimation d'incertitude en utilisant des techniques issues de l'inversion Bayésienne. Enfin, on cherche à résoudre des problèmes plus exigeants, à savoir un problème transitoire en temps, un cas non-linéaire, et on donne des éléments pour effectuer des résolutions sur des géométries ayant un très grand nombre de degrés de liberté en s'aidant de la décomposition de domaine.Pour ce qui est du problème d'identification de fissures par la méthode de l'écart à la réciprocité, on commence par proposer et tester numériquement différents moyens de stabiliser la résolution (utilisation de fonctions-tests différentes, minimisation des gradients a posteriori ou régularisation de Tikhonov). Puis on présente une autre variante de la méthode de l'écart à la réciprocité, qui est applicable aux cas pour lesquels les mesures sont incomplètes. Cette méthode, basée sur une approche de Petrov-Galerkine, est confrontée entre autres à un cas expérimental. Enfin, on s'intéresse à certaines idées permettant d'étendre la méthode de l'écart à la réciprocité à l'identification de fissures non planes. / The goal of this thesis is to study and to develop some methods in order to solve two types of identification problems in the framework of elliptical equations. As those problems are known to be particularly unstable, the proposed methods are accompanied with regularization procedures, that ensure that the obtained solutions keep a physical meaning.Firstly, we study the resolution of the Cauchy problem (boundary conditions identification) by the Steklov-Poincaré method. We start by proposing some improvements based on the used Krylov solver, especially by introducing a regularization method that consists in truncating the Ritz values decomposition of the operator in question. We study afterwards the estimation of uncertainties by the mean of techniques stemming from Bayesian inversion. Finally, we aim at solving more demanding problems, namely a time-transient problem, a non-linear case, and we give some elements to carry out resolutions on geometries that have a very high number of degrees of freedom, with help of domain decomposition.As for the problem of crack identification by the reciprocity gap method, we firstly propose and numerically test some ways to stabilize the resolution (use of different test-functions, a posteriori minimization of the gradients or Tikhonov regularization). Then we present an other variant of the reciprocity gap method, that is applicable on cases for which the measurements are incomplete. This method, based on a Petrov-Galerkin approach, is confronted, among others, with an experimental case. Finally, we investigate some ideas that allow to extend the reciprocity gap method for the identification of non-plane cracks.

Les expériences sismiques actives sont largement utilisées pour caractériser la structure de la subsurface. Les méthodes dites d’analyse de vitesse par migration ont pour but la détermination d’un macro-modèle de vitesse, lisse, et contrôlant la cinématique de propagation des ondes. Le modèle est estimé par des critères de cohérence d’image ou de focalisation d’image. Les images de réflectivité obtenues par les techniques de migration classiques sont cependant contaminées par des artefacts, altérant la qualité de la remise à jour du macro-modèle. Des résultats récents proposent de coupler l’inversion asymptotique, qui donne des images beaucoup plus propres en pratique, avec l’analyse de vitesse pour la version offset en profondeur. Cette approche cependant demande des capacités de calcul et de mémoire importantes et ne peut actuellement être étendue en 3D.Dans ce travail, je propose de développer le couplage entre l’analyse de vitesse et la migration plus conventionnelle par point de tir. La nouvelle approche permet de prendre en compte des modèles de vitesse complexes, comme par exemple en présence d’anomalies de vitesses plus lentes ou de réflectivités discontinues. C’est une alternative avantageuse en termes d’implémentation et de coût numérique par rapport à la version profondeur. Je propose aussi d’étendre l’analyse de vitesse par inversion au domaine des données pour les cas par point de tir. J’établis un lien entre les méthodes formulées dans les domaines données et images. Les méthodologies sont développées et analysées sur des données synthétiques 2D. / Active seismic experiments are widely used to characterize the structure of the subsurface. Migration Velocity Analysis techniques aim at recovering the background velocity model controlling the kinematics of wave propagation. The first step consists of obtaining the reflectivity images by migrating observed data in a given macro velocity model. The estimated model is then updated, assessing the quality of the background velocity model through the image coherency or focusing criteria. Classical migration techniques, however, do not provide a sufficiently accurate reflectivity image, leading to incorrect velocity updates. Recent investigations propose to couple the asymptotic inversion, which can remove migration artifacts in practice, to velocity analysis in the subsurface-offset domain for better robustness. This approach requires large memory and cannot be currently extended to 3D. In this thesis, I propose to transpose the strategy to the more conventional common-shot migration based velocity analysis. I analyze how the approach can deal with complex models, in particular with the presence of low velocity anomaly zones or discontinuous reflectivities. Additionally, it requires less memory than its counterpart in the subsurface-offset domain. I also propose to extend Inversion Velocity Analysis to the data-domain, leading to a more linearized inverse problem than classic waveform inversion. I establish formal links between data-fitting principle and image coherency criteria by comparing the new approach to other reflection-based waveform inversion techniques. The methodologies are developed and analyzed on 2D synthetic data sets.

Ce travail de thèse traite du développement d'une méthode de télédiagnostique de l'état de santé des suspensions des trains à grande vitesse à partir de mesures de la réponse dynamique du train en circulation par des accéléromètres embarqués. Un train en circulation est un système dynamique dont l'excitation provient des irrégularités de la géométrie de la voie ferrée. Ses éléments de suspension jouent un rôle fondamental de sécurité et de confort. La réponse dynamique du train étant dépendante des caractéristiques mécaniques des éléments de suspension, il est possible d'obtenir en inverse des informations sur l'état de ces éléments à partir de mesures accélérométriques embarquées. Connaître l'état de santé réel des suspensions permettrait d'améliorer la maintenance des trains. D’un point de vue mathématique, la méthode de télédiagnostique proposée consiste à résoudre un problème statistique inverse. Elle s'appuie sur un modèle numérique de dynamique ferroviaire et prend en compte l'incertitude de modèle ainsi que les erreurs de mesures. Les paramètres mécaniques associés aux éléments de suspension sont identifiés par calibration Bayésienne à partir de mesures simultanées des entrées (les irrégularités de la géométrie de la voie) et sorties (la réponse dynamique du train) du système. La calibration Bayésienne classique implique le calcul de la fonction de vraisemblance à partir du modèle stochastique de réponse et des données expérimentales. Le modèle numérique étant numériquement coûteux d'une part, ses entrées et sorties étant fonctionnelles d'autre part, une méthode de calibration Bayésienne originale est proposée. Elle utilise un métamodèle par processus Gaussien de la fonction de vraisemblance. Cette thèse présente comment un métamodèle aléatoire peut être utilisé pour estimer la loi de probabilité des paramètres du modèle. La méthode proposée permet la prise en compte du nouveau type d'incertitude induit par l'utilisation d'un métamodèle. Cette prise en compte est nécessaire pour une estimation correcte de la précision de la calibration. La nouvelle méthode de calibration Bayésienne a été testée sur le cas applicatif ferroviaire, et a produit des résultats concluants. La validation a été faite par expériences numériques. Par ailleurs, l'évolution à long terme des paramètres mécaniques de suspensions a été étudiée à partir de mesures réelles de la réponse dynamique du train / The work presented here deals with the development of a health-state monitoring method for high-speed train suspensions using in-service measurements of the train dynamical response by embedded acceleration sensors. A rolling train is a dynamical system excited by the track-geometry irregularities. The suspension elements play a key role for the ride safety and comfort. The train dynamical response being dependent on the suspensions mechanical characteristics, information about the suspensions state can be inferred from acceleration measurements in the train by embedded sensors. This information about the actual suspensions state would allow for providing a more efficient train maintenance. Mathematically, the proposed monitoring solution consists in solving a statistical inverse problem. It is based on a train-dynamics computational model, and takes into account the model uncertainty and the measurement errors. A Bayesian calibration approach is adopted to identify the probability distribution of the mechanical parameters of the suspension elements from joint measurements of the system input (the track-geometry irregularities) and output (the train dynamical response).Classical Bayesian calibration implies the computation of the likelihood function using the stochastic model of the system output and experimental data. To cope with the fact that each run of the computational model is numerically expensive, and because of the functional nature of the system input and output, a novel Bayesian calibration method using a Gaussian-process surrogate model of the likelihood function is proposed. This thesis presents how such a random surrogate model can be used to estimate the probability distribution of the model parameters. The proposed method allows for taking into account the new type of uncertainty induced by the use of a surrogate model, which is necessary to correctly assess the calibration accuracy. The novel Bayesian calibration method has been tested on the railway application and has achieved conclusive results. Numerical experiments were used for validation. The long-term evolution of the suspension mechanical parameters has been studied using actual measurements of the train dynamical response