La geoquímica de los sedimentos de llanuras de inundación como guía para la valoración ambiental en España

En la presente tesis se estudian los contenidos geoquímicos de los sedimentos de llanura de inundación en diversas cuencas fluviales seleccionadas con el objetivo de contribuir a un mejor conocimiento de sus condiciones medioambientales. En cada cuenca, se ha muestreado un perfil vertical de llanura de inundación, dividiéndolo en tramos que generalmente corresponden a diferentes episodios de inundación. Estos sedimentos se depositan durante los episodios de crecidas una vez que la corriente sobrepasa los límites del canal. Se caracterizan normalmente por tener un tamaño de grano fino y una estructura en capas horizontales que corresponden a los sucesivos episodios de inundación. Muestran dos ventajas principales con respecto a otros medios de muestreo en geoquímica como suelos o sedimentos de corriente: • Pueden almacenar sedimento antiguo así como actual, con lo que se puede estudiar la historia geoquímica de una zona específica en un mismo punto de muestreo (perfil vertical). • Los sedimentos de llanuras de inundación son capaces de caracterizar grandes áreas de drenaje. El origen de los sedimentos es más diverso que en el sedimento de corriente, debido a las mayores áreas de donde proceden las aguas de las avenidas. Las cuencas han sido seleccionadas según las actividades antropogénicas que en ellas se llevan a cabo, en concreto, actividades urbanas e industriales, minería y agricultura. Así mismo, se han estudiado, como referencia, dos cuencas donde no se espera encontrar ningún tipo de actividad contaminante. Una vez hecha la selección, los sedimentos aluviales de las cuencas se han estudiado cuidadosamente para asegurar que no existen depósitos de acreción lateral en el punto seleccionado. Posteriormente se ha procedido al muestreo del perfil vertical. Las muestras han sido analizadas mediante ICP-MS (ataque total) e INAA para conocer los contenidos totales de los elementos traza y mayoritarios. Los análisis de la fracción extraíble se han llevado a cabo mediante ICP-MS (ataque con agua regia). Así mismo, algunas muestras seleccionadas han sido sometidas a una extracción secuencial para un estudio más detallado. La presencia de materia orgánica ha sido estimada mediante el análisis de Carbono Orgánico Total (TOC). Finalmente, se ha llevado a cabo un análisis de isótopos de Pb en muestras escogidas en los perfiles, con el objetivo de hacer una evaluación ambiental. Los contenidos metálicos aumentan hacia la superficie en algunos de los perfiles, mientras en otros muestran una distribución muy constante exceptuando algún nivel específico con un aumento de los contenidos de la mayoría de los metales. Ha sido posible determinar la influencia de las actividades antropogénicas en algunos de los perfiles. Aquellos que pertenecen a cuencas mineras, urbanas o industrializadas muestran generalmente altos contenidos en elementos metálicos. Es el caso de los perfiles muestreados en los ríos Odiel y Tinto, Besaya, Besós y Manzanares. Algunos de estos perfiles pueden incluso correlacionarse con periodos de tiempo en los que ha tenido lugar una actividad antropogénica más intensa. Los perfiles que mejor se correlacionan con la actividad antropogénica de la cuenca son el perfil de Rivas en el río Manzanares (Madrid), que refleja un crecimiento de la contaminación producida por las actividades urbana e industrial en las últimas décadas, y el río Tinto, que muestra un crecimiento llamativo de los contenidos en su mayoría metálicos que puede estar relacionado con el incremento de la actividad minera que tuvo lugar hace aproximadamente 125 años. El análisis de los isótopos de Pb ha resultado ser una herramienta útil en la evaluación ambiental de estos sedimentos. Con este estudio y mediante la comparación con fuentes naturales y antropogénicas, ha sido posible diferenciar las muestras afectadas por diferentes fuentes de plomo, así como detectar las más afectadas antropogénicamente. ABSTRACT The geochemical composition of overbank sediments of some selected river basins is studied in this thesis in order to contribute to a better knowledge of the environmental conditions surrounding them. In each basin a vertical overbank profile has been sampled, dividing it into stretches that usually correspond to different flood events. The overbank sediments are those deposited during a flood event once the flow spills over the channel banks. They are usually characterized by a very fine grain size and a structure of horizontal layers, which correspond to successive flood events. These sediments show two main advantages regarding other sampling media in geochemistry, like soils or stream sediments: • They can store sediment deposited in the past as well as in current times, so that the history of a specific location can be studied at the very same point (vertical profile). • The overbank sediments are able to characterize a large drainage area. The origin of the sediment is wider than in the stream sediments due to the larger areas where the flood water comes from. The basins have been selected depending on the anthropogenic activities developed in them, namely, urban and industrial activities, mining activities and agricultural activities. As well, two pristine basins have been studied as a reference. Afterwards, the alluvial sediments in the basins have been carefully studied in order to sample a vertical profile and make sure that lateral accretion materials are not present in the profile. The samples have been analysed by ICP-MS (total digestion) and INAA to know the total contents of trace and major elements. Analysis of the mobile fraction has been carried out by ICP-MS (aqua regia); as well some of the samples have been subjected to sequential extraction for a more detailed study. The presence of organic matter has been estimated by the analysis of the Total Organic Carbon (TOC). Finally, a lead isotope analysis of some of the samples in the profiles was carried out in order to make an environmental assessment. Metal contents grow towards the surface in some of the profiles, while others show a very steady distribution, except for some of them with a growth of most of the metals in a specific level. XI It has been possible to determine the influence of the anthropogenic activities in some of the profiles. The ones that belong to mining and urban or industrialized basins show generally high contents of metal elements. This is the case of the profiles sampled in the Odiel and Tinto Rivers, the Besaya River, the Besós River and the Manzanares River. Some of these profiles can even correlate with the periods of time when a more intense activity in their respective basins has taken place. The profiles which best correlate with the anthropogenic activity are the Rivas profile in the Manzanares River, which reflects a growth of the pollution produced by urban and industrial activities in the city of Madrid in the last decades and the Tinto profile, which shows a very dramatic growth of the elemental contents (mostly metals) which can be related to the increase of the mining activities that took place in the last 125 years. The analysis of lead isotopes has turned out to be a powerful tool in the environmental assessment in this type of sediments. With this study and through the comparison with natural and anthropogenic sources, it has been possible to determine samples affected by different sources of lead and to detect the most anthropogenicaly affected ones.

Estrategias de monitorización de CO2 y otros gases en los estudios de análogos naturales

