Ce travail de recherche vise à caractériser les propriétés hydrodynamiques des sols des Cévennes, une zone de moyenne montagne du sud de la France connue pour être sujette à des crues à cinétique rapide, appelées « crues éclair ». Il s'agit donc de déterminer les paramètres relatifs à la conductivité hydraulique et à la rétention en eau du modèle de Mualem-van Genuchten, leur variabilité spatiale à l'échelle du versant ou du petit bassin, et les facteurs qui expliquent cette variabilité spatiale. Il s'agit également d'évaluer la perméabilité du sous-sol constitué d'une roche plus ou moins altérée, et d'en déduire des informations sur les processus de saturation des sols pouvant être à l'origine de ruissellement.Les propriétés hydrodynamiques des sols et du sous-sol ont été estimées par modélisation inverse de la teneur en eau sols, à partir d'une station pilote sur le site de Sumène et d'un réseau de 50 stations de mesure répartis sur 4 versants granitiques (site de Valescure) et 2 versants schisteux (site de Tourgueille). La station pilote comprend 3 sondes de mesure installées à 20, 40, 60 cm de profondeur, et les stations installées sur les versants comprennent chacune 2 sondes de mesure installées à des profondeurs variables allant de 15 à 45 cm. Les teneurs en eau du sol ainsi que les précipitations ont été mesurées en continu au pas de temps 15 min sur les différents sites. La procédure de modélisation inverse est basée sur l'algorithme génétique multi-objectif NSGA-II. Cette procédure a été appliquée pour chaque station, en considérant que le sol d'épaisseur variable constitué de la superposition de 2 horizons, surmontant un troisième horizon figurant la roche altérée sous-jacente. Au total, quinze paramètres doivent être calibrés pour chaque station, et sont estimés à partir de périodes sélectionnées correspondant aux épisodes de pluie les plus importants. Les processus d'évapotranspiration ont été considérés comme négligeables au cours de ces épisodes. L'analyse des propriétés de rétention et de conductivité obtenues a mis en évidence la très forte conductivité hydraulique des sols étudiés, de l'ordre de 1000 à 2000 mm/h, et leur faible rétention en eau. Les teneurs en eau à saturation estimées s'échelonnent de 0.30 à 0.60 cm3.cm-3, traduisant une grande hétérogénéité locale. On retrouve cette forte variabilité pour l'estimation des profondeurs de sol qui varient de 31 à 120 cm. Pour l'horizon représentant le sous-sol, les perméabilités sont très contrastées, avec des valeurs de conductivités hydrauliques variant de quelques unités à plus d'une centaine de mm/h. Les stations présentant les sols les moins profondes et des sous-sols moins perméables génèrent des zones saturées qui se développent à la base du profil de sol, comme ça a été observé et simulé lors des épisodes pluvieux les plus intenses. Si les propriétés des sols apparaissent très variables d'une station à une autre et d'un versant à un autre, les résultats n'ont pas montré d'organisation particulière de cette variabilité spatiale. La position de la station sur le versant n'apparaît pas être déterminante pour les porosités, les profondeurs, les perméabilités et les rétentions en eau des sols ; les distributions de ces propriétés, calculées à l'échelle des différents versants, ne semblent pas reliées aux grands traits du paysage : géologie, orientation, végétation. / This research work aims to characterize soil hydrodynamic properties of the Cévennes area, a low mountain region known to be prone to flash floods. We thus chose to estimate hydraulic conductivity and water retention properties of the Mualem –van Genuchten model and their spatial variability at the hillslope scale and between multiple hillslope profiles. The objectives are also to evaluate the weathered bedrock permeability, and consequently to enhance our understanding of the soil saturation processes.The soil and bedrock hydrodynamic properties were estimated by inverse modelling of soil moisture from a pilot station and a network of 50 stations distributed on 4 granitic hillslope and 2 schist hillslope. The pilot station was instrumented with 3 soil moisture sensors located at 20, 40 and 60 cm deep whereas the hillslope stations were instrumented with 2 soil moisture sensors located variably from 15 to 45 cm deep. Both soil moisture and precipitation were recorded with a 15 min time step for every station. The inverse modelling procedure is based on the multi-objective genetic algorithm NSGA-II. This procedure was used for every station considering a variably deep soil composed by 2 layers surmounting a 100 cm deep third layer representing the weathered bedrock. Fifteen parameters were calibrated for every station, and were estimated from small selected rainfall periods of 1 to 15 days corresponding to the major rainfall events during the monitoring period. Evapotranspiration has been considered as negligible during those events. The analysis of those retention and conductivity properties shows a very important hydraulic conductivity for the studied soils, ranging from 1000 to 2000 mm/h, and their low retention capacity. The soil saturated water content varies from 0.30 to 0.60 cm3.cm-3, which is a consequence of the important soil heterogeneity at the hillslope scale. We also find that heterogeneity on the estimation of soil depth which varies from 31 to 120 cm. For the weathered bedrock layer, we also found contrasted permeability varying from a few units to a hundred of mm/h. Stations with the shallowest soils and the less permeable bedrock presenting a saturation process at the soil/bedrock contact, as monitored and simulated during the most intense rainfalls. If the estimated soil hydrodynamic properties exhibit important variations at the hillslope scale from a station to another, results showed no specific spatial organization of this variability. The station's position on the hillslope is not decisive to estimate saturated soil moisture, soil depth, soil permeability or water retention capacity. The distribution of those properties, calculated for the 6 studied hillslope is not either correlated to the landscape general characteristic: geology, orientation, land cover.

