Sistemas de Detecção de Intrusões Análise de intrusões na rede multimédia da Unipiaget

Os sistemas computacionais ligados através da rede precisam, cada vez mais, de medidas e mecanismos de protecção que os permitem realizar as suas tarefas de acordo com os objectivos traçados pelos indivíduos ou organizações que os utilizam. Neste sentido, os Intrusion Detection System (IDS) ou Sistemas de Detecção de Intrusão desempenham um papel preponderante na monitorização contínua de actividades de um sistema pretendendo, essencialmente, garantir a autenticidade, confidencialidade, integridade e, sobretudo, a disponibilidade dos seus recursos. O presente trabalho apresenta um estudo sobre os Sistemas de Detecção de Intrusões, apresentando ainda um caso de estudo sobre a Rede Multimédia da Universidade Jean Piaget de Cabo Verde, estudo esse que incide sobre uma ferramenta Open Source, o Snort

Anàlisi geomètrica del Diaclasat de les Guilleries en el sector de la Vall del Ter entre els embassaments del Pasteral i Susqueda (Girona)

The Guilleries are a small and mountainous area located in the north-westem part of the Catalonian Coastal Ranges where metamorphic and igneous Paleozoic rocks are exposed. After the main hercynian folding this area was affected by a brittle deformation that is mainly manifested by the intrusion of a very large number of dykes of granodiorite and the development of a complex joint system. Trends of dykes indicate that their intrusion was related to a SE-NW extension, whose estimated value is 40% on an average. This extension seems to stand, although without any associated igneous event, with the development of NE-SW directed joints which make the main set. Five families more were developed later, onegently-dipping and fou upright; the latter trending roughly SE-NW, ENE-WSW, ESE-WNW and N-S. AU the joint sets appear in the metasedimentary Paleozoic rocks and in the hercynian intrusive bodies. Concerning the ages, joints that belong to the NE-SW and SE-NW directed sets and also those slightly dipping have been attributed to the late-hercynian times and all the other are considered to be later

Late Pan-African magmatism in the Himalaya: new geochronological and geochemical data from the Ordovician Tso Morari metagranites (Ladakh, NW India)

Two granitic plutons, the Tso Morari gneiss and the Rupshu metagranite, crop out in the Tso Morari area. The Polokongka La granite, classically interpreted as a young intrusion in the Tso Morari gneiss, has been recognized as the undeformed facies of the latter. Conventional isotope dilution U-Pb zircon dating on single-grain and small multi-grain fractions yielded magmatic ages of 479 +/- 2 Ma for the Tso Morari gneiss and the Polokongka La granite, and 482.5 +/- 1 Ma for the Rupshu granite. There is a great difference in zircon morphology between the Tso Morari gneiss (peraluminous type) and the Rupshu granite (alkaline type). This difference is confirmed by whole-rock chemistry. The Tso Morari gneiss is a typical deformed S-type granite, resulting from crustal anatexis. On the other hand, the Rupshu granite is an essentially metaluminous alkali-calcic intrusion derived from a different source material. Data compilation from other Himalayan Cambro-Ordovician granites reveals huge and widespread magmatic activity all along and beyond the northern Indian plate between 570 and 450 Ma, with a peak at 500-480 Ma. A major, continental-scale tectonic event is required to generate such a large magmatic belt; it has been tentatively compared to the Variscan post-orogenic extensional regime of Western Europe, as a late evolution stage of a Pan-African orogenic event.

Incremental growth of the Patagonian Torres del Paine laccolith over 90 k.y.

The Miocene Paine Granite in the Torres del Paine Intrusive Complex, southern Chile, is an extraordinary example of an upper crustal mafic and granitic intrusion. The granite intruded as a series of three sheets, each one underplating the previous sheet along the top of the basal Paine Mafic Complex. High-precision U/Pb geochronology on single zircons using isotope dilution-thermal ionization mass spectrometry yields distinct ages of 12.59 +/- 0.02 Ma and 12.50 +/- 0.02 Ma, respectively, for the first and last sheet of the laccolith. This age relationship is consistent with field observations. The zircon ages define a time frame of 90 +/- 40 k.y. for the emplacement of a >2000-m-thick granite laccollith. These precise U-Pb zircon ages permit identification of the pulses in a 20 k.y. range. The data obtained for the Paine Granite fill the gap between 100 k.y. and 100-1000 yr pulses described in the literature for crustal magma chambers.