La Universidad Politécnica de Madrid (UPM) y la Università degli Studi di Firenze (UniFi), bajo la coordinación técnica de AMPHOS21, participan desde 2009 en el proyecto de investigación “Estrategias de Monitorización de CO2 y otros gases en el estudio de Análogos Naturales”, financiado por la Fundación Ciudad de la Energía (CIUDEN) en el marco del Proyecto Compostilla OXYCFB300 (http://www.compostillaproject.eu), del Programa “European Energy Program for Recovery - EEPR”. El objetivo principal del proyecto fue el desarrollo y puesta a punto de metodologías de monitorización superficiales para su aplicación en el seguimiento y control de los emplazamientos donde se realice el almacenamiento geológico de CO2, analizando técnicas que permitan detectar y cuantificar las posibles fugas de CO2 a la atmósfera. Los trabajos se realizaron tanto en análogos naturales (españoles e italianos) como en la Planta de Desarrollo Tecnológico de Almacenamiento de CO2 de Hontomín. Las técnicas analizadas se centran en la medición de gases y aguas superficiales (de escorrentía y manantiales). En cuanto a la medición de gases se analizó el flujo de CO2 que emana desde el suelo a la atmósfera y la aplicabilidad de trazadores naturales (como el radón) para la detección e identificación de las fugas de CO2. En cuanto al análisis químico de las aguas se analizaron los datos geoquímicos e isotópicos y los gases disueltos en las aguas de los alrededores de la PDT de Hontomín, con objeto de determinar qué parámetros son los más apropiados para la detección de una posible migración del CO2 inyectado, o de la salmuera, a los ambientes superficiales. Las medidas de flujo de CO2 se realizaron con la técnica de la cámara de acúmulo. A pesar de ser una técnica desarrollada y aplicada en diferentes ámbitos científicos se estimó necesario adaptar un protocolo de medida y de análisis de datos a las características específicas de los proyectos de captura y almacenamiento de CO2 (CAC). Donde los flujos de CO2 esperados son bajos y en caso de producirse una fuga habrá que detectar pequeñas variaciones en los valores flujo con un “ruido” en la señal alto, debido a actividad biológica en el suelo. La medida de flujo de CO2 mediante la técnica de la cámara de acúmulo se puede realizar sin limpiar la superficie donde se coloca la cámara o limpiando y esperando al reequilibrio del flujo después de la distorsión al sistema. Sin embargo, los resultados obtenidos después de limpiar y esperar muestran menor dispersión, lo que nos indica que este procedimiento es el mejor para la monitorización de los complejos de almacenamiento geológico de CO2. El protocolo de medida resultante, utilizado para la obtención de la línea base de flujo de CO2 en Hontomín, sigue los siguiente pasos: a) con una espátula se prepara el punto de medición limpiando y retirando el recubrimiento vegetal o la primera capa compacta de suelo, b) se espera un tiempo para la realización de la medida de flujo, facilitando el reequilibrio del flujo del gas tras la alteración provocada en el suelo y c) se realiza la medida de flujo de CO2. Una vez realizada la medición de flujo de CO2, y detectada si existen zonas de anomalías, se debe estimar la cantidad de CO2 que se está escapando a la atmósfera (emanación total), con el objetivo de cuantificar la posible fuga. Existen un amplio rango de metodologías para realizar dicha estimación, siendo necesario entender cuáles son las más apropiadas para obtener el valor más representativo del sistema. En esta tesis se comparan seis técnicas estadísticas: media aritmética, estimador insegado de la media (aplicando la función de Sichel), remuestreo con reemplazamiento (bootstrap), separación en diferentes poblaciones mediante métodos gráficos y métodos basados en criterios de máxima verosimilitud, y la simulación Gaussiana secuencial. Para este análisis se realizaron ocho campañas de muestreo, tanto en la Planta de Desarrollo Tecnológico de Hontomón como en análogos naturales (italianos y españoles). Los resultados muestran que la simulación Gaussiana secuencial suele ser el método más preciso para realizar el cálculo, sin embargo, existen ocasiones donde otros métodos son más apropiados. Como consecuencia, se desarrolla un procedimiento de actuación para seleccionar el método que proporcione el mejor estimador. Este procedimiento consiste, en primer lugar, en realizar un análisis variográfico. Si existe una autocorrelación entre los datos, modelizada mediante el variograma, la mejor técnica para calcular la emanación total y su intervalo de confianza es la simulación Gaussiana secuencial (sGs). Si los datos son independientes se debe comprobar la distribución muestral, aplicando la media aritmética o el estimador insesgado de la media (Sichel) para datos normales o lognormales respectivamente. Cuando los datos no son normales o corresponden a una mezcla de poblaciones la mejor técnica de estimación es la de remuestreo con reemplazamiento (bootstrap). Siguiendo este procedimiento el máximo valor del intervalo de confianza estuvo en el orden del ±20/25%, con la mayoría de valores comprendidos entre ±3,5% y ±8%. La identificación de las diferentes poblaciones muestrales en los datos de flujo de CO2 puede ayudar a interpretar los resultados obtenidos, toda vez que esta distribución se ve afectada por la presencia de varios procesos geoquímicos como, por ejemplo, una fuente geológica o biológica del CO2. Así pues, este análisis puede ser una herramienta útil en el programa de monitorización, donde el principal objetivo es demostrar que no hay fugas desde el reservorio a la atmósfera y, si ocurren, detectarlas y cuantificarlas. Los resultados obtenidos muestran que el mejor proceso para realizar la separación de poblaciones está basado en criterios de máxima verosimilitud. Los procedimientos gráficos, aunque existen pautas para realizarlos, tienen un cierto grado de subjetividad en la interpretación de manera que los resultados son menos reproducibles. Durante el desarrollo de la tesis se analizó, en análogos naturales, la relación existente entre el CO2 y los isótopos del radón (222Rn y 220Rn), detectándose en todas las zonas de emisión de CO2 una relación positiva entre los valores de concentración de 222Rn en aire del suelo y el flujo de CO2. Comparando la concentración de 220Rn con el flujo de CO2 la relación no es tan clara, mientras que en algunos casos aumenta en otros se detecta una disminución, hecho que parece estar relacionado con la profundidad de origen del radón. Estos resultados confirmarían la posible aplicación de los isótopos del radón como trazadores del origen de los gases y su aplicación en la detección de fugas. Con respecto a la determinación de la línea base de flujo CO2 en la PDT de Hontomín, se realizaron mediciones con la cámara de acúmulo en las proximidades de los sondeos petrolíferos, perforados en los ochenta y denominados H-1, H-2, H-3 y H-4, en la zona donde se instalarán el sondeo de inyección (H-I) y el de monitorización (H-A) y en las proximidades de la falla sur. Desde noviembre de 2009 a abril de 2011 se realizaron siete campañas de muestreo, adquiriéndose más de 4.000 registros de flujo de CO2 con los que se determinó la línea base y su variación estacional. Los valores obtenidos fueron bajos (valores medios entre 5 y 13 g•m-2•d-1), detectándose pocos valores anómalos, principalmente en las proximidades del sondeo H-2. Sin embargo, estos valores no se pudieron asociar a una fuente profunda del CO2 y seguramente estuvieran más relacionados con procesos biológicos, como la respiración del suelo. No se detectaron valores anómalos cerca del sistema de fracturación (falla Ubierna), toda vez que en esta zona los valores de flujo son tan bajos como en el resto de puntos de muestreo. En este sentido, los valores de flujo de CO2 aparentemente están controlados por la actividad biológica, corroborado al obtenerse los menores valores durante los meses de otoño-invierno e ir aumentando en los periodos cálidos. Se calcularon dos grupos de valores de referencia, el primer grupo (UCL50) es 5 g•m-2•d-1 en las zonas no aradas en los meses de otoño-invierno y 3,5 y 12 g•m-2•d-1 en primavera-verano para zonas aradas y no aradas, respectivamente. El segundo grupo (UCL99) corresponde a 26 g•m-2•d- 1 durante los meses de otoño-invierno en las zonas no aradas y 34 y 42 g•m-2•d-1 para los meses de primavera-verano en zonas aradas y no aradas, respectivamente. Flujos mayores a estos valores de referencia podrían ser indicativos de una posible fuga durante la inyección y posterior a la misma. Los primeros datos geoquímicos e isotópicos de las aguas superficiales (de escorrentía y de manantiales) en el área de Hontomín–Huermeces fueron analizados. Los datos sugieren que las aguas estudiadas están relacionadas con aguas meteóricas con un circuito hidrogeológico superficial, caracterizadas por valores de TDS relativamente bajos (menor a 800 mg/L) y una fácie hidrogeoquímica de Ca2+(Mg2+)-HCO3 −. Algunas aguas de manantiales se caracterizan por concentraciones elevadas de NO3 − (concentraciones de hasta 123 mg/l), lo que sugiere una contaminación antropogénica. Se obtuvieron concentraciones anómalas de of Cl−, SO4 2−, As, B y Ba en dos manantiales cercanos a los sondeos petrolíferos y en el rio Ubierna, estos componentes son probablemente indicadores de una posible mezcla entre los acuíferos profundos y superficiales. El estudio de los gases disueltos en las aguas también evidencia el circuito superficial de las aguas. Estando, por lo general, dominado por la componente atmosférica (N2, O2 y Ar). Sin embargo, en algunos casos el gas predominante fue el CO2 (con concentraciones que llegan al 63% v/v), aunque los valores isotópicos del carbono (<-17,7 ‰) muestran que lo más probable es que esté relacionado con un origen biológico. Los datos geoquímicos e isotópicos de las aguas superficiales obtenidos en la zona de Hontomín se pueden considerar como el valor de fondo con el que comparar durante la fase operacional, la clausura y posterior a la clausura. En este sentido, la composición de los elementos mayoritarios y traza, la composición isotópica del carbono del CO2 disuelto y del TDIC (Carbono inorgánico disuelto) y algunos elementos traza se pueden considerar como parámetros adecuados para detectar la migración del CO2 a los ambientes superficiales. ABSTRACT Since 2009, a group made up of Universidad Politécnica de Madrid (UPM; Spain) and Università degli Studi Firenze (UniFi; Italy) has been taking part in a joint project called “Strategies for Monitoring CO2 and other Gases in Natural analogues”. The group was coordinated by AMPHOS XXI, a private company established in Barcelona. The Project was financially supported by Fundación Ciudad de la Energía (CIUDEN; Spain) as a part of the EC-funded OXYCFB300 project (European Energy Program for Recovery -EEPR-; www.compostillaproject.eu). The main objectives of the project were aimed to develop and optimize analytical methodologies to be applied at the surface to Monitor and Verify the feasibility of geologically stored carbon dioxide. These techniques were oriented to detect and quantify possible CO2 leakages to the atmosphere. Several investigations were made in natural analogues from Spain and Italy and in the Tecnchnological Development Plant for CO2 injection al Hontomín (Burgos, Spain). The studying techniques were mainly focused on the measurements of diffuse soil gases and surface and shallow waters. The soil-gas measurements included the determination of CO2 flux and the application to natural trace gases (e.g. radon) that may help to detect any CO2 leakage. As far as the water chemistry is concerned, geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of the PDT of Hontomín were analyzed to determine the most suitable parameters to trace the migration of the injected CO2 into the near-surface environments. The accumulation chamber method was used to measure the diffuse emission of CO2 at the soil-atmosphere interface. Although this technique has widely been applied in different scientific areas, it was considered of the utmost importance to adapt the optimum methodology for measuring the CO2 soil flux and estimating the total CO2 output to the specific features of the site where CO2 is to be stored shortly. During the pre-injection phase CO2 fluxes are expected to be relatively low where in the intra- and post-injection phases, if leakages are to be occurring, small variation in CO2 flux might be detected when the CO2 “noise” is overcoming the biological activity of the soil (soil respiration). CO2 flux measurements by the accumulation chamber method could be performed without vegetation clearance or after vegetation clearance. However, the results obtained after clearance show less dispersion and this suggests that this procedure appears to be more suitable for monitoring CO2 Storage sites. The measurement protocol, applied for the determination of the CO2 flux baseline at Hontomín, has included the following steps: a) cleaning and removal of both the vegetal cover and top 2 cm of soil, b) waiting to reduce flux perturbation due to the soil removal and c) measuring the CO2 flux. Once completing the CO2 flux measurements and detected whether there were anomalies zones, the total CO2 output was estimated to quantify the amount of CO2 released to the atmosphere in each of the studied areas. There is a wide range of methodologies for the estimation of the CO2 output, which were applied to understand which one was the most representative. In this study six statistical methods are presented: arithmetic mean, minimum variances unbiased estimator, bootstrap resample, partitioning of data into different populations with a graphical and a maximum likelihood procedures, and sequential Gaussian simulation. Eight campaigns were carried out in the Hontomín CO2 Storage Technology Development Plant and in natural CO2 analogues. The results show that sequential Gaussian simulation is the most accurate method to estimate the total CO2 output and the confidential interval. Nevertheless, a variety of statistic methods were also used. As a consequence, an application procedure for selecting the most realistic method was developed. The first step to estimate the total emanation rate was the variogram analysis. If the relation among the data can be explained with the variogram, the best technique to calculate the total CO2 output and its confidence interval is the sequential Gaussian simulation method (sGs). If the data are independent, their distribution is to be analyzed. For normal and log-normal distribution the proper methods are the arithmetic mean and minimum variances unbiased estimator, respectively. If the data are not normal (log-normal) or are a mixture of different populations the best approach is the bootstrap resampling. According to these steps, the maximum confidence interval was about ±20/25%, with most of values between ±3.5% and ±8%. Partitioning of CO2 flux data into different populations may help to interpret the data as their distribution can be affected by different geochemical processes, e.g. geological or biological sources of CO2. Consequently, it may be an important tool in a monitoring CCS program, where the main goal is to demonstrate that there are not leakages from the reservoir to the atmosphere and, if occurring, to be able to detect and quantify it. Results show that the partitioning of populations is better performed by maximum likelihood criteria, since graphical procedures have a degree of subjectivity in the interpretation and results may not be reproducible. The relationship between CO2 flux and radon isotopes (222Rn and 220Rn) was studied in natural analogues. In all emissions zones, a positive relation between 222Rn and CO2 was observed. However, the relationship between activity of 220Rn and CO2 flux is not clear. In some cases the 220Rn activity indeed increased with the CO2 flux in other measurements a decrease was recognized. We can speculate that this effect was possibly related to the route (deep or shallow) of the radon source. These results may confirm the possible use of the radon isotopes as tracers for the gas origin and their application in the detection of leakages. With respect to the CO2 flux baseline at the TDP of Hontomín, soil flux measurements in the vicinity of oil boreholes, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber. Seven surveys were carried out from November 2009 to summer 2011. More than 4,000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were relatively low (from 5 to 13 g•m-2•day-1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g•m-2•d-1 for non-ploughed areas in autumn-winter seasons and 3.5 and 12 g•m-2•d-1 for in ploughed and non-ploughed areas, respectively, in spring-summer time, and UCL99 of 26 g•m-2•d-1 for autumn-winter in not-ploughed areas and 34 and 42 g•m-2•d-1 for spring-summer in ploughed and not-ploughed areas, respectively, were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project. The first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín–Huermeces (Burgos, Spain) are presented and discussed. The chemical and features of the spring waters suggest that they are related to a shallow hydrogeological system as the concentration of the Total Dissolved Solids approaches 800 mg/L with a Ca2+(Mg2+)-HCO3 − composition, similar to that of the surface waters. Some spring waters are characterized by relatively high concentrations of NO3 − (up to 123 mg/L), unequivocally suggesting an anthropogenic source. Anomalous concentrations of Cl−, SO4 2−, As, B and Ba were measured in two springs, discharging a few hundred meters from the oil wells, and in the Rio Ubierna. These contents are possibly indicative of mixing processes between deep and shallow aquifers. The chemistry of the dissolved gases also evidences the shallow circuits of the Hontomín– Huermeces, mainly characterized by an atmospheric source as highlighted by the contents of N2, O2, Ar and their relative ratios. Nevertheless, significant concentrations (up to 63% by vol.) of isotopically negative CO2 (<−17.7‰ V-PDB) were found in some water samples, likely related to a biogenic source. The geochemical and isotopic data of the surface and spring waters in the surroundings of Hontomín can be considered as background values when intra- and post-injection monitoring programs will be carried out. In this respect, main and minor solutes, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Development of a methodology for nuclear data uncertainty propagation on isotopic evolution calculations for advanced nuclear systems

Una apropiada evaluación de los márgenes de seguridad de una instalación nuclear, por ejemplo, una central nuclear, tiene en cuenta todas las incertidumbres que afectan a los cálculos de diseño, funcionanmiento y respuesta ante accidentes de dicha instalación. Una fuente de incertidumbre son los datos nucleares, que afectan a los cálculos neutrónicos, de quemado de combustible o activación de materiales. Estos cálculos permiten la evaluación de las funciones respuesta esenciales para el funcionamiento correcto durante operación, y también durante accidente. Ejemplos de esas respuestas son el factor de multiplicación neutrónica o el calor residual después del disparo del reactor. Por tanto, es necesario evaluar el impacto de dichas incertidumbres en estos cálculos. Para poder realizar los cálculos de propagación de incertidumbres, es necesario implementar metodologías que sean capaces de evaluar el impacto de las incertidumbres de estos datos nucleares. Pero también es necesario conocer los datos de incertidumbres disponibles para ser capaces de manejarlos. Actualmente, se están invirtiendo grandes esfuerzos en mejorar la capacidad de analizar, manejar y producir datos de incertidumbres, en especial para isótopos importantes en reactores avanzados. A su vez, nuevos programas/códigos están siendo desarrollados e implementados para poder usar dichos datos y analizar su impacto. Todos estos puntos son parte de los objetivos del proyecto europeo ANDES, el cual ha dado el marco de trabajo para el desarrollo de esta tesis doctoral. Por tanto, primero se ha llevado a cabo una revisión del estado del arte de los datos nucleares y sus incertidumbres, centrándose en los tres tipos de datos: de decaimiento, de rendimientos de fisión y de secciones eficaces. A su vez, se ha realizado una revisión del estado del arte de las metodologías para la propagación de incertidumbre de estos datos nucleares. Dentro del Departamento de Ingeniería Nuclear (DIN) se propuso una metodología para la propagación de incertidumbres en cálculos de evolución isotópica, el Método Híbrido. Esta metodología se ha tomado como punto de partida para esta tesis, implementando y desarrollando dicha metodología, así como extendiendo sus capacidades. Se han analizado sus ventajas, inconvenientes y limitaciones. El Método Híbrido se utiliza en conjunto con el código de evolución isotópica ACAB, y se basa en el muestreo por Monte Carlo de los datos nucleares con incertidumbre. En esta metodología, se presentan diferentes aproximaciones según la estructura de grupos de energía de las secciones eficaces: en un grupo, en un grupo con muestreo correlacionado y en multigrupos. Se han desarrollado diferentes secuencias para usar distintas librerías de datos nucleares almacenadas en diferentes formatos: ENDF-6 (para las librerías evaluadas), COVERX (para las librerías en multigrupos de SCALE) y EAF (para las librerías de activación). Gracias a la revisión del estado del arte de los datos nucleares de los rendimientos de fisión se ha identificado la falta de una información sobre sus incertidumbres, en concreto, de matrices de covarianza completas. Además, visto el renovado interés por parte de la comunidad internacional, a través del grupo de trabajo internacional de cooperación para evaluación de datos nucleares (WPEC) dedicado a la evaluación de las necesidades de mejora de datos nucleares mediante el subgrupo 37 (SG37), se ha llevado a cabo una revisión de las metodologías para generar datos de covarianza. Se ha seleccionando la actualización Bayesiana/GLS para su implementación, y de esta forma, dar una respuesta a dicha falta de matrices completas para rendimientos de fisión. Una vez que el Método Híbrido ha sido implementado, desarrollado y extendido, junto con la capacidad de generar matrices de covarianza completas para los rendimientos de fisión, se han estudiado diferentes aplicaciones nucleares. Primero, se estudia el calor residual tras un pulso de fisión, debido a su importancia para cualquier evento después de la parada/disparo del reactor. Además, se trata de un ejercicio claro para ver la importancia de las incertidumbres de datos de decaimiento y de rendimientos de fisión junto con las nuevas matrices completas de covarianza. Se han estudiado dos ciclos de combustible de reactores avanzados: el de la instalación europea para transmutación industrial (EFIT) y el del reactor rápido de sodio europeo (ESFR), en los cuales se han analizado el impacto de las incertidumbres de los datos nucleares en la composición isotópica, calor residual y radiotoxicidad. Se han utilizado diferentes librerías de datos nucleares en los estudios antreriores, comparando de esta forma el impacto de sus incertidumbres. A su vez, mediante dichos estudios, se han comparando las distintas aproximaciones del Método Híbrido y otras metodologías para la porpagación de incertidumbres de datos nucleares: Total Monte Carlo (TMC), desarrollada en NRG por A.J. Koning y D. Rochman, y NUDUNA, desarrollada en AREVA GmbH por O. Buss y A. Hoefer. Estas comparaciones demostrarán las ventajas del Método Híbrido, además de revelar sus limitaciones y su rango de aplicación. ABSTRACT For an adequate assessment of safety margins of nuclear facilities, e.g. nuclear power plants, it is necessary to consider all possible uncertainties that affect their design, performance and possible accidents. Nuclear data are a source of uncertainty that are involved in neutronics, fuel depletion and activation calculations. These calculations can predict critical response functions during operation and in the event of accident, such as decay heat and neutron multiplication factor. Thus, the impact of nuclear data uncertainties on these response functions needs to be addressed for a proper evaluation of the safety margins. Methodologies for performing uncertainty propagation calculations need to be implemented in order to analyse the impact of nuclear data uncertainties. Nevertheless, it is necessary to understand the current status of nuclear data and their uncertainties, in order to be able to handle this type of data. Great eórts are underway to enhance the European capability to analyse/process/produce covariance data, especially for isotopes which are of importance for advanced reactors. At the same time, new methodologies/codes are being developed and implemented for using and evaluating the impact of uncertainty data. These were the objectives of the European ANDES (Accurate Nuclear Data for nuclear Energy Sustainability) project, which provided a framework for the development of this PhD Thesis. Accordingly, first a review of the state-of-the-art of nuclear data and their uncertainties is conducted, focusing on the three kinds of data: decay, fission yields and cross sections. A review of the current methodologies for propagating nuclear data uncertainties is also performed. The Nuclear Engineering Department of UPM has proposed a methodology for propagating uncertainties in depletion calculations, the Hybrid Method, which has been taken as the starting point of this thesis. This methodology has been implemented, developed and extended, and its advantages, drawbacks and limitations have been analysed. It is used in conjunction with the ACAB depletion code, and is based on Monte Carlo sampling of variables with uncertainties. Different approaches are presented depending on cross section energy-structure: one-group, one-group with correlated sampling and multi-group. Differences and applicability criteria are presented. Sequences have been developed for using different nuclear data libraries in different storing-formats: ENDF-6 (for evaluated libraries) and COVERX (for multi-group libraries of SCALE), as well as EAF format (for activation libraries). A revision of the state-of-the-art of fission yield data shows inconsistencies in uncertainty data, specifically with regard to complete covariance matrices. Furthermore, the international community has expressed a renewed interest in the issue through the Working Party on International Nuclear Data Evaluation Co-operation (WPEC) with the Subgroup (SG37), which is dedicated to assessing the need to have complete nuclear data. This gives rise to this review of the state-of-the-art of methodologies for generating covariance data for fission yields. Bayesian/generalised least square (GLS) updating sequence has been selected and implemented to answer to this need. Once the Hybrid Method has been implemented, developed and extended, along with fission yield covariance generation capability, different applications are studied. The Fission Pulse Decay Heat problem is tackled first because of its importance during events after shutdown and because it is a clean exercise for showing the impact and importance of decay and fission yield data uncertainties in conjunction with the new covariance data. Two fuel cycles of advanced reactors are studied: the European Facility for Industrial Transmutation (EFIT) and the European Sodium Fast Reactor (ESFR), and response function uncertainties such as isotopic composition, decay heat and radiotoxicity are addressed. Different nuclear data libraries are used and compared. These applications serve as frameworks for comparing the different approaches of the Hybrid Method, and also for comparing with other methodologies: Total Monte Carlo (TMC), developed at NRG by A.J. Koning and D. Rochman, and NUDUNA, developed at AREVA GmbH by O. Buss and A. Hoefer. These comparisons reveal the advantages, limitations and the range of application of the Hybrid Method.

Modelización y validación experimental del transporte de tritio en metal líquido y otros materiales para reactores de fusión nuclear magnética e inercial

La presente tesis se centra en el estudio de los fenómenos de transporte de los isótopos de hidrógeno, y más concretamente del tritio, en materiales de interés para los reactores de fusión nuclear. Los futuros reactores de fusión nuclear necesitarán una Planta de Tritio, con una envoltura regeneradora (breeding blanket) y unos sistemas auxiliares claves para su diseño. Por lo tanto su desarrollo y cualificación son cruciales para demostrar que los reactores de fusión son una opción viable como futura fuente de energía. Se han resaltado los diferentes retos de la difusión y retención de estas especies ligeras para cada sistema de la Planta de Tritio, y se han identificado las necesidades experimentales y paramétricas para abordar las simulaciones de difusión, como factores de transporte como la difusividad, absorción/desorción, solubilidad y atrapamiento. Se han estudiado los fenómenos de transporte y parámetros del T en el metal líquido LiPb, componente del breeding blanket tanto para una planta de fusión magnética como inercial. Para ello se han utilizado dos experimentos con características diversas, uno de ellos se ha llevado a cabo en un reactor de alto flujo, y por lo tanto, en condiciones de irradiación, y el otro sin irradiación. Los métodos de simulación numérica aplicados se han adaptado a los experimentos para las mediciones y para estudiar el régimen de transporte. En el estudio de estos experimentos se ha obtenido un valor para algunos de los parámetros claves en el transporte y gestión del tritio en el reactor. Finalmente se realiza un cálculo de la acumulación y difusión de tritio en una primera pared de tungsteno para un reactor de fusión inercial. En concreto para el proyecto de fusión por láser europeo, HiPER (para sus fases 4a y 4b). Se ha estudiado: la implantación de los isótopos de H y He en la pared de W tras una reacción de fusión por iluminación directa con un láser de 48MJ; el efecto en el transporte de T de los picos de temperatura en el W debido a la frecuencia de los eventos de fusión; el régimen de transporte en la primera pared. Se han identificado la naturaleza de las trampas más importantes para el T y se ha propuesto un modelo avanzado para la difusión con trampas. ABSTRACT The present thesis focuses into study the transport phenomenons of hydrogen isotopes, more specifically tritium, in materials of interest for nuclear fusion reactors. The future nuclear reactors will be provided of a Tritium Plant, with its breeding blanket and its auxiliary systems, all of them essential components for the plant. Therefore a reliable development and coalification are key issues to prove the viability of the nuclear fusion reactors as an energy source. The currently challenges for the diffusion and accumulation of these light species for each system of the TP has been studied. Experimental and theoretical needs have been identified and analyzed, specially from the viewpoint of the parameters. To achieve reliable simulations of tritium transport, parameters as diffusivity, absorption/desorption, solubility and trapping must be reliables. Transport phenomenon and parameters of T in liquid metal have been studied. Lead lithium is a key component of the breeding blanket, either in magnetic or inertial fusion confinement. Having this aim in mind, two experiments with different characteristics have been used; one of them has been realized in a high flux reactor, and hence, under irradiation conditions. The other one has been realized without radiation. The mathematical methods for the simulation have been adapted to the experiments, for the measures and also to study the transport behavior. A value for some key parameters for tritium management has been obtained in these studies. Finally, tritium accumulation and diffusion in a W first wall of an inertial nuclear fusion reactor has been assessed. A diffusion model of the implanted H, D, T and He species for the two initial phases of the proposed European laser fusion Project HiPER (namely, phase 4a and phase 4b) has been implemented using Tritium Migration Analysis Program, TMAP7. The effect of the prompt and working temperatures and the operational pulsing modes on the diffusion are studied. The nature of tritium traps in W and their performance has been analyzed and discussed.

Las mineralizaciones auríferas de Galicia Occidental en la Zona de Cizalla Busto-Limideiro (A Coruña)

Las mineralizaciones de oro de Galicia Occidental ligadas a grandes estructuras comienzan a ser interpretadas de manera más rigurosa, al mismo tiempo que la continuada elevación de las cotizaciones del oro han hecho más atractiva la reactivación de áreas hasta ahora casi olvidadas. El propósito de esta investigación participa de las razones anteriores. La zona de cizalla Busto Limideiro se inscribe en una de estas estructuras, la cizalla Santa Comba-Punta Langosteira, con la que se relacionan un importante número de yacimientos auríferos entre los que destacan Corcoesto y Monte Piñor. Por su situación geológica, estos indicios se enmarcan en el Macizo Ibérico, en la Zona de Galicia Tras-os-Montes. Este trabajo aporta más conocimiento sobre un modelo general del oro de los yacimientos del noroeste español, integrándolo en los yacimientos de oro orogénico del Varisco europeo. Este estudio ha supuesto un importante esfuerzo investigador centrado sobre de los procesos genéticos que han dado lugar a las mineralizaciones auríferas. Los estudios petrográficos, metalogénicos, geoquímicos, microtemométricos y el análisis de isótopos han permitido establecer un modelo general y otros más específicos sobre la existencia y condiciones de formación en la concentración del oro mineralógico. De esta forma se ha caracterizado su mineralogía, quimismo, edad, y naturaleza de los paleofluidos que han intervinieron en su formación. De los resultados obtenidos se destaca que estas mineralizaciones se encuentran fuertemente ligadas a las estructuras y transformaciones físicas y químicas de las rocas afectadas por la banda de cizalla. Derivada de esta situación, se comprende la disposición espacial, a diversas escalas, con estructuras claramente en la transición del dominio dúctil a frágil. Además, resulta evidente, a escala macro y microscópica, el emplazamiento preferente del oro en estructuras de deformación miloníticas. Los estudios anteriores posibilitan formalizar un modelo de génesis para el yacimiento de Monte Piñor que resulta compatible con el establecido por Groves (1998), e incluirlas en un mismo contexto, como yacimientos de oro orogénico mesotermales, válido para el conjunto de todos los de la zona de cizalla Busto Limideiro. El estudio no omite la consideración geoestadística de este yacimiento, justificando la elección de direcciones de anisotropía explicadas por consideraciones de carácter estructural. Con ello, buscando el lado práctico desde el punto de vista minero, esta investigación propone líneas generales de aplicación de actuación con sentido de exploración, contribuyendo a la puesta en valor del conocimiento de las mineralizaciones de oro del noroeste español. ABSTRACT Old and nearly forgotten gold mines from Western Galicia are currently attractive to be recovered due to both, the new strict interpretation of gold mineralization linked to Variscan large structures and the increasing gold prices. The shear zone Busto Limiderio is relevant as a part of a Variscan large structure, Santa Comba-Punta Langosteira, which is related with a high number of gold deposits as Corcoesto and Monte Piñor. For its geological situation, these signs are part of the Iberian Massif in the area of Galicia Tras-os-Montes. Hence, the main aim of this study is to provide additional and important knowledge of Busto Limiderio about the general model of gold deposits in the Spanish northwest, integrating them in the orogenic gold deposits of European Variscan, enabling to analyze the economic recovery possibility of this area. This thesis has been focused on the genetic processes that have led to the gold mineralization. Petrographic studies, metalogenic, geochemical, microthermometric and isotope analysis have allowed to establish a general and specific models about the existence and training conditions in the concentration of mineralogical gold. Consequently, the mineralogy, chemical characteristics, age, and nature of paleofluids that have participated in its formation have been characterized. The results showed that Busto Limiderio gold mineralization are strongly linked to the structures and physical and chemical transformations of affected rocks by the shear band. The mineralization spatial distribution at various scales is due to the link with the large structures in the ductile-to-brittle transition domain. Furthermore, this factor determine at macro and microscopic scale, the preferred location of the gold structures mylonitic deformation. Previous studies have permitted to establish a model of genesis for the site of Monte Piñor which is compatible with the established by Groves (1998), and also to include them in the same context as mesothermal orogenic gold deposits, valid throughout all of Busto Limideiro shear zone. This study does not omit consideration of the geostatistical deposit, justifying the choice of directions of anisotropy explained by structural considerations. This thesis concludes with general action guidelines for mining exploration, contributing to the enhancement of knowledge of Spanish gold mineralization northwest.

Sais de diazônio: síntese e eletrorredução

Sais de diazônio são um classe de compostos amplamente usados em química orgânica. Sua aplicação abrange uma gama de sínteses desde corantes até reações de hetero-acoplamento para produção de fármacos, mas pouco é conhecido de sua redução eletroquímica para fins sintéticos. As metodologias empregadas na redução de sais de diazônio geralmente envolvem o uso de metais ou compostos capazes de transferir elétrons como Pd, Cu ou tetratiafulvaleno. Neste trabalho é descrita a redução eletroquímica de dois sais de diazônio: tetrafluoroborato de 2-(2-propen-1-ilóxi)benzenodiazônio (1) e tetrafluoroborato de 2-(2-propen-1-iltio)benzenodiazônio (2) usando três eletrodos: Pt, Hg e pó de grafite. Quando foi feita a eletrólise de (1) utilizando cátodo de Hg vários produtos foram formados envolvendo uma reação de ciclização intramolecular, porém não conseguimos separá-los pelos métodos cromatográficos. A eletrólise de (2) em condições experimentais similares conduziu a uma mistura complexa de produtos provavelmente devido a uma interação do Hg com o átomo de enxofre do substrato e seus produtos de redução. Usando o cátodo de Pt e sal (1) a reação não foi eficiente pois ocorria uma queda brusca da corrente, provavelmente devido ao bloqueio da superfície do eletrodo. Concernente ao eletrodo de pó de grafite, apenas alguns experimentos preliminares foram feitos, portanto uma análise de seu desemprenho é prematura.

Caracterização microbiológica, química e presença de poluentes orgânicos em amostras de lodo de esgoto de São Paulo

Objetivou-se com este trabalho avaliar o potencial agrícola do lodo de esgoto produzido no estado de São Paulo, bem como, verificar a possibilidade de interação entre a composição química e a abundância relativa de bactérias no lodo. Foram realizadas coletas de amostra de lodo de esgoto em 19 estações de tratamento de esgoto, em três épocas distintas. Nas amostras provenientes das três épocas foram determinados as concentrações dos 16 hidrocarbonetos policíclicos aromáticos (HPAs) listados como prioritários no monitoramento ambiental pela USEPA (acenafteno, acenaftileno, antraceno, benzo(a)antraceno, benzo(a)pireno, benzo(b)fluoranteno, benzo(ghi)perileno, benzo(k)fluoranteno, criseno, dibenzo(a,h)antraceno, fenantreno, fluoranteno, fluoreno, indeno(1,2,3-cd)pireno, naftaleno e pireno). Nas amostras da segunda época de coleta, além da presença de HPAs, determinou-se as concentrações de poluentes orgânicos emergentes (hormônios, produtos farmacêuticos e produtos de uso industrial), realizou-se a caracterização completa segundo a Resolução CONAMA 375/2006 (umidade, pH, N-Kjeldahl e inorgânico, carbono orgânico, cálcio, potássio, fósforo, magnésio, enxofre, boro, cobre, ferro, níquel, manganês, molibdênio, selênio, zinco, alumínio, arsênio, bário, cádmio, cromo, chumbo, mercúrio e sódio) e a caracterização da comunidade bacteriana através de metodologia independente de cultivo (sequenciamento illumina). Os macronutrientes em maiores concentrações no lodo de esgoto são: N > Ca > S > P > Mg > K. Os elementos inorgânicos Ni e Zn apresentaram concentração superior à máxima permitida para utilização agrícola pela resolução Conama 375/2006 em 1 e 3 amostras, respectivamente. A substância inorgânica que mais limita o enquadramento do lodo de esgoto como adubo orgânico (Instrução Normativa 27/2006) é o Hg. Os compostos benzilparabeno, bisfenol AF (BPAF), ácido perfluorooctanoico (PFOA) e tetrabromobisfenol A (TBBPA) não foram detectados. Por outro lado, cimetidina, metilparabeno, bisfenol A (BPA) e triclocarban foram detectados nas 19 amostras avaliadas. O composto presente em maior concentração é o triclocarban. As concentrações de hidrocarbonetos policíclicos aromáticos são baixas, de acordo com a norma Europeia. Os filos Proteobacteria e Bacteroidetes estão presentes em maior abundância relativa. Existe uma comunidade bacteriana núcleo nas estações de tratamento de esgoto do estado de São Paulo, composta por 81 gêneros, presentes nas 19 ETEs avaliadas, dos quais, os que estão em maior abundância relativa são Treponema, Clostridium, Propionibacterium, Syntrophus e Desulfobulbus. A elevação do pH a valores próximos de 12 reduz a diversidade microbiana. Considerando a abundância relativa e a composição química do lodo de esgoto, as estações podem ser agrupadas em três grupos distintos, sendo que um deles é influenciado principalmente pelos teores de Ca, Zn e Cu, o outro pelos teores de Fe e S e o terceiro grupo que foi influenciado pelos demais fatores avaliados.

Climatic evolution of eastern South America during the last deglaciation: a paleoceanographic approach

Para dar suporte ao atual debate sobre as consequências climáticas da liberação antropogênica de CO2 na atmosfera, o refinamento do conhecimento sobre mudanças climáticas e oceanográficas no passado é necessário. A Circulação de Revolvimento Meridional do Atlântico (CRMA) tem papel fundamental na oceanografia e clima das áreas sob influência do Oceano Atlântico, controlando diretamente a estratificação e distribuição de massas d\'água, a quantidade de calor transportada pelo oceano e os ciclo e armazenamento de compostos químicos, como o CO2 em mar profundo. A formação e circulação da Água Intermediária Antártica (AIA), envolvida no transporte de calor e sal para o giro subtropical do Hemisfério Sul e nas teleconexões climáticas entre altas e baixas latitudes, é componente importante do ramo superior da CRMA. A reconstrução de propriedades de massas de água intermediárias é, portanto, importante para a compreensão dos sistemas de retroalimentação entre oceano-clima. No entanto, informações quanto a evolução da AIA continuam limitadas. Oscilações da CRMA e consequentes mudanças na distribuição de calor tem implicações importantes para o clima Sul Americano, influenciando a disponibilidade de umidade para o Sistema de Monções Sul Americano (SMSA), via temperatura da superfície marinha e posicionamento da Zona de Convergência Intertropical. Neste trabalho nós reconstruímos a assinatura isotópica da AIA durante os estágios isotópicos marinhos 2 e 3 (41-12 cal ka AP) usando isótopos de carbono e oxigênio de foraminíferos bentônicos (gêneros Cibicidoides e Uvigerina) de um testemunho de sedimentos marinhos datados por radiocarbono (1100 m de profundidade e a 20°S na costa do Brasil). Concluímos que propriedades físicas e químicas da AIA mudaram durante os estadiais Heinrich 3 e 4, provavelmente como consequência de enfraquecimento da CRMA durante estes períodos. Também reconstruímos as condições continentais do leste brasileiro entre o último máximo glacial e a deglaciação (23-12 cal ka AP) baseadas em razões Ti/Ca de nosso testemunho de sedimentos marinhos como indicadoras de aporte terrígeno do Rio Doce. A maior parte da chuva que cai na Bacia do Rio Doce está relacionada a atividade do SMAS. Nosso registro de Ti/Ca em conjunto com \'\'delta\' POT.18\'O de espeleotemas da Caverna Lapa Sem Fim, também no leste do Brasil, sugere diminuição marcante da chuva durante o interestadial Bølling-Allerød, provavelmente relacionada a enfraquecimento do SMAS. Ademais comparamos as razões de Ti/Ca com dados de saída da rodada SYNTRACE do modelo climático CCSM3 com forçantes transientes para a última deglaciação. Os registros geoquímicos e a saída do modelo mostram resultados consistentes entre si e sugerem que o leste da América do Sul passou pelo seu período mais seco de toda a última deglaciação durante o interestadial Bølling-Allerød, provavelmente relacionado ao enfraquecimento do SMAS.

Metodo Rb-Sr em rochas sedimentares: aplicação para as Bacias do Paraná e do Amazonas

Estudo radiométrico em sedimentos com posição estratigráfica conhecida, das Bacias do Paraná e do Amazonas, permite estabelecer critérios de seleção de amostras e procedimentos experimentais adequados para a obtenção de idades significativas. Foram efetuadas 120 determinações Rb-Sr e 44 K-Ar. As idades K-Ar foram empregadas essencialmente para auxiliar a interpretação dos dados Rb-Sr. As interpretações Rb-Sr foram efetuadas mediante gráficos com linhas isócronas. Na avaliação estatística dos dados, o método usual dos mínimos quadrados, revelou-se deficiente na estimação dos parâmetros, em alguns exemplos testados. Na presente investigação, a melhor isócrona em cada caso foi definida levando-se em conta uma ponderação adequada para os pontos e uma correlação entre os erros específica e conveniente. A formação Ponta Grossa foi estudada em 8 amostras provenientes de testemunhos de 5 sondagens da Petrobrás S. A. na Bacia do Paraná. Os resultados obtidos em pelo menos 2 deles são concordantes com a idade estratigráfica. Nos outros 2 poços, apesar dos poucos dados, as isócronas mostram uma possível concordância, indicando que a diagênese teria ocorrido logo após a sedimentação marinha. A formação Rio Bonito, estudada em 8 amostras do poço TV-4-SC, apresentou 3 isócronas aproximadamente paralelas, com idades mais ou menos semelhantes e concordantes com a idade estratigráfica. Tal comportamento indicaria uma homogeneização isotópica mesmo entre as frações grosseiras. Dentre as amostras estudadas na Bacia do Paraná, os sedimentos da Formação Itararé foram os únicos que não puderam ser interpretados adequadamente. Isto evidenciou que as isócronas Rb-Sr devem ser restritas a rochas de apenas um determinado tipo litológico, quando examinamos Formações como a Itataré, de ambientes variados (fluvial, lacustre, glacial, marinho). A formação Trombetas foi estudada em 7 amostras de 2 poços localizados no Médio Amazonas. Ambas as isócronas obtidas, indicando idade ordoviciana-siluriana permitem supor que houve apenas uma homogeneização isotópica parcial após a deposição. Os dados podem ser considerados concordantes se forem levados em conta os erros experimentais das isócronas. As Formações Maecuru e Ereré foram estudadas em 6 amostras do poço MS-4-AM. A litologia desfavorável das 4 amostras da Formação Ereré (Membro Ariramba), levaram o autor a definir uma "isócrona mínima", cuja idade revelou-se próxima da admitida estratigraficamente. A Formação Maecuru (Membro Jatapu), estudada em 2 arenitos arcozianos, apesar do material não ser considerado satisfatório para datações, evidenciaram uma isócrona de referência cuja idade é compatível com a situação estratigráfica. A Formação Curuá foi analisada em 3 amostras do poço NA-1-PA. Tanto as rochas totais como as frações situaram-se sobre uma isócrona de referência cuja idade calculada apresentou concordância, com a idade estratigráfica, dentro do erro experimental. A boa correlação linear verificada leva a admitir uma diagênese precoce, acompanhada de equilíbrio entre os isótopos de Sr. A Formação Itaituba foi investigada em 7 amostras de 2 poços, situados um de cada lado do Alto de Purus. Os folhetos evidenciaram grande dispersão dos pontos sobre o diagrama Rb87/Sr86 x Sr87/Sr86, devida a teores variáveis de minerais detríticos difíceis de serem identificados petrograficamente. Novamente foi traçada uma "isócrona mínima" da qual participaram materiais calcíferos. As idades identicas, bem como a concordância com a idade estratigráfica da formação, parecem demonstrar a validade da técnica empregada. Os dados do presente trabalho indicam que rochas sedimentares podem ser datadas pelo método Rb-Sr, desde que sejam obedecidos alguns critérios importantes de seleção do material. Além disso devem ser empregadas técnicas apropriadas, tais como separação ) granulométrica de frações menores que 2 ou 4 μ, ou lixiaviação com HCl. As isócronas a serem traçadas, as quais indicariam a época da diagênese, devem incluir material de litologia semelhante, de um só ambiente de formação. As análises K-Ar podem servir como dados auxiliares, principalmente para avaliar a quantidade de material detrítico existente no sistema.