Les cycles océaniques du carbone et des principaux nutriments sont mal connus car ils sont affectés par de nombreux puits et sources physiques, chimiques ou biologiques difficiles à estimer par des mesures directes. Une manière de mieux contraindre ces processus importants est d’utiliser l’information contenue dans des traceurs plus simples : les proxies. Le radium 228 (228Ra), émis par les plateaux continentaux, est utilisé comme proxy des flux d’eau et d’éléments minéraux de la côte vers l’océan ouvert. Il permet en particulier d’estimer les flux d’eau souterraine ou SGDs (Submarine Groundwater Discharge). Le thorium 234(234Th), insoluble, permet quant à lui de contraindre la dynamique des particules par lesquelles il est adsorbé. Il est régulièrement utilisé pour estimer la pompe biologique du carbone (PBC), c’est-à-dire le flux de carbone de la surface vers l’océan profond.Au cours de cette thèse, un modèle numérique à une résolution de 2° a été construit pour chacun de ces deux radio-isotopes, en s’appuyant sur la circulation du modèle NEMO-OPA et les champs de particules du modèle PISCES. Plusieurs paramètres inconnus des modèles ont été contraints par des observations dans le cadre d’une méthode inverse.La modélisation inverse du 228Ra a permis d’estimer les flux de 228Ra venant de 38 régions côtières. En revanche, l’estimation des SGDs est imprécise, car les SGDs sont difficiles à distinguer d’une autre source de 228Ra: la diffusion par les sédiments.La modélisation inverse du 234Th a permis d’estimer les coefficients de partage du 234Th, qui représentent l’affinité de différents types de particules pour cet isotope. Elle a aussi permis d’estimer les erreurs associées à quelques simplifications courantes dans les études de la PBC fondées sur le 234Th. / The oceanic cycles of carbon and the main nutrients are poorly known since they are affected by many physical, chemical or biological sources and sinks that are difficult to estimate by direct measurements.One way to better constrain these important processes is to use the information contained in more simple tracers called "proxies". As radium 228 (228Ra) flows from the continental shelves, it is used as a proxy of water and mineral elements fluxes from the coast to the open ocean. In particular, it is often used to estimate the SGD (Submarine Groundwater Discharge). For its part, thorium 234 (234Th), an insoluble radio-isotope, is used to constrain the dynamics of the solid particles onto which it is adsorbed. The carbon flux from the surface to the deep ocean, called "biological carbon pump" (BCP), is often estimated by a 234Th-based method.During this PhD, a numerical model with a resolution of 2°, based on the circulation of the NEMO-OPA model and the particle fields of the PISCES model, was built for each of the two radioisotopes.Several unknown model parameters were constrained by observations using an inverse technique.The inverse modeling of 228Ra was used to constrain 228Ra fluxes from 38 coastal regions.However, the SGD fluxes are poorly constrained by this method, because SGD can be confused with another source of 228Ra: diffusion from sediments.The inverse modeling of 234Th produced estimates of partition coefficients, representing the affinity of different particle types for this isotope. It was also used to estimate the errors associated with some common simplifications made in 234Th-based BCP studies.

