Cette thèse porte sur l'analyse mathématique de modèles de réaction-dispersion de la forme [delta]tu=D(u) +f(x,u). L'objectif est de comprendre l'influence du terme de réaction f, de l'opérateur de dispersion D, et de la donnée initiale u0 sur la propagation des solutions de ces équations. Nous nous sommes intéressés principalement à deux types d'équations de réaction-dispersion : les équations de réaction-diffusion où l'opérateur de dispersion différentielle est D=[delta]2z et les équations intégro-différentielles pour lesquelles D est un opérateur de convolution, D(u)=J* u-u. Dans le cadre des équations de réaction-diffusion en milieu homogène, nous proposons une nouvelle approche plus intuitive concernant les notions de fronts progressifs tirés et poussés. Cette nouvelle caractérisation nous a permis de mieux comprendre d'une part les mécanismes de propagation des fronts et d'autre part l'influence de l'effet Allee, correspondant à une diminution de la fertilité à faible densité, lors d'une colonisation. Ces résultats ont des conséquences importantes en génétique des populations. Dans le cadre des équations de réaction-diffusion en milieu hétérogène, nous avons montré sur un exemple précis comment la fragmentation du milieu modifie la vitesse de propagation des solutions. Enfin, dans le cadre des équations intégro-différentielles, nous avons montré que la nature sur- ou sous-exponentielle du noyau de dispersion J modifie totalement la vitesse de propagation. / This thesis deals with the mathematical analysis of reaction-dispersion models of the form [delta]tu=D(u) +f(x,u). We investigate the influence of the reaction term f, the dispersal operator D and the initial datum u0 on the propagation of the solutions of these reaction-dispersion equations. We mainly focus on two types of equations: reaction-diffusion equations (D=[delta]2z and integro-differential equations (D is a convolution operator, D(u)=J* u-u). We first investigate the homogeneous reaction-diffusion equations. We provide a new and intuitive explanation of the notions of pushed and pulled traveling waves. This approach allows us to understand the inside dynamics the traveling fronts and the impact of the Allee effect, that is a low fertility at low density, during a colonisation. Our results also have important consequences in population genetics. In the more general and realistic framework of heterogeneous reaction-diffusion equations, we exhibit examples where the fragmentation of the media modifies the spreading speed of the solution. Finally, we investigate integro-differential equations and prove that emph{fat-tailed} dispersal kernels J, that is kernels which decay slower than any exponentially decaying function at infinity, lead to acceleration of the level sets of the solution u.

Durant la dernière décade, les radars millimétriques en bande W (75 - 110 GHz) pour les applications civiles que ce soit dans le domaine de l'aide à la conduite ou de la sécurité. La maturité de ces systèmes et les exigences accrues en termes d'application, orientent actuellement les recherches vers l'insertion de fonctions permettant l'identification. Ainsi, des systèmes d'imagerie radar ont été développés, notamment à l'aide d'imagerie qualitative (SAR). Les premiers résultats sont très prometteurs, cependant, afin de reconstruire les propriétés électromagnétiques des objets, il faut travailler de manière quantitative. De nombreux travaux ont déjà été conduits en ondes centimétriques, mais aucun système d'imagerie quantitative n'existe, à notre connaissance, en gamme millimétrique. L'objectif du travail présenté dans ce manuscrit est de poser les bases d'un système d'imagerie quantitative en gamme millimétrique et de le comparer à l'imagerie radar de systèmes développés en collaboration avec l'Université d'Ulm (Allemagne). L'ensemble des résultats obtenus valide le processus développé pour d'imagerie quantitative. Les recherches doivent être poursuivies. D'une part le système de mesure doit évoluer vers un vrai système multi-incidences/multivues. D'autre part, le cas 2D-TE doit être implémenté afin de pouvoir traiter un objet 2D quelconque dans n'importe quelle polarisation. Enfin, les mesures à partir de systèmes radar réels doivent être poursuivies, en particulier pour rendre exploitables les mesures des coefficients de transmission. Ces dernières sont indispensables si l'on veut un jour appliquer les algorithmes d'inversion à des mesures issues de systèmes radar. / In the last decade, microwave and millimeter-wave systems have gained importance in civil and security applications. Due to an increasing maturity and availability of circuits and components, these systems are getting more compact while being less expensive. Furthermore, quantitative imaging has been conducted at lower frequencies using computational intensive inverse problem algorithms. Due to the ill-posed character of the inverse problem, these algorithms are, in general, very sensitive to noise: the key to their successful application to experimental data is the precision of the measurement system. Only a few research teams investigate systems for imaging in the W-band. In this manuscript such a system is presented, designed to provide scattered field data to quantitative reconstruction algorithms. This manuscript is divided into six chapters. Chapter 2 describes the theory to compute numerically the scattered fields of known objects. In Chapter 3, the W-band measurement setup in the anechoic chamber is shown. Preliminary measurement results are analyzed. Relying on the measurement results, the error sources are studied and corrected by post-processing. The final results are used for the qualitative reconstruction of all three targets of interest and to image quantitatively the small cylinder. The reconstructed images are compared in detail in Chapter 4. Close range imaging has been investigated using a vector analyzer and a radar system. This is described in Chapter 5, based on a future application, which is the detection of FOD on airport runways. The conclusion is addressed in Chapter 6 and some future investigations are discussed.

