Le travail présenté dans cette thèse porte sur la modélisation et la commande non linéaire le système du forage pétrolier. La problématique principale consiste à proposer un modèle de forage non linéaire et à élaborer un ensemble de stratégies de commande avancées pour améliorer les performances du trépan.<br />A l'heure actuelle, les commandes des systèmes de forage pétroliers sont réalisées à l'aide de régulateurs classiques PI ou PID. Leurs performances se dégradent notablement lors de vibration de torsion causée par le « stick- slip », phénomène qui peut provoquer des couples non linéaires et donc endommager le trépan et le train de tiges. Par conséquent, il est nécessaire de disposer de systèmes de commande plus performants afin de mieux répondre aux exigences et aux performances demandées.<br />En dépit des diverses lois de commandes linéaires tentatives pour commander ce système mais qui ne garantissent pas la stabilité vis-à-vis des incertitudes. Ceci nous a conduits à proposer une modèle non linéaire en prenant en compte la totalité de l'information du phénomène du stick-slip, ce qui a conduit à une modélisation globale du forage rotary approchant donc au mieux le comportement réel de la plate-forme pétrolière. Pour cela, nous avons suggéré deux approches pour améliorer les performances du système de forage rotary:<br />• Commande par la technique de linéarisation par bouclage.<br />• Développement d'une commande par la technique des modes glissant associé à la commande linéarisation par Bouclage.<br />Comme le système de forage est généralement, est un système soumis aux incertitudes paramétriques, nous avons proposé, de plus, une nouvelle loi de commande non linéaire par Backstepping permettant de mieux améliorer la réponse de la vitesse de rotation du trépan.<br />Ces trois stratégies de commandes apparaissent la première fois dans ce domaine et la approche de Backsteping s'avère la meilleure. Les résultats de simulations obtenus dans les trois cas sont concluants et la synthèse par Backstepping montré la supériorité de, où la question de robustesse a pu être considérée.

L'objet de cette thèse est l'étude de quelques aspects de <br />la notion d'hyperbolicité, plus particulièrement de la <br />relation qui existe entre celle-ci et la nature bien posée <br />d'un problème de Cauchy obtenu à partir d'un système<br />d'équations aux dérivées partielles issu de la mécanique <br />des fluides ou la réalisation de la simulation numérique <br />d'un tel problème.<br /><br />Dans un premier temps, nous rappelons en quoi la notion de<br />linéarisation d'un système d'équations aux dérivées <br />partielles semble naturelle à l'étude de ce système et <br />comment, de l'étude de ces problèmes linéarisés, plus <br />précisément de leur nature bien posée c'est-à-dire de leur <br />stabilité, découle la notion d'hyperbolicité.<br /><br />Nous étudions ensuite le cas particulier d'un modèle à <br />quatre équations pour un écoulement bifluide comportant des <br />termes de diffusion pour les équations de quantité de <br />mouvement. Nous montrons alors que, bien que, pour ce <br />système, l'ajout des termes de diffusion n'entraîne pas <br />l'hyperbolicité du modèle obtenu, les problèmes de Cauchy <br />construits à partir de la linéarisation de ce système, <br />autour d'un état constant, sont désormais bien posés.<br /><br />Enfin, nous considérons le cas d'un modèle à cinq équations <br />pour un écoulement bifluide. Ce modèle ne nécessite pas de <br />loi de fermeture algébrique (équations d'état ou lois <br />tabulées) mais comporte une équation aux dérivées <br />partielles portant sur la pression. Le système ainsi <br />obtenu n'est pas hyperbolique mais les valeurs propres de <br />l'opérateur d'advection sont toutes réelles. La simulation<br />numérique d'un écoulement régi par ce modèle, pour le cas <br />test du robinet de Ransom, ne fait néanmoins pas apparaître <br />les instabilités numériques que la nature mal posée du<br />linéarisé nous faisait craindre et qui sont présentes dans <br />les simulations réalisées à partir du modèle isentropique <br />classique à quatre équations.