Estudo da substituição da fluorita por alumina ou sodalita e de cal por resíduo de mármore em escórias sintéticas dessulfurantes.

A siderurgia vem sofrendo transformações que buscam inovação e matérias-primas alternativas. Dentro deste contexto, o uso de resíduos industriais para a formação de escórias sintéticas é tido como alternativa na busca de novos materiais e rotas de reaproveitamento de resíduos. Portanto, este trabalho teve como objetivo estudar o uso de escórias sintéticas na etapa de dessulfuração do ferro-gusa, aço e ferro fundido. Assim como, propor a utilização da sodalita e da alumina em substituição à fluorita e o resíduo de mármore em substituição à cal convencional. Inicialmente, o resíduo foi caracterizado utilizando as seguintes técnicas: análise química, análise granulométrica, área de superfície específica, difração de raios-X, microscopia eletrônica de varredura (MEV) e análise de espectroscopia por energia dispersiva (EDS). Os resultados da caracterização mostraram que aproximadamente 90% das partículas do resíduo de mármore estão abaixo de 100m e sua área superficial foi de 0,24m²/g. Através da difração de raios-X foi observado que o resíduo é composto por CaCO3, MgCO3 e SiO2. Na sequência, foram feitas simulações com o software Thermo-Calc para obter dados termodinâmicos das fases presentes nas misturas e compará-los com os resultados experimentais. Além disso, também foram calculados dados de capacidade de sulfeto (Cs), partição de enxofre (Ls) e basicidade ótica () das misturas iniciais. Posteriormente, foram realizados os ensaios experimentais em escala laboratorial para ferro-gusa, ferro fundido e aço, respectivamente nas temperaturas de 1400°C, 1550°C e 1600°C. Nos ensaios de dessulfuração do aço e do ferro-gusa, utilizou-se um rotor de alumina com o objetivo de favorecer a agitação no metal e aumentar a remoção de enxofre. Na etapa de dessulfuração do ferro-gusa, constatou-se que a fase sólida de CaO é a responsável pela remoção de enxofre e que a presença das fases silicato tricálcio e aluminato tricálcio (3CaO.SiO2 e 3CaO.Al2O3) limitam a reação, sendo maiores suas concentrações nas escórias que utilizaram o resíduo de mármore e sodalita, devido a presença de SiO2 e Al2O3 nestas matérias-primas. Já para o aço e o ferro fundido, que foram estudados com escórias à base de CaO e Al2O3, observou-se que o aumento da fase líquida favoreceu a dessulfuração. Verificou-se que a dessulfuração no ferro fundido foi por escória de topo e no aço por um processo misto, onde a fase líquida e fase sólida participaram da dessulfuração.