Le sujet de l'étude concerne la caractérisation thermique à hautes températures de matériaux anisotropes dont la diffusivité thermique varie selon la direction considérée. Cette mesure de la diffusivité est permise par l'observation des variations transitoires de température d'un matériau soumis à un flux de chaleur de type impulsionnel. L’excitation provient d’un Laser et la mesure de température est réalisée par thermographie infrarouge sur la face opposée à l'excitation thermique. Le champ de température ainsi obtenu permet de déterminer les trois diffusivités du matériau selon ses directions d'anisotropie. En effet, grâce à des transformations intégrales du champ de température, il est possible d'obtenir un modèle théorique décrivant les variations de température au sein du matériau. Les estimations des diffusivités s'obtiennent alors par la minimisation de la somme des écarts quadratiques entre les modèles théoriques et leurs équivalents expérimentaux. Il s'agit de problèmes d'optimisation non linéaire et les estimations sont réalisées dans le domaine des fréquences spatiales et dans le temps grâce à une inversion numérique de Laplace. Basée sur des dispositifs optiques, cette méthode est non intrusive et grâce aux modèles analytiques les mesures sont rapides et précises même à haute température. La méthode ainsi que le nouveau banc expérimental mis en place rendent possible la mesure des trois diffusivités en une unique expérience pour des excitations de forme quelconque en espace et non nécessairement Dirac en temps / The study concerns the thermal characterization at high temperatures of anisotropic materials whose thermal diffusivity varies according to the direction considered. This measurement of diffusivity is allowed by the observation of the transient variations of temperature of a material subjected to a heat pulse source. The excitation is performed by a Laser and the temperature measurement is carried out by infrared thermography on the opposite face of the thermal excitation. The temperature field thus obtained makes it possible to determine the three diffusivities of the material according to its directions of anisotropy. Indeed, thanks to integral transforms of the temperature field, it is possible to obtain a theoretical model describing the temperature variations within the material. The estimates of diffusivities are then obtained by the minimization of the sum of squared residuals between the theoretical models and their experimental equivalents. These are problems of nonlinear optimization and the estimations are carried out in the spatial frequency domain and in time thanks to a numerical inversion of Laplace. Based on optical devices this method is non-intrusive and thanks to the use of analytical models the estimations are fast and accurate even at high temperatures. The method and the new experimental facility make it possible to estimate the three thermal diffusivities in a single experiment and this for excitations of any shape in space and not necessarily Dirac’s delta function in time