La modélisation inverse de la dispersion atmosphérique consiste à reconstruire les caractéristiques d’une source (quantité de polluants rejetée, position) à partir de mesures de concentration dans l’air, en utilisant un modèle direct de dispersion et un algorithme d’inversion. Nous avons utilisé dans cette étude deux modèles directs de dispersion atmosphérique SIRANE (Soulhac, 2000; Soulhac et al., 2011) et SIRANERISK (Cierco et Soulhac, 2009a; Lamaison et al., 2011a, 2011b). Il s’agit de deux modèles opérationnels de « réseau des rues », basés sur le calcul du bilan de masse à différents niveaux du réseau. Leur concept permet de décrire correctement les différents phénomènes physiques de dispersion et de transport de la pollution atmosphérique dans des réseaux urbains complexes. L’étude de validation de ces deux modèles directs de dispersion a été effectuée après avoir évalué la fiabilité des paramétrages adoptés pour simuler les échanges verticaux entre la canopée et l'atmosphère, les transferts aux intersections de rues et la canalisation de l’écoulement à l’intérieur du réseau de rues. Pour cela, nous avons utilisé des mesures en soufflerie effectuées dans plusieurs configurations académiques. Nous avons développé au cours de cette thèse un système de modélisation inverse de dispersion atmosphérique (nommé ReWind) qui consiste à déterminer les caractéristiques d’une source de polluant (débit, position) à partir des concentrations mesurées, en résolvant numériquement le système matriciel linéaire qui relie le vecteur des débits au vecteur des concentrations. La fiabilité des résultats et l’optimisation des temps de calcul d’inversion sont assurées par le couplage de plusieurs méthodes mathématiques de résolution et d’optimisation, bien adaptées pour traiter le cas des problèmes mal posés. L’étude de sensibilité de cet algorithme d’inversion à certains paramètres d’entrée (comme les conditions météorologiques, les positions des récepteurs,…) a été effectuée en utilisant des observations synthétiques (fictives) fournies par le modèle direct de dispersion atmosphérique. La spécificité des travaux entrepris dans le cadre de ce travail a consisté à appliquer ReWind dans des configurations complexes de quartier urbain, et à utiliser toute la variabilité turbulente des mesures expérimentales obtenues en soufflerie pour qualifier ses performances à reconstruire les paramètres sources dans des conditions représentatives de situations de crise en milieu urbain ou industriel. L’application de l’approche inverse en utilisant des signaux instantanés de concentration mesurés en soufflerie plutôt que des valeurs moyennes, a montré que le modèle ReWind fournit des résultats d’inversion qui sont globalement satisfaisants et particulièrement encourageants en termes de reproduction de la quantité de masse totale de polluant rejetée dans l’atmosphère. Cependant, l’algorithme présente quelques difficultés pour estimer à la fois le débit et la position de la source dans certains cas. En effet, les résultats de l’inversion sont assez influencés par le critère de recherche (d’optimisation), le nombre de récepteurs impactés par le panache, la qualité des observations et la fiabilité du modèle direct de dispersion atmosphérique. / The aim of this study is to develop an inverse atmospheric dispersion model for crisis management in urban areas and industrial sites. The inverse modes allows for the reconstruction of the characteristics of a pollutant source (emission rate, position) from concentration measurements, by combining a direct dispersion model and an inversion algorithm, and assuming as known both site topography and meteorological conditions. The direct models used in these study, named SIRANE and SIRANERISK, are both operational "street network" models. These are based on the decomposition of the urban atmosphere into two sub-domains: the urban boundary layer and the urban canopy, represented as a series of interconnected boxes. Parametric laws govern the mass exchanges between the boxes under the assumption that the pollutant dispersion within the canopy can be fully simulated by modelling three main bulk transfer phenomena: channelling along street axes, transfers at street intersections and vertical exchange between a street canyon and the overlying atmosphere. The first part of this study is devoted to a detailed validation of these direct models in order to test the parameterisations implemented in them. This is achieved by comparing their outputs with wind tunnel experiments of the dispersion of steady and unsteady pollutant releases in idealised urban geometries. In the second part we use these models and experiments to test the performances of an inversion algorithm, named REWind. The specificity of this work is twofold. The first concerns the application of the inversion algorithm - using as input data instantaneous concentration signals registered at fixed receptors and not only time-averaged or ensemble averaged concentrations. - in urban like geometries, using an operational urban dispersion model as direct model. The application of the inverse approach by using instantaneous concentration signals rather than the averaged concentrations showed that the ReWind model generally provides reliable estimates of the total pollutant mass discharged at the source. However, the algorithm has some difficulties in estimating both emission rate and position of the source. We also show that the performances of the inversion algorithm are significantly influenced by the cost function used to the optimization, the number of receptors and the parameterizations adopted in the direct atmospheric dispersion model.

