Terrestrial exposure of oilfield flowline additives diminish soil structural stability and remediative microbial function

Onshore oil production pipelines are major installations in the petroleum industry, stretching many thousands of kilometres worldwide which also contain flowline additives. The current study focuses on the effect of the flowline additives on soil physico-chemical and biological properties and quantified the impact using resilience and resistance indices. Our findings are the first to highlight deleterious effect of flowline additives by altering some fundamental soil properties, including a complete loss of structural integrity of the impacted soil and a reduced capacity to degrade hydrocarbons mainly due to: (i) phosphonate salts (in scale inhibitor) prevented accumulation of scale in pipelines but also disrupted soil physical structure; (ii) glutaraldehyde (in biocides) which repressed microbial activity in the pipeline and reduced hydrocarbon degradation in soil upon environmental exposure; (iii) the combinatory effects of these two chemicals synergistically caused severe soil structural collapse and disruption of microbial degradation of petroleum hydrocarbons.

Impact of the Arctic Ocean Atlantic water layer on Siberian shelf hydrography

This paper examines the role of the Arctic Ocean Atlantic water (AW) in modifying the Laptev Sea shelf bottom hydrography on the basis of historical records from 1932 to 2008, field observations carried out in April–May 2008, and 2002–2009 cross‐slope measurements. A climatology of bottom hydrography demonstrates warming that extends offshore from the 30–50 m depth contour. Bottom layer temperature‐time series constructed from historical records links the Laptev Sea outer shelf to the AW boundary current transporting warm and saline water from the North Atlantic. The AW warming of the mid‐1990s and the mid‐2000s is consistent with outer shelf bottom temperature variability. For April–May 2008 we observed on‐shelf near‐bottom warm and saline water intrusions up to the 20 m isobath. These intrusions are typically about 0.2°C warmer and 1–1.5 practical salinity units saltier than ambient water. The 2002–2009 cross‐slope observations are suggestive for the continental slope upward heat flux from the AW to the overlying low‐halocline water (LHW). The lateral on‐shelf wind‐driven transport of the LHW then results in the bottom layer thermohaline anomalies recorded over the Laptev Sea shelf. We also found that polynya‐induced vertical mixing may act as a drainage of the bottom layer, permitting a relatively small portion of the AW heat to be directly released to the atmosphere. Finally, we see no significant warming (up until now) over the Laptev Sea shelf deeper than 10–15 m in the historical record. Future climate change, however, may bring more intrusions of Atlantic‐modified waters with potentially warmer temperature onto the shelf, which could have a critical impact on the stability of offshore submarine permafrost.

Absolute Thellier paleointensities from Ponta Grossa dikes (southern Brazil) and the early Cretaceous geomagnetic field strength

We report a detailed rock magnetic and Thellier paleointensity study from similar to 130.5 Ma Ponta Grossa Dike Swarms in Southern Brazil. Twenty-nine samples from seven cooling units were pre-selected for paleointensity experiments based on their low viscosity index, stable remanent magnetization and close to reversible continuous thermomagnetic curves. 19 samples characterized by negative pTRM tests, Arai concave- up curves or positive pTRM tests with NRM loss uncorrelated with TRM acquisition were rejected. High quality reliable paleointensity determinations are determined from detailed evaluation criteria, with 10 samples belonging to three dikes passing the tests. The site-mean paleointensity values obtained in this study range from 25.6 +/- 4.3 to 11.3 +/- 2.1 mu T and the corresponding VDM`s range from 5.7 +/- 0.9 to 2.5 +/- 0.5 (10(22) Am(2)). These data yield a VDM mean value of 4.1 +/- 1.6 x 10(22) Am(2). Significant variability of Earth`s magnetic field strength is observed for Ponta Grossa Dikes with the mean value being significantly lower as compared to the mean VDM obtained from the nearby Parana Magmatic Province. The paleointensities for the Ponta Grossa Dikes are in agreement with absolute paleointensities retrieved from the submarine basaltic glasses from 130 to 120 Ma. It seems that a relatively low field prevailed just before the Cretaceous Normal Superchron.

Boron isotope composition of tourmalinite and vein tourmalines associated with gold mineralization, Serra do Itaberaba Group, central Ribeira Belt, SE Brazil