The Palaeozoic metamorphic evolution of the Alpine External Massifs

The pre-Mesozoic metamorphic pattern of the External Massifs, composed of subunits of different metamorphic histories, resulted from the telescoping of Variscan, Ordovician and older metamorphic and structural textures and formations. During an early period, the future External Massifs were part of a peri-Gondwanian microplate evolving as an active margin. Precambrian to lower Palaeozoic igneous and sedimentary protoliths were reworked during an Ordovician subduction cycle (eclogites, granulites) preceding Ordovician anatexis and intrusion of Ordovician granitoids. Little is known about the time period when the microcontinent containing the future External Massifs followed a migration path leading to collision with Laurussia. Corresponding rock-series have not been identified. This might be because they have been eroded or transformed by migmatisation or because they remain hidden in the monocyclic areas. Besides the transformations which originated during the Ordovician subduction cycle, strong metamorphic transformations resulted from Variscan collision when many areas underwent amphibolite facies transformations and migmatisation. The different subunits composing the External Massifs and their corresponding P-T evolution are the expression of different levels in a nappe pile, which may have formed before Visean erosion and cooling. The presence of durbachitic magmatic rocks may be the expression of a large scale Early Variscan upwelling line which formed after Variscan lithospheric subduction. Late Variscan wrench fault tectonics and crustal thinning accompanied by high thermal gradients triggered several pulses of granite intrusions.

EEG Patterns and Imaging Correlations in Encephalopathy: Encephalopathy Part II.

SUMMARY:: The EEG patterns seen with encephalopathies can be correlated to cerebral imaging findings including head computerized tomography and MRI. Background slowing without slow-wave intrusion is seen with acute and chronic cortical impairments that spare subcortical white matter. Subcortical/white matter structural abnormalities or hydrocephalus may produce projected slow-wave activity, while clinical entities involving both cortical and subcortical regions (diffuse cerebral abnormalities) engender both background slowing and slow-wave activity. Triphasic waves are seen with hepatic and renal insufficiency or medication toxicities (e.g., lithium, baclofen) in the absence of a significant cerebral imaging abnormality, Conversely, subcortical/white matter abnormalities may facilitate the appearance of triphasic waves without significant hepatic, renal, or toxic comorbidities. More specific syndromes, such as Jakob-Creutzfeldt disease, autoimmune limbic encephalitis, autoimmune corticosteroid-responsive encephalopathy with thyroid autoimmunity, sepsis-associated encephalopathy, and acute disseminated encephalomyelitis, have imaging/EEG changes that are variable but which may include slowing and epileptiform activity. This overview highlighting EEG-imaging correlations may help the treating physician in the diagnosis, and hence the appropriate treatment, of patients with encephalopathy.

Synorogenic extension: Quantitative constraints on the age and displacement of the Zanskar shear zone (northwest Himalaya)

subsequent extension-induced exhumation. Geochronological dating of various Structural, thermobarometric, and geochronological data place limits on the age and tectonic displacement along the Zanskar shear zone, a major north-dipping synorogenic extensional structure separating the high-grade metamorphic sequence of the High Himalayan Crystalline Sequence from the overlying low-grade sedimentary rocks of the Tethyan Himalaya, A complete Barrovian metamorphic succession, from kyanite to biotite zone mineral assemblages, occurs within the I-km-thick Zanskar shear zone. Thermobarometric data indicate a difference In equilibration depths of 12 +/- 3 km between the lower kyanite zone and the garnet zone, which is Interpreted as a minimum estimate for the finite vertical displacement accommodated by the Zanskar shear zone. For the present-day dip of the structure (20 degrees), a simple geometrical model shows that a net slip of 35 +/- 9 km is required to regroup these samples to the same structural level. Because the kyanite to garnet zone rocks represent only part of the Zanskar shear zone, and because its original dip may have been less than the present-day dip, these estimates fur the finite displacement represent minimum values. Field relations and petrographic data suggest that migmatization and associated leucogranite intrusion in the footwall of the Zanskar shear zone occurred as a continuous profess starting at the Barrovian metamorphic peak and lasting throughout the subsequent extension-induced exhumation. Geochronological dataing of various leucogranitic plutons and dikes in the Zanskar shear zone footwall indicates that the main ductile shearing along the structure ended by 19.8 Ma and that extension most likely initiated shortly before 22.2 Ma.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