Les eaux souterraines constituent un des principaux facteurs de contrôle influençant l'activité des glissements de terrain. L'étude hydrogéologique constitue donc une étape incontournable dans la compréhension de leur fonctionnement, en vue de prédire leur déplacement ou d'adapter les solutions d'assainissement. Cependant, les aquifères qui se développent au sein des glissements de terrain sont généralement complexes, et la compréhension de l'hydrogéologie des grands glissements reste encore aujourd'hui une problématique majeure. Dans un contexte climatique extrême de type tropical humide, les glissements du cirque de Salazie (île de La Réunion) constituent des cas d'étude particulièrement porteurs pour l'apport de nouvelles connaissances sur l'hydrogéologie et le fonctionnement des mouvements de terrain de grande ampleur. Pour caractériser ces glissements, une étude pluridisciplinaire associant géologie, hydrochimie, modélisation numérique (globale et déterministe), est mise en œuvre. Sur le secteur nord du cirque de Salazie, l'étude géologique réalisée permet de reconstruire la géométrie de la base d'une instabilité de très grande ampleur, dont le volume est estimé à 370 Mm3. Le glissement de Grand Ilet (175 Mm3) en est la partie la plus active. Les formations volcano-détritiques (brèches) constituant le glissement sont des dépôts d'une ancienne avalanche de débris (> 350 Ma), réactivés localement. Cette étude met également en évidence un régime de déformation extensif associée à l'activité du glissement, favorisant le développement de zones décomprimées subverticales. Au sein des brèches constituant la zone instable, un aquifère continu est identifié. Les écoulements souterrains y sont contrôlés par la géométrie de la base du glissement. Par ailleurs, la nature du couvert végétal en surface, le climat, les propriétés des sols et l'épaisseur de la zone non saturée font que seuls les épisodes pluvieux intenses dont le cumul dépasse 80 mm/ évènement génèrent une recharge de l'aquifère des brèches. Lors de cette recharge, les zones décomprimées favorisent des transferts rapides (< 1 jour) à travers la zone non-saturée, malgré des épaisseurs pouvant atteindre localement 160 mètres. Une infiltration plus lente percole à travers la matrice poreuse des brèches. Pour l'étude de la dynamique des grands glissements dans le cirque de Salazie, l'application de l'outil de modélisation inverse met en évidence que les variations de vitesses de déplacement sont directement reliée aux processus de recharge et de tarissement de la nappe des brèches. Ces modèles s'avèrent particulièrement performants pour prédire les déplacements des grands glissements et comprendre leur fonctionnement. / The landslide activity is commonly controlled by the variation of hydraulic head inside the instable mass. Thus, the hydrogeological study of landslides is an essential step to predict landslide dynamic, and for the remediation choices. However, the aquifers developed in landslide are generally complexes. Actually, the comprehension of landslide hydrogeological functioning is a major problematic, especially for the deep-seated landslides.Under humid tropical climate, the Salazie landslides (Reunion Island) allow to implement an interesting study to characterized deep-seated landslide hydrogeology and functioning. In this study, we performed a multidisciplinary approach, combining geology, hydrochemistry and numerical modeling (global and deterministic);The geological study allows the construction of the deep seated landslide geological model (Grand Ilet sector). The Grand Ilet landslide, corresponding to 175 Mm3, is the most active part of deep seated instability whose the total volume is estimated at 370 Mm3. The volcano-detritic lithologies (i.e basic breccia) constitute the main geological formation of the unstable mass. This breccia is a present-day reactivation of an old destabilization (> 350 Ma). Furthermore, the extensive deformation generated by the landslide activity allows the formation of decompressed zones.Inside breccias, a continuous aquifer is identified. Here, the groundwater flows are controlled by the geometry of the landslide base. Natural land cover, soils properties, unsaturated zone thickness and warm temperature limit the groundwater recharge. Only the intense rainfall episode (80 mm/event) can generated the recharge of landslide aquifer. During these events, rapid transfers circulating inside the decompressed zones have a significant effect on hydraulic charge variations. A slow component is infiltrated inside the porous medium of breccias.The inverse modeling methods with bimodal Gaussian-Exponential transfer function is applied to study the Salazie landslides dynamic. The results show that landslide speed variations are directly controlled by groundwater level variations during the hydrologic cycle (recharge and recession stage). Consequently, the inverse model is a powerful tool for predicting deep-seated landslide movements and for studying how they function.