Important concentrations of tourmaline occur as gold-bearing stratiform tourmalinites and in mineralized quartz-tourmaline veins at the Tapera Grande and Quartzito gold prospects in the Mesoproterozoic Serra do Itaberaba Group, central Ribeira Belt (Sao Paulo State, SE Brazil). The main rock types in both prospects constitute the volcanic-sedimentary Morro da Pedra Preta Formation, which formed in a submarine back-arc setting. At Tapera Grande, the volcanic-sedimentary sequence is composed of metabasic and metavolcaniclastic rocks, graphitic and sulfur-rich metapelites, banded iron formation, metandesite, metarhyolite, calcsilicates, tourmalinites and metahydrothermalites derived from mafic and felsic rocks. The Mesoproterozoic rocks at Quartzito prospect are lithologically similar but they have been affected by Neoproterozoic faulting and shearing and by the emplacement of granitic rocks, resulting in the formation of tourmaline-rich quartz-carbonate veins with gold and base metal mineralization. We conducted a chemical and B-isotope study of tourmalines in order to better understand the origin of the stratiform tourmalinites in the Morro da Pedra Preta Formation and their relationship with gold mineralization. The overall range of delta(11)B values obtained for the tourmalinite and vein tourmalines is between - 15%. and -5 parts per thousand, with the tourmalinites failing at the low end of this range (-15 to -8 parts per thousand). Such values are typical for continental crust and inconsistent with a primary marine boron signature as expected from the submarine-exhalative model for the gold prospects. We conclude from this that tourmaline formed or recrystallized from crustal fluids related to the amphibolite-grade metamorphism which affected the Serra do Itaberaba Group and that gold deposition occurred syn- to post-peak metamorphism by phase immiscibility, as attested by fluid inclusions in Tapera Grande tourmalinite tourmaline and quartz. The vein-hosted tourmalines at Quartzito have isotopically variable boron signatures, with heavier delta(11)B values of -5 parts per thousand to -8 parts per thousand for acicular green tourmalines and lighter values (-15 parts per thousand to -7 parts per thousand for light blue, Ti-firee tourmaline from quartz-carbonate veins). We attribute the heavier boron to fluids derived from the volcano-sedimentary rocks of marine affinity whereas the lighter boron was contributed by crustal fluids related to the granitoids or metasediments in the continental crust. (c) 2009 Elsevier B.V. All rights reserved.

Sr and Nd isotopic geochemistry of the early ultramafic-mafic rocks of the Mako bimodal volcanic belt of the Kedougou-Kenieba inlier (Senegal)

The Mako bimodal volcanic belt of the Kedougou-Kenieba inlier is composed of volcanic basalts and peridotites interbedded by quartzites and limestones intruded by different generations of granitoids. The early volcanic episode of the belt is constituted of submarine basalts with peridotite similar to those of the oceanic abyssal plains. It is intruded by the Badon Kakadian TTG-granitic batholite dated around 2200 Ma. The second volcanic phase is constituted of basaltic, andesitic, and felsitic flows exhibit structures of aerial volcanic rocks. It is intruded by granites dated between 2160 and 2070 Ma. The general pattern of trace element variation of submarine volcanic rocks is consistent with those of basalts from oceanic plateaus which are the modern equivalent of the Archean greenstones belts. The Nd and Sr isotopic systematics typical of juvenile material indicates that the source of these igneous rocks is derived from a depleted mantle source. These results are consistent with the idea of a major accretion within the West African Craton occurring at about 2.1 Ga and corresponding to an important process of mantle-oceanic lithosphere differentiation.

Geologia e geoquimica do complexo cambaizinho, sao gabriel - RS