The Variscan evolution in the External massifs of the Alps and place in their Variscan framework

In the general discussion on the Variscan evolution of central Europe the pre-Mesozoic basement of the Alps is, in many cases, only included with hesitation. Relatively well-preserved from Alpine metamorphism, the Alpine External massifs can serve as an excellent example of evolution of the Variscan basement, including the earliest Gondwana-derived microcontinents with Cadomian relics. Testifying to the evolution at the Gondwana margin, at least since the Cambrian, such pieces took part in the birth of the Rheic Ocean. After the separation of Avalonia, the remaining Gondwana border was continuously transformed through crustal extension with contemporaneous separation of continental blocks composing future Pangea, but the opening of Palaeotethys had only a reduced significance since the Devonian. The Variscan evolution in the External domain is characterised by an early HP-evolution with subsequent granulitic decompression melts. During Visean crustal shortening, the areas of future formation of migmatites and intrusion of monzodioritic magmas in a general strike-slip regime, were probably in a lower plate situation, whereas the so called monometamorphic areas may have been in an upper plate position of the nappe pile. During the Latest Carboniferous, the emplacement of the youngest granites was associated with the strike-slip faulting and crustal extension at lower crustal levels, whereas, at the surface, detrital sediments accumulated in intramontaneous transtensional basins on a strongly eroded surface. To cite this article: J.R von Raumer et al., C. R. Geoscience 341 (2009). (C) 2008 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Flank stability and processes off the western Canary Islands: a review from El Hierro and La Palma

The morphological characterisation of the western submarine island flanks of El Hierro and La Palma differentiates four type-zones that may give new insights into the evolution of oceanic island slopes. The different type-zones result from the interplay between constructive volcanic processes, hemipelagic settling and volcano collapses. The latter results in massive debris avalanche deposits, which form large volcaniclastic aprons. In most cases, the headwall scarps are clearly exposed on the emerged part of the islands. The events that occurred in the youngest and westernmost islands of El Hierro and La Palma have vertical runouts exceeding 6,000 m and volumes that can reach several hundred km3. The landslide frequency for the entire Canaries is one major event per 90 ka. Triggering mechanisms are closely related to magmatic processes. The increase in the shear stress is directly linked with the forceful intrusion of magma along ridge-rift systems, while in the western Canary Islands it seems that the main process reducing shear resistance may be related to the rise in pore pressure due to hydrothermal circulation.

[Illustrations de Histoires ou contes du temps passé avec des moralités] / [Non identifié] ; Charles Perrault, aut. du texte

