Aplicação dos métodos eletromagnéticos VLF e HLEM na prospecção hidrogeológica no município de São Domingos do Araguaia-PA

Inserido no convênio de cooperação firmado entre a Fundação Nacional de Saúde e a Universidade Federal do Pará, através do Departamento de Geofísica e do Curso de Pós-Graduação em Geofísica, este trabalho tem como objetivo o mapeamento geológico com a aplicação dos métodos eletromagnéticos VLF e HLEM, visando o estudo hidrogeológico do Município de São Domingos do Araguaia no Sudeste do Estado do Pará. Também, estudamos o desempenho quantitativo do método VLF, no que diz respeito a prospecção de água subterrânea. Primeiramente, analisamos a resposta VLF, através dos parâmetros da elipse de polarização (inclinação e elipsidade), de diferentes situações geológicas possíveis de serem encontradas no local de estudo e voltadas a hidrogeologia. Também, avaliamos a influência que os parâmetro físicos do modelo exercem na anomalia, utilizando para estas finalidades modelagem numérica por elemento finito. A análise quantitativa foi realizada através de ajuste, por tentativa e erro, com curvas teóricas obtidas da simulação de modelos idealizados na interpretação qualitativa, a qual fizemos com o auxilio da técnica de filtragem de Fraser e com perfís do método HLEM. Os resultados mostraram a existência de feições estruturais de "trends" N-S e NE-SW, que correlacionadas à geologia regional são falhas normais e de transferências, respectivamente. Com base nesses resultados, foram selecionados alguns locais com maiores probabilidades de sucesso na captação de água subterrânea através de poços.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Modelamento de perfis de indução

Neste trabalho foi implementado um algoritmo que permite o modelamento das respostas das sondas de indução (2C40-6FF40-ILD) em meios estratificados. O algoritmo é baseado na solução exata do problema de Green para a obtenção do campo eletromagnético em meios estratificados com a fonte no interior dos estratos. O código computacional desenvolvido neste trabalho para meios horizontalmente estratificados tem o objetivo de auxiliar o desenvolvimento de códigos mais avançados, uma vez que é difícil o acesso para a pesquisa dos códigos citados na literatura. Os códigos disponíveis são comerciais, mas estes têm o objetivo único de produção, além de serem, bastante dispendiosos. O código também foi utilizado para a obtenção da resposta da sonda ILD em reservatórios do tipo laminado, onde a sonda tem baixa resolução. E verificamos, a exemplo de outros autores, que a condutividade obtida pela sonda pode ser expressa pela soma das condutividades dos constituintes litológicos do reservatório ponderada pelos seus percentuais em volume. Esta relação foi definida para reservatórios ainda não analisados na literatura, tipo arenito moderadamente silicificado/arenito/folhelho, onde a lâmina de interesse tem resistividade superior a resistividade de uma das lâminas, e inferior a resistividade da outra. Foi analisado também reservatórios tipo arenito moderadamente silicificado/arenito, onde a lâmina de interesse tem a menor resistividade. Com esta relação pode-se obter para estes reservatórios laminados estimativas mais precisas de Rt e consequentemente determinar as saturações água/óleo e produtibilidade, também de forma mais precisa.

Análise dos impactos das diferentes formas de ocupação da superfície sobre as condições meteorológicas na região de Santarém, Pará

