Ce travail a pour objectif de développer et d'appliquer une méthode originale permettant de simuler l'écoulement dans un milieu poreux fracturé. Cette méthode repose sur une approche multicouches double continuum permettant de séparer le comportement des différents aquifères présents sur un site. La résolution des écoulements, basée sur la méthode des Eléments Finis de Crouzeix-Raviart, est associée à une méthode inverse (minimisation de type Quasi-Newton combinée à la méthode de l'état adjoint) et à une paramétrisation multi-échelle.La méthode est appliquée dans un premier temps sur l'aquifère fracturé du site expérimental de Poitiers. Les résultats montrent une bonne restitution du comportement de l'aquifère et aboutissent à des champs de transmissivité plus réguliers par rapport à ceux de l'approche simple continuum. L'application finale est réalisée sur le site de Cadarache (taille plus importante et données d'entrée moins denses). Le calage des deux aquifères présents sur le site est satisfaisant et montre que ceux-ci se comportent globalement de façon indépendante. Ce calage pourra être amélioré localement grâce à données de recharge plus fines.

Le concept de radio logicielle (SDR) est une solution pertinente pour concevoir des équipements multi-standards. Une façon de réaliser de tels équipements est d'identifier les fonctions et opérateurs communs entre les standards. Cette approche s'appelle la paramétrisation et est divisée en deux catégories : l'approche pragmatique qui est une version pratique pour créer et développer des opérateurs communs à partir d'opérateurs existants, et l'approche théorique dont l'objectif est de réaliser une exploration graphique d'un équipement multi-standards selon différents niveaux de granularité, accompagnée d'un problème d'optimisation. C'est cette dernière approche qui a constitué le sujet de base de cette thèse. Ainsi, une fonction de coût doit être optimisée afin de sélectionner les opérateurs communs entre les différentes normes, ce qui permet de proposer une configuration optimale à partir de laquelle sont déduits les opérateurs communs. Dans notre travail, nous avons dans un premier temps modélisé théoriquement la structure graphique d'un système multi-standards par un hypergraphe orienté. En outre, nous avons fourni une expression mathématique alternative de la fonction de coût suggérée, en utilisant des définitions propres à la théorie des graphes. Ensuite, nous avons montré que le problème d'optimisation associé était un problème NP sous une certaine contrainte, ce qui a entraîné une preuve d'exclusion de certaines configurations dont les coûts ne peuvent être minimaux. Ceci a constitué la deuxième contribution de cette thèse. Enfin, nous avons proposé un nouvel algorithme permettant de résoudre le problème d'optimisation donné, et dont l'intérêt est de donner une solution optimale du problème au lieu d'une solution approchée fournie par les méthodes heuristiques classiques. Un programme associé à cet algorithme a été développé en langage C, puis appliqué à plusieurs exemples de cas génériques afin d'en étudier les performances.