Mapeamentos geológicos feitos pelo autor na região do Arroio Cambaizinho - são Gabriel/RS, resultaram na definição do Complexo Cambaizinho, representado pelas seqüências meta-sedimentar e máfica-ultramáfica, intimamente associadas ao longo de toda a extensão da associação supracrustal. A seqüência meta-se dimentar é constitu1da por gnaisses quartzo-feldspáticos dominan tes, anfibolitos bandados e quartzitos subordinados, derivados de sedimentos areno-pelítico-carbonatados estruturados, de forma ritmica em ambiente subaquoso. Niveis composicionais de ocorrência restrita contendo estaurolita, definem o grau metamórfico (médio) para esta região. Intercalações de serpentinitos, xistos magnesianos variados e anfibolitos a granulação fina, na forma de camadas e/ou lentes interestratificadas nos meta-sedimentos indicam suas derivações a partir de derrames e/ou intrusões , ígneas de pequena profundidade de composição básica-ultrabásica. Este juntamente com níveis de sedimentos qu1micos intercalados e corpos de gabros, constituem a seqüência máfica-ultramáfica. O complexo, representa o segmento norte de um cintu rão supracrustal polideformado de forma geométrica aproximadamente linear, com orientação NNE, que se extende desde a localidade de Passo do Ivo, situado mais a sul, até a região objeto deste trabalho. Quatro fases de deformação dúcteis foram das para a área, estando as duas primeiras (Dl e D2) identifica associadas xx aos eventos metamórficos regionais Ml e M2. O metamorfismo mais antigo (Ml), assinalado por paragêneses diagnósticas em metapelitos alcançou o fácies anfibolito (zona da estaurolita), estando representado em outras litologias pela ocorrência de olivina metamórfica em paragênese com tremolita e/ou talco (meta-serpentinitos) e hornblenda mais oligoclásio/andesina em meta-básicas. O M2, mais jovem, atingiu o fácies xistos verdes, cujas assembléias mineralógicas se associam à foliação S2, de distribuição ir regular ao longo do cinturão. As condições fisicas de Ml foram de média P/T, similares às do metamorfismo Dalradiano. Intrusões graniticas na forma de lâminas (corpos ta bulares) durante a segunda fase de deformação D2, datados pelo mé todo Rb/Sr em 661 :: 29 Ma e agrupados sob a denominação de Granatóides Sanga do Jobim, fornecem idades mínimas para o complexo. Os vários grupos composiciconais da seqüência máfic~ -ultramáfica, individualizados com base em critérios petrográf~ cos e conteúdo de elementos maiores correspondem a: serpentinitos e olivina-talco ultramafitos (cumulados komatilticos); xis tos magnesianos à talco e clorita e anfibólio xistos (komatiitos); clorita-hornblenda xistos (basaltos komatiiticos), litos e metagabros (basaltos e gabros tolelticos). anfibo Estes vários tipos litológicos foram originados através de diferentes graus de fusão parcial do manto como indicado pelo hiato composicional de MgO (11 à 17%) e os diferentes padrões de ETR existentes entre os anfibolitos/metagabros (toleitos) e os serpentinitos/xistos magnesianos (komatiltos).As variações composicionais no interior de cada grupo, foram controladas pelo fracionamento (acumulação ou extração) de olivina e pouco ortopiroxênio (serpentinitos e olivina-talco ultramafitos)clinopiroxênios (clorita e anfibólios xistos, clorita hornblenda xistos), clinopiroxênio e plagioclásio (anfibolitos e metagabros). As abundâncias e os padrões de ERTL (elementos terras raras leves) enriquecidos, juntamente com os baixos valores das razões A1203/Ti02 e CaO/Ti02 das amostras de xistos magnesiinos das camadas A e B sugerem derivações deste material a partir de baixas percentagens de fusão de um manto enriquecido em mentos incompatíveis. As anomalias negativas de Ce e Eu na ele maio ria das rochas da seqüência máfica-ultramáfica indicam que protólitos ígneossofreram alterações em ambiente submarino.

Subsídios para a gestão dos planos de aplicação de vinhaça (PAV): um estudo de caso da região de Piracicaba