O objetivo desta dissertação foi obter informações referentes ao uso e ocupação do solo na região de Santarém, em diferentes anos das últimas décadas, para melhor representar os efeitos causados pelas modificações das propriedades da superfície sobre as condições atmosféricas simuladas por modelos numéricos de tempo e clima. As superfícies continentais caracterizam-se por causar efeitos substanciais sobre a atmosfera e, consequentemente, influir na qualidade das previsões de tempo e de clima. Por outro lado, o desmatamento contribui com as mudanças climáticas, por eliminar grandes quantidades de gases de efeito estufa para a atmosfera. Estas atividades também causam efeitos na saúde publica, na agricultura, nos recursos florestais, nos recursos faunísticos e nos recursos hídricos. Além disso, a substituição da superfície natural por pastagem ou agricultura altera as propriedades térmicas e radiativas da superfície, gerando modificações nas condições atmosféricas locais, regionais e globais. Neste trabalho foram analisados períodos representativos de possíveis mudanças climáticas na região, identificados a partir do tratamento e analise estatística de dados climatológicos de estações meteorológicas de superfície, bem como a evolução temporal e quantitativa do desmatamento na região de estudo com os dados do Projeto PRODES Monitoramento da Floresta Amazônica Brasileira por Satélite). Para avaliar os efeitos atmosféricos das mudanças no uso e ocupação do solo utilizou-se como base o mapa de vegetação do IBGE, e a inclusão da classe “desmatamento” ao mesmo em diferentes períodos analisados (anos de 1997e 2009) trabalhadas no software Arc. Gis. 9.2. Foram criados arquivos de dados de tipos de superfície compatíveis com a leitura do modelo BRAMS, que foi então utilizado para simular os diferentes efeitos desses mapas temáticos de uso e ocupação do solo na atmosfera local. Os resultados indicam uma tendência de aumento da precipitação média anual e da frequência média de dias com precipitação, diminuição da temperatura média das máximas e aumento da temperatura média das mínimas ao longo dos anos na região de Santarém. A área de estudo, até o ano de 1997, registrou um desmatamento de 19,44% e até o ano de 2009 passou para 25,54%. As simulações com os arquivos gerados de uso e ocupação do solo para 1997 e 2009 apresentaram poucas variações para os diferentes mapas temáticos em suas variáveis (temperatura, umidade e fluxos de calor sensível e latente), quando considerado os valores médios da área total simulada. Porém, quando se considera pequenas áreas localizadas somente sobre as regiões que sofreram maiores modificações, observam-se maiores influências com o aumento do desmatamento.

Análise do escoamento turbulento por um queimador industrial a gás utilizando dinâmica dos fluidos computacional

Pós-graduação em Engenharia Mecânica - FEG

Estudo da efetividade da atenuação natural de compostos BTEX em área contaminada por querosene de aviação

Pós-graduação em Geociências e Meio Ambiente - IGCE

Índices para a representação da paisagem como apoio para levantamento pedológico em ambiente de geoprocessamento