Dans cette thèse, nous proposons des solutions pour la conception de systèmes de contrôle actif de vibration robustes (AVC). Le manuscrit de thèse comporte deux grandes parties.Dans la première, les problèmes d'incertitude paramétrique dans les systèmes de contrôle actif de vibration sont étudiés. En plus des incertitudes sur la fréquence des perturbations, nous avons trouvé que la présence de zéros complexes peu amortis soulevait des problèmes de conception difficiles, même pour des systèmes et des modèles parfaitement connus. Dans ce contexte, nous avons proposé des solutions pour le problème linéaire. Une procédure améliorée d'identification en boucle fermée a été développée pour réduire les incertitudes dans l'identification de ces zéros. Pour traiter les incertitudes sur la perturbation, l'adaptation de la fréquence est de toute façon incontournable.La seconde partie est consacrée au développement et/ou à l'amélioration de deux algorithmes, désormais classiques, de compensation par feedback adaptatif direct, fondés sur la paramétrisation de Youla-Kučera. Le premier résulte de l'amélioration d'un précédent travail (Landau et al., 2005) ; les contributions concernent la synthèse du contrôleur central robuste et l'utilisation optionnelle de la surparamétrisation du filtre Q-FIR (réponse à temps fini) avec pour effet de minimiser l'effet « waterbed » sur la fonction de sensibilité de sortie. Le second algorithme présente une structure hybride directe/indirecte qui utilise un filtre Q-IIR (à temps de réponse infini). Les améliorations sont dues principalement au dénominateur du filtre, obtenu à partir d'une estimation de la perturbation. Cette solution permet également de simplifier la conception du contrôleur central.Les algorithmes ont été testés, comparés et validés sur un procédé réel du laboratoire Gipsa-lab, dans le cadre d'un benchmark international. / In this thesis, solutions for the design of robust Active Vibration Control (AVC) systems are presented. The thesis report is composed of two main parts.In the first part of the thesis uncertainties issues in Active Vibration Control systems are examined. In addition of the uncertainties on the frequency of the disturbances it has been found that the presence of low damped complex zeros raise difficult design problems even if plant and models are perfectly known. Solutions for the linear control in this context have been proposed. In order to reduce the uncertainties in the identification of low complex zeros and improved closed loop identification procedure has been developed. To handle the uncertainties on the disturbance frequency adaptation has any way to be used.The second part is concerned with the further development and/or the improvement of the now classical direct adaptive feedback compensation algorithms using Youla Kucera controller parametrization. Two new solutions have been proposed in this context. The first one results from the improvement of a previous work (Landau et al., 2005). The contributions are a new robust central controller design to the optional use of over parameterization of the Q-FIR filter which aims to ensure a small waterbed effect for the output sensitivity function and therefore reducing the unwanted amplification. The second algorithm presents a mixed direct/indirect structure which uses a Q-IIR filter. The improvements are mainly the effect of the Q filter denominator, which is obtained from a disturbance identification. This solution in addition drastically simplifies the design of the central controller.The algorithms have been tested, compared and validated on an international benchmark setup available at the Control System Department of GIPSA-Lab, Grenoble, France.

Ce travail porte sur le développement d'algorithmes et l'application de méthodes numériques pour la résolution des problèmes inverses d'estimation de paramètres, d'identification de conditions aux limites et d'identification de sources dans un milieu poreux. Ces outils seront d'une grande utilité pour l'aide à la gestion des ressources en eaux souterraines et à leur préservation quant aux dégradations. L'objectif de cette thèse est de résoudre ces problèmes inverses en se basant sur différentes approches : Une résolution basée sur l'optimisation de forme topologique qui est la recherche de la géométrie d'un objet qui soit optimale vis à vis d'un critère donné, et ce sans aucun a priori sur sa topologie, c'est-à-dire sur le nombre de "trous" qu'il peut contenir. Sachant que ces trous représentent les puits recherchés. Pour ce faire, nous avons adopté la méthode du gradient topologique, qui consiste à étudier le comportement d'une fonction objectif lors de la création d'un petit trou à l'intérieur du domaine. Une résolution basée sur la minimisation d'une fonctionnelle d'erreur en énergie en utilisant des données surabondantes sur une partie de la frontière du domaine afin de compléter les données sur toute la frontière du domaine et de déterminer les positions, les débits et le nombre de puits existants à l'intérieur du domaine. Une résolution par le couplage de la méthode de paramétrisation adaptative qui a l'avantage de minimiser le nombre des inconnus de paramètres permettant d'interpréter au mieux les données disponibles et la méthode du gradient topologique. Ce couplage nous permet à la fois d'identifier les zones géologiques, de déterminer les valeurs de la transmissivité hydraulique dans chaque zone et de localiser les positions des puits. / This work deals with the development of algorithms and application of numerical methods for solving inverse problems of parameters estimation, identification of boundary conditions and localisation of sources in porous media. These tools will be usefull in the management of groundwater resources and their preservation as to damage. The objective of this thesis is to solve the inverse problem based on different approaches: A resolution based on topological shape optimization is to find an optimal design without any priori assumption about its topology, that is, about the number of holes it may contain. Knowing that these holes represent the searched wells. To do this, we have adopted the method of topological gradient, which is to study the behavior of an objective function when creating a small hole inside the domain. A resolution based on the minimization of a constitutive law gap functional by using overspecified data on a part of the boundary of the domain to complete the data on all the boundary of the domain and determine the positions, the flows and the number of existing wells inside the domain. A resolution by the coupling of the adaptive parameterization method which has the advantage to minimize the number of the unknowns of parameters allowing to interpret at best the available data and the method of the topological gradient. This coupling allows us at the same time to identify the geological zones, to determine the values of the hydraulic transmissivity in every zone and to locate wells' positions.