La tomographie de résistivité électrique (ERT) est une méthode courante géophysique de terrain, souvent utilisée pour détecter et l’évolution suivre les panaches de polluants en zone saturée. L’ERT est cependant une méthode intégratrice dont la fiabilité des modèles est confronte aux problèmes de non unicité des solutions du problème inverse. Ces contraintes limitent l’interprétation des modèles a un aspect qualitatif de la distribution des contrastes de résistivité modélises en 2D ou 3D, résultant du choix des paramètres d’inversion et de l’association de paramètres du milieu non identifiables a l’échelle du volume poreux.Cette thèse propose de tester la faisabilité de la méthode pour quantifier les paramètres de transport de polluants et de solutés miscibles au contact des eaux souterraines, ainsi que la sensibilité des paramètres d’inversion les plus influents sur la modélisation.Les tests expérimentaux sont réalisés en laboratoire sur des empilements 2D de billes de verre sphériques (de l’ordre de la 100èn de μm) dans un réservoir en plexiglas transparent. Deux réseaux verticaux de 21 électrodes sont disposes sur les bords latéraux du réservoir pour effectuer le suivi ERT du traceur sale (NaCl dissout) a partir de 210 points de mesures en dipôle-dipôle transverse acquis toutes les 5 minutes afin d’optimiser la résolution temporelle. Le dispositif est également dispose face a un panneau lumineux permettant de réalise un suivi vidéo simultané du colorant.L’analyse vidéo révèle une propagation plus rapide du colorant sur les bords latéraux qui reste somme toute négligeable pour les débits a la pompe les plus faibles. En revanche les mesures ERT sont fortement perturbées par les effets résistant de la cellule plexiglas qui se répercutent sur les modèles. La normalisation des mesures de résistivité apparente à partir d’une série acquise à l’état initial permet de les atténuer fortement.La modélisation est particulièrement sensible au choix du maillage, aux normes appliques (L1 ou L2) sur les données et les paramètres, et au facteur d’amortissement _. Des valeurs trop élevées de _ et du facteur d’acceptance tendent à lisser les contrastes au niveau du front de dispersion et augmentent l’impact des effets des bords horizontaux sur D et _. A l’inverse, une modélisation contrainte par de faibles valeurs de α et du facteur d’acceptante donne des résultats plus proches l’analyse vidéo, mais produit des effets de bosses à l’avant et à l’arrière du front.La vitesse interstitiel u est indépendante du choix des paramètres d’inversion pour l’ERT. Pour les deux méthodes u est toujours inférieure au débit impose par la pompe, dont le décalage est exprimé par le facteur retard Rf . Les effets de retard résultent de l’adsorption du Na+ sur les surfaces des billes de verre chargées négativement qui retarde le front de dispersion du suivi ERT. Pour le suivi vidéo, la taille importante de la molécule du colorant favorise son piégeage dans les zones ou la perméabilité est plus faible, en plus d’une éventuelle affinité avec la surface solide. Les contrastes de conductivité et la stabilité de l’interface créent par la différence de densité entre les fluides testes ici n’ont pas d’influences significatives sur la dispersion qui est dominée par le débit impose a la pompe. Les estimations du coefficient de dispersion D en fonction du nombre de Péclet sont cohérentes avec la courbe théorique de Bachmat (1968). Cependant la dispersivité α augmente pour les vitesses d’écoulement les plus élevées. Les premières expérimentations de terrain réalises en 2D sur des sables de Fontainebleau présentent l’avantage de s’affranchir des effets de bords inhérents au laboratoire. En revanche la recalibration des données normalisées par la loi de Archie est plus complexe puisqu’il est nécessaire de tenir compte de l’état de saturation de la résistivité des fluides initialement présents. De plus l’erreur importante sur les modèles ne permet pas de déduire une estimation fiable des paramètres de transport u, α (ou D), et Rf . / Electrical resistivity tomography is a common geophysical method often used to detect and follow plumepollutants in aquifers. However ERT is an integrative method whose reliability of the models is faced tothe non-unicity of the inverse problem solutions. These constraints limit the interpretation to a qualitativeview of the resistivity contrasts modelled in 2D or 3D, resulting of the chosen inverse parameters and thecombination of several hydrodynamic paramaters related to the poral network.The purpose of this thesis was to test the abilities of the ERT imaging to quantify solute transport parametersin miscible displacement occurred in groundwater and the sensitivity of inverse parameters most affectingthe modelled dispersion front.Laboratory experiments are conducted on glass beads poured into a transparent plexiglas container. Twovertical lines of 21 stainless steel electrodes are fixed on the lateral sides of the container to perform the ERTmonitoring, of the NaCl dissolved in the tracer, from a sequence of 210 quadripole measurements acquiredin transverse dipole-dipole each 5 minutes. A light panel is placed behind the experimental device and avideo follow up of the dyed part of the tracer is acquired from the other side.Video analysis reveal a faster propagation of the dye in contact of the vertical edge, which is negligible forthe lowest flow rates imposed by the pump. In contrast, ERT mesurements are strongly disturbed by theresistant edges of the plexiglas container which affect the resulting models. Normalisation of the apparentresisitivity measurements acquired at the experimental stage and by the Archie’s law strongly tones downthese resistive artefacts.ERT modelling is here particularly sensitive to the grid mesh, the norm (L1 or L2 ) applied on data andparameters, and the damping factor λ. High values of λ and the cutoff factor tend to smooth the resistivityconstrasts in the area of the mixing front and increase the weight of the horizontal edge effects on D andα. While results from inverse modelling constraint by low λ and cutoff factors are much closer to the videoanalysis but with enhanced side slope effects at the rear and the front of the mixing area.The interstitial velocity u is independant of the chosen inverse parameters. For both methods u is alwaysinferior to the flow rate provided by the pump, whose the gap is expressed as a retardation factor Rf . Thisretardation is due to adsorption of Na+ on the beads surfaces, which contributes to delay the dispersionfront followed by ERT. The retardation expressed by the video analysis can be due to the important sizeof the molecule of the dye which is easily slowed down in lower permeability areas, added to an eventualaffinity to the solid surface.The ranges of fluid conductivity contrasts and stability of the interfaces tested here have no influences onthe dispersion which is dominated by the flow velocity u. Estimations of the dispersion coefficient D asfunction of the Péclet is consistent with the theoretical curve of Bachmat (1968) and Bijeljic & al (2004).Field experiments are first conducted in 2D on homogeneous unsaturated sand which is considered as aninfinite half-space. However, data normalisation is much more complicated since the saturation state andthe initial fluid conductivities need to be estimated to calibrate the Archie’s law. Because of the 3D tracerinflitration, the RMS error of 2D-ERT models highlights that the inversion process is not enough constraintby data which does not allow to quantify the transport parameters. 3D experiments were then adaptedto detect and follow plumes of saline tracers injected in the centre of the electrode device. From 3D ERTmeasurements we are able to produce reliable models in order to estimate such transport parameters as themean flow velocity, and transverse and longitudinal dispersivities.