Actuellement, le principal obstacle à la mise en œuvre de la FWI élastique en trois dimensions sur des cas d'étude réalistes réside dans le coût de calcul associé aux taches de modélisation sismique. Pour surmonter cette difficulté, je propose deux contributions. Tout d'abord, je propose de calculer le gradient de la fonctionnelle avec la méthode de l'état adjoint à partir d'une forme symétrisée des équations de l'élastodynamique formulées sous forme d'un système du premier ordre en vitesse-contrainte. Cette formulation auto-adjointe des équations de l'élastodynamique permet de calculer les champs incidents et adjoints intervenant dans l'expression du gradient avec un seul opérateur de modélisation numérique. Le gradient ainsi calculé facilite également l'interfaçage de plusieurs outils de modélisation avec l'algorithme d'inversion. Deuxièmement, j'explore dans cette thèse dans quelle mesure les encodages des sources avec des algorithmes d'optimisation du second-ordre de quasi-Newton et de Newton tronqué permettait de réduire encore le coût de la FWI. Finalement, le problème d'optimisation associé à la FWI est mal posé, nécessitant ainsi d'ajouter des contraintes de régularisation à la fonctionnelle à minimiser. Je montre ici comment une régularisation fondée sur la variation totale du modèle fournissait une représentation adéquate des modèles du sous-sol en préservant le caractère discontinu des interfaces lithologiques. Pour améliorer les images du sous-sol, je propose un algorithme de débruitage fondé sur une variation totale locale au sein duquel j'incorpore l'information structurale fournie par une image migrée pour préserver les structures de faible dimension. / Currently, the main limitation to perform 3D elastic full waveform inversion on a production level is the computational cost it represents. With this in mind, we provide two contributions. First, we develop a self adjoint formulation of the isotropic first order velocity-stress elastic equations that allow to implement only one forward modeling operator in the gradient computation. Second, we combine Newton and quasi-Newton optimization methods with source encoding techniques to see to what extent the computational cost could be further reduced. Finally, the optimization process associated to FWI is ill posed and requires regularization constraints. I show that the total variation of the model as a regularization term provides and adequate description of earth models, preserving the discontinuous character of the lithological layers. To improve the quality of the images, we propose a local total variation denoising algorithm based on the incorporation of the information provided by a migrated image.