Dans cette thèse, nous présentons des solutions pour la conception des systèmes de contrôle actif de vibrations. Dans la première partie, des méthodes de contrôle par action anticipatrice (feedforward) sont développées. Celles-ci sont dédiées à la suppression des perturbations bande large en utilisant une image de la perturbation mesurée par un deuxième capteur, en amont de la variable de performance à minimiser. Les algorithmes présentés dans cette mémoire sont conçus pour réaliser de bonnes performances et maintenir la stabilité du système en présence du couplage positif interne qui apparaît entre le signal de commande et l'image de la perturbation. Les principales contributions de cette partie sont l'assouplissement de la condition de « Stricte Positivité Réelle » (SPR) par l'utilisation des algorithmes d'adaptation « Intégrale + Proportionnelle » et le développement de compensateurs à action anticipatrice (feedforward) sur la base de la paramétrisation Youla-Kučera. La deuxième partie de la thèse concerne le rejet des perturbations bande étroite par contre-réaction adaptative (feedback). Une méthode d'adaptation indirecte est proposée pour le rejet de plusieurs perturbations bande étroite en utilisant des filtres Stop-bande et la paramétrisation Youla-Kučera. Cette méthode utilise des Filtres Adaptatifs à Encoche en cascade pour estimer les fréquences de perturbations sinusoïdales puis des Filtres Stop-bande pour introduire des atténuations aux fréquences estimées. Les algorithmes sont vérifiés et validés sur un dispositif expérimental disponible au sein du département Automatique du laboratoire GIPSA-Lab de Grenoble. / In this thesis, solutions for the design of robust Active Vibration Control (AVC) systems are presented. The thesis report is composed of two parts. In the first one, feedforward adaptive methods are developed. They are dedicated to the suppression of large band disturbances and use a measurement, correlated with the disturbance, obtained upstream from the performance variable by the use of a second transducer. The algorithms presented in this thesis are designed to achieve good performances and to maintain system stability in the presence of the internal feedback coupling which appears between the control signal and the image of the disturbance. The main contributions in this part are the relaxation of the Strictly Positive Real (SPR) condition appearing in the stability analysis of the algorithms by use of “Integral + Proportional” adaptation algorithms and the development of feedforward compensators for noise or vibration reduction based on the Youla-Kučera parameterization. The second part of this thesis is concerned with the negative feedback rejection of narrow band disturbances. An indirect adaptation method for the rejection of multiple narrow band disturbances using Band-Stop Filters (BSF) and the Youla-Kučera parameterization is presented. This method uses cascaded Adaptive Notch Filters (ANF) to estimate the frequencies of the disturbances' sinusoids and then, Band-stop Filters are used to shape the output sensitivity function independently, reducing the effect of each narrow band signal in the disturbance. The algorithms are verified and validated on an experimental setup available at the Control Systems Department of GIPSA-Lab, Grenoble, France.