Le travail réalisé dans cette thèse s'inscrit dans le cadre de l'étude de la commande optimale des systèmes non linéaires à temps discret par des méthodes de type programmation dynamique. Une formulation particulière de l'équation fonctionnelle considérée ici nous a amené à examiner le problème de commandabilité que nous avons encore qualifié de problème d'inversion entrée-états. L'étude comporte deux phases. La première phase consiste à analyser la propriété de commandabilité des systèmes. Plusieurs approches sont considérées: composition de fonction, géométrie différentielle et algèbre différentielle. On propose ensuite une condition nécessaire de commandabilité basée sur l'utilisation du principe du minimum. Finalement, une application directe de cette condition nécessaire de commandabilité est développée pour les systèmes bilinéaires. Une fois la propriété de commandabilité du système vérifiée, la seconde phase consiste à résoudre le problème d'inversion. On propose tout d'abord une technique algébrique conduisant à la détermination d'une séquence de commandes à partir de la matrice de commandabilité. On propose ensuite une approche basée sur une formulation du problème d'inversion en terme d'un problème de commande optimale. On considère également les techniques de linéarisations exactes qu'on aborde selon deux voies. On présente dans la première voie une approche de l'algèbre différentielle. Dans la deuxième voie, on propose une approche de géométrique différentielle dans laquelle on définit une séquence de distributions basée sur l'utilisation de champs vecteurs associés au système, ainsi qu'une sortie virtuelle linéaire. La transformation et le bouclage linéarisant sont calculés à partir de ces distributions. Finalement, ces résultats sur l'accessibilité et la linéarisation exacte sont présentés sur des applications pratiques dans le cadre de la commande optimale par programmation dynamique.

Les modulations multiporteuses apparaissent aujourd'hui comme une technologie éprouvée pour la transmission de données à haut-débits sur des canaux pouvant être très perturbés. L'OFDM (Orthogonal Frequency Division Multiplexing) a d'ailleurs été choisie dans plusieurs normes de télécommunications (ADSL, Wi-Max, IEEE 802.11a/g/n, LTE, DVB,...). Cependant un des inconvénients de ce type de modulation est la forte variation de la puissance instantanée à transmettre. Cette propriété rend ces modulations très sensibles aux non-linéarités des composants analogiques, en particulier celles de l'amplificateur de puissance à l'émission. Or l'amplificateur de puissance est un élément déterminant dans une chaîne de communication dans la mesure où il a une influence prépondérante sur le bilan global de la transmission en termes de puissance, de rendement et de distorsion. Plus l'on souhaite que l'impact de ses non linéarités soit faible et plus son rendement est faible, et inversement. Il est donc nécessaire d'effectuer un compromis linéarité/rendement.L'objectif de la thèse est d'éviter cette détérioration du rendement tout en conservant de bonnes performances de linéarité, de surcroit pour des signaux OFDM. Pour ce faire nous proposons d'utiliser conjointement des méthodes de linéarisation (prédistorsion numérique en bande de base) et d'amélioration du rendement (envelope tracking) de l'amplificateur de puissance ainsi qu'une méthode de réduction de la dynamique du signal (active constellation extension). La prédistorsion numérique classique échouant aux fortes puissances, nous proposons une méthode d'amélioration de cette technique à ces puissances. Nos résultats sont validés par des mesures sur un amplificateur de puissance 50W. Nous proposons également une association des méthodes permettant d'améliorer simultanément les performances en terme de linéarité hors bande et de rendement en minimisant les dégradations des performances de taux d'erreur binaire.