Utilização de resíduo de mármore e de óxido de boro em escórias sintéticas dessulfurantes para aços.

O mercado atual exige das indústrias siderúrgicas aços de melhor qualidade produzidos por meio de processos que causem menor impacto ao meio ambiente. Dessa forma, este trabalho teve como objetivo reciclar o resíduo de mármore gerado na indústria de rochas ornamentais, que possui em sua composição óxido de cálcio (CaO) e óxido de magnésio (MgO). O CaO é suficiente para substituir a cal nas escórias e o MgO contribui para a diminuição do desgaste dos refratários, através do emprego do resíduo no processo produtivo do aço. Além disso, foi realizada a substituição da fluorita por óxido de boro como fluxante na composição de misturas dessulfurantes. O resíduo de mármore foi caracterizado utilizando as seguintes técnicas: análise química via EDXFR, análise granulométrica via espalhamento de luz, área de superfície específica pelo método BET, difração de raios-X, microscopia eletrônica de varredura (MEV) e análise de micro-regiões por EDS. Visando verificar a eficiência na dessulfuração, foram formuladas misturas sintéticas utilizando a cal convencional ou resíduo de mármore, e a fluorita ou o óxido de boro. As misturas foram formuladas com o auxílio dos programas de termodinâmica computacional, Thermo-Calc e FactSage. Estas misturas foram adicionadas no aço fundido a temperatura de 1600°C sob atmosfera de argônio e agitadas por meio de um rotor de alumina. Amostras de metal foram retiradas para verificar a variação do teor de enxofre durante o experimento. O resíduo de mármore caracterizado, apresentou em sua composição 40% de CaO e 14% de MgO, na forma dos carbonatos CaCO3 e MgCO3. Obteve uma perda de massa de 42,1%, na forma de CO2 a temperatura de 780°C. Os experimentos mostraram que, as misturas testadas apresentaram, na maioria dos casos, eficiência de dessulfuração acima de 60%.