La caractérisation expérimentale de sources acoustiques est l'une des étapes essentielles pour la réduction des nuisances sonores produites par les machines industrielles. L'objectif de la thèse est de mettre au point une procédure complète visant à localiser et à quantifier des sources acoustiques stationnaires ou non sur un maillage surfacique par la rétro-propagation d'un champ de pression mesuré par un réseau de microphones. Ce problème inverse est délicat à résoudre puisqu'il est généralement mal-conditionné et sujet à de nombreuses sources d'erreurs. Dans ce contexte, il est capital de s'appuyer sur une description réaliste du modèle de propagation acoustique direct. Dans le domaine fréquentiel, la méthode des sources équivalentes a été adaptée au problème de l'imagerie acoustique dans le but d'estimer les fonctions de transfert entre les sources et l'antenne, en prenant en compte le phénomène de diffraction des ondes autour de l'objet d'intérêt. Dans le domaine temporel, la propagation est modélisée comme un produit de convolution entre la source et une réponse impulsionnelle décrite dans le domaine temps-nombre d'onde. Le caractère sous-déterminé du problème acoustique inverse implique d'utiliser toutes les connaissances a priori disponibles sur le champ sources. Il a donc semblé pertinent d'employer une approche bayésienne pour résoudre ce problème. Des informations a priori disponibles sur les sources acoustiques ont été mises en équation et il a été montré que la prise en compte de leur parcimonie spatiale ou de leur rayonnement omnidirectionnel pouvait améliorer significativement les résultats. Dans les hypothèses formulées, la solution du problème inverse s'écrit sous la forme régularisée de Tikhonov. Le paramètre de régularisation a été estimé par une approche bayésienne empirique. Sa supériorité par rapport aux méthodes communément utilisées dans la littérature a été démontrée au travers d'études numériques et expérimentales. En présence de fortes variabilités du rapport signal à bruit au cours du temps, il a été montré qu'il est nécessaire de mettre à jour sa valeur afin d'obtenir une solution satisfaisante. Finalement, l'introduction d'une variable manquante au problème reflétant la méconnaissance partielle du modèle de propagation a permis, sous certaines conditions, d'améliorer l'estimation de l'amplitude complexe des sources en présence d'erreurs de modèle. Les développements proposés ont permis de caractériser, in situ, la puissance acoustique rayonnée par composant d'un groupe motopropulseur automobile par la méthode de la focalisation bayésienne dans le cadre du projet Ecobex. Le champ acoustique cyclo-stationnaire généré par un ventilateur automobile a finalement été analysé par la méthode d'holographie acoustique de champ proche temps réel. / Experimental characterization of acoustic sources is one of the essential steps for reducing noise produced by industrial machinery. The aim of the thesis is to develop a complete procedure to localize and quantify both stationary and non-stationary sound sources radiating on a surface mesh by the back-propagation of a pressure field measured by a microphone array. The inverse problem is difficult to solve because it is generally ill-conditioned and subject to many sources of error. In this context, it is crucial to rely on a realistic description of the direct sound propagation model. In the frequency domain, the equivalent source method has been adapted to the acoustic imaging problem in order to estimate the transfer functions between the source and the antenna, taking into account the wave scattering. In the time domain, the propagation is modeled as a convolution product between the source and an impulse response described in the time-wavenumber domain. It seemed appropriate to use a Bayesian approach to use all the available knowledge about sources to solve this problem. A priori information available about the acoustic sources have been equated and it has been shown that taking into account their spatial sparsity or their omnidirectional radiation could significantly improve the results. In the assumptions made, the inverse problem solution is written in the regularized Tikhonov form. The regularization parameter has been estimated by an empirical Bayesian approach. Its superiority over methods commonly used in the literature has been demonstrated through numerical and experimental studies. In the presence of high variability of the signal to noise ratio over time, it has been shown that it is necessary to update its value to obtain a satisfactory solution. Finally, the introduction of a missing variable to the problem reflecting the partial ignorance of the propagation model could improve, under certain conditions, the estimation of the complex amplitude of the sources in the presence of model errors. The proposed developments have been applied to the estimation of the sound power emitted by an automotive power train using the Bayesian focusing method in the framework of the Ecobex project. The cyclo-stationary acoustic field generated by a fan motor was finally analyzed by the real-time near-field acoustic holography method.