On s'intéresse dans ce travail à la description des enrichissements des disques de Siegel d'une fraction rationnelle f. Dans un premier temps nous étudions les enrichissements qui sont définis sur un ouvert de la grande orbite d'un disque de Siegel donné. Ce sont nécessairement des applications qui commutent à f là où les compositions ont un sens. Ce sont donc des applications linéaires en coordonnées linéarisantes. Le résultat principal de ce travail est que l'on peut obtenir toutes les applications linéaires en coordonnées linéarisantes définies sur un sous-disque du disque de Siegel de f. Pour démontrer ce résultat nous utilisons la compacité des applications linéarisantes normalisées, le théorème des fonctions implicites dans l'espace des fractions rationnelles de degré fixé et une étude du comportement du rayon d'univalence des applications linéarisantes. Nous identifions également les approches donnant lieu à des enrichissements définis ou à valeurs dans le disque de Siegel tout entier (enrichissements maximaux). Au passage nous généralisons aux limites avec ordre de contact fini par rapport au cercle unité un théorème de JC.Yoccoz sur le comportement du rayon d'univalence pour la famille quadratique lorsque le paramètre converge vers un nombre complexe de module un et d'argument un nombre de Brjuno.Ensuite, nous nous intéressons au cas où f a plusieurs cycles de disques de Siegel. Nous utilisons le théorème de transversalité d'A.Epstein pour décrire les enrichissements de f dans ce cas là. La linéarisabilité de f et la convergence des applications linéarisantes permet de transférer le problème de la description des enrichissements de Siegel de f à un problème de limite géométrique de sous-semigroupes de l'ensemble des nombres complexes non-nuls engendrés par un élément. Nous donnons dans ce travail un modèle topologique de l'adhérence de cet ensemble de sous-semigroupes. Nous déduisons de ces résultats une interprétation en terme de convergence géométrique de dynamiques de polynômes quadratiques et une description des points d'accumulation, pour la topologie de Hausdorff sur les compacts non-vides, des ensembles de Julia lorsque le paramètre tend vers un paramètre de Siegel. / In this work we are interested in giving the description of Siegel discs enrichments of a rational map f. We first study the case of enrichments that are defined on an open subset of the grand orbit of a given Siegel disc. These maps commute with f where it makes sense. Thus they are linear in linearizing coordinates. The main result of this work is that we can obtain all linear maps in linearizing coordinates that are defined in a subdisc of the Siegel disc. For this we use the compactness of the set of normalized linearizing maps, the implicit functions theorem in the space of rational maps with fixed degree and a study on the behaviour on the univalent radius of the linearizing maps. We identify approaches giving enrichments that are defined or take values on the whole Siegel disc (maximal enrichments). We generalize to finite order of contact approaches with respect to the unit circle a theorem of JC.Yoccoz on the behaviour of the univalent radius for the quadratic family when the parameter converges to a complex number of modulus one with argument a Brjuno number.We then focus on the case where f has more than one Siegel disc. We make use of A.Epstein's transversality theorem to describe Siegel enrichments of f in this case. The linearisability of f and the convergence of the linearizing maps reduces the problem of Siegel enrichments description to a geometric limit problem on one generated closed sub-semigroups ofthe set of non zero complex numbers. We give in this work a topological model fot the closure of this set of sub-semigroups.We deduce from these results an interpretation in terms of geometric convergence of quadratic polynomial dynamics and we describe the accumulation points (for the Hausdorff topology on non empty compact subsets) of Julia sets when the parameter converges to a Siegel parameter.