Depuis les années 80 la cartographie optique est devenu un outil important pour l'étude et la compréhension des arythmies cardiaques. Cette expérience permet la visualisation de flux de fluorescence à la surface du tissu ; fluorescence qui est directement liée au potentiel transmembranaire. Dans les observations en surface se cachent des informations sur la distribution en trois dimensions de ce potentiel. Nous souhaitons exploiter ces informations surfaciques afin de reconstruire le front de dépolarisation dans l'épaisseur. Pour cela nous avons développé une méthode basée sur la résolution d'un problème inverse. Le modèle direct est composée de deux équations de diffusion et d'une paramétrisation du front de dépolarisation. La résolution du problème inverse permet l'identification des caractéristiques du front. La méthode a été testée sur des données in silico avec différentes manières de caractériser le front (sphère qui croît au cours du temps, équation eikonale). Les résultats obtenus sont très satisfaisants et comparés à une méthode développée par Khait et al. [1]. Le passage à l'étude sur données expérimentales a mis en évidence un problème au niveau du modèle. Nous détaillons ici les pistes explorées pour améliorer le modèle : illumination constante, paramètres optiques, précision de l'approximation de diffusion. Plusieurs problèmes inverses sont considérés dans ce manuscrit, ce qui implique plusieurs fonctionnelles à minimiser et plusieurs gradients associés. Pour chaque cas, le calcul du gradient est explicité, le plus souvent par la méthode de l'adjoint. La méthode développée a aussi été appliquée à des données autres que la cartographie optique cardiaque. / Since the 80's optical mapping has become an important tool for the study and the understanding of cardiac arythmias. This experiment allows the visualization of fluorescence fluxes through tissue surface. The fluorescence is directly related to the transmembrane potential. Information about its three-dimension distribution is hidden in the data on the surfaces. Our aim is to exploit this surface measurements to reconstruct the depolarization front in the thickness. For that purpose we developed a method based on the resolution of an inverse problem. The forward problem is made of two diffusion equations and the parametrization of the wavefront. The inverse problem resolution enables the identification of the front characteristics. The method has been tested on in silico data with different ways to parameter the front (expanding sphere, eikonal equation). The obtained results are very satisfying, and compared to a method derived by Khait et al. [1]. Moving to experimental data put in light an incoherence in the model. We detail the possible causes we explored to improve the model : constant illumination, optical parameters, accuracy of the diffusion approximation. Several inverse problems are considered in this manuscript, that involves several cost functions and associated gradients. For each case, the calculation of the gradient is explicit, often with the gradient method. The presented method was also applied on data other than cardiac optical mapping.