Il ne fait aucun doute que le coût relativement faible des échographes, la procédure rapide et la capacité d’imager des tissus biologiques mous ont permis à l’échographie de devenir l’un des outils de diagnostic médical les plus courants. Malheureusement, l’imagerie ultrasonore présente encore des inconvénients par rapport à d’autres techniques d’imagerie médicale, principalement en termes de qualité d’image et de détails fournis. Dans le souci d’améliorer la qualité d’image, le prix à payer est généralement la baisse de la cadence d’acquisition. Ce compromis profondément enraciné entre la qualité d'image fournie et le temps d'acquisition est peut-être l'un des plus difficiles verrous de la recherche dans le domaine de l’échographie et son dépassement pourrait mener à des meilleures performances des applications existantes et même à des nouvelles utilisations de l’échographie. Ce travail de thèse porte sur l’amélioration du compromis précédemment indiqué entre la qualité d’image et la cadence d’acquisition. En combinant des concepts tels que l’imagerie par ondes planes, les systèmes à entrées multiples / sorties multiples et les problèmes inverses, ce travail vise à acquérir simultanément des images ultrasonores du tissu insonifié, ce qui permet d’augmenter la cadence d’acquisition sans réduire la qualité de l’image fournie. À travers cette étude, nous avons mis au point un modèle mathématique permettant de modéliser la propagation des ondes ultrasonores dans les tissus mous. Ce modèle a été utilisé pour passer en revue un grand nombre de schémas d’acquisition en échographie déjà existants et pour exposer leurs avantages et leurs inconvénients. Nous avons proposé de surmonter le compromis qualité / cadence d’acquisition d'images en utilisant des ondes ultrasonores codées temporellement émises simultanément, et le modèle direct généré a permis l'utilisation des différentes approches de problèmes inverses afin de reconstruire la réponse impulsionnelle pulse / écho du milieu insonifié et donc de son image. De plus, nous avons encore amélioré le modèle direct, ce qui nous a permis de relier directement les échos rétrodiffusés à la position / magnitude des diffuseurs à l'intérieur du milieu imagé. Les résultats fournies par les approches de problèmes inverses basées sur ce modèle amélioré nous ont mis devant une méthode de pointe qui non seulement améliorent d’un facteur conséquent la qualité de l'image en termes de résolution et de cohérence du speckle, mais permettent également d'améliorer la cadence d'acquisition d’images. / It is beyond doubt that the relative low cost of ultrasound scanners, the quick procedure and the ability to image soft biological tissues helped ultrasound imaging to become one of the most common medical diagnostic tools. Unfortunately, ultrasound still has some drawbacks when compared to other medical imaging techniques mainly in terms of the provided image quality and details. In the quest for an improved image quality, usually, the price to pay is the drop in the frame acquisition rate. This deep rooted trade-off between the provided image quality and the acquisition time is perhaps one of the most challenging in today’s ultrasound research and its overcoming could lead to diagnostic improvements in the already existing ultrasound applications and even pave the way towards novel uses of echography. This study addresses the previously stated trade-off. Through a mix of such concepts as plane wave imaging, multiple-input /multiple-output systems and inverse problems, this work aims at acquiring ultrasound images of the insonified tissue simultaneously, thus providing an increased frame rate while not degrading the image quality. Through this study we came up with a mathematical model that allows modelling the ultrasound wave propagation inside soft tissues. This model was used to review a great number of existing ultrasound acquisition schemes and to expose their advantages and drawbacks. We proposed to overcome the image quality / frame rate trade-off by using temporally encoded ultrasound waves emitted simultaneously, and the generated direct model enabled the use of different inverse problem approaches in order to reconstruct the pulse-echo impulse response of the insonified medium and thus its image. Moreover, we further improved the direct model, which allowed us to directly link the backscattered echoes to the position / magnitude of the scatterers inside the imaged medium. The results yielded by the inverse problem approaches based on the former model put us face to face with state of the art method that not only increase the image quality several times in terms resolution and speckle coherence but also provide a boost in frame acquisition rate.

Die vorliegende Arbeit hat nichtparametrische Schätzmethoden für diskret beobachtete Lévyprozesse zum Gegenstand. Ein Lévyprozess mit endlichen zweiten Momenten und endlicher Variation auf Kompakta wird niederfrequent beobachtet. Die Sprungdynamik wird vollständig durch das endliche signierte Maß my(dx):= x ny(dx) beschrieben. Ein lineares Funktional von my soll nichtparametrisch geschätzt werden. Im ersten Teil werden Kernschätzer konstruiert und obere Schranken für das korrespondierende Risiko bewiesen. Daraus werden Konvergenzraten unter Glattheitsannahmen an das Lévymaß hergeleitet. Für Spezialfälle werden untere Schranken bewiesen und daraus Minimax-Optimalität gefolgert. Der Schwerpunkt liegt auf dem Problem der datengetriebenen Wahl des Glättungsparameters, das im zweiten Teil untersucht wird. Da die nichtparametrische Schätzung für Lévyprozesse starke strukturelle Ähnlichkeiten mit Dichtedekonvolutionsproblemen mit unbekannter Fehlerdichte aufweist, werden beide Problemstellungen parallel diskutiert und die Methoden allgemein sowohl für Lévyprozesse als auch für Dichtedekonvolution entwickelt. Es werden Methoden der Modellwahl durch Penalisierung angewandt. Während das Prinzip der Modellwahl im üblichen Fall darauf beruht, dass die Fluktuation stochastischer Terme durch Penalisierung mit einer deterministischen Größe beschränkt werden kann, ist die Varianz im hier betrachteten Fall unbekannt und der Strafterm somit stochastisch. Das Hauptaugenmerk der Arbeit liegt darauf, Strategien zum Umgang mit dem stochastischen Strafterm zu entwickeln. Dabei ist ein modifizierter Schätzer für die charakteristische Funktion im Nenner zentral, der es erlaubt, die punktweise Kontrolle der Abweichung dieses Objects von seiner Zielgröße auf die gesamte reelle Achse zu erweitern. Für die Beweistechnik sind insbesondere Talagrand-Konzentrationsungleichungen für empirische Prozesse relevant. / This thesis deals with nonparametric estimation methods for discretely observed Lévy processes. A Lévy process X having finite variation on compact sets and finite second moments is observed at low frequency. The jump dynamics is fully described by the finite signed measure my(dx)=x ny(dx). The goal is to estimate, nonparametrically, some linear functional of my. In the first part, kernel estimators are constructed and upper bounds on the corresponding risk are provided. From this, rates of convergence are derived, under regularity assumptions on the Lévy measure. For particular cases, minimax lower bounds are proved. The rates of convergence are thus shown to be minimax optimal. The focus lies on the data driven choice of the smoothing parameter, which is being considered in the second part. Since nonparametric estimation methods for Lévy processes have strong structural similarities with with nonparametric density deconvolution with unknown error density, both fields are discussed in parallel and the concepts are developed in generality, for Lévy processes as well as for density deconvolution. The choice of the bandwidth is realized, using techniques of model selection via penalization. The principle of model selection via penalization usually relies on the fact that the fluctuation of certain stochastic quantities can be controlled by penalizing with a deterministic term. Contrarily to this, the variance is unknown in the setting investigated here and the penalty term is hence itself a stochastic quantity. It is the main concern of this thesis to develop strategies to dealing with the stochastic penalty term. The most important step in this direction will be a modified estimator of the unknown characteristic function in the denominator, which allows to make the pointwise control of this object uniform on the real line. The main technical tools involved in the arguments are concentration inequalities of Talagrand type for empirical processes.