L'observateur est un outil essentiel pour accéder à l'information, non mesurable directement, d'un système dynamique. Dans le cas des systèmes non linéaires, il y a une grande carence concernant la synthèse d'observateur. Motivée par l'absence d'une solution générique, cette thèse élargit la classe des systèmes non linéaires pour lesquels on peut appliquer les observateurs connus. Dans l'approche adoptée, le système non linéaire est transformé à travers un changement de coordonnées sous forme normale d'observabilité qui a une structure adéquate à la synthèse d'observateurs. Les difficultés liées aux conditions d'existence d'un changement de coordonnées sont mises en évidence et des solutions sont proposées. Par conséquent, la classe des systèmes non linéaires qui peuvent se mettre sous une forme normale d'observabilité est élargie. Dans un premier temps, nous avons proposé une forme normale d'observabilité étendue dépendante de la sortie en augmentant l'espace d'état par des variables auxiliaires. Ainsi, nous avons établi les conditions nécessaires et suffisantes à l'existence d’un changement de coordonnées permettant d’obtenir une telle forme. En outre, nous avons proposé, pour certains modèles, des procédures heuristiques pour la mise sous forme normale d'observabilité étendue dépendante de la sortie. Dans un deuxième temps, nous avons traité la mise sous forme normale d'observabilité d'un système non linéaire partiellement observable. Enfin, nous avons abordé la transformation d'une classe spéciale de systèmes non linéaires sous la forme normale d'observabilité réduite. L'efficacité et l'intérêt des méthodes développées sont établis au travers de plusieurs applications. / Observer is an essential means to access to no-measurable information of a dynamical system. In the case of nonlinear systems, there is a great deficiency concerning the observer design theory. Motivated by the lack of a generic solution to observer design, this thesis enlarges the class of nonlinear systems which admit a standard observer. Using a geometrical approach, the considered nonlinear systems are transformed through a change of coordinates into observer form, which has an adequate structure to the observer design. The difficulties related to the conditions on the existence of such a change of coordinates are highlighted and solutions are proposed. Therefore, the class of nonlinear systems which can be transformed under an observer form is expanded. Firstly, we proposed an extended output depending observer form which does not preserve the size of the original state space and we established the sufficient conditions for the existence of a change of coordinates enabling to construct the proposed form. Moreover, we proposed a heuristic procedure to construct the extended output depending observer form of some models. Secondly, we dealt with the observer form for some partially observable nonlinear systems. Finally, we discussed the transformation into the reduced observer form for a class of nonlinear systems. The efficiency and interest of the developed methods is established through several applications.