Dans cette recherche, nous tentons d'évaluer la pertinence d'une approche de simulation inverse des facteurs physiques d'ambiances dans le processus de conception en architecture et en aménagement urbain. Dans cette perspective, un outil de simulation inverse est développé et interfacé en sur-couche du logiciel de CAO AutoCAD largement utilisé dans les agences d'architectures. Il agrége les méthodes développées au Cerma ces dernières années sur la simulation inverse, de l'ensoleillement d'une part, et de la visibilité d'autre part. Dans le projet architectural ou urbain, la simulation directe solaire ou visuelle permet l'analyse d'une situation à un instant donné ou pour un point de vue donné a posteriori (après la conception d'un projet), tandis que notre démarche inverse consiste à cadrer la génération de solutions architecturales à partir des intentions d'ambiances du concepteur, et ce pour des périodes solaires et des zones de contraintes définies au début d'un projet. Les modélisations inverses de volumes de contraintes solaires et visuelles – ensembles de solutions à une contrainte donnée – ou d'uniques solutions d'ouvertures ou d'écrans sont des réponses optimales aux problèmes spatio-temporels de l'ensoleillement et de la visibilité. Ces volumes englobent l'ensemble des rayons solaires ou visuels joignant une base de contrainte (une zone à ensoleiller ou à protéger du soleil, ou un ensemble de positions potentielles d'un observateur) à une cible de contrainte (période solaire sur la "voûte solaire" ou cible visuelle du type monument, repère...). Ces volumes de contraintes sont aisément manipulables dans le modeleur géométrique. Comme les autres entités qui composent une scène numérique, il est possible de leur faire subir des opérations booléennes, des sections, des translations... Dans un premier temps, nous présentons un état de l'art de cette recherche sur les différents outils et méthodes de simulation inverse 2D puis 3D de l'ensoleillement et de la visibilité. Dans un second temps, nous développons et présentons l'outil SVR ainsi que ses applications. Enfin, nous livrons un état de l'art de quelques observations et méthodes d'observation des concepteurs à l'œuvre et nous évaluons notre logiciel par l'observation des raisonnements et comportements de onze architectes et étudiants en architecture concevant avec notre outil durant des phases d'esquisses et d'avant-projet. Pour cette évaluation, nous mettons en œuvre un protocole d'observation simultané ("think aloud protocol").