Durant l’épisode de canicule-sécheresse de 2003, les peuplements de douglas et de hêtres en Bourgogne ont été lourdement affectés, présentant des symptômes de dépérissement et de surmortalité. Cet épisode semble être la première occurrence d’aléas climatiques attendus dans un futur proche et remet en question leur pérennité en Bourgogne puisque leur vulnérabilité au climat est attribuable à l'amplitude et au cumul des contraintes hydriques exercées durant leur cycle de végétation. Dans le contexte du changement climatique et en réponses aux demandes des gestionnaires forestiers qui s’appuient partiellement sur une cartographie de l’évolution des contraintes climatiques jusqu’à la fin de ce siècle, ce travail explore la capacité d’une chaîne de modélisation hydroclimatique haute résolution couplant le modèle de climat régional WRF alimenté par les réanalyses ERA-Interim au modèle de bilan hydrique Biljou© ˆ simuler des indices de déficit hydrique pour ces deux essences. La première partie de ce travail propose une analyse de la capacité du modèle WRF à simuler chacune des variables atmosphériques de surface qui sont utilisées en entrée du modèle du bilan hydrique. L’analyse de la capacité du modèle à simuler ces variables repose (i) sur une approche comparative directe entre les données simulées par WRF et les observations enregistrées par le réseau de stations Météo-France et les réanalyses SAFRAN à l’échelle de la région, de la station, et du peuplement forestier, (ii) sur une approche indirecte utilisant l’évapotranspiration potentielle (ETP) et la relation entre les indices de croissance radiale et les indices de déficit hydrique calculés par le modèle d’impact pour les deux essences. Les résultats montrent une amélioration significative des données ERA-Interim par le modèle WRF pour chacune des variables ainsi qu’une capacité certaine à les spatialiser à haute résolution. Toutefois, la bonne reproduction de l’ETP par WRF, combinée à la faible corrélation entre la moyenne annuelle des indices de déficit hydrique estimés avec les données WRF et la moyenne annuelle des indices de croissance radiale montrent que les difficultés de WRF à simuler le déficit hydrique sont principalement imputables à ses biais de précipitations. La seconde partie propose l’application d’une post-correction statistique aux données de précipitations WRF. Bien que cette méthode améliore significativement la distribution spatiale des précipitations, leurs variabilités saisonnière et interannuelle et surtout les cumuls précipités, les données post-corrigées ne permettent pas de reproduire un indice de déficit hydrique suffisamment proche de celui estimé à partir des observations ou des analyses SAFRAN. Deux nouvelles simulations résolvant explicitement les processus convectifs et utilisant un guidage spectral ont permis de montrer à partir de deux années types que cette déficience est imputable à l’incapacité de la méthode de correction à résoudre les différences de timing de la variabilité climatique transitoire simulée par WRF. Deux types d’erreurs de modélisation climatique, survenant indépendamment, sont donc d'une importance primordiale pour les études d'impact: (i) la chronologie des événements pluvieux ; (ii) la distribution statistique des précipitations quotidiennes. La combinaison de ces deux éléments contrôle le nombre de jours franchissant le seuil de 40% de réserve relative en eau du sol et indirectement l’intensité des indices de déficit hydrique. / During the 2003 drought and heat wave event, douglas-fir and common beech stands in Burgundy have been heavily affected, and presented symptoms of dieback and mortality. This event seems to be the first occurrence of expected climatic changes in the near future and questions their sustainability in Burgundy since their climate vulnerability is mainly due to the amplitude and accumulated water constraints exercised during their growing cycle. In the context of climate change and in order to provide information to forest managers who partly rely on a mapping of the climatic constraints until the end of this century, this work explores the ability of a high resolution hydroclimatic modelling chain, coupling the regional climate model WRF to the daily lumped water balance model Biljou© in order to simulate soil water deficit indices for these two species. The first part of this paper analyzes the capacity of WRF model to simulate each surface atmospheric variable used as input for the water balance computation. The analysis of model's ability to simulate these variables is based on (i) a direct and comparative approach between WRF simulated data and observations recorded by the Météo-France stations network and SAFRAN reanalyses across the whole region, over stations and forest stands, (ii) on an indirect approach using the potential evapotranspiration and soil water deficit index calculated by Biljou©. Results show a significant improvement upon the ERA-Interim data for each variable and a strong ability to produce reliable data at high resolution. However, the WRF capability to estimate a realist potential evapotranspiration, combined to the the low correlation between the average annual soil water deficit and radial growth indexes, show that the WRF deficiencies in simulating water deficit are mainly attributable to its precipitation biases. The second part proposes to apply a statistical post-correction to the WRF precipitation data. Although this method significantly improves the spatial distribution of precipitation, their seasonal and interannual variability and precipitation amounts, post-corrected data do not produce a water deficit index sufficiently close to those ones estimated from observations or SAFRAN reanalysis. Two new simulations explicitly solving convective processes and using a spectral nudging have shown that this deficiency is mainly attributable to the inability of the correction method to solve timing differences of the transient climate variability simulated by WRF. This work showed that two types of climate modeling errors occurring independently, are major issues for impact studies: (i) the timing of precipitations events ; (ii) the statistical distribution of daily precipitation. Combined together, they control the number of days crossing the 40% threshold of relative extractable water and indirectly the soil water deficit index intensity.