Classificação das amostras do ensaio de Baumann através do processamento digital de imagens.

O presente trabalho apresenta uma alternativa ao processo de classificação do defeito da segregação central em amostras de aço, utilizando as imagens digitais que são geradas durante o ensaio de Baumann. O algoritmo proposto tem como objetivo agregar as técnicas de processamento digital de imagens e o conhecimento dos especialistas sobre o defeito da segregação central, visando a classificação do defeito de referência. O algoritmo implementado inclui a identificação e a segmentação da linha segregada por meio da aplicação da transformada de Hough e limiar adaptativo. Adicionalmente, o algoritmo apresenta uma proposta para o mapeamento dos atributos da segregação central nos diferentes graus de severidade do defeito, em função dos critérios de continuidade e intensidade. O mapeamento foi realizado por meio da análise das características individuais, como comprimento, largura e área, dos elementos segmentados que compõem a linha segregada. A avaliação do desempenho do algoritmo foi realizada em dois momentos específicos, de acordo com sua fase de implementação. Para a realização da avaliação, foram analisadas 255 imagens de amostras reais, oriundas de duas usinas siderúrgicas, distribuídas nos diferentes graus de severidade. Os resultados da primeira fase de implementação mostram que a identificação da linha segregada apresenta acurácia de 93%. As classificações oriundas do mapeamento realizado para as classes de criticidade do defeito, na segunda fase de implementação, apresentam acurácia de 92% para o critério de continuidade e 68% para o critério de intensidade.