En ingénierie sismique, il est admis que le comportement d’une structure soumise à de forts séismes soit caractérisé par des boucles d’hystérésis qui peuvent être amples ou étroites selon le type de structure impactée. La prise en compte de ce type de comportement non-linéaire dans un calcul temporel présente des difficultés liées à l’identification des paramètres, au coût numérique élevé, au risque de non-convergence. Dans ce contexte, la méthode de linéarisation équivalente, a été introduite en géotechnique dès les années 70. Elle reste peu utilisée dans le domaine des structures malgré les efforts de nombreux auteurs. Ce travail de thèse a pour objet l’étude du comportement linéaire équivalent dans le contexte des méthodes simplifiées d'évaluation de la réponse non-linéaire d'une structure en ingénierie sismique. Nous passons en revue les critères de linéarisation adoptés par les différentes méthodes qui recherchent l’équivalence (1) du déplacement maximum ou (2) de la quantité d’énergie dissipée ou (3) de la force de rappel. Nos analyses montrent que ces trois critères ne sont pas pertinents et/ou efficaces, conduisant à des méthodes peu robustes qui conduisent dans certains cas à des résultats inexplicables. Nous montrons le rôle important, négligé par toutes les méthodes disponibles, du contenu fréquentiel respectif des signaux et du système dans la détermination de la ductilité appelée. Sur cette constatation, nous introduisons une nouvelle méthode de linéarisation équivalente basée sur la fonction de transfert. Nous utilisons cette méthode pour explorer un plan d’expérience numérique dans lequel nous calculons les caractéristiques de fréquence et d’amortissement équivalents en fonction de la ductilité appelée pour différente configurations caractérisées par (a) le rapport entre fréquence de l’oscillateur et fréquence centrale du signal excitateur, (b) la pente d’écrouissage et (c) le modèle de comportement qui varie continument de élastoplastique à endommageant. Nous proposons deux nouvelles approches du comportement linéaire équivalent. La première, visant à améliorer la procédure statique non-linéaire de l’ATC40, utilise la rigidité sécante et le déplacement maximal. Elle fait intervenir une estimation de l’amortissement différente de celle de l’ATC40. Sa pertinence est établie par le fait qu’elle permet d’évaluer avec exactitude le déplacement maximal de systèmes canoniques non-linéaires. La seconde consiste à restituer la dynamique de la réponse d'un oscillateur non-linéaire au travers de la fonction de transfert. Sa pertinence est démontrée au travers des critères d’Anderson, avec notamment un critère relatif au spectre transféré. La détermination du comportement linéaire équivalent par fonction de transfert est validée sur des structures réelles au travers des essais sur voiles en béton armé (SAFE) et sur systèmes des tuyauteries (BARC et EPRI) / In earthquake engineering, it is common that the behaviour of a structure undergoing a strong motion is characterized by wise or narrow hysteresis loops depending on the type of behaviour of the structure. Considering this non-linear behaviour in a transient calculation requires a huge need of resources in terms of calculation time and memory. In this context, the method of equivalent linearization, consisting in the evaluation of the non-linear response of the structure has been introduced by geotechnical engineers In the 1970s. Despite efforts of many authors, this method is still not used in structural field. The goal of this research is to examine the linear equivalent behaviour in the context of the simplified method of evaluating the non-linear response of a structure in earthquake engineering. We review the criteria of equivalence adopted by many methods searching for the equivalence of (1) the maximum of displacement or (2) quantity of dissipated energy or (3) the restore force. Our argumentative analyses carry out that these three criteria are not pertinent and/or efficient. This leads, in some cases, to some unexplained results. We show the important role, which is mostly neglected in existing method, of frequency content while evaluating the ductile demand. Based on this recognition, we introduce a new method of equivalent linearization based on the transfer function. We use this method in order to explore a numerical experimental plan in which we calculate the equivalent characteristics (frequency and damping) versus the ductile demand for different configuration characterized by (a) the ratio between the frequency of the oscillator and the central frequency of the input signal, (b) the hardening and (c) the behaviour which covers the elastoplastic and damaged ones. We propose two new approaches of the linear equivalent behaviour. The first one, aiming to improve the non-linear static procedure of ATC40, use the secant stiffness and the maximal displacement. This approach consists in an estimation of damping which is different to ATC40. Its pertinence is established by evaluating with accuracy the maximal displacement of the canonical non-linear systems. The second approach consists in restitution of the dynamic of the response of a non-linear oscillator by using the transfer function. The pertinence of this proposition is shown through the criteria of Anderson, especially in terms of transferred motion. In this effect, the linear equivalent behaviour based on the transfer function allows to cope the transferred motion through the non-linear oscillator without performing the non-linear transient calculation. The validation of the linear equivalent behaviour based on the transfer function has been examined on real structures through some experimental tests such as the reinforced concrete wall (SAFE) or piping systems (BARC and EPRI)

Cette thèse porte sur les techniques d'analyse formelle de systèmes hybrides à dynamiques continues non linéaire. Ses contributions portent sur les algorithmes d'atteignabilité et sur les problèmatiques liées à la representation des ensembles atteignables. This thesis deals with formal analysis of hybrid system with non linear continous dynamic. It contributes to the fields of reachability analysis algorithm and the set representation. / In this thesis, we presented our contributions to the formal analysis of dynamical systems. We focused on the problem of efficiently computing an accurate approximation of the reachable sets under nonlinear dynamics given by differential equations. Our aim was also to design scalable methods which can handle large systems. The first contribution of this thesis concerns the dynamic hybridization technique for a large class of nonlinear systems. We focused on the hybridization domain construction such that the linear interpolation realized in this domain ensures a desired error between the original system trajectories and those computed with the approximated system. We propose a construction method which tends to maximize the domain volume which reduce the number of creation of new domains during the analysis. The second research direction that we followed concerns a subclass of nonlinear dynamical systems which are the polynomial systems. Our results for the reachability analysis of these systems are based on the Bernstein expansion properties. We approximate an initial reachability computation (which requires solving polynomial optimization problems) with an accurate over-approximation (which requires solving linear optimization problems). The last theoretical contribution concerns the reachability analysis of linear systems with polyhedral input which often result from approximation of nonlinear systems. We proposed a technique to refine