Comprend : [Frontispice : Devant une cheminée, une vieille dame fait le récit de contes à trois enfants.] [Cote : BnF-Impr. microfilm R.24829] ; [Pl. en reg. p. 1 : La rencontre dans la forêt du Petit Chaperon rouge et du Loup.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Le Petit Chaperon rouge, p. 1 : Le Loup dévorant Mère-grand.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Les Fées, p. 5 : Jeune fille donnant à boire à une fée déguisée en vieille femme.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte La Barbe bleue, p. 10 : Les deux frères de la dernière épouse de Barbe bleue viennent à la rescousse de celle-ci et terrassent d'un coup d'épée le mari assassin.] [Cote : BnF-Impr. microfilm R. 24829] ; [Bandeau au titre du conte La Belle au bois dormant, p. 18 : Le jeune prince découvre la princesse endormie.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Le Maître Chat ou Le Chat botté, p. 32 : Le Chat botté montrant le château du Marquis de Carabas à son maître accompagné de la Princesse du royaume et d'une escorte royale.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Cendrillon ou La petite pantoufle de verre, p. 40 : Cendrillon quitte précipitamment la salle de bal.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Riquet à la Houpe, p. 51 : Dans la forêt, le prince Riquet à la Houpe demande sa main à la Princesse. En arrière plan les gens du prince s'activent à la préparation du banquet de noce.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte Le Petit Poucet, p. 62 : Le Petit Poucet retire les bottes de sept lieues au géant endormi.] [Cote : BnF-Impr. microfilm R.24829] ; [Bandeau au titre du conte L'Adroite princesse ou les Aventures de Finette, p. 77 : Intrusion de Riche-Cautelle, prince indigne, dans les appartements des trois princesses.] [Cote : BnF-Impr. microfilm R.24829]

Lifetime of an ocean island volcano feeder zone: constraints from U-Pb dating on coexisting zircon and baddeleyite, and 40Ar/39Ar age determinations, Fuerteventura, Canary Islands

High-precision isotope dilution - thermal ionization mass spectrometry (ID-TIMS) U-Pb zircon and baddeleyite ages from the PX1 vertically layered mafic intrusion Fuerteventura, Canary Islands, indicate initiation of magma crystallization at 22.10 +/- 0.07 Ma. The magmatic activity lasted a minimum of 0.52 Ma. Ar-40/Ar-39 amphibole dating yielded ages from 21.9 +/- 0.6 to 21.8 +/- 0.3, identical within errors to the U-Pb ages, despite the expected 1% theoretical bias between Ar-40/Ar-39 and U-Pb dates. This overlap could result from (i) rapid cooling of the intrusion (i. e., less than the 0.3 to 0.6 Ma 40Ar/39Ar age uncertainties) from closure temperatures (T-c) of zircon (699-988 degrees C) to amphibole (500-600 degrees C); (ii) lead loss affecting the youngest zircons; or (iii) excess argon shifting the plateau ages towards older values. The combination of the Ar-40/Ar-39 and U/Pb datasets implies that the maximum amount of time PX1 intrusion took to cool below amphibole T-c is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 Ma. Age disparities among coexisting baddeleyite and zircon (22.10 +/- 0.07/0.08/0.15 Ma and 21.58 +/- 0.15/0.16/0.31 Ma) in a gabbro sample from the pluton margin suggest complex genetic relationships between phases. Baddeleyite is found preserved in plagioclase cores and crystallized early from low silica activity magma. Zircon crystallized later in a higher silica activity environment and is found in secondary scapolite and is found close to calcite veins, in secondary scapolite that recrystallised from plagioclase. close to calcite veins. Oxygen isotope delta O-18 values of altered plagioclase are high (+7.7), indicating interaction with fluids derived from host-rock carbonatites. The coexistence of baddeleyite and zircon is ascribed to interaction of the PX1 gabbro with CO2-rich carbonatite-derived fluids released during contact metamorphism.

Multiscale descriptions of density-driven flow instabilities in porous media

Les instabilités engendrées par des gradients de densité interviennent dans une variété d'écoulements. Un exemple est celui de la séquestration géologique du dioxyde de carbone en milieux poreux. Ce gaz est injecté à haute pression dans des aquifères salines et profondes. La différence de densité entre la saumure saturée en CO2 dissous et la saumure environnante induit des courants favorables qui le transportent vers les couches géologiques profondes. Les gradients de densité peuvent aussi être la cause du transport indésirable de matières toxiques, ce qui peut éventuellement conduire à la pollution des sols et des eaux. La gamme d'échelles intervenant dans ce type de phénomènes est très large. Elle s'étend de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères à laquelle interviennent les phénomènes à temps long. Une reproduction fiable de la physique par la simulation numérique demeure donc un défi en raison du caractère multi-échelles aussi bien au niveau spatial et temporel de ces phénomènes. Il requiert donc le développement d'algorithmes performants et l'utilisation d'outils de calculs modernes. En conjugaison avec les méthodes de résolution itératives, les méthodes multi-échelles permettent de résoudre les grands systèmes d'équations algébriques de manière efficace. Ces méthodes ont été introduites comme méthodes d'upscaling et de downscaling pour la simulation d'écoulements en milieux poreux afin de traiter de fortes hétérogénéités du champ de perméabilité. Le principe repose sur l'utilisation parallèle de deux maillages, le premier est choisi en fonction de la résolution du champ de perméabilité (grille fine), alors que le second (grille grossière) est utilisé pour approximer le problème fin à moindre coût. La qualité de la solution multi-échelles peut être améliorée de manière itérative pour empêcher des erreurs trop importantes si le champ de perméabilité est complexe. Les méthodes adaptatives qui restreignent les procédures de mise à jour aux régions à forts gradients permettent de limiter les coûts de calculs additionnels. Dans le cas d'instabilités induites par des gradients de densité, l'échelle des phénomènes varie au cours du temps. En conséquence, des méthodes multi-échelles adaptatives sont requises pour tenir compte de cette dynamique. L'objectif de cette thèse est de développer des algorithmes multi-échelles adaptatifs et efficaces pour la simulation des instabilités induites par des gradients de densité. Pour cela, nous nous basons sur la méthode des volumes finis multi-échelles (MsFV) qui offre l'avantage de résoudre les phénomènes de transport tout en conservant la masse de manière exacte. Dans la première partie, nous pouvons démontrer que les approximations de la méthode MsFV engendrent des phénomènes de digitation non-physiques dont la suppression requiert des opérations de correction itératives. Les coûts de calculs additionnels de ces opérations peuvent toutefois être compensés par des méthodes adaptatives. Nous proposons aussi l'utilisation de la méthode MsFV comme méthode de downscaling: la grille grossière étant utilisée dans les zones où l'écoulement est relativement homogène alors que la grille plus fine est utilisée pour résoudre les forts gradients. Dans la seconde partie, la méthode multi-échelle est étendue à un nombre arbitraire de niveaux. Nous prouvons que la méthode généralisée est performante pour la résolution de grands systèmes d'équations algébriques. Dans la dernière partie, nous focalisons notre étude sur les échelles qui déterminent l'évolution des instabilités engendrées par des gradients de densité. L'identification de la structure locale ainsi que globale de l'écoulement permet de procéder à un upscaling des instabilités à temps long alors que les structures à petite échelle sont conservées lors du déclenchement de l'instabilité. Les résultats présentés dans ce travail permettent d'étendre les connaissances des méthodes MsFV et offrent des formulations multi-échelles efficaces pour la simulation des instabilités engendrées par des gradients de densité. - Density-driven instabilities in porous media are of interest for a wide range of applications, for instance, for geological sequestration of CO2, during which CO2 is injected at high pressure into deep saline aquifers. Due to the density difference between the C02-saturated brine and the surrounding brine, a downward migration of CO2 into deeper regions, where the risk of leakage is reduced, takes place. Similarly, undesired spontaneous mobilization of potentially hazardous substances that might endanger groundwater quality can be triggered by density differences. Over the last years, these effects have been investigated with the help of numerical groundwater models. Major challenges in simulating density-driven instabilities arise from the different scales of interest involved, i.e., the scale at which instabilities are triggered and the aquifer scale over which long-term processes take place. An accurate numerical reproduction is possible, only if the finest scale is captured. For large aquifers, this leads to problems with a large number of unknowns. Advanced numerical methods are required to efficiently solve these problems with today's available computational resources. Beside efficient iterative solvers, multiscale methods are available to solve large numerical systems. Originally, multiscale methods have been developed as upscaling-downscaling techniques to resolve strong permeability contrasts. In this case, two static grids are used: one is chosen with respect to the resolution of the permeability field (fine grid); the other (coarse grid) is used to approximate the fine-scale problem at low computational costs. The quality of the multiscale solution can be iteratively improved to avoid large errors in case of complex permeability structures. Adaptive formulations, which restrict the iterative update to domains with large gradients, enable limiting the additional computational costs of the iterations. In case of density-driven instabilities, additional spatial scales appear which change with time. Flexible adaptive methods are required to account for these emerging dynamic scales. The objective of this work is to develop an adaptive multiscale formulation for the efficient and accurate simulation of density-driven instabilities. We consider the Multiscale Finite-Volume (MsFV) method, which is well suited for simulations including the solution of transport problems as it guarantees a conservative velocity field. In the first part of this thesis, we investigate the applicability of the standard MsFV method to density- driven flow problems. We demonstrate that approximations in MsFV may trigger unphysical fingers and iterative corrections are necessary. Adaptive formulations (e.g., limiting a refined solution to domains with large concentration gradients where fingers form) can be used to balance the extra costs. We also propose to use the MsFV method as downscaling technique: the coarse discretization is used in areas without significant change in the flow field whereas the problem is refined in the zones of interest. This enables accounting for the dynamic change in scales of density-driven instabilities. In the second part of the thesis the MsFV algorithm, which originally employs one coarse level, is extended to an arbitrary number of coarse levels. We prove that this keeps the MsFV method efficient for problems with a large number of unknowns. In the last part of this thesis, we focus on the scales that control the evolution of density fingers. The identification of local and global flow patterns allows a coarse description at late times while conserving fine-scale details during onset stage. Results presented in this work advance the understanding of the Multiscale Finite-Volume method and offer efficient dynamic multiscale formulations to simulate density-driven instabilities. - Les nappes phréatiques caractérisées par des structures poreuses et des fractures très perméables représentent un intérêt particulier pour les hydrogéologues et ingénieurs environnementaux. Dans ces milieux, une large variété d'écoulements peut être observée. Les plus communs sont le transport de contaminants par les eaux souterraines, le transport réactif ou l'écoulement simultané de plusieurs phases non miscibles, comme le pétrole et l'eau. L'échelle qui caractérise ces écoulements est définie par l'interaction de l'hétérogénéité géologique et des processus physiques. Un fluide au repos dans l'espace interstitiel d'un milieu poreux peut être déstabilisé par des gradients de densité. Ils peuvent être induits par des changements locaux de température ou par dissolution d'un composé chimique. Les instabilités engendrées par des gradients de densité revêtent un intérêt particulier puisque qu'elles peuvent éventuellement compromettre la qualité des eaux. Un exemple frappant est la salinisation de l'eau douce dans les nappes phréatiques par pénétration d'eau salée plus dense dans les régions profondes. Dans le cas des écoulements gouvernés par les gradients de densité, les échelles caractéristiques de l'écoulement s'étendent de l'échelle poreuse où les phénomènes de croissance des instabilités s'opèrent, jusqu'à l'échelle des aquifères sur laquelle interviennent les phénomènes à temps long. Etant donné que les investigations in-situ sont pratiquement impossibles, les modèles numériques sont utilisés pour prédire et évaluer les risques liés aux instabilités engendrées par les gradients de densité. Une description correcte de ces phénomènes repose sur la description de toutes les échelles de l'écoulement dont la gamme peut s'étendre sur huit à dix ordres de grandeur dans le cas de grands aquifères. Il en résulte des problèmes numériques de grande taille qui sont très couteux à résoudre. Des schémas numériques sophistiqués sont donc nécessaires pour effectuer des simulations précises d'instabilités hydro-dynamiques à grande échelle. Dans ce travail, nous présentons différentes méthodes numériques qui permettent de simuler efficacement et avec précision les instabilités dues aux gradients de densité. Ces nouvelles méthodes sont basées sur les volumes finis multi-échelles. L'idée est de projeter le problème original à une échelle plus grande où il est moins coûteux à résoudre puis de relever la solution grossière vers l'échelle de départ. Cette technique est particulièrement adaptée pour résoudre des problèmes où une large gamme d'échelle intervient et évolue de manière spatio-temporelle. Ceci permet de réduire les coûts de calculs en limitant la description détaillée du problème aux régions qui contiennent un front de concentration mobile. Les aboutissements sont illustrés par la simulation de phénomènes tels que l'intrusion d'eau salée ou la séquestration de dioxyde de carbone.

Magmatic-dominated fluid evolution in the Jurassic Nambija gold skarn deposits (southeastern Ecuador)

The Jurassic (approximately 145 Ma) Nambija oxidized gold skarns are hosted by the Triassic volcanosedimentary Piuntza unit in the sub-Andean zone of southeastern Ecuador. The skarns consist dominantly of granditic garnet (Ad(20-98)) with subordinate pyroxene (Di(46-92)Hd(17-42)Jo(0-19)) and epidote and are spatially associated with porphyritic quartz-diorite to granodiorite intrusions. Endoskarn is developed at the intrusion margins and grades inwards into a potassic alteration zone. Exoskarn has an outer K- and Na-enriched zone in the volcanosedimentary unit. Gold mineralization is associated with the weakly developed retrograde alteration of the exoskarn and occurs mainly in sulfide-poor vugs and milky quartz veins and veinlets in association with hematite. Fluid inclusion data for the main part of the prograde stage indicate the coexistence of high-temperature (500A degrees C to > 600A degrees C), high-salinity (up to 65 wt.% eq. NaCl), and moderate- to low-salinity aqueous-carbonic fluids interpreted to have been trapped at pressures around 100-120 MPa, corresponding to about 4-km depth. Lower-temperature (510-300A degrees C) and moderate- to low-salinity (23-2 wt.% eq. NaCl) aqueous fluids are recorded in garnet and epidote of the end of the prograde stage. The microthermometric data (Th from 513A degrees C to 318A degrees C and salinity from 1.0 to 23 wt.% eq. NaCl) and delta(18)O values between 6.2aEuro degrees and 11.5aEuro degrees for gold-bearing milky quartz from the retrograde stage suggest that the ore-forming fluid was dominantly magmatic. Pressures during the early retrograde stage were in the range of 50-100 MPa, in line with the evidence for CO(2) effervescence and probable local boiling. The dominance of magmatic low-saline to moderately saline oxidizing fluids during the retrograde stage is consistent with the depth of the skarn system, which could have delayed the ingression of external fluids until relatively low temperatures were reached. The resulting low water-to-rock ratios explain the weak retrograde alteration and the compositional variability of chlorite, essentially controlled by host rock compositions. Gold was precipitated at this stage as a result of cooling and pH increase related to CO(2) effervescence, which both result in destabilization of gold-bearing chloride complexes. Significant ingression of external fluids took place after gold deposition only, as recorded by delta(18)O values of 0.4aEuro degrees to 6.2aEuro degrees for fluids depositing quartz (below 350A degrees C) in sulfide-rich barren veins. Low-temperature (< 300A degrees C) meteoric fluids (delta(18)O(water) between -10.0aEuro degrees and -2.0aEuro degrees) are responsible for the precipitation of late comb quartz and calcite in cavities and veins and indicate mixing with cooler fluids of higher salinities (about 100A degrees C and 25 wt.% eq. NaCl). The latter are similar to low-temperature fluids (202-74.5A degrees C) with delta(18)O values of -0.5aEuro degrees to 3.1aEuro degrees and salinities in the range of 21.1 to 17.3 wt.% eq. CaCl(2), trapped in calcite of late veins and interpreted as basinal brines. Nambija represents a deep equivalent of the oxidized gold skarn class, the presence of CO(2) in the fluids being partly a consequence of the relatively deep setting at about 4-km depth. As in other Au-bearing skarn deposits, not only the prograde stage but also the gold-precipitating retrograde stage is dominated by fluids of magmatic origin.

Evaluation of Dense Bridge Floor Concrete Using High Range Water Reducer, May 1983

The objectives of this research project are: (1) To determine the feasibility of proportioning, mixing, placing and finishing a dense portland cement concrete in a bridge floor using conventional mixing, placing and finishing equipment. (2) To determine the economics, longevity, maintenance performance and protective qualities of a dense portland cement concrete bridge floor when using a high rangewater reducing admixture. The purpose of a high range water reducing admixture is to produce a dense, high quality concrete at a low water-cement ratio witj adequate workability. A low water-cement ratio contributes greatly to increased strength. The normal 7 day strength of untreated concrete would be expected i n 3 days using a superplasticizer. A dense concrete also has the desirable properties of excellent durability and reduced permeability. It is felt that a higher quality, denser, higher strength portland cement concrete can be produced and placed, using conventional equipment, by the addition of a high range water reducing admixture. Such a dense concrete, w i t h a water/cement ratio of approximately 0.30 to 0.35, would be expected to be much less permeable and thus retard the intrusion of chloride. With care and attention given to obtaining the design cover over steel (2% inches clear), it i s hoped that protection for the design life of the structure will be obtained. Evaluation of this experimental concrete bridge floor included chloride content and delamination testing of the concrete floor five years after construction. A comparitive evaluation o f a control section o f concrete without the water reducing admixture was conducted. Other items o f comparison include workability during construction, strength, density, water-cement ratio and chloride penetration.