La caracterisation de l'agencement spatial des proprietes hydrauliques est essentielle pour predire les ecoulements et le transport des solutes dans les milieux heterogenes. Les methodes de tomographie hydraulique, principalement developpees pour estimer les proprietes des milieux poreux, n'ont qu'une faible r'esolution spatiale qui ne reflete pas la vraie heterogeneite des distributions de fractures des milieux fractures. Le principal objectif de cette these est de developper une nouvelle methode d'inversion specifique pour imager les proprietés hydrauliques et de transport des milieux fractures a l'echelle du site. Pour atteindre ces objectifs, des experiences in situ ainsi qu'une nouvelle approche de modelisation inverse sont proposees, notamment en utilisant la temperature comme marqueur des ecoulements. Nous proposons tout d'abord la tomographie d'ecoulement bas'ee sur des tests s'equentiels de debimetrie entre puits, comme une nouvelle approche pour caracteriser la connectivit'e des fractures ainsi que leur transmissivite. A partir de simulations numeriques reproduisant des cas d'etudes synth'etiques, nous montrons que l'approche par tomographie r'eduit significativement l'incertitude sur les parametres estimes, et fournit une caracterisation detaillee du reseau de fracture sans requerir a l'utilisation d'obturateurs hydrauliques. Nous montrons ensuite comment les mesures de temperature peuvent etre utilisees pour quantifier les ecoulements dans les milieux fractur'es. Le grand int'erˆet d'utiliser la temperature est d'obtenir facilement et de facon continue en puits des profils de temp'erature. En utilisant un mod'ele numerique d'ecoulement et de transfert de chaleur a l'echelle du puits, une methode d'inversion pour estimer les vitesses d'ecoulement dans le puits 'a partir des donnes de temperature est proposee. Nous couplons ensuite les deux approches presentees precedemment dans une nouvelle approche experimentale consistant en des enregistrements sequentiels de temperature dans un puits dans des conditions de pompage entre puits. L'application de cette approche de tomographie en temperature sur le site de Stanger Brune montre des resultats encourageants pour l'identification du reseau global de connectivite et des zones d'ecoulement principales. Enfin, nous discutons de l'interet d'utiliser la chaleur comme traceur par rapport 'a l'utilisation de traceurs classiques. Nous montrons que realiser des tests de tracage thermiques en milieu fracture fournit des contraintes supplementaires importantes sur les propri'et'es de transport du milieu.