ESTUDIO DE LOS MODOS DE VIDA E HISTORIA BIOLÓGICA DE LAS POBLACIONES HUMANAS DE LA REGIÓN CENTRAL DE ARGENTINA DESDE EL REGISTRO BIOARQUEOLÓGICO

Se pretende aportar al estudio de la estructura, historia biológica y estilos de vida de las poblaciones que habitaron la región central de Argentina durante el Holoceno, desde una perspectiva que combina los aportes teóricos y metodológicos de la Genética del paisaje y la Bioarqueología. Interesa a) identificar barreras de diferenciación morfológica entre poblaciones, b) poner a prueba modelos poblacionales para explicar la variación observada e identificar las variables que contribuyan a dicha diferenciación, c) evaluar la congruencia de los resultados obtenidos, d) reconstruir los patrones de movilidad residencial de las poblaciones, e) estudiar sus patrones dietarios considerando diferencias temporales y espaciales, f) identificar indicadores de diversos tipos de estrés (nutricional, funcional), así como traumas, g) estudiar las historias tafonómicas del registro bioarqueológico regional, y h) proponer un modelo para explicar el poblamiento y la evolución local de las poblaciones que habitaron esta región, a partir de la información arqueológica y bioantropológica. Para el análisis de los patrones espaciales de variación biológica se trabajará a partir del registro de rasgos epigenéticos craneales, medidas lineales y datos obtenidos a partir de morfometría geométrica sobre fotografías en 2D sobre muestras arqueológicas procedentes de esta región y de otras regiones geográficas de la Argentina. Para el análisis de la estructura de la población se trabajará a partir del cálculo de la matriz R para datos morfológicos y sus estimaciones derivadas (distancia D², Fst, coordenadas principales) y la aplicación del modelo de Harpending y Ward. Desde la genética del paisaje, se realizarán análisis de autocorrelación espacial, barreras genéticas y análisis geoestadísticos (kriging). Para el estudio de los modos de vida a partir del registro bioarqueológico se relevarán patologías dento-alveolares y alteraciones vinculadas con la salud bucal tales como desgaste dental –a nivel micro y macroscópico- caries, abscesos, pérdidas dentales antemortem, cálculos, hipoplasias, marcadores esqueletales de salud y lesiones traumáticas. Se analizarán isótopos estables (δ13C, δ15N, 86Sr y 87Sr) en restos óseos humanos de diversos sitios arqueológicos con el objetivo de reconstruir patrones dietarios y analizar la movilidad residencial y migración de las poblaciones. Paralelamente, se establecerán procedimientos de control tafonómico de los restos óseos, y se harán análisis específicos para estudiar las historias tafonómicas y evaluar el grado de integridad de los contextos de depositación y de las colecciones en general. Estimamos que el análisis de los patrones espaciales y temporales de variabilidad morfológica craneofacial, así como el estudio de las dietas a partir de información isotópica y bioarqueológica, de las migraciones y la movilidad residencial de las poblaciones a partir de isótopos de estroncio, la reconstrucción de comportamientos y actividades cotidianas a partir de marcadores de estrés músculo-esqueletal, en un marco cronológico y espacial constituye un aporte novedoso y eficaz que permitirá incrementar de manera substancial la información sobre la evolución de las poblaciones originarias del centro del territorio argentino