Concomitantemente à produção de etanol, a partir da fermentação de diferentes matérias-primas (cana-de-açúcar, beterraba, milho, trigo, batata, mandioca), produz-se a vinhaça. A vinhaça da cana-de-açúcar é gerada na proporção média de 12 litros para cada litro de etanol. O Brasil é o maior produtor mundial de cana-de-açúcar (570 milhões de toneladas, em 2009) com a produção de 27 bilhões de litros/ano de etanol, no ciclo 2009/2010, basicamente para fins carburantes, e, portanto, a quantidade de vinhaça produzida é da ordem de 320 bilhões de litros/ano. Dentre as possíveis soluções para a destinação da vinhaça, tais como: concentração, tratamento químico e biológico ou produção de biomassa, é sua aplicação no solo a forma mais usual de disposição. No entanto, sua aplicação no solo não pode ser excessiva nem indiscriminada, sob pena de comprometer o ambiente e a rentabilidade agrícola e industrial da unidade sucroalcooleira. A necessidade de medidas de controle sobre a aplicação de vinhaça no solo do Estado de São Paulo, que concentra a produção de 60% do etanol produzido no Brasil, levou à elaboração da Norma Técnica P4.231, em 2005. A constatação do cumprimento dos itens desta Norma é uma das ações do agente fiscalizador, no monitoramento da geração e destinação da vinhaça pelos empreendimentos produtores de etanol. Os itens 5.7.1 e 5.7.2, desta Norma, tornam obrigatório o encaminhamento à CETESB (Companhia Ambiental do Estado de São Paulo), para fins de acompanhamento e fiscalização, até o dia 02 de abril de cada ano, o PAV (Plano de Aplicação de Vinhaça). Com o intuito de proporcionar melhor entendimento da Norma P4.231, contribuir para otimização de sua aplicabilidade e melhoria e perceber a realidade prática da aplicação da vinhaça, foram analisados PAVs protocolados na Agência Ambiental de Piracicaba, licenças concedidas, processos de licenciamento e realizadas vistorias a campo. Conclui-se que a Norma P4.231 representa um avanço no gerenciamento do uso da vinhaça no Estado de São Paulo, por disciplinar a disposição de vinhaça no solo, tornando obrigatórios: demarcações de áreas protegidas e núcleos populacionais, caracterização de solo e da vinhaça, doses máximas a serem aplicadas, estudos da geologia e hidrogeologia locais, monitoramento das águas subterrâneas, impermeabilização de tanques e dutos. Estabelece critérios a serem obedecidos por lei, e todos são condutas de boas práticas de proteção ao meio, que repercutem em maior rentabilidade agrícola e industrial. A Norma P4.231 é bem elaborada, complexa e extensa, o que a torna de difícil execução. Para facilitar sua aplicabilidade, sugere-se o estabelecimento de um cronograma de prioridades: impermeabilização dos tanques e dutos, drenos testemunha com funcionamento otimizado e boa caracterização do que está sendo aplicado, se vinhaça pura ou se misturada a águas residuárias. Das oito empresas deste estudo, apenas quatro (as maiores) vêm protocolando os PAVs com regularidade, desde 2005. As informações apresentadas foram incompletas e, em alguns casos, precárias. Mesmo com lacunas, os dados fornecidos, desde o início da obrigatoriedade do PAV, foram de utilidade para esboçar o perfil da atividade sucroalcooleira na região estudada.

Previsão de desvios em arquiteturas multitarefas simultâneas

A exploração do paralelismo no nível de instrução (ILP) em arquiteturas superescalares é limitada fortemente pelas dependências de controle, as quais são ocasionadas pelas instruções de desvio, e pelas dependências de dados. As arquiteturas SMT (Simultaneous MultiThreaded) buscam explorar um novo nível de paralelismo, denominado paralelismo no nível de tarefa (TLP), para buscar e executar instruções de diversas tarefas ao mesmo tempo. Com isso, enquanto uma tarefa está bloqueada por dependências de controle e de dados, outras tarefas podem continuar executando, mascarando assim as latências de previsões incorretas e de acessos à memória, usando mais eficientemente as unidades funcionais e demais recursos disponíveis. Contudo, o projeto dessas arquiteturas continua a esbarrar nos mesmos problemas associados ao uso de técnicas utilizadas para a exploração de ILP, como a previsão de devios. Além disso, essas arquiteturas trazem novos desafios, como a determinação da maneira mais eficiente de distribuição/compartilhamento de recursos entre as tarefas. Nesse trabalho será apresentada uma topologia para as tabelas de previsão de desvios em arquiteturas multitarefas simultâneas. Além disso, serão desenvolvidas duas análises complementares acerca de previsão de desvios: o impacto da taxa de acertos da previsão de desvios em arquiteturas com pipelines profundos e o impacto da taxa de acerto na previsão do alvo de um desvio. Entre as principais contribuições do trabalho pode-se citar a definição de uma estrutura particionada para as tabelas de previsão em arquiteturas SMT, aliando desempenho a um menor custo de implementação em uma arquitetura real. Além disso, é mostrado que a taxa de acerto da previsão de desvios tem um grande impacto no desempenho das arquiteturas SMT com pipelines profundos, bem como nas causas de bloqueio do estágio de busca quando utiliza-se cache de instruções bloqueantes.