La gestion et la dépollution de sites contaminés peuvent être complexe et demandent un investissement important pour localiser les sources de contaminations, zones émettant les flux de polluants les plus importants. Les travaux réalisés proposent une stratégie pour localiser les sources de pollution à partir de mesures in situ de flux massiques et de modélisation inverse. Ainsi, dans le cadre de l’étude, un outil innovant a d’abord été développé afin de mesurer la vitesse des eaux souterraines dans un puits. L’outil appelé DVT (Direct Velocity Tool) a permis de répondre aux contraintes imposées par les outils existants et de mesurer des vitesses d’écoulement très lentes. Des essais en laboratoire et des tests en site réels ont été réalisés et comparés à d’autres outils de mesure. Le DVT permet aussi indirectement de définir la portion de source conduisant au flux de polluant maximum, en le combinant avec une mesure locale de concentration. L’étude présente ensuite l’utilisation de la modélisation inverse pour localiser une source de contaminant et d’estimer les paramètres définissant les caractéristiques du domaine. Pour cela, l'étude s'est faite sur deux cas synthétiques. Pour adapter les méthodes à une véritable gestion de sites pollués, une stratégie itérative est développée en imposant un ajout limité de nouvelles observations à chaque phase de modélisation, basée sur l’approche de type Data Worth. Les résultats de la position de la source sur les deux cas synthétiques ont permis d’évaluer la méthode mise en place et de juger son applicabilité à une problématique réelle. Cette stratégie de localisation de source est par la suite testée sur un site réel à partir (i) de mesures in situ de flux massiques avec les vitesses au DVT et les concentrations et (ii) la modélisation. Les essais ont permis de cibler les forages à mettre en place sur site aidant à localiser la source. Néanmoins, en analysant plus précisément les résultats, le champ de conductivité hydraulique estimé par l'optimisation ne correspond pas à la réalité. De plus, les flux massiques de contaminants ainsi que le ratio des polluants du site, mettent en valeur deux panaches distincts. Une phase finale de modélisation a donc été lancée afin d'estimer (i) la présence potentielle de deux sources et (ii) la chimie de la zone étudiée. Les résultats de la stratégie sont comparés aux mesures geoprobe qui a pu confirmer la présence d’une des deux sources identifiées. / Contaminated sites management and remediation can be complex and require a significant investment to locate the contaminant source, which delivers the higher pollutant mass fluxes. The study proposes a strategy for contaminant source localisation using in situ measurement and inverse modelling. First, an innovative tool was developed to measure groundwater velocity in a well. The developed tool called DVT (Direct Velocity Tool) made it possible to measure a low Darcy flux. Laboratory and field tests were performed with the DVT and compared to other velocity measurement tools. By combining the DVT with a local concentration measurement, it is possible to calculate the mass fluxes passing through wells. Then the thesis present the inverse modeling used for source localisation and parameters estimation. The study was done on two synthetics cases using the non-linear optimisation method. To adapt the method to a real management of polluted sites, an iterative strategy is developed by imposing a limited addition of new observations to each modeling phase. This strategy is base on the Data Worth approach. Source localisation results on the two synthetic cases made it possible to judge the method applicability to a real site problem. The source localisation strategy is then applied to a real site with (i) mass flux measurement with velocities (DVT) and concentrations and (ii) inverse modeling. The modeling phases made it possible to locate the new wells and helped the source localisation. Nevertheless, by analysing the results more precisely, the hydraulic conductivity field estimated by the optimisation did not correspond to reality. In addition, contaminant mass fluxes highlightes two distinct zones of flux. By analysing the pollutant ratio of the site, it appears that two plumes are potentially present. Thus, another inverse modeling phase has been tested (i) to locate the two potential sources and (ii) to estimate the chemistry of the site. Results of the strategy were compared to the geoprobe campaign which confirmed the second source location.

Dans cette thèse, plusieurs approches de modélisation de composites renforcés par des fibres sont proposées. Le matériau étudié est le béton fibré, et dans ce modèle, on tient compte de l’influence de trois constituants : le béton, les fibres, et la liaison entre eux. Le comportement du béton est analysé avec un modèle d’endommagement, les fibres d'acier sont considérées comme élastiques linéaires, et le comportement sur l'interface est décrit avec une loi de glissement avec l’extraction complète de la fibre. Une approche multi-échelle pour coupler tous les constituants est proposée, dans laquelle le calcul à l'échelle macro est effectué en utilisant la procédure de solution operator-split. Cette approche partitionnée divise le calcul en deux phases, globale et locale, dans lesquelles différents mécanismes de rupture sont traités séparément, ce qui est conforme au comportement du composite observé expérimentalement. L'identification des paramètres est effectuée en minimisant l'erreur entre les valeurs calculées et mesurées. Les modèles proposés sont validés par des exemples numériques. / In this thesis, several approaches for modeling fiber-reinforced composites are proposed. The material under consideration is fiber-reinforced concrete, which is composed of a few constituents: concrete, short steel fibers, and the interface between them. The behavior of concrete is described by a damage model with localized failure, fibers are taken to be linear elastic, and the behavior of the interface is modeled with a bond-slip pull-out law. A multi-scale approach for coupling all the constituents is proposed, where the macro-scale computation is carried out using the operator-split solution procedure. This partitioned approach divides the computation in two phases, global and local, where different failure mechanisms are treated separately, which is in accordance with the experimentally observed composite behavior. An inverse model for fiber-reinforced concrete is presented, where the stochastic caracterization of the fibers is known from their distribution inside the domain. Parameter identification is performed by minimizing the error between the computed and measured values. The proposed models are validated through numerical examples.