Recursos hídricos superficiais vs subterrâneos na Várzea da Nazaré : implicações na qualidade da água

Tese de doutoramento, Geologia (Hidrogeologia), Universidade de Lisboa, Faculdade de Ciências, 2016

Obtenção de Mo2C com adição de Co por reação gás-sólido em reator de leito fixo

Sulfur compounds emissions have been, on the late years, subject to more severe environmental laws due to its impact on the environment (causing the acid rain phenomena) and on human health. It has also been object of much attention from the refiners worldwide due to its relationship with equipment’s life, which is decreased by corrosion, and also with products’ quality, as the later may have its color, smell and stability altered by the presence of such compounds. Sulfur removal can be carried out by hydrotreating (HDT) which is a catalytic process. Catalysts for HDS are traditionally based on Co(Ni)-Mo(W)/Al2O3. However, in face of the increased contaminants’ content on crude oil, and stricter legislation on emissions, the development of new, more active and efficient catalysts is pressing. Carbides of refractory material have been identified as potential materials for this use. The addition of a second metal to carbides may enhance catalytic activities by increasing the density of active sites. In the present thesis Mo2C with Co addition was produced in a fixed bed reactor via gas-solid reaction of CH4 (5%) and H2(95%) with a precursor made of a mix of ammonium heptamolybdate [(NH4)6[Mo7O24].4H2O] and cobalt nitrate[Co(NO3)2.6H2O] at stoichiometric amounts. Precursors’ where analyzed by XRF, XRD, SEM and TG/DTA. Carboreduction reactions were carried out at 700 and 750°C with two cobalt compositions (2,5 and 5%). Reaction’s products were characterized by XRF, XRD, SEM, TOC, BET and laser granulometry. It was possible to obtain Mo2C with 2,5 and 5% cobalt addition as a single phase at 750°C with nanoscale crystallite sizes. At 700°C, however, both MoO2 and Mo2C phases were found by XRD. No Co containing phases were found by XRD. XRF, however, confirmed the intended Co content added. SEM images confirmed XRD data. The increase on Co content promoted a more severe agglomeration of the produced powder. The same effect was noted when the reaction temperature was increased. The powder synthesized at 750°C with 2,5% Co addition TOC analysis indicated the complete conversion from oxide material to carbide, with a 8,9% free carbon production. The powder produced at this temperature with 5% Co addition was only partially converted (86%)