Au cours de la production d'hydrocarbures, l'extraction de fluides fait décroître la pression de pore dans les réservoirs (« depletion »). Ceci induit un changement du champ de contraintes qui résulte en une augmentation des contraintes effectives appliquées sur le réservoir. Les mesures in situ montrent que les variations de contraintes peuvent être décrites par un paramètre appelé chemin de chargement (stress path), défini comme le rapport entre la variation de contrainte effective horizontale et la variation de contrainte effective verticale par rapport aux conditions initiales dans le réservoir. La compaction induite par la production d'hydrocarbures peut avoir de graves conséquences dans le cas de roches faiblement consolidées car elle induit des variations des propriétés pétrophysiques des roches in situ, notamment de la perméabilité, un des paramètres les plus importants pour estimer la performance d'un réservoir mais aussi un des plus difficiles à mesurer. Pour compliquer encore les choses, la perméabilité est souvent anisotrope dans les réservoirs avec de forts contrastes entre la perméabilité horizontale kh et la perméabilité verticale kv.L'objectif de cette étude est de comprendre l'influence des chemins de chargement sur le comportement mécanique et les évolutions couplées de perméabilité pour un sable quartzeux (Sable de la Durance, DS) et un grès faiblement consolidé (grès d'Otter Sherwood, OSS, qui constitue la roche réservoir du champ pétrolier de Wytch Farm en Angleterre). Nos résultats montrent que le grès peu consolidé présente un comportement mécanique similaire à celui de roches consolidées. Au contraire, le sable présente un comportement différent, avec une transition plus graduelle entre les régimes de déformation qui nécessite d'utiliser un critère basé sur l'évolution du rayon de courbure des courbes contraintes-déformations pour déterminer les contraintes limites : cette méthode a été validée par une étude d'analyse des émissions acoustiques pour caractériser l'endommagement. Les domaines de déformation élastique et plastique ont été bien définis et les contraintes limites ont été comparées aux prédictions du modèle Cam-Clay modifié et du modèle d'enveloppe limite normalisée. Les perméabilités horizontale et verticale ont été mesurées sous contraintes. Pour analyser l'influence des effets de bord dans les essais mécaniques, les perméabilités mesurées soit classiquement sur toute la longueur de l'échantillon, soit entre deux points intermédiaires ont été comparées. Pour l'écoulement horizontal, les facteurs géométriques et facteurs d'anisotropie ont été déterminés par des simulations numériques en éléments finis afin de pouvoir déterminer les vraies valeurs de perméabilité horizontale. L'évolution de la perméabilité suit l'évolution de la déformation des matériaux et est contrôlée aussi bien par la déformation volumique que par la déformation en cisaillement. A partir de nos mesures il est possible de séparer l'effet de la pression moyenne de l'effet de la contrainte déviatorique sur l'évolution de la perméabilité en construisant des cartes d'isoperméabilités dans l'espace des contraintes. Enfin une modélisation élasto-plastique a été réalisée pour prédire le comportement hydro-mécanique du grès faiblement consolidé. L'approche utilisée permet de prédire de manière satisfaisante l'évolution de la perméabilité avec les contraintes, à partir d'une loi exponentielle fonction de la déformation effective. Au contraire, pour le sable de la Durance le lien entre l'évolution de la perméabilité et la déformation est loin d'être évidente, notamment aux faibles contraintes où la réduction de perméabilité est très rapide. Pour mieux comprendre ces évolutions de perméabilité, une analyse de l'endommagement a été réalisée par des mesures sur échantillons et en utilisant des techniques d'imagerie à plusieurs échelles. / During hydrocarbon production, the extraction of fluid induces a decrease of pore pressure called depletion. This depletion causes a change in the stress field that results in an increased stress on the rock by enhancement of the effective stress in the reservoir. In situ measurements show that the stress variations can be described by the so-called stress path parameter, defined as the ratio of the change in effective horizontal stress by the change in effective overburden stress from initial reservoir conditions. This production induced compaction can have severe consequences in the case of poorly consolidated reservoirs. Compaction induces variations of petrophysical properties of in situ rocks and particularly permeability variations, one of the most important parameters controlling reservoir performance. Yet it is one of the most difficult property to measure. To complicate matters further, permeability anisotropy is often found in reservoirs. Therefore the horizontal permeability kh, may be different from the vertical permeability kv.The aim of this study is to understand the influence of stress paths on the mechanical behavior and coupled permeability evolutions of a Quartz sand (Durance Sand, DS) and a weakly consolidated sandstone (Otter Sherwood Sandtone, OSS which is the reservoir rock of the Wytch Farm oil field, UK). We found that the weakly consolidated rock presents a mechanical behavior similar to that of consolidated rocks. However, the sand shows a different behavior, with a gradual transition regime which requires the use of a curvature criterion to peak yield stresses on the stress-strain evolution plot; this criterion has been validated on the basis of Acoustic Emission analysis. The elastic and plastic deformation regimes are well identified and the determined yield stresses are fitted using the modified Cam-Clay and Elliptic Cap models for all observed onsets of plastic yielding. Both vertical and horizontal permeability have been measured during loading. To analyze the influence of end effects during loading in the triaxial cell, permeabilities measured over the mid-section and over the total core length were compared. For the horizontal flow, the geometrical and anisotropy factors were determined using Finite Element simulations in order to calculate the correct horizontal permeability. Permeability evolution follows closely the material deformation and is controlled by both volumetric and shear strains. It is possible to infer the effect of the mean pressure and/or the deviatoric stress on the permeability evolution by building isopermeability maps in the stress space. Finally, an application of elasto-plastic modeling to predict the hydromechanical behavior of the weakly consolidated rock is presented. This approach allows a satisfying prediction of the permeability evolution with stresses, using an exponential function of an effective strain. Reversely for DS, the link between strain and permeability is not obvious as permeability reduction is pronounced at early stage of loading. To understand these permeability evolutions, a damage analysis has been performed using core analysis measurements and multi-scale imaging.

Le changement climatique pourrait avoir des conséquences importantes pour les ressources en eaux et les infrastructures hydrauliques. Le dimensionnement et le fonctionnement des réservoirs devraient ainsi être modifiés en prenant en compte des scénarios de changement climatique. Un modèle de dimensionnement coût-bénéfice d'un réservoir à l'échelle du bassin versant est appliqué en Grèce afin d'évaluer le coût de l'incertitude sur le climat futur et les dommages du changement climatique. Dans cette étude de cas, le coût de l'erreur est faible, et l'adaptation n'est pas efficace. Une méthodologie sur toute la région méditerranéenne, avec une modélisation générique à l'échelle des bassins versant est ensuite développée. Les réseaux de réservoirs et les liens réservoirs-demandes sont reconstruits et le fonctionnement coordonné des réseaux de réservoirs est déterminé, en utilisant uniquement des données disponibles à l'échelle globale. La méthodologie de reconstruction des liens est appliquée à l'irrigation et validée qualitativement sur l'Algérie. Le changement de fiabilité, avec adaptation des règles opérationnelles, sous changement climatique, semble être plus influencé par les changements de ruissellement que par les changements de demande. Les changements obtenus pour le Nil, l'Europe et le Moyen Orient ne sont pas très marqués, alors que les pays d'Afrique du Nord sous influence méditerranéenne comme la Tunisie ou l'Algérie voient une diminution importante de la fiabilité des apports d'eau pour l'irrigation / Climate change could affect water resources and hydraulic infrastructures seriously. Dimensioning and operation of reservoirs should therefore be modified according to climatic change scenarios. To assess the effect of climate uncertainty on reservoir volumes dimensioning using cost-benefit analysis, a model of reservoir dimensioning at the river basin scale is applied in Greece. For the case study, there is no cost of error and adaptation seems to be inefficient. A methodology at the scale of the Mediterranean region with a generic modeling at the river basin level is developed. Reservoirs networks and reservoirs-demands links are reconstructed and coordinated reservoirs networks operation is determined, using only globally available data. The link reconstruction methodology is applied on irrigation demand and validated qualitatively on Algeria. Change in reliability with adaptation of reservoir operating rules under climate change over the Mediterranean region is then assessed. Reliability changes seem to be more influenced by inflow changes than by demand changes. They are not important for the Nile basin and the European and Middle East sub-regions, while in North African countries with a more pronounced Mediterranean influence, like Tunisia or Algeria, reliability decrease can be significant

Les modèles génétiques de réservoirs sont construits par la simulation des principaux processus de sédimentation dans le temps. En particulier, les modèles tridimensionnels de systèmes chenalisés méandriformes peuvent être construits à partir de trois processus principaux : la migration du chenal, l’aggradation du système et les avulsions, comme c’est réalisé dans le logiciel Flumy pour les environnements fluviatiles. Pour une utilisation opérationnelle, par exemple la simulation d'écoulements, ces simulations doivent être conditionnées aux données d'exploration disponibles (diagraphie de puits, sismique, …). Le travail présenté ici, basé largement sur des développements antérieurs, se concentre sur le conditionnement du modèle Flumy aux données de puits.Deux questions principales ont été examinées au cours de cette thèse. La première concerne la reproduction des données connues aux emplacements des puits. Cela se fait actuellement par une procédure de "conditionnement dynamique" qui consiste à adapter les processus du modèle pendant le déroulement de la simulation. Par exemple, le dépôt de sable aux emplacements de puits est favorisé, lorsque cela est souhaité, par une adaptation des processus de migration ou d'avulsion. Cependant, la manière dont les processus sont adaptés peut générer des effets indésirables et réduire le réalisme du modèle. Une étude approfondie a été réalisée afin d'identifier et d'analyser les impacts indésirables du conditionnement dynamique. Les impacts ont été observés à la fois à l'emplacement des puits et dans tout le modèle de blocs. Des développements ont été réalisés pour améliorer les algorithmes existants.La deuxième question concerne la détermination des paramètres d’entrée du modèle, qui doivent être cohérents avec les données de puits. Un outil spécial est intégré à Flumy - le "Non-Expert User Calculator" (Nexus) - qui permet de définir les paramètres de simulation à partir de trois paramètres clés : la proportion de sable, la profondeur maximale du chenal et l’extension latérale des corps sableux. Cependant, les réservoirs naturels comprennent souvent plusieurs unités stratigraphiques ayant leurs propres caractéristiques géologiques. L'identification de telles unités dans le domaine étudié est d'une importance primordiale avant de lancer une simulation conditionnelle avec des paramètres cohérents pour chaque unité. Une nouvelle méthode de détermination des unités stratigraphiques optimales à partir des données de puits est proposée. Elle est basée sur la Classification géostatistique hiérarchique appliquée à la courbe de proportion verticale (VPC) globale des puits. Les unités stratigraphiques ont pu être détectées à partir d'exemples de données synthétiques et de données de terrain, même lorsque la VPC globale des puits n'était pas visuellement représentative. / Process-based reservoir models are generated by the simulation of the main sedimentation processes in time. In particular, three-dimensional models of meandering channelized systems can be constructed from three main processes: migration of the channel, aggradation of the system and avulsions, as it is performed in Flumy software for fluvial environments. For an operational use, for instance flow simulation, these simulations need to be conditioned to available exploration data (well logging, seismic, …). The work presented here, largely based on previous developments, focuses on the conditioning of the Flumy model to well data.Two main questions have been considered during this thesis. The major one concerns the reproduction of known data at well locations. This is currently done by a "dynamic conditioning" procedure which consists in adapting the model processes while the simulation is running. For instance, the deposition of sand at well locations is favored, when desired, by an adaptation of migration or avulsion processes. However, the way the processes are adapted may generate undesirable effects and could reduce the model realism. A thorough study has been conducted in order to identify and analyze undesirable impacts of the dynamic conditioning. Such impacts were observed to be present both at the location of wells and throughout the block model. Developments have been made in order to improve the existing algorithms.The second question is related to the determination of the input model parameters, which should be consistent with the well data. A special tool is integrated in Flumy – the Non Expert User calculator (Nexus) – which permits to define the simulation parameters set from three key parameters: the sand proportion, the channel maximum depth and the sandbodies lateral extension. However, natural reservoirs often consist in several stratigraphic units with their own geological characteristics. The identification of such units within the studied domain is of prime importance before running a conditional simulation, with consistent parameters for each unit. A new method for determining optimal stratigraphic units from well data is proposed. It is based on the Hierarchical Geostatistical Clustering applied to the well global Vertical Proportion Curve (VPC). Stratigraphic units could be detected from synthetic and field data cases, even when the global well VPC was not visually representative.

Les Acinetobacter sont des bactéries opportunistes impliquées dans les infections nosocomiales.Le but de ce travail était d’étudier leur épidémiologie et écologie au Liban.Tout d’abord, nous avons analysé 119 souches d’A.baumannii isolées de plusieurs hôpitaux. 76.5 % étaient résistantes aux carbapénèmes et le gène OXA-23 était le plus fréquemment trouvé. Le typage par Multilocus sequence typing a montré que le clone international II était majoritairement détecté. L’électrophorèse en champ pulsé a révélé que 72.6% des souches appartenant au ST2 ont été classées dans un même cluster qui semble être prédominant à Beirut et Tripoli. Ensuite, les réservoirs extrahospitaliers ont été investigués sur 2361 prélèvements collectés au Liban. Au total, 171 souches ont été isolées dans l’environnement, les produits alimentaires ainsi que chez l’homme et les animaux. La majorité de ces souches, globalement sensibles aux antibiotiques, était des Acinetobacter non baumannii, Seuls 15 A.baumannii, de 14 STs différents dont 10 nouveaux ont été isolés. Enfin, nous avons conduit une étude taxonomique approfondie sur plusieurs souches d’Acinetobacter non identifiées au rang d’espèce et retrouvées dans notre étude. Nous avons ainsi caractérisé une nouvelle espèce, nommée « Acinetobacter lebanonensis ».Ce travail a montré que le Liban était un pays à forte endémie d’A.baumannii résistants aux carbapénèmes. Nous n’avons toutefois pas mis en évidence de lien entre les souches cliniques et extrahospitalières, les clones correspondants étant globalement différents. D’autres études sont nécessaires pour élucider l’origine des souches multi-résistantes émergeant dans les hôpitaux. / Acinetobacter spp are opportunistic bacteria widely involved in nosocomial infections. The aim of this work was to study the epidemiology and the ecology of these bacteria in Lebanon. First, we have analyzed 119 clinical strains of A.baumannii. 76.5% of them were resistant to carbapenems and the production of OXA-23 was the main mechanism. Multi-locus sequence typing revealed the predominance of international clone II. Pulsed field gel electrophoresis showed that 72.6% of strains belonging to ST2 were classified in the same cluster which appeared to be predominant in Beirut and Tripoli. On the other hands, Acinetobacter reservoirs were investigated on 2361 samples collected in Lebanon. A total number of 171 strains have been isolated in the environment, food, humans and animals. The majority of these strains was identified as non baumannii Acinetobacter and was susceptible to antibiotics. Besides, typing of A.baumannii revealed the presence of 14 STs including 10 new ones. Finally, we have described a novel species called “Acinetobacter lebanonensis” by conducting a taxonomic study on several strains isolated in Lebanon and other countries. Although the data may be limited, this work has shown the endemic situation of carbapenem resistant A.baumannii circulating in the Lebanese hospitals while the extra hospital ones were different. However, further studies are needed to elucidate the origin of these emerging multidrug resistant strains.

En exploration pétrolière, la technique d'imagerie la plus couramment utilisée reste l'imagerie par sismique réflexion. Elle permet d'obtenir une première représentation des structures géologiques mais aussi une estimation du coefficient de réflexion des interfaces du sous-sol. Les variations de ce coefficient dépendent des paramètres pétrophysiques des milieux situés de part et d'autre de la discontinuité donnant naissance à la réflexion (densité et vitesses des ondes de propagation et de cisaillement). Une analyse de ces variations permet ainsi d'évaluer les paramètres élastiques qui sont des données indispensables et complémentaires à l'interprétation des réservoirs. Les traitements classiques consistent à étudier ces variations directement sur des données enregistrées en fonction du déport et du temps de propagation. L'analyse s'effectue par inversion d'une approximation linéaire des équations de Zoeppritz régissant les variations théoriques de l'amplitude. Nous proposons ici d'étudier les variations de l'amplitude après la phase d'imagerie, directement dans le domaine profondeur. Pour réaliser cette analyse, nous proposons un traitement spécifique après l'imagerie et avant l'analyse de l'amplitude pour corriger les impacts liés à la qualité du modèle de vitesse, à l'interférence des ondelettes entre deux réflexions proches et à l'étirement de l'ondelette. Ce traitement a aussi permis l'obtention d'un squelette du sous-sol conduisant à une première interprétation en profondeur des réflecteurs majeurs.

A l'Est du Bassin de Paris, les carbonates du Dogger présentent habituellement de faibles porosités et perméabilités, en relation notamment avec une intense cimentation sparitique et d'importants phénomènes de pression-solution (stylolitisation). Il existe néanmoins quelques exceptions à cette tendance : des niveaux réservoirs sont recensés au toit de séquences de dépôt réputées émersives. Pour comprendre la répartition des propriétés pétrophysiques du Dogger calcaire, ce travail de thèse basé sur une étude sédimentologique, diagénétique, géochimique et pétrophysique, se propose d'examiner les influences de la diagenèse précoce puis de la diagenèse tardive sur les propriétés réservoirs dans ces carbonates. L'analyse sédimentologique d'un grand nombre de coupes et de 4 forages carottés répartis sur un profil Langres-Nancy révèle 18 faciès attribuables à 5 environnements de dépôts caractéristiques d'un profil de rampe carbonatée. Ces 4 intervalles sont marqués par des producteurs carbonatés et par des morphologies de rampe différents. Les δ18O et δ13C des fragments biogènes et l'évolution des cortèges argileux suggèrent que des changements environnementaux tels que la température, l'humidité et la pCO2 ont, en partie, contrôlé l'existence de ces intervalles bio-sédimentaires. Huit phases minérales mésogénétiques comprenant 3 générations de calcite, de la pyrite, de la silice, 2 générations de dolomite et de la fluorine se sont succédées depuis le Jurassique moyen, réduisant la porosité intergranulaire. En moyenne 70 % du volume de ces ciments sont formés par 2 zones successives de sparites de blocage (Bc1 et Bc2). Leurs signatures isotopiques indiquent des circulations de fluides parents météoriques ou mixtes probablement liées aux émersions du Crétacé inférieur. Les analyses pétrophysiques ont confirmé les faibles porosités et perméabilités générales de l'Est du Bassin de Paris (φ<5% and k<0.5mD). Toutefois, ces analyses pétrophysiques ont mis en exergue la présence de deux niveaux poreux (Niveau Poreux 1 et Niveau Poreux 2) situés dans la partie supérieure du Dogger calcaire. D'une dizaine de mètres d'épaisseur, ces niveaux poreux ont une porosité dépassant 15%. Leur perméabilité est variable, elle est élevée (k=100-700mD) lorsque les niveaux sont macroporeux et correspondent à des grainstones de shoals oolitiques, elle est faible (k=0,1-1mD) lorsque les faciès sont boueux (dépôts de lagon protégé) et dominés par de la microporosité inter-cristalline. L'architecture stratigraphique issue de l'étude sédimentologique permet de supposer que celle-ci a joué un rôle sur le positionnement des niveaux poreux au sommet du Dogger. En effet, le modèle stratigraphique montre la présence d'un écran imperméable marneux positionné latéralement aux calcaires du sommet du Dogger. Cet écran marneux a pu perturber ou empêcher les recharges en eaux météoriques, pouvant expliquer le sous-développement des ciments de blocage dans ces niveaux poreux du sommet du Dogger. L'étude intégrée sédimentologique et diagénétique montre que l'existence des niveaux poreux réservoirs (NP1 et NP2) semble tributaire d'une histoire sédimento-diagénétique particulière qui peut être résumée en 3 actes : (1) lithification et stabilisation minéralogique précoce sous des surfaces d'émersion (SB7 et SB8), rendant le sédiment résistant aux phénomènes de compaction et pression-solution ; (2) faibles recharges météoriques au Crétacé en raison de l'écran latéral formé par les Marnes à Rhynchonelles (Bathonien), empêchant le développement de ciments de blocage dans le sommet du Dogger (3) faibles recharges hydrothermales et cimentations associées durant l'Oligocène. Ces 3 facteurs ne se sont conjugués qu'au niveau de NP1 et de NP2 (niveaux réservoirs), sur une épaisseur cumulée de moins de 20 m, à comparer aux 210 m totalement compactés et/ou cimentés formant le reste du Dogger (niveaux non-réservoirs).

La rupture dynamique des réservoirs métalliques pressurisés résulte de la compétition entre la fissuration de la paroi et la sortie du fluide qui en est à l'origine. Pour comprendre ce phénomène, deux études expérimentales ont été mises en œuvre.<br />· La première a consisté à soumettre des plaques métalliques à des chargements sous forme d'explosions, et à examiner leur mode de fissuration. Les vitesses de fissuration semblent indépendantes du chargement, tandis que la taille des fragments augmente avec la violence de l'explosion.<br />· La deuxième a consisté à faire exploser des enceintes d'environ vingt litres. Elle a permis de confirmer une partie des tendances observées précédemment.

La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture "Trans-Mexicain Volcanique Belt " (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide.

Les réservoirs carbonatés sont exploités pour leur réserve d’eau potable, de ressource géothermique ou le stockage géologique du dioxyde de carbone. Ces réservoirs sont difficiles à caractériser à cause d’une histoire diagénétique souvent complexe. Cette thèse offre de nouvelles perspectives sur la caractérisation des propriétés pétrophysiques et élastiques des calcaires urgoniens de Provence. Une approche intégrée et multi-échelle est proposée pour caractériser les propriétés pétrophysiques et élastiques des carbonates. Cette étude est basée sur des mesures de vitesse d’ondes P (Vp) et S (Vs) à l’échelle du laboratoire (centimètre–décimètre) et du terrain (mètre–décamètre). En laboratoire, les Vp et Vs ainsi que l’anisotropie sont mesurées sur des plugs et sur des carottes, en utilisant différentes fréquences centrales ultrasonores. Sur le terrain, l’approche consiste à mesurer les Vp et Vs entre deux puits distant de 2 m sur une profondeur de 14 m. Les mesures sont ensuite interprétées en fonction de la géologie observée aux échelles macro- et microscopiques. Les principaux résultats montrent que les Vp et Vs moyennes sont indépendantes de l’échelle de mesure, car elles sont dictées par la porosité. L’anisotropie causée par les fractures (15%) et les hétérogénéités (5%) se manifestent par une variabilité autour des vitesses moyennes. L’approche adoptée ici a permis de définir les interactions entre les propriétés de la matrice, les hétérogénéités, les fractures et les propriétés élastiques des roches carbonatées. Elle a montré que les propriétés élastiques tout comme les structures géologiques varient en fonction de l’échelle. / Carbonate reservoirs are also exploited for water production, geothermal energy, and carbon geological storage. Their Geophysical characterization remains challenging because of complex diagenetic history. This work offers new insights into the characterization of petrophysical, and elastic properties of the Urgonian limestones in the Provence region. An integrated multi-scale approach is proposed to characterize carbonate rocks petrophysical and elastic properties. This study relies on P- and S-wave velocity (Vp and Vs) measurements carried out at laboratory (centimeter–decimeter) and field (meter–decameter) scales. Laboratory scale Vp, Vs, and anisotropy are measured on plugs and cores, while on the field they are measured between two boreholes (crosshole) over a distance of 2 m and 14 m depth. The measurements are then compared to the geology from the macro- to the microscopic scale. The main results show that the average Vp and Vs are porosity related, and are independent from scale. Anisotropy caused by fractures (15%) and heterogeneities (5%) is responsible for variations around the mean velocities. The approach adopted during this work has enabled to scope out the interplay between matrix properties, heterogeneity, fracturing, and elastic properties in carbonate rocks. It has shown that the elastic properties evolve with scale as well as the geological structures.

La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture “Trans-Mexicain Volcanique Belt » (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide. / La qualité globale des eaux continentales (lacs et réservoirs) continue à se détériorer dans de nombreuses régions du Mexique. Le réservoir Cointzio, situé dans la partie sud du Plateau Central Mexicain, sur la ceinture “Trans-Mexicain Volcanique Belt » (TMVB), ne fait pas exception. Ce réservoir tropical turbide est un système monomictique chaud (surface = 6 km2, capacité = 66 Mm3, temps de séjour ~ 1 an). Il est stratégique pour l'approvisionnement en eau potable de la ville de Morelia, capitale de l'état du Michoacán, et pour l'irrigation en aval pendant la saison sèche. Le réservoir est menacé par l'accumulation des sédiments et des nutriments provenant d'eaux non traitées dans le bassin versant amont. La haute teneur en particules d'argile très fines et le manque d'installations de traitement de l'eau conduisent à de graves épisodes d'eutrophisation (jusqu'à 70 µg chl. a L-1) ainsi qu'à des périodes saisonnières de forte turbidité (profondeur de Secchi < 30 cm) et d'anoxie (de Mai à Octobre). Sur la base de mesures intensives de terrain réalisées en 2009 (échantillonnage dans le bassin versant, sédiments déposés, profils verticaux de l'eau, bilan entrées et sorties) nous avons présenté une étude intégrée du fonctionnement hydrodynamique et biogéochimique du réservoir Cointzio. Les profils verticaux de température (T), de concentration en sédiment en suspension (TSS), d'oxygène dissous (DO), de chlorophylle a, de carbone et de nutriments ont été utilisés pour évaluer les cycles internes au réservoir. Pour compléter l'analyse des données de terrain, nous avons examiné la capacité de différents modèles numériques 1DV (Aquasim modèle biogéochimique couplé avec le modèle de mélange k-ε) à (i) reproduire les principaux cycles biogéochimiques dans le réservoir Cointzio et (ii) préciser les scénarios de réduction des éléments nutritifs (P and N) pour limiter l'eutrophisation dans les prochaines décennies. Le modèle k-ε s'est avéré le plus robuste pour reproduire les conditions hydrodynamiques, en particulier la stratification modérée liée à la très grande turbidité et au régime soutenu de vent thermique. Le modèle Aquasim s'est avéré pertinent pour reproduire les principaux cycles de l'oxygène dissous DO, des nutriments et de la chlorophylle a au cours de l'année 2009. Les différentes simulations ont souligné l'impact négatif à long terme du réchauffement climatique. A la fin du siècle (2090), une augmentation de la température de l'air atteignant 4.4° C a été prédite à partir de modèles de circulation mondiale. Couplé avec une année hydrologique sèche, une telle situation pourrait conduire à des conditions anoxiques sévères et à des blooms importants de chlorophylle a (jusqu'à de 94 µg L-1). Diverses simulations ont montré que la réduction drastique de l'apport de nutriments (à 90%) serait nécessaire pour réduire de façon significative les teneurs en chlorophylle a. Si de telles mesures d'atténuation sont adoptées, le pic maximal de chlorophylle a se stabiliserait à 40 µg L-1, au lieu de 94 mg L-1, après une période de cinq ans d'efforts. À notre connaissance, cette étude fournit la première application numérique de modèles k-ε et AQUASIM pour simuler les niveaux d'eutrophisation élevés dans un réservoir tropical très turbide.