Le présent travail de thèse fait suite à une première étude où une stratégie d’identification des paramètres et formes des lois de zones cohésives a été élaborée pour des matériaux homogènes. L’extension au cas de matériaux présentant des hétérogénéités nécessite d’accéder localement aux champs de contraintes.Ainsi, l’objectif principal de cette étude est de mettre au point une méthode de caractérisation locale des propriétés mécaniques et des contraintes. Cette méthode est basée sur l’erreur en relation de comportement combinée à l’exploitation de la richesse des mesures de champs cinématiques planes et plus particulièrement des champs de déformations, obtenus par dérivation numérique des champs de déplacements. Cette mesure cinématique est réalisée par une technique de corrélation d’images numériques enrichie.La méthode d’identification est basée sur la minimisation itérative d’une norme énergétique faisant intervenir le tenseur élastoplastique sécant. Différentes simulations numériques ont illustré la capacité de la procédure à identifier localement des champs de propriétés hétérogènes et sa robustesse et sa stabilité vis-à-vis du bruit de mesure, du choix du jeu de paramètres d’initialisation de l’algorithme et de la finesse du maillage.Pour finir, des essais plans avec différentes géométries d’éprouvettes ont été effectués et un essai a été mis au point pour obtenir de manière maîtrisée un état initial très hétérogène. Les résultats d’identification élastoplastique multilinéaire ont montré la capacité de la méthode à identifier les lois de comportements locales sur ce matériau hétérogène. / The present work follows a first approach where a strategy for identifying the shape and the parameters of cohesive-zone laws has been developed for homogeneous materials. The extension of this method to heterogeneous material requires the knowledge of the local stress state.The study aims at developing a local characterization method for mechanical properties and stresses. This method is based on the constitutive equation gap principles and relies on the knowledge of mechanical kinematic fields and particularly of the strain fields. These fields are obtained by the numerical differentiation of displacement fields measured by digital image correlation.This identification method is based on the iterative minimization of an energy norm involving the secant elastoplastic tensor. Various numerical simulations were used to illustrate the performance of the procedure for locally identifying heterogeneous property fields, and to characterize its robustness and its stability with respect to noise to the values of the algorithm initialization parameter and to the mesh refinement.Finally, various experimental tests with different specimen geometries were performed and a test has been developed to obtain a controlled heterogeneous initial state. The multilinear elastoplastic identification results showed the ability of the method to identify the local behavior properties on heterogeneous materials.