Desenvolvimento de tubulação em compósito polimérico revestida externamente com poliuretano de alta densidade

Pipelines for the transport of crude oil from the production wells to the collecting stations are named production lines . These pipes are subjected to chemical and electrochemical corrosion according to the environment and the type of petroleum transported. Some of these lines, depending upon the composition of the fluid produced, may leak within less than one year of operation due to internal corrosion. This work aims at the development of composite pipes with an external protecting layer of high density polyurethane for use in production lines of onshore oil wells, meeting operational requirements. The pipes were manufactured using glass fibers, epoxy resin, polyester resin, quartz sand and high density polyurethane. The pipes were produced by filament winding with the deposition of high density polyurethane on the external surface and threaded ends (API 15 HR/PM-VII). Three types of pipes were manufactured: glass/epoxy, glass/epoxy with an external polyurethane layer and glass/epoxy with an intermediate layer of glass fiber, polyester, sand and with an external polyurethane layer. The three samples were characterized by Scanning Electronic Microscopy (SEM) and for the determination of constituent content. In addition, the following tests were conducted: hydrostatic test, instant rupture, shorttime failure pressure, Gardner impact, transverse stiffness and axial tension. Field tests were conducted in Mossoró RN (BRAZIL), where 1,677 meters of piping were used. The tests results of the three types of pipes were compared in two events: after two months from manufacturing of the samples and after nine months of field application. The results indicate that the glass/epoxy pipes with an intermediate layer of fiber glass composite, polyester e sand and with an external layer of high density polyurethane showed superior properties as compared to the other two and met the requirements of pressure class, axial tensile strength, transverse stiffness, impact and environmental conditions, for onshore applications as production lines

Implementação de um Módulo de Supervisão para um Sistema de Detecção de Vazamentos em Dutos de Petróleo

The transport of fluids through pipes is used in the oil industry, being the pipelines an important link in the logistics flow of fluids. However, the pipelines suffer deterioration in their walls caused by several factors which may cause loss of fluids to the environment, justifying the investment in techniques and methods of leak detection to minimize fluid loss and environmental damage. This work presents the development of a supervisory module in order to inform to the operator the leakage in the pipeline monitored in the shortest time possible, in order that the operator log procedure that entails the end of the leak. This module is a component of a system designed to detect leaks in oil pipelines using sonic technology, wavelets and neural networks. The plant used in the development and testing of the module presented here was the system of tanks of LAMP, and its LAN, as monitoring network. The proposal consists of, basically, two stages. Initially, assess the performance of the communication infrastructure of the supervisory module. Later, simulate leaks so that the DSP sends information to the supervisory performs the calculation of the location of leaks and indicate to which sensor the leak is closer, and using the system of tanks of LAMP, capture the pressure in the pipeline monitored by piezoresistive sensors, this information being processed by the DSP and sent to the supervisory to be presented to the user in real time

Avaliação da corrosão em dutos por técnica gravimétrica e de resistência elétrica

With the increasing of demand for natural gas and the consequent growth of the pipeline networks, besides the importance of transport and transfer of oil products by pipeline, and when it comes to product quality and integrity of the pipeline there is an important role regarding to the monitoring internal corrosion of the pipe. This study aims to assess corrosion in three pipeline that operate with different products, using gravimetric techniques and electrical resistance. Chemical analysis of residues originated in the pipeline helps to identify the mechanism corrosive process. The internal monitoring of the corrosion in the pipelines was carried out between 2009 and 2010 using coupon weight loss and electrical resistance probe. Physico-chemical techniques of diffraction and fluorescence X-rays were used to characterize the products of corrosion of the pipelines. The corrosion rate by weight loss was analyzed for every pipeline, only those ones that has revealed corrosive attack were analyzed located corrosion rate. The corrosion potential was classified as low to pipeline gas and ranged from low to severe for oil pipelines and the pipeline derivatives. Corrosion products were identified as iron carbonate, iron oxide and iron sulfide

Simulador de escoamento multifásico em poços de petróleo (SEMPP)