La liquéfaction des sols saturés lors des séismes est l’un des problèmes les plus importants auxquels sont confrontés les ingénieurs. Il n’y a guère eu de séisme majeur sans au moins quelques cas de liquéfaction. Des tassements, des basculements de bâtiments, des écoulements latéraux, des cônes de liquéfaction et des instabilités de pentes, sont certaines de ses manifestations. La conception sismique des centrales nucléaires et autres installations critiques comprend systématiquement une évaluation du risque de liquéfaction.Dans ce cadre, des méthodes de nature entièrement empirique sont couramment utilisées en ingénierie. Ces approches procurent des marges à la conception et des limitations d’utilisation. Pour exploiter ces marges en situation de réévaluation, on a recours à des calculs transitoires non linéaires avancés dans lesquels on doit modéliser finement la loi de comportement du sol pour mettre en évidence les montées de pression interstitielle.Ces derniers calculs sont coûteux en termes de temps et de compétences numériques. L’objectif de ce travail de recherche est notamment de réduire les conservatismes en vigueur lors de l’utilisation de la méthode simplifiée sans pour autant mettre en œuvre d’emblée les méthodes les plus sophistiqués. On propose pour cela une méthode de complexité intermédiaire qui élargit l’applicabilité des modèles semi-empiriques pour une analyse plus fine du risque sismique.Dans un premier temps, en partant d’un calcul linéaire équivalent conventionnel, une nouvelle approche pour la prise en compte de la montée de pression interstitielle est proposée sous le nom « X-ELM ». Le modèle de comportement employé est basé sur la relation entre la pression interstitielle et les déformations volumiques plastiques. La nouvelle approche «X-ELM » est utilisée pour modéliser la réponse des sols pour le séisme de Tōhoku (Mw=9.0) à la ville d’Urayasu au Japon. Le modèle a été appliqué sur douze profils de sols différents. L’étude de ces cas rend possible la validation du modèle par comparaison des résultats des calculs aux observations in situ. Le modèle peut donc être considéré comme un outil fiable pour la prédiction de déclenchement de liquéfaction des sols saturés.Ensuite, un outil de prédiction rapide a été conçu en se basant sur des approximations de processus aléatoire, sur les propriétés mécaniques de base du sol et sur les caractéristiques du chargement sismique. Outre son ampleur, une caractéristique importante du signal sismique d’entrée est sa durée qui peut conduire à de fortes non linéarités et à un état de liquéfaction étendu. En considérant donc la durée de phase forte, le spectre de réponse, la fréquence propre du modèle et les caractéristiques de densification du sol, l’outil de prédiction proposé procure des estimations rapides du taux de montée de pression interstitielle et du tassement pré-liquéfaction sans devoir exécuter des calculs transitoires.Enfin, un modèle 2D de barrage est étudié, en examinant l’influence de la montée de pression interstitielle et celle des déformations de cisaillement sur la réponse sismique de l’ouvrage. Un calcul linéaire équivalent adapté aux situations bidimensionnelles est élaboré et le prédicteur est employé pour évaluer la montée de pression interstitielle. On compare les résultats de la simulation aux observations in situ, piézométriques et accélérométriques.En conclusion, ce travail de recherche fournit des méthodes et outils de calculs numériques performants et accessibles aux ingénieurs pour l’évaluation sismique des profils de sols et des ouvrages en terre tels que digues ou barrages / The liquefaction of saturated soils during earthquakes is one of the most important problems facing engineers. There has hardly been a major earthquake without at least some cases of liquefaction. Settlements, tilting of buildings, lateral flows, sand boilings and slope instabilities have been some of its manifestations. The seismic design of nuclear power plants and other critical facilities systematically includes a liquefaction risk assessment.In this context, fully empirical methods are commonly used in engineering. These approaches provide design margins and limitations of use. To exploit these margins in a re-evaluation situation, we use advanced nonlinear transient calculations in which the soil behavior must be finely modeled to highlight the pore-water pressure build-up. These last calculations are expensive in terms of time and numerical skills. The objective of this research work is to reduce the conservatisms in force when using the simplified method without necessarily implementing the most sophisticated methods from the outset. To this end, we propose a method of intermediate complexity that broadens the applicability of semi-empirical models for a more detailed analysis of seismic risk.First, starting from a conventional equivalent linear calculation, a new approach for taking into account excess pore pressure is proposed under the name "X-ELM". The behavioral model employed is based on the relationship between pore pressure and plastic volumetric deformations. The new "X-ELM" approach is used to model soil response in the city of Urayasu,Japan during the Tohoku earthquake (Mw = 9.0). The model has been applied to twelve different soil profiles. The study of these cases makes possible the validation of themodel by comparing the results of the calculations with the observations in situ. The model can therefore be considered as a reliable tool for the prediction of liquefaction triggering of saturated soils.Then, a prediction tool was designed based on random process approximations, the basic mechanical properties of the soil and the characteristics of the seismic loading. In addition to its magnitude, an important feature of the input seismic signal is its duration which can lead to strong nonlinearities and an extended liquefaction state. Considering the strong phase duration, the response spectrum, the natural frequency of the model and the characteristics of soil compaction, this tool provides fast estimations of the rate of pore pressure build-up and pre-liquefaction settlement without having to perform transient calculations.Finally, a 2D dam model is studied, by examining the influence of excess pore pressure and that of the shear strains on the seismic response of the structure. An equivalent linear computation adapted to two-dimensional situations is elaborated and the predictor is used to evaluate pore pressure increase. Simulation results are compared with in situ, piezometric and accelerometric observations.In conclusion, this research work provides methods and tools of numerical computation that are efficient and accessible to engineers for the seismic evaluation of soil profiles and earth structures such as dikes or dams

La tomographie est un ensemble de techniques visant à reconstruirel’intérieur d’un objet sans toucher l’objet lui même comme dans le casd’un scanner. Les principes théoriques de la tomographie ont été énoncéspar Radon en 1917. On peut assimiler l’objet à reconstruire à une image,matrice, etc.Le problème de reconstruction tomographique consiste à estimer l’objet àpartir d’un ensemble de projections obtenues par mesures expérimentalesautour de l’objet à reconstruire. La tomographie discrète étudie le cas où lenombre de projections est limité et l’objet est défini de façon discrète. Leschamps d’applications de la tomographie discrète sont nombreux et variés.Citons par exemple les applications de type non destructif comme l’imageriemédicale. Il existe d’autres applications de la tomographie discrète, commeles problèmes d’emplois du temps.La tomographie discrète peut être considérée comme un problème d’optimisationcombinatoire car le domaine de reconstruction est discret et le nombrede projections est fini. La programmation mathématique en nombres entiersconstitue un outil pour traiter les problèmes d’optimisation combinatoire.L’objectif de cette thèse est d’étudier et d’utiliser les techniques d’optimisationcombinatoire pour résoudre les problèmes de tomographie. / The tomographic imaging problem deals with reconstructing an objectfrom a data called a projections and collected by illuminating the objectfrom many different directions. A projection means the information derivedfrom the transmitted energies, when an object is illuminated from a particularangle. The solution to the problem of how to reconstruct an object fromits projections dates to 1917 by Radon. The tomographic reconstructingis applicable in many interesting contexts such as nondestructive testing,image processing, electron microscopy, data security, industrial tomographyand material sciences.Discete tomography (DT) deals with the reconstruction of discret objectfrom limited number of projections. The projections are the sums along fewangles of the object to be reconstruct. One of the main problems in DTis the reconstruction of binary matrices from two projections. In general,the reconstruction of binary matrices from a small number of projections isundetermined and the number of solutions can be very large. Moreover, theprojections data and the prior knowledge about the object to reconstructare not sufficient to determine a unique solution. So DT is usually reducedto an optimization problem to select the best solution in a certain sense.In this thesis, we deal with the tomographic reconstruction of binaryand colored images. In particular, research objectives are to derive thecombinatorial optimization techniques in discrete tomography problems.