Ce travail a consisté à développer un modèle de canopée (CIM), qui pourrait servir d’interface entre des modèles méso-échelles de calcul du climat urbain et des modèles micro-échelles de besoin énergétique du bâtiment. Le développement est présenté en conditions atmosphériques variées, avec et sans obstacles, en s’appuyant sur les théories précédemment proposées. Il a été, par exemple, montré que, pour être en cohérence avec la théorie de similitude de Monin-Obukhov, un terme correctif devait être rajouté au terme de flottabilité de la T.K.E. CIM a aussi été couplé au modèle méso-échelle WRF. Une méthodologie a été proposée pour profiter de leurs avantages respectifs (un plus résolu, l’autre intégrant des termes de transports horizontaux) et pour assurer la cohérence de leurs résultats. Ces derniers ont montré que ce système, en plus d’être plus précis que le modèle WRF à la même résolution, permettait, par l’intermédiaire de CIM, de fournir des profils plus résolus près de la surface. / This study consisted in the development of a canopy model (CIM), which could be use as an interface between meso-scale models used to simulate urban climate and micro-scale models used to evaluate building energy use. The development is based on previously proposed theories and is presented in different atmospheric conditions, with and without obstable. It has been shown, for example, that to be in coherence with the Monin-Obukhov Similarity Theory, that a correction term has to be added to the buoyancy term of the T.K.E. CIM has also been coupled with the meteorological meso-scale model WRF. A methodology was proposed to take advantage of both models (one being more resolved, the other one integrating horizontal transport terms) and to ensure a coherence of the results. Besides being more precise than the WRF model at the same resolution, this system allows, through CIM, to provide high resolved vertical profiles near the surface.

Les bâtiments représentent 40 pourcents de la consommation finale d'énergie. Ils sont ainsi le fer de lance des politiques de réduction des dépenses énergétiques. Récemment, des systèmes de modèles climatiques qui incluent un modèle atmosphérique régional et des paramétrisations urbaines sophistiquées ont été développés. Ils considèrent la complexité de l’îlot de chaleur urbain et ses interactions avec les besoins énergétiques des bâtiments. Dans quelle mesure ces systèmes constituent-ils des outils d’aide à la décision pour les autorités locales ? Cette étude menée sur le territoire de l'Eurodistrict (Strasbourg - Kehl) en 2010, puis en 2030, à l’aide du système de modèles de climat WRF/ARW-BEP+BEM a démontré que si le système de modèles estimait de manière fiable les besoins en chauffage des bâtiments, ces derniers étaient davantage sensibles aux caractéristiques intrinsèques des bâtiments qu'aux formes urbaines et à l'îlot de chaleur urbain induit par ces formes. / Buildings represent 40 percent of the end-use energy. Thus, they constitute a key point of the energy saving policies. Recently, climate modeling systems that include a mesoscale atmospheric model, sophisticated urban parameterizations have been developed to account for the complexity of the urban climate and its interactions with the building energy loads. This study aims to assess the capability of such climate modeling systems to provide climate and energy guidelines to urban planners. For this, we used the research collaborative WRF/ARW-BEP+BEM climate modeling system and performed sensitivity tests considering the territory of the Eurodistrict in 2010, and then in 2030. The results reveal that the climate modeling system achieves estimating the building energy needs over the study area, but also indicate that the building energy needs are more sensitive to the building intrinsic properties and occupant behavior than to the urban forms and their induced urban heat island.