Dans les puits de production d’eau, de pétrole, de gaz et de chaleur géothermique, ou dans les puits d’accès à un stockage d’hydrocarbures, il est précieux de connaître la perméabilité de la formation ou de sa couverture en fonction de la profondeur, soit pour améliorer le modèle de réservoir, soit pour choisir les zones dans lesquelles procéder à des opérations spéciales.On propose une technique qui consiste à balayer la hauteur du découvert par une interface entre deux liquides de viscosités très contrastées. Le débit total qui pénètre la formation à chaque instant est ainsi une fonction de la position de l’interface et de l’historique des pressions dans le puits. On doit alors résoudre un problème inverse : rechercher la perméabilité fonction de la profondeur à partir de l’historique des débits dans le temps. Dans la pratique, le puits est équipé d’un tube central. Le balayage est effectué par injection d’un liquide à pression d’entrée constante dans le tube central et soutirage d’un autre liquide par l’espace annulaire. On mesure les débits d’injection et de soutirage dont la différence est le débit qui entre dans la formation.Pour valider et améliorer cette technique, on a d’abord utilisé une maquette simulant un découvert multi-couches disponible au LMS. On a exploité aussi des essais en place réalisés dans la couverture peu perméable d’un stockage souterrain de gaz. Dans ces essais, un liquide visqueux placé dans le découvert était déplacé par un liquide moins visqueux (méthode dite « opening »). Les couches plus perméables étaient correctement identifiées (Manivannan et al. 2017), mais une estimation quantitative était un défi en raison des phénomènes transitoires qui affectent le voisinage immédiat des puits. De plus, le rayon investigué dans le massif était petit.La thèse a relevé ces défis en proposant un essai légèrement différent et une nouvelle technique d’interprétation. Les essais avec une maquette modifiée ont montré la supériorité d’une méthode « closing » dans laquelle le puits est d’abord rempli du liquide le moins visqueux. On ménage une période de stabilisation avant l’injection du liquide visqueux pour réduire les effets transitoires ; elle permet aussi d’estimer la perméabilité moyenne et l’influence de la zone endommagée à la paroi (le « skin »).Puis on conduit l’essai proprement dit. L’historique des débits mesurés en tête de puits constitue le profil d’injection dont on déduit le profil de perméabilité.. Cette estimation suppose un écoulement monophasique dans chaque couche et la même « skin » pour toute la formation. Les incertitudes principales portent sur les pressions de formation et les variations possibles du « skin ». Elles sont estimées au moyen d’un calcul analytique. On a vérifié sur la maquette que les profils de perméabilité estimés présentent une bonne concordance avec les perméabilités mesurées avant les essais.On a réalisé un essai sur un sondage de 1750 m de long atteignant une couche de sel dont on a correctement estimé la perméabilité moyenne pendant la période de stabilisation. Toutefois elle était si faible (4.0E-21 m²) que l’utilisation de deux fluides n’a pas permis de faire une différence entre les diverses parties du puits. / In wells producing water, oil, gas or geothermal energy, or in access wells to hydrocarbon storage, it is critical to evaluate the permeability of the formation as a function of depth, to improve the reservoir model, and also to identify the zones where additional investigation or special completions are especially useful.A new technique is proposed, consisting of scanning the open hole (uncased section of the wellbore) with an interface between two fluids with a large viscosity contrast. The injection rate into the formation depends on interface location and well pressure history. An inverse problem should be solved: estimate permeability as a function of depth from the evolution of flow rates with time. The wells are usually equipped with a central tube. The scanning is done by injecting a liquid in the central tube at constant wellhead pressure. Injection and withdrawal rates are measured at the wellhead; the difference between these two rates is the formation injection rate.To validate and improve this technique, we used a laboratory model mimicking a multi-layer formation, already available at LMS. We also made use of in-situ tests performed on an ultra-low permeable cap rock above an underground gas storage reservoir. In these tests, a viscous fluid contained in the open hole was displaced by a less-viscous fluid (a method called opening WTLog). The more permeable layers were correctly identified (Manivannan et al. 2017), but a quantitative estimation was challenging due to transient phenomena in the vicinity of the wellbore (near-wellbore zone). In addition, the investigation radius was small.These challenges are addressed by proposing a slightly modified test procedure and a new interpretation workflow. Laboratory tests with a modified test setup showed the advantages of the ‘closing’ method in which the well is filled with a less-viscous fluid at the start of the test. We also added a stabilization period before the injection of viscous fluid to minimize the transient effects; this period is also used to estimate the average permeability of the open hole and the effect of near-wellbore damage (skin).Then the test proper is performed (closing WTLog). The injection profile of the less-viscous fluid is computed from the wellhead flow rate history. A permeability profile is estimated from the injection profile. The permeability estimation considers a monophasic flow in each layer and the same skin value for all the formation layers. Major uncertainties in the permeability estimates are caused by formation pressures and heterogeneities in skin values; they are estimated using an analytical formula. We have verified on the laboratory setup that the estimated permeability profiles are well correlated to the permeabilities measured before the tests.An attempt was made to perform a WTLog in a 1750-m long wellbore opening in a salt formation. The first phase was successful and the average permeability was correctly assessed. However, this permeability was so small (4.0E-21 m² or 4 nD) that the gauges and the flowmeters were not accurate enough to allow a clear distinction between the permeabilities of the various parts of the open hole.