O mapeamento de solos tem ganhado destaque dentro da comunidade científica, pois, à medida em que a preocupação com o meio ambiente aumenta, ocorre a necessidade de se entender cada vez mais a distribuição dos solos na paisagem, tal como seu potencial e limitações de uso. Desta forma, o principal objetivo do trabalho foi aplicar índices de representação da paisagem com o apoio de geoprocessamento, suporte na delimitação dos diferentes compartimentos da paisagem. Foram utilizados índices primários Altitude above channel network (AACN) e secundários Channel network base level (CNBL, Multiresolution índex of valley bottom flatness (MRVBF) e Wetness índex (ITW), cujo estudo foi a Fazenda Experimental Canguiri, no município de Pinhais, região Metropolitana de Curitiba. Para correlacionar os atributos químicos e granulométricos por grupo de amostragem, totalizando 17 pontos (Sugamosto, 2002) foi gerada, no Software Statistica, uma matriz de correlação linear simples (Pearson) com os índices de representação da paisagem, os quais foram não apenas utilizados na análise de agrupamentos mas também eficientes no apoio do mapeamento dos solos, a nível de subordem do Sistema Brasileiro de Classificação de Solo.

Aplicação de parâmetros convectivos para um caso de linha de instabilidade

Este trabalho visa avaliar a aplicação de parâmetros convectivos (PC) para um caso de linha de instabilidade (LI) ocorrida sobre a Região Metropolitana de São Paulo (RMSP), com o intuito de demonstrar como o uso da modelagem numérica, utilizando o modelo BRAMS (Brazilian developments on the Regional Atmospheric Modeling System), pode contribuir principalmente com a obtenção de perfis atmosféricos, necessários para uma boa previsão das tempestades.

Impacto da temperatura da superfície do mar na simulação da Zona de Convergência do Atlântico Sul

Durante o verão austral, diversas regiões do Brasil são afetadas por precipitação intensa, geralmente associada à Zona de Convergência do Atlântico Sul (ZCAS). O objetivo geral deste trabalho foi investigar a influência da resolução espacial e temporal dos dados de temperatura da superfície do mar (TSM) na simulação da precipitação associada à ZCAS. Foram realizadas simulações com o modelo BRAMS (Brazilian developments on the Regional Atmospheric Modeling System) para dois eventos de ZCAS, ocorridos em 1998 (El Niño) e 1999 (La Niña). A quantidade de precipitação acumulada na parte oceânica da ZCAS foi maior nos experimentos com TSM mais quente. Índices estatísticos foram utilizados para verificação do desempenho do modelo na simulação de precipitação nas regiões que compõem a ZCAS (oceânica, costeira e amazônica), com diferentes dados de TSM. A resolução espaço-temporal dos dados de TSM influencia de forma pouco significativa na representação da ZCAS pelo modelo BRAMS. O modelo é mais eficiente em identificar a ocorrência/não ocorrência de chuva do que em localizar núcleos mais intensos e seu desempenho foi superior (inferior) na região amazônica (oceânica) da ZCAS.

Estudo das forças capilar, viscosa e gravitacional através de grupos adimensionais na simulação da drenagem gravitacional assistida por gás (GAGD)

The distribution and mobilization of fluid in a porous medium depend on the capillary, gravity, and viscous forces. In oil field, the processes of enhanced oil recovery involve change and importance of these forces to increase the oil recovery factor. In the case of gas assisted gravity drainage (GAGD) process is important to understand the physical mechanisms to mobilize oil through the interaction of these forces. For this reason, several authors have developed physical models in laboratory and core floods of GAGD to study the performance of these forces through dimensionless groups. These models showed conclusive results. However, numerical simulation models have not been used for this type of study. Therefore, the objective of this work is to study the performance of capillary, viscous and gravity forces on GAGD process and its influence on the oil recovery factor through a 2D numerical simulation model. To analyze the interplay of these forces, dimensionless groups reported in the literature have been used such as Capillary Number (Nc), Bond number (Nb) and Gravity Number (Ng). This was done to determine the effectiveness of each force related to the other one. A comparison of the results obtained from the numerical simulation was also carried out with the results reported in the literature. The results showed that before breakthrough time, the lower is the injection flow rate, oil recovery is increased by capillary force, and after breakthrough time, the higher is the injection flow rate, oil recovery is increased by gravity force. A good relationship was found between the results obtained in this research with those published in the literature. The simulation results indicated that before the gas breakthrough, higher oil recoveries were obtained at lower Nc and Nb and, after the gas breakthrough, higher oil recoveries were obtained at lower Ng. The numerical models are consistent with the reported results in the literature

Transporte longitudinal de sedimento na zona costeira de Natal/RN

The longshore sediment transport (LST) is determinant for the occurrence of morphological changes in coastal environments. Understanding their movement mechanisms and transport is an essential source of information for the project design and coastal management plans. This study aims to characterize, initially, the active hydrodynamic circulation in the study area, comprised of four beach sectors from the south coast of Natal, assessing the average annual LST obtained through three proven equations (CERC, Kamphuis and Bayram et al.), defining the best formulation for the study area in question, and analyze the seasonal variability and the decadal transport evolution. The coastal area selected for this work constitutes one of the main tourist corridors in the city, but has suffered serious damage resulting from associated effects of hydrodynamic forcings and their disorderly occupation. As a tool was used the Coastal Modelling System of Brazil (SMC-Brazil), which presents integrated a series of numerical models and a database, properly calibrated and validated for use in developing projects along the Brazilian coastline. The LST rates were obtained for 15 beach profiles distributed throughout the study area. Their extensions take into account the depth of closure calculated by Harllermeier equation, and regarding the physical properties of the sediment, typical values of sandy beaches were adopted, except for the average diameter, which was calculated through an optimization algorithm based on equilibrium profile formulation proposed by Dean. Overall, the results showed an intensification of hydrodynamic forcings under extreme sea wave conditions, especially along the headlands exist in the region. Among the analyzed equations, Bayram et al. was the most suitable for this type of application, with a predominant transport in the south-north direction and the highest rates within the order of 700.000 m3 /year to 2.000.000 m3 /year. The seasonal analysis also indicated a longitudinal transport predominance in the south to north, with the highest rates associated with the fall and winter seasons. In these periods are observed erosive beach states, which indicate a direct relationship between the sediment dynamics and the occurrence of more energetic sea states. Regarding the decadal evolution of transportation, it was found a decrease in transport rate from the 50’s to the 70’s, followed by an increase until the 2000’s, coinciding with the beginning of urbanization process in some stretches of the studied coastline.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

This work approaches the forced air cooling of strawberry by numerical simulation. The mathematical model that was used describes the process of heat transfer, based on the Fourier's law, in spherical coordinates and simplified to describe the one-dimensional process. For the resolution of the equation expressed for the mathematical model, an algorithm was developed based on the explicit scheme of the numerical method of the finite differences and implemented in the scientific computation program MATLAB 6.1. The validation of the mathematical model was made by the comparison between theoretical and experimental data, where strawberries had been cooled with forced air. The results showed to be possible the determination of the convective heat transfer coefficient by fitting the numerical and experimental data. The methodology of the numerical simulations was showed like a promising tool in the support of the decision to use or to develop equipment in the area of cooling process with forced air of spherical fruits.

Aproximação numérica – analítica para a modelagem da conversão termoquímica de combustíveis sólidos

Um algoritmo numérico foi criado para apresentar a solução da conversão termoquímica de um combustível sólido. O mesmo foi criado de forma a ser flexível e dependente do mecanismo de reação a ser representado. Para tanto, um sistema das equações características desse tipo de problema foi resolvido através de um método iterativo unido a matemática simbólica. Em função de não linearidades nas equações e por se tratar de pequenas partículas, será aplicado o método de Newton para reduzir o sistema de equações diferenciais parciais (EDP’s) para um sistema de equações diferenciais ordinárias (EDO’s). Tal processo redução é baseado na união desse método iterativo à diferenciação numérica, pois consegue incorporar nas EDO’s resultantes funções analíticas. O modelo reduzido será solucionado numericamente usando-se a técnica do gradiente bi-conjugado (BCG). Tal modelo promete ter taxa de convergência alta, se utilizando de um número baixo de iterações, além de apresentar alta velocidade na apresentação das soluções do novo sistema linear gerado. Além disso, o algoritmo se mostra independente do tamanho da malha constituidora. Para a validação, a massa normalizada será calculada e comparada com valores experimentais de termogravimetria encontrados na literatura, , e um teste com um mecanismo simplificado de reação será realizado.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

Este trabalho abordou o resfriamento rápido com ar forçado de morango via simulação numérica. Para tanto, foi empregado o modelo matemático que descreve o processo de transferência de calor, com base na lei de Fourier, escrito em coordenadas esféricas e simplificado para descrever o processo unidimensional. A resolução da equação expressa pelo modelo matemático deu-se por meio da implementação de um algoritmo, fundamentado no esquema explícito do método numérico das diferenças finitas, executado no ambiente de computação científica MATLAB 6.1. A validação do modelo matemático foi realizada a partir da comparação de dados teóricos com dados obtidos num experimento, no qual morangos foram resfriados com ar forçado. Os resultados mostraram que esse tipo de investigação para a determinação do coeficiente de transferência de calor por convecção é promissora como ferramenta no suporte à decisão do uso ou desenvolvimento de equipamentos na área de resfriamento rápido de frutos esféricos com ar forçado.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

Este trabalho abordou o resfriamento rápido com ar forçado de morango via simulação numérica. Para tanto, foi empregado o modelo matemático que descreve o processo de transferência de calor, com base na lei de Fourier, escrito em coordenadas esféricas e simplificado para descrever o processo unidimensional. A resolução da equação expressa pelo modelo matemático deu-se por meio da implementação de um algoritmo, fundamentado no esquema explícito do método numérico das diferenças finitas, executado no ambiente de computação científica MATLAB 6.1. A validação do modelo matemático foi realizada a partir da comparação de dados teóricos com dados obtidos num experimento, no qual morangos foram resfriados com ar forçado. Os resultados mostraram que esse tipo de investigação para a determinação do coeficiente de transferência de calor por convecção é promissora como ferramenta no suporte à decisão do uso ou desenvolvimento de equipamentos na área de resfriamento rápido de frutos esféricos com ar forçado.