The multiphase flow occurrence in the oil and gas industry is common throughout fluid path, production, transportation and refining. The multiphase flow is defined as flow simultaneously composed of two or more phases with different properties and immiscible. An important computational tool for the design, planning and optimization production systems is multiphase flow simulation in pipelines and porous media, usually made by multiphase flow commercial simulators. The main purpose of the multiphase flow simulators is predicting pressure and temperature at any point at the production system. This work proposes the development of a multiphase flow simulator able to predict the dynamic pressure and temperature gradient in vertical, directional and horizontal wells. The prediction of pressure and temperature profiles was made by numerical integration using marching algorithm with empirical correlations and mechanistic model to predict pressure gradient. The development of this tool involved set of routines implemented through software programming Embarcadero C++ Builder® 2010 version, which allowed the creation of executable file compatible with Microsoft Windows® operating systems. The simulator validation was conduct by computational experiments and comparison the results with the PIPESIM®. In general, the developed simulator achieved excellent results compared with those obtained by PIPESIM and can be used as a tool to assist production systems development

Ecogenômica de archaea e monitoramento de comunidades de procariotos redutores de sulfato: aplicações na indústria de petróleo e gás

The human activities responsible for the ambient degradation in the modern world are diverse. The industrial activities are preponderant in the question of the impact consequences for brazilian ecosystems. Amongst the human activities, the petroliferous industry in operation in Potiguar Petroliferous Basin (PPB) displays the constant risk of ambient impacts in the integrant cities, not only for the human populations and the environment, but also it reaches the native microorganisms of Caatinga ground and in the mangrove sediment. Not hindering, the elaboration of strategies of bioremediation for impacted areas pass through the knowledge of microbiota and its relations with the environment. Moreover, in the microorganism groups associated to oil, are emphasized the sulfate-reducing prokaryotes (SRP) that, in its anaerobic metabolism, these organisms participate of the sulfate reduction, discharging H2S, causing ambient risks and causing the corrosion of surfaces, as pipelines and tanks, resulting in damages for the industry. Some ancestries of PRS integrate the Archaea domain, group of microorganisms whose sequenced genomes present predominance of extremophilic adaptations, including surrounding with oil presence. This work has two correlated objectives: i) the detection and monitoring of the gene dsrB, gift in sulfate-reducing prokaryotes, through DGGE analysis in samples of mDNA of a mangrove sediment and semiarid soil, both in the BPP; ii) to relate genomic characteristics to the ecological aspects of Archaea through in silico studies, standing out the importance to the oil and gas industry. The results of the first work suggest that the petrodegraders communities of SRP persist after the contamination with oil in mangrove sediment and in semiarid soil. Comparing the populations of both sites, it reveals that there are variations in the size and composition during one year of experiments. In the second work, functional and structural factors are the probable cause to the pressure in maintenance of the conservation of the sequences in the multiple copies of the 16S rDNA gene. Is verified also the discrepancy established between total content GC and content GC of the same gene. Such results relating ribosomal genes and the ambient factors are important for metagenomic evaluations using PCR-DGGE. The knowledge of microbiota associated to the oil can contribute for a better destination of resources by the petroliferous industry and the development of bioremediation strategies. Likewise, search to lead to the best agreement of the performance of native microbiota in biogeochemical cycles in Potiguar Petroliferous Basin ecosystem

Influence of a Brazilian sewage outfall on the toxicity and contamination of adjacent sediments

The submarine sewage outfall of Santos (SSOS) is situated in the Santos Bay (São Paulo, Brazil) and is potentially a significant source of contaminants to the adjacent marine ecosystem. The present study aimed to assess the influence of SSOS on the sediment toxicity and contamination at Santos Bay. At the disposal site, sediments tended to be finer, organically richer and exhibited higher levels of surfactants and metals, sometimes exceeding the Threshold Effect Level values. The SSOS influence was more evident toward the East, where the sediments exhibited higher levels of TOC, total S and metals during the summer 2000 sampling campaign. Sediment toxicity to amphipods was consistently detected in four of the five stations studied. Amphipod survival tended to correlate negatively to Hg, total N and % mud. This work provides evidence that the SSOS discharge affects the quality of sediments from Santos Bay, and that control procedures are warranted. (c) 2005 Elsevier Ltd. All rights reserved.

ESR dating of pleistocene mammal teeth and its implications for the biostratigraphy and geological evolution of the coastal plain, Rio Grande do Sul, southern Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)