L'apprentissage de dictionnaire pour la représentation parcimonieuse est bien connu dans le cadre de la résolution de problèmes inverses. Les méthodes d'optimisation et les approches paramétriques ont été particulièrement explorées. Ces méthodes rencontrent certaines limitations, notamment liées au choix de paramètres. En général, la taille de dictionnaire doit être fixée à l'avance et une connaissance des niveaux de bruit et éventuellement de parcimonie sont aussi nécessaires. Les contributions méthodologies de cette thèse concernent l'apprentissage conjoint du dictionnaire et de ces paramètres, notamment pour les problèmes inverses en traitement d'image. Nous étudions et proposons la méthode IBP-DL (Indien Buffet Process for Dictionary Learning) en utilisant une approche bayésienne non paramétrique. Une introduction sur les approches bayésiennes non paramétriques est présentée. Le processus de Dirichlet et son dérivé, le processus du restaurant chinois, ainsi que le processus Bêta et son dérivé, le processus du buffet indien, sont décrits. Le modèle proposé pour l'apprentissage de dictionnaire s'appuie sur un a priori de type Buffet Indien qui permet d'apprendre un dictionnaire de taille adaptative. Nous détaillons la méthode de Monte-Carlo proposée pour l'inférence. Le niveau de bruit et celui de la parcimonie sont aussi échantillonnés, de sorte qu'aucun réglage de paramètres n'est nécessaire en pratique. Des expériences numériques illustrent les performances de l'approche pour les problèmes du débruitage, de l'inpainting et de l'acquisition compressée. Les résultats sont comparés avec l'état de l'art.Le code source en Matlab et en C est mis à disposition. / Dictionary learning for sparse representation has been widely advocated for solving inverse problems. Optimization methods and parametric approaches towards dictionary learning have been particularly explored. These methods meet some limitations, particularly related to the choice of parameters. In general, the dictionary size is fixed in advance, and sparsity or noise level may also be needed. In this thesis, we show how to perform jointly dictionary and parameter learning, with an emphasis on image processing. We propose and study the Indian Buffet Process for Dictionary Learning (IBP-DL) method, using a bayesian nonparametric approach.A primer on bayesian nonparametrics is first presented. Dirichlet and Beta processes and their respective derivatives, the Chinese restaurant and Indian Buffet processes are described. The proposed model for dictionary learning relies on an Indian Buffet prior, which permits to learn an adaptive size dictionary. The Monte-Carlo method for inference is detailed. Noise and sparsity levels are also inferred, so that in practice no parameter tuning is required. Numerical experiments illustrate the performances of the approach in different settings: image denoising, inpainting and compressed sensing. Results are compared with state-of-the art methods is made. Matlab and C sources are available for sake of reproducibility.