En viabilité hivernale, la prévision de l’état de surface des infrastructures s’avère indispensable, et permet une anticipation, une meilleure coordination et une efficacité d’intervention des services d’exploitation. La majorité des pays dispose de modèles de prévision de la température de surface d'infrastructures et des routes en particulier (TSR). La complexité de ces outils d’aide à la décision est croissante, pour servir au mieux les usagers et l’exploitant. Le microclimat urbain influence le bilan énergétique de surface selon différents processus : radiatifs, aérodynamiques et hydrologiques. Néanmoins, d’autres processus physiques anthropiques influencent cette TSR, tel que le trafic. Des travaux ont été menés par le passé concernant l’apport énergétique du trafic dans le bilan thermique de la ville. Celui-ci a fait l’objet d’études sur les périodes estivales et les îlots de chaleur urbains associés. Cependant, dans les cas de dégradations des conditions hivernales de circulation, ces apports énergétiques ont été intégrés de façon marginale dans la modélisation des paramètres de surface de la route. L’absence de cette contribution du trafic dans la modélisation du bilan énergétique de surface explique, dans une certaine limite, la prévision imparfaite de l’état de surface de la route. La bibliographie recense plusieurs études conduites afin d’identifier et de quantifier ces effets du trafic. Elles n'ont pas ou peu traité la perte ou le gain d'énergie causé par le passage des véhicules sur le bilan énergétique de surface, ou sur la modélisation de la TSR. Dans la présente étude, deux approches ont été proposées pour paramétrer le trafic dans le modèle numérique Town Energy Balance (TEB), l'une globale et la seconde détaillée. Leur analyse comparée indique que la seconde a significativement amélioré les résultats de la modélisation de la TSR. Les apports thermiques du trafic ont augmenté la TSR de 2 à 4°C pour la rapprocher des mesures expérimentales (écart de 0.5 à 1°C en moyenne). Elle est le résultat de l’effet cumulatif des différents processus physiques associés au trafic, et varie en fonction de ses paramètres (densité, vitesse de circulation, fluidité, etc.). Une étude de sensibilité a été menée afin d’apprécier les processus physiques responsables de l’amélioration de la TSR. Les résultats indiquent que l’effet turbulent, la densité de flux radiatif ainsi que la densité de flux due aux frottements des pneumatiques contribuent le plus à l'augmentation la TSR. Néanmoins la contribution énergétique de chacun de ces processus dépend à la fois des paramètres du trafic et des conditions météorologiques. Cette étude a présenté également une analyse préliminaire de l’influence de la lame d’eau présente en surface sur la TSR. Cette dernière décroit en fonction de l’épaisseur de la lame d’eau. Les facteurs anthropiques tels que le trafic et l’adjonction de fondants routiers sur la lame d’eau présente en surface sont décrits et discutés, et une paramétrisation proposée en perspectives / In winter maintenance, forecasting the infrastructure surface status is mandatory, to allow anticipation, better coordination and efficiency of services. The majority of countries have forecast models of the infrastructure surface temperature and especially roads one (RST). The complexity of these decision tools is growing, to best serve the users and managers. The urban microclimate influences the surface energy balance according to various processes: radiative, aerodynamic and hydrologic. However, other anthropogenic physical processes influence this RST, such as traffic. Studies have been carried out in the past about the traffic heat input in the town heat balance. These were conducted on the summer periods and associated urban heat islands. However, in case of traffic in adverse winter conditions, these energy contributions were marginally integrated into the modeling of the road surface parameters. The absence of this traffic's contribution in the surface energy balance modeling explains, to a given limit, the imperfect forecasting of road surface status. The literature identifies several studies conducted to identify and to quantify these effects of traffic. They have insufficiently or not treated the loss or gain on energy caused by circulating vehicles on the surface energy balance, or on modeling the RST. In this study, two approaches have been proposed to parameterize the traffic in the Town Energy Balance (TEB) numerical model, this first one being overall and the second one detailed. Their comparative analysis indicates that the second significantly improved the results of the RST modeling. The traffic heat inputs increased RST by 2 to 4°C, results being then closer to experimental measurements (average difference of 0.5 to 1°C). It is the result of the cumulative effect of the various traffic physical processes, and varies according to its parameters (density, vehicle velocity, fluidity, etc.). A sensitivity analysis was conducted to assess the physical processes responsible for the improvement of the RST. The results indicate that the turbulent effect, the radiative heat flux and flux density due to tires friction represent the greatest contribution to RST increase. Nevertheless the energy contribution of each of these processes depends both traffic parameters and weather conditions. This study also presented a preliminary analysis of the influence of a water layer over the surface on the RST. The latter decreases as a function of the thickness of the water layer. Anthropogenic factors such as traffic and the addition of de-icing products into the water layer present on the surface are described and discussed, and a parameterization proposed as a perspective