La propagation sismique dans les milieux poreux est classiquement associée à des phénomènes de dispersion et d’atténuation des ondes sous l’effet des mouvements fluides. Dans certaines conditions, celle-ci peut également être associée à une conversion d’énergie sismique en énergie électromagnétique dite « sismoélectrique ». La théorie des phénomènes sismoélectriques combinant la théorie de l’électrocinétique à la poroélasticité de Biot, repose en grande partie sur les développements de Pride (1994). Sur la base de ces développements théoriques, Pride et Haartsen (1996) relient le champ électrique cosismique à l’accélération sismique qui le génère par une fonction de transfert. Nous proposons une étude quantitative des couplages sismoélectriques en vue de valider la théorie de Pride et sa généralisation en milieu non saturé. Dans ce but, nous avons développé une expérience en laboratoire sur un sable de quartz non-consolidé, menée dans la gamme du kilohertz sur la base d’un dispositif d’acquisition électrique modulable. Deux méthodes de traitement des signaux sont proposées, l’une temporelle, l’autre spectrale, permettant d’obtenir une analyse complète des vitesses de phase, atténuations et fonctions de transfert. Les expériences réalisées se sont focalisées sur l’étude du rôle de la conductivité du fluide et de la saturation en eau dans le phénomène sismoélectrique cosismique. Une étude time-lapse a ainsi pu être réalisée dans des situations de changements de salinité et de teneur en eau. Dans tous les contextes, l’étude quantitative des rapports d’amplitudes des champs sismoélectriques et sismiques E/ü montre une bonne corrélation avec les prédictions théoriques. Par ailleurs, l’étude des variations de saturation dans une gamme allant de la saturation résiduelle en eau (Sw = 0.3) à la saturation totale, montre que: i) les atténuations et fonctions de transfert ont des comportements reliés à la distribution des fluides qui influencent fortement les propriétés mécaniques du milieu ; ii) une inversion de polarité du champ sismoélectrique peut être observée dans le cas très particulier des milieux non consolidés. / Seismic propagation within porous media is usually associated with wave attenuation and dispersion phenomena related to fluid flow. Under certain circumstances, it may also be correlated to a conversion of seismic into electromagnetic energy known as “seismoelectric”. The understanding of seismoelectric phenomena, combining the theory of electrokinetic to Biot’s poroelasticity, relies mainly on the formulation by Pride (1994). On basis of these theoretical developments, Pride and Haartsen (1996) defined a transfer function expressing the link between the coseismic seismoelectric field and the seismic acceleration at its origin. We propose a quantitative analysis of coseismic seismoelectric couplings with the purpose of validating Pride’s theory and generalizing it to partially saturated media. With this aim in view we developed a laboratory experiment involving an adjustable device for electric acquisitions, conducted within the kilohertz range on unconsolidated quartz sand. Experimental data were subsequently processed in both time and frequency domains, enabling a full analysis that embraces phase velocities, attenuations and transfer functions. The conducted experiments focused on the impact of fluid conductivity and water saturation with regard to the coseismic seismoelectric phenomenon. Time-lapse monitoring were accordingly run under varying salinity or water content. In all scenarios, the quantitative analysis of the electric-to-seismic amplitude ratio E/ü appeared in good agreement with theoretical projections. Moreover, investigations of saturation variations, ranging from the residual water saturation (Sw = 0.3) to full saturation, showed that: i) the behavior of attenuations and transfer functions are directly related to fluid distribution, that greatly impacts the mechanical properties of the medium; ii) in the very peculiar case of unconsolidated media, polarity inversion of the coseismic seismoelectric field may be experienced.

Un nanocalorimètre permet l'analyse thermique de très faibles quantités d'échantillons (quelques nanogrammes ou picogrammes), ainsi que l'étude de films minces dont l'épaisseur varie de quelques nanomètres à plus d'un micron. Les vitesses de chauffe et de refroidissement sont nettement plus élevées que celles réalisées avec une DSC classique : les vitesses de chauffe peuvent atteindre 103 à 106 K/s, par conséquent, les mesures réalisées avec ce type d'appareil sont très rapides (quelques millisecondes). Du fait de sa sensibilité élevée, des vitesses de chauffe rapides atteintes et de l'acquisition rapide des données, le nanocalorimètre peut être utilisé pour la caractérisation et la détection de quelques nanogrammes de matériaux énergétiques. Les objectifs de cette thèse sont d'une part de mettre au point une procédure de calibration des capteurs nanocalorimétriques (calibration de la température, de la puissance, de la masse du microcristal à analyser) et d'autre part de caractériser des matériaux énergétiques afin de pouvoir effectuer des analyses quantitatives en vue d'applications pour la détection d'explosifs. Les matériaux énergétiques étudiés sont des films de nitrocellulose, des cristaux de penthrite, d'hexogène (principal constituant du C4), de 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12- hexaazaisowurtzitane (Cl20), et des nano-cristaux d'explosifs. Le travail réalisé a montré qu'il est possible de rapidement différencier des explosifs par leurs températures de fusion, de décomposition et d'évaporation. Il est aussi possible de déterminer des paramètres cinétiques d'un cristal isolé d'explosif. / Calorimetry is one of the main techniques of thermal analysis. Most of physical or chemical modifications of material are associated with thermal effects whereby heat is absorbed (i.e., melting) or released (i.e., thermal decomposition). Typically, calorimetric experiments are performed with Differential Scanning Calorimetry (DSC), which measures the heat flux absorbed or released by the sample following the same temperature program as a reference material. In these experiments, measurements are typically carried out on a few milligrams of sample. However, for many applications one has to handle nanograms or even picograms of sample. One of such applications is relevant to studies of materials which can release a significant amount of energy during their decomposition (energetic materials). Calorimetry able to handle nanograms of sample could find potential applications in the field of explosives detection. Nanocalorimetry allows to heat small amounts of sample (a few nanograms to a few hundred picograms) at extreme heating rates, i.e. up to one million °C/s. The temperature increase can initiate several phenomena in energetic materials, therefore the calorimetry could be an appropriate technique to characterize and to detect energetic materials. The energetic materials used in this study are nitrocellulose (NC), hexogen (RDX), 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-ltexaazaisowurtzitane (CL-20) and penthrite (PETN). The manuscript presents our results on the nanocalorimeter calibration, on the thermal behaviour of the explosives studied with nanocalorimetry and also includes an evaluation of nanocalorimetry as a tool for explosives detection.

La magnésite (MgCO3) est la forme anhydre la plus stable d'une série de carbonates de magnésium qui présentent différents degrés d'hydratation. Malgré sa rareté dans les environnements naturels, elle constitue une phase minérale fondamentale pour le piégeage minéral permanent du CO2. La détermination expérimentale des vitesses de précipitation et de dissolution de la magnésite dans des conditions représentatives de la séquestration géologique est donc fondamentale pour l'estimation du potentiel de séquestration de CO2 par cette phase dans les basaltes et dans les roches ultrabasiques et pour l'optimisation des procédés de stockage du CO2. Nous avons mesuré les vitesses de précipitation de la magnésite en utilisant des réacteurs à circulation et des réacteurs fermés, en fonction de la température (100 ≤ T ≤ 200 °C), de la composition de la solution aqueuse et de la pression partielle de CO2 (de 0 à 30 bar). Les vitesses mésurées sont indépendantes de la force ionique de la solution pour 0.1 M < I < 1.1 M, mais elles diminuent significativement avec l'augmentation de l'activité des ions CO32- pour des pH supérieures à 8. Les vitesses mesurées dans les réacteurs à circulation sont cohérentes avec le modèle de coordination chimique surfacique de Pokrovsky et al. (1999) selon lequel les vitesses de précipitation de la magnésite sont proportionnelles à la concentration des sites surfaciques >MgOH2+. L'étude des vitesses de cristallisation conduite par microscopie à force atomique hydrothermale (HAFM) a montré un bon accord entre les vitesses déduites de mesures microscopiques et les vitesses macroscopiques et a aussi démontré que la précipitation de la magnésite s'effectue selon un mécanisme de croissance spirale. Suivant les observations effectuées par AFM, ce mécanisme contrôle la vitesse de croissance de la magnésite dans un grand intervalle de température et d'indice de saturation (15≤ Ω ≤ 200 pour 80 ≤ T < 120 °C). En raison de l'inhibition de la précipitation de la magnésite par les ions carbonates, il est recommandé d'opérer sous des pressions partielles de CO2 assez élevées, ce qui présente en outre l'avantage d'accélérer la cinétique de dissolution des silicates magnésiens, grâce à l'acidification de la solution par le CO2. La détermination des vitesses de dissolution de la magnésite dans des réacteurs à circulation à 150 et 200 °C et en milieu neutre à alcalin nous a permis d'améliorer le modèle de complexation de surface et d'étendre son application aux températures considérées. La diminution des vitesses de dissolution observée de 150 à 200 °C peut être expliquée par l'augmentation, en fonction de la température, de la carbonatation et de l'hydrolyse des sites >MgOH2+ qui contrôlent la vitesse de dissolution de la magnésite. Des températures supérieures à 100 °C qui entrainent une diminution de la vitesse de dissolution de la magnésite et des autres carbonates sont donc favorables au stockage de CO2 sous forme dissoute dans les aquifères profonds riches en minéraux carbonatés. L'utilisation d'une cellule à électrodes d'hydrogène (HECC) nous a permis de préciser les données cinétiques à proximité de l'équilibre grâce à la détermination précise du produit de solubilité de la magnésite en fonction de la température (50-200 °C). De plus, ces mesures nous ont permis de générer les propriétés thermodynamiques de la magnésite et de les comparer à celles obtenues par mesures calorimétriques et par équilibres de phases. Les résultats de cette étude représentent une importante contribution à la compréhension des cinétiques de réaction des minéraux carbonatés dans les systèmes hydrothermaux et permettent de proposer une base de données essentielle pour la quantification des réactions de dissolution/précipitation des carbonates dans les systèmes complexes. En outre, ce travail fournit des contraintes cinétiques pour la modélisation géochimique des processus de séquestration du CO2 et sera utile à l'évaluation de l'impact et des risques liés au stockage de CO2 à long terme.

The notion of visibility is central in computer graphics and computational geometry. A visibility computation amounts to determining the shape, or merely the existence of the set of segments in space linking two specified objects without crossing any other. We examine the problem of maintaining the visibility of a moving viewpoint with the aid of a complex describing all the visibility relationships in a scene, at once: the visibility complex. <br />We take two different approaches to this problem. In one, we show how one can maintain the so-called visibility polyhedron of the moving point, in an exact fashion. We apply a variant of the algorithm to the construction of the visibility complex of a set of disjoint polytopes. <br />The second approach is motivated by the need to render complex 3D scenes at interactive rates. We propose an algorithm that decomposes a 3D scene into simple cells related together with simple visibility relationships represented as a graph. The scene is drawn as the graph is traversed, allowing pruning to speed up rendering. We also provide an efficient algorithm for real-time rendering of hard shadows.

Mon travail de recherche a couvert ces dernières années un spectre assez large de modélisation, analyse numérique et simulation pour des problèmes physiques et biologiques, de la mécanique à l'échelle moléculaire ou particulaire, niveau dit " microscopique ", à la diffusion non-linéaire, niveau " macroscopique ", en passant par des équations cinétiques décrivant la distribution en vitesse de particules, niveau " mésoscopique ". Le point commun de ces travaux est l'étude de comportements asymptotiques et la recherche d'explications de phénomènes observables macroscopiques par des descriptions micro ou mésoscopiques à l'aide d'outils numériques. Les applications auxquelles on s'intéresse ici sont, pour la partie physique, liées à la thermodynamique couplée ou non avec du transfert radiatif ou une dynamique particulaire raréfiée et, pour la partie biologie-chimie, à des problèmes de propagation d'information par des mécanismes de transport ou de diffusion, ainsi qu'à la recherche de formation de motifs et à l'étude d'extinction de populations. Les équations aux dérivées partielles étudiées proviennent de modèles déterministes ou probabilistes et se classent dans les catégories de transport et de diffusion évolutifs. L'apparition, lors de l'adimensionnement des problèmes, de petits paramètres qui augmentent l'influence de certains des phénomènes caractéristiques dans la solution peut entraîner des difficultés importantes lors du traitement numérique, ce qui impose le recours à des solutions nouvelles permettant de recouvrer au minimum le comportement macroscopique prédit par les observations et par l'analyse mathématique.

La conversion photovoltaïque présente de nombreux avantages. Actuellement, les technologiesbasées sur l'élaboration de wafers de silicium cristallins dominent le marché, mais sont responsablesde pertes de matières importantes, très néfastes au coût de production des cellules. Le défi à releverest donc la réalisation de matériaux bas coûts en silicium par un procédé de métallurgie des poudres.Cependant, le frittage du silicium est dominé par des mécanismes de grossissement de grains quirendent la densification difficile par frittage naturel. Dans la littérature, l'identification de cesmécanismes est sujette à controverse. En particulier, le rôle de la couche d'oxyde natif (SiO2) à lasurface des particules de silicium reste inexploré. Dans ce manuscrit, l'influence de l'atmosphère surla réduction de cette couche de silice au cours du frittage est étudiée par analysethermogravimétrique. Les cinétiques de réduction sont en accord avec un modèle thermochimiqueprenant en compte, les quantités d'oxygène initialement présentes dans poudre, la pression partielleen espèces oxydantes autour de l'échantillon et l'évolution de la porosité du fritté. Pour la premièrefois, des données expérimentales permettent de montrer que la couche de silice inhibe legrossissement de grain. Des nouveaux procédés, basés sur un contrôle de l'atmosphère enmonoxyde de silicium (SiO(g)) autour de l'échantillon, sont alors proposés afin de maitriser la stabilitéde cette couche. Bien que la couche d'oxyde retarde les cinétiques de diffusion en volume, sonmaintien à des températures de 1300 - 1400 °C permet d'améliorer significativement la densification.Dans ces conditions, le comportement au frittage du silicium peut être séparé en deux étapes,clairement mises en évidences par la présence de deux pics de retrait sur les courbes de dilatométrie.Ce résultat est inhabituel compte tenu de l'aspect monophasé du matériau étudié. Cependant, il peutêtre expliqué à l'aide d'un modèle cinétique de frittage, basé sur des simplifications géométriques enaccord avec l'évolution microstructurale du matériau.

Cette thèse propose une étude mathématique du comportement des solutions autour d'états stationnaires pour des systèmes cinétiques gravitationnels de type Vlasov. Les trois premières parties présentent des résultats théoriques. Tout d'abord, par une ap- proche variationnelle, on construit des états stationnaires pour le système de Vlasov-Manev et on montre leur stabilité orbitale. Ensuite, on prouve l'existence de solutions autosimi- laires explosant en temps fini autour d'un état stationnaire pour le système dit de "Vlasov- Manev pur". Enfin on démontre la stabilité orbitale d'une large classe d'états stationnaires pour le système de Vlasov-Poisson relativiste. Ces résultats s'appuient sur de nouvelles méthodes utilisant la rigidité du flot. Celles-ci permettent notamment d'obtenir la séparation d'états stationnaires en évitant l'étude d'équations d'Euler-Lagrange non locales, de résoudre un problème variationnel avec une infinité de contraintes et de prouver la stabilité orbitale de solutions stationnaires non nécessairement obtenues de manière variationnelle. Dans la quatrième et dernière partie, nous étudions numériquement l'équation de Vlasov-Poisson en coordonnées radiales. Après avoir choisi un système de variables adéquates, nous présentons des schémas numériques de différences finies conservant la masse et le Hamiltonien du système.

Saccharothrix algeriensis est une bactérie filamenteuse capable de produire différentes molécules de la famille des dithiolopyrrolones présentant des activités antibiotiques et anticancéreuses. Ce travail porte sur l’analyse macroscopique et la modélisation du comportement dynamique de cette bactérie en vue de concevoir et d’optimiser un procédé de bioproduction de thiolutine. Les expériences menées en fermenteur de deux litres ont permis de décrire le comportement macroscopique de Sa. algeriensis. Sur substrat mixte (extrait de levure et glucose), la bactérie présente une croissance diauxique. La production de thiolutine, découplée de la croissance, est provoquée par l’épuisement du substrat préférentiel (extrait de levure). L’analyse des bilans carbone, azote et rédox nous a conduit à considérer deux phénomènes physico-chimiques qui n’étaient pas pris en compte dans le modèle proposé par Strub (2008) : l’entrainement de l’ammonium dans la phase gazeuse et l’adsorption des ions ammonium et des acides aminés sur la biomasse. Les données réconciliées sur la base des bilans élémentaires sont utilisées comme support pour la construction d’un modèle stœchio-cinétique. Différentes approches de modélisation des cinétiques et de représentation des bascules métaboliques ont été testées afin de retenir un schéma stœchio-cinétique apparent cohérent avec les données réconciliées. Cette étape de discrimination a été réalisée à l’aide de la toolbox ExOptim, développée en partenariat avec AgroParisTech. Cet outil, alliant une méthode d’optimisation robuste tout en limitant les temps de calcul, permet de prendre en compte les erreurs de mesure dès la mesure initiale et de reparamétrer indépendamment les paramètres du modèle et les données expérimentales. Les cinétiques réactionnelles sont représentées de manière satisfaisante par l’association du modèle de Monod et la loi logistique. Les seize paramètres du modèle sont identifiés sur sept expériences indépendantes. Le modèle finalement obtenu est ensuite testé sur cinq nouvelles expériences aux conditions initiales différentes. Ce modèle permet de prédire raisonnablement les cinétiques de croissance et de production de thiolutine. Toutefois, son utilisation reste limitée à des conditions initiales proches de celles de cette étude. / Saccharothrix algeriensis is a filamentous bacterium that produces several dithiolopyrrolone compounds with antibiotic and antitumor properties. This study focuses on the macroscopic analysis and modelling of the dynamic behaviour of this bacterium, in order to design and optimize the bioproduction process of thiolutin. Experiments were conducted in a two-litres bioreactor to characterise the macroscopic behaviour of Sa. algeriensis. On a medium containing several substrates (yeast extract and glucose), the bacterium shows a diauxic growth. Thiolutin production, which appears to be dissociated from growth, is triggered by the preferential substrate (yeast extract) depletion. The analysis of the carbon, nitrogen and redox balances led us to consider two physicochemical phenomena, which were not taken into consideration in the model proposed by Strub (2008): ammonium stripping in gaseous phase and adsorption of ammonium ions and amino acids on biomass. The data, reconciled on the basis of elementary balances, were used for the construction of a stoichio-kinetic model. Different modelling of kinetic and metabolic switches were tested in order to select a stoichio-kinetic model consistent with the reconciled data. This discrimination step was carried out using the ExOptim toolbox, developed in partnership with AgroParisTech. This toolbox combines a robust optimisation method and a reduced computation time. The experimental biases are taken into account from the first measurement and model parameters and experimental data reconditioned independently. The reaction kinetics are well represented by association of the Monod model and logistic laws. The sixteen model parameters were estimated from seven experiments. Finally, the model was validated by five additional experiments using different initial conditions. The kinetics of growth and thiolutin production are reasonably well predicted. However, the model is limited to initial conditions close to those of this study.

Le travail de thèse est préliminaire à une étude consacrée à l'élaboration et la validation d'une méthodologie permettant de prédire l'avancement des réactions de Maillard lors de la cuisson de la génoise. L'objectif du travail est de développer des outils expérimentaux pour suivre des marqueurs réactionnels et des variables de procédés de manière synchrone. Dans un premier temps, nous avons donc proposé un schéma réactionnel à la lumière de la connaissance expérimentale des composés volatils générés au cours de la cuisson de la génoise et des données de la littérature. Ensuite, un four pilote a été spécialement conçu et instrumenté. Deux dispositifs d'échantillonnage non perturbateur pour prélever les génoises et les vapeurs de cuisson en continu (au cours de la cuisson) ont été développés et validés. Le dispositif d'échantillonnage des vapeurs de cuisson est un système qui assure le piégeage des composés volatils libérés au cours de la cuisson de la génoise. Aprés avoir stabilisé l'environnement thermique de l'extraction et optimisé la méthode pour les marqueurs sélectionnés, nous avons examiné la répétabilité, la linéarité et la robustesse de la méthode et capacité d'une telle méthode à quantifier certains composés d'arôme libérés au cours de la cuisson de la génoise. À cet effet, les coefficients de proportionnalité entre les concentrations dans la phase gazeuse et dans la fibre SPME ont été déterminés. Nous nous sommes enfin intéressés à l'acquisition des données cinétiques des marqueurs avec les outils validés afin d'observer le schéma réactionnel sélectionné. Les essais de cuissons ont été menés selon un plan de cuisson à trois températures d'air, 140°C, 170°C et 200°C. La dégradation des réactifs et la génération des produits de réactions semblent être étroitement liés aux variables de procédés étudiées (température, teneur en eau, humidité de l'air). De plus le suivi de l'évolution des précurseurs au cours de la cuisson de la génoise montre que bien que la réaction de Maillard reste prépondérante, d'autres réactions sont impliquées dans leurs dégradations en fin de cuisson, probablement la réaction de caramélisation et l'isomérisation des sucres. / The thesis is a preliminary study on the development and validation of a methodology to predict the progress of Maillard reactions during the baking of sponge cake. The aim of this work is to develop experimental tools to track markers and reaction process variables synchronously. Initially, we proposed a reaction scheme based on experimental knowledge of the volatile compounds generated during sponge cake baking and literature data. Then, a pilot oven was designed and instrumented. Two non-invasive sampling devices to collect the sponge cakes and baking vapors continuously (during baking) have been developed and validated. The volatiles sampling device is a system that ensures the trapping of volatile compounds released during the cake baking. After we have stabilized the thermal environment of the extraction and optimized the method for the selected markers, we examined the repeatability, linearity and robustness of the method and capacity of such a method to quantify some aroma compounds released during the cake baking. Fot this reason, the proportionality coefficients between the concentrations in the gas phase and the SPME fiber were determined. Finally, we were interested in the acquisition of kinetic data markers with validated tools to observe the selected reaction scheme. The baking experiments were conducted using at three air temperatures, 140 ° C, 170 ° C and 200 ° C. The degradation of the reactants and the generation of reaction products appear to be closely related to studied process variables (temperature, moisture, humidity). Also tracking the evolution of precursors during the baking of the sponge cakes shows that although the Maillard reaction is still predominant, other reactions are involved in their degradation after baking, probably the caramelization and the isomerization of sugars reactions.

La Nouvelle-Calédonie détient d’importantes réserves de nickel latéritique et est devenue, en 2017, le cinquième producteur mondial de Ni. Ces dépôts sont habituellement considérés comme résultant d’altération latéritique intense de la péridotite, qui constitue la principale source de nickel. Ainsi, le principal modèle conceptuel de la formation des minerais de nickel latéritique en Nouvelle-Calédonie est un modèle per descensum où la plupart des éléments (Mg, Ni et Si) ont été lessivés depuis la surface, en particulier lors du développement du sol latéritique. Le nickel est ensuite reprécipité, soit dans la goethite de la latérite fine, soit au niveau de la saprolite, sous forme de goethite et de silicates Mg-Ni, dont des talc-like ou kérolite. Les observations minéralogiques et structurales récentes ainsi que les données minières ont cependant mis en évidence de nombreux types d’hétérogénéités dans les concentrations, et la distribution des porteurs de Ni. Comprendre les facteurs la mobilité de cet élément, ses mécanismes de piégeage ainsi que les paramètres chimiques et hydrodynamiques à l’origine de ce piégeage, est essentiel afin de prévoir la distribution du nickel dans les profils latéritiques en Nouvelle-Calédonie, et constitue l’objectif de cette thèse. Ce travail est basé sur le développement (i) d’un modèle 1D s’intéressant en particulier au comportement géochimique du nickel lors de l’altération de l’ophiolite, sa comparaison avec les observations in situ et une compréhension détaillée de la mobilité des éléments traces pendant le processus, et (ii) d’un modèle 2D hydro-géochimique couplé avec l’hydrodynamique complexe des profils latéritiques, améliorant ainsi la connaissance du contrôle structural sur la redistribution et la minéralisation du nickel. Tandis que les simulations 1D permettent de mieux comprendre les aspects chimiques contrôlant les processus de rétention du nickel au sein d’un profil, le modèle 2D se révèle être un outil puissant pour la compréhension de la formation des dépôts locaux les plus riches en nickel. Les résultats du modèle 2D montrent une remobilisation du nickel depuis les horizons supérieurs puis sa reprécipitation sous forme de silicates dans la saprolite. Le nickel remobilisé provient principalement de la zone saprolitique à cause de la dissolution des silicates de nickel formés précédemment ainsi que de l’olivine résiduelle de cette zone. Ce modèle a également révélé que l’horizon latéritique (et en particulier les oxy- hydroxydes de nickel) avait un faible impact dans la remobilisation du nickel. L’infiltration latérale de l’eau contenant le nickel dissout issu des formations surincombantes est à l’origine de la formation des zones les plus riches dans les parties inférieures du profil. Cette redistribution est entièrement contrôlée par l’hydrodynamique locale, la topographie ainsi que l’orientation et la position des fractures. Les modèles présentés permettent d’expliquer les processus de formation des minerais de nickel latéritique saprolitique, améliorant ainsi la compréhension des paramètres contrôlant la mobilité des éléments traces dans un environnement ultramafique. Ceci donne une nouvelle clé de distribution du nickel dans les profils actuels, qui peut devenir un outil pour la prospection minière, et la recherche de nouvelles ressources exploitables / New Caledonia hosts significant lateritic nickel reserves, and presently became the fifth largest Ni producer in the world. These deposits are generally thought to be closely as- sociated with the intense chemical and mechanical weathering of peridotite bedrock that is a principal source of nickel. Thus, the main ore genesis model for Ni ores in New Caledonia is based on a single per descensum model where most elements (Mg, Ni, and Si) are leached from the surface, particularly, during lateritic soil development. Nickel is then concentrated either in the fine-grained laterite where goethite is the main Ni bearer, the so-called ’lateritic ore’, or below the laterites, in the saprolite level, where nickel occurs as goethite and several types of Mg-Ni silicates, in particular kerolite. Recent mineralogical and structural observations together with mining data have revealed a lot of different types of heterogeneities associated with the distribution and mineralogy of Ni bearing minerals. Therefore, in depth investigations of Ni mobility, its retardation processes along with its governing chemical and hydrodynamic parameters are of big importance for understanding and subsequent prediction of Ni distribution in profiles of New Caledonia. Such an investigation is an objective of the present work. The concept is based on the development of i) a powerful 1D model with particular emphasis on Ni geochemical behavior during ophiolite weathering, its comparison with in situ observations, and detailed understanding of trace elements mobility, and ii) 2D hydro-geochemical model coupled with complex hydrodynamics, that would additionally provide new insight into the structural control on Ni redistribution and mineralization. While the 1D simulations provide a remarkable result for understanding the chemical features that drive Ni retention processes in a profile, 2D model appears to be a powerful tool for understanding how local Ni-enrichments may form. The results of this model show the reactivation of Ni from upper horizons and its concentration in neo-formed silicates in bottom of the saprolite. The reactivated Ni comes mostly from the saprolite horizon due to the redissolution of previously formed Ni-bearing silicates and still persisting in this olivine zone. Modeling has revealed minor contribution of the laterite horizon (Ni-oxi-hydroxides) into the Ni remobilization. The lateral infiltration of water with remobilized Ni from areas such as topographic highs to downstream slope areas leads to the formation of richest deposits in this lower part of profile. The manner of redistribution is fully governed by the topographic slopes, orientation and position of the fractures. Presented models appear to be of importance in attempt of explanation of Ni mineralization processes, revealing the main keys to understanding the control of trace elements mobility in ultramafic environment. The latter gives new insights into the Ni distribution in present day profiles and, therefore, may greatly help in mineral prospecting and forecasting the distribution of future resources