Geoestatística, simulação estocástica e sensoriamento remoto na estimativa de produção do café conilon

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelação estocástica do depósito mineral do Zambujal (Mina de Neves-Corvo): contribuição da densidade como indicador morfológico

Este estudo tem como objetivo construir um modelo estocástico de alta resolução da morfologia e dos teores em metal, do depósito mineral do Zambujal, Mina de Neves-Corvo. O depósito do Zambujal é um corpo vulcanogénico que se localiza no setor Português, na parte mais a Sul da Faixa Piritosa Ibérica, com sulfuretos maciços no topo e fissurais na base. Para construir um modelo estocástico deste depósito, onde os teores em metal evidenciam zonamento, é importante ter em conta a proporção local de sulfuretos na matriz rochosa. Como esta variável não é quantificada em laboratório, propõe-se a utilização da densidade das amostras como um indicador indireto da proporção de minérios na matriz rochosa. Conhecida esta variável, a modelação dos teores em metal pode ser feita para os chamados teores relativos (teores em metal na fração de sulfuretos). As principais etapas da metodologia proposta são: (a) para cada amostra analisada no laboratório, estimação de soluções para a variável proporção de sulfuretos na matriz rochosa ????(????), tendo em conta a paragénese principal do depósito, os teores e a densidade da rocha; (b) construção de um modelo morfológico 3D de baixa resolução com duas regiões, minérios maciços e minérios fissurais, por digitalização de limites em perfis, interpolação de superfícies e conversão para o modelo de blocos do depósito; (c) construção de um modelo morfológico 3D de alta resolução para todo o depósito da variável ????(????) por Simulação Sequencial Direta (SSD), tendo como informação condicionante os histogramas regionais de ????(????) para minérios maciços e minérios fissurais e as soluções obtidas em (a) para ????(????) na localização das sondagens; (d) SSD dos teores relativos em cobre, zinco e prata; (e) discussão dos resultados e quantificação de recursos. Os resultados foram validados por comparação com os equivalentes obtidos numa estimação por krigagem normal dos teores em metal e mostraram ser da mesma ordem de grandeza. O conjunto de imagens simuladas das variáveis ????(????) e teores permite quantificar a incerteza do conhecimento do depósito relativamente à informação disponível.

O estudo geoestatístico de sondagens spt para geração de mapas auxiliares em obras de engenharia

A sondagem com SPT ("Standard Penetration Test") é o método de investigação com mais freqüência empregado na caracterização geotécnica de solos. Devido a relevância desse tipo de ensaio, desenvolveu-se um procedimento para geração de mapas de SPT com base em métodos geoestatísticos. Esses mapas devem conter a estimativa do índice de penetração (SPT) em diferentes posições no terreno em estudo, partindo de ensaios pré-existentes, objetivando auxiliar de maneira mais efetiva o planejamento de obras de engenharia. O trabalho foi baseado em um banco de dados com 141 sondagens realizadas na cidade de Passo Fundo/RS, compreendendo uma área total de 4,16km2. O banco de dados compreende sondagens de uma única empresa (o que minimiza os erros decorrentes da própria sondagem SPT) e uma única unidade geotécnica. Essas informações foram utilizadas para a avaliação do comportamento espacial dos parâmetros modelados e, posteriormente, construção de mapas de SPT. As estimativas foram efetuadas por krigagem ordinária em blocos com dimensões de 100mx100m, em quatro níveis distintos definidos de acordo com a variação média do índice de penetração com a profundidade Os mapas resultantes da krigagem ordinária revelaram-se satisfatórios para o fim a que se destina esse trabalho. Para a avaliação desses mapas, é relevante que se quantifiquem as incertezas associadas a essas estimativas. A metodologia para tal é baseada em Simulação Geoestatística, onde foram efetuadas simulações seqüenciais Gaussianas (ssG), em blocos de mesma dimensão que a utilizada na krigagem. Os resultados positivos mostraram ser adequado o emprego da simulação seqüencial Gaussiana, para geração de mapas auxiliares de valores de SPT a serem empregados no planejamento de investigações geotécnicas complementares.

Análise de incerteza associada à determinação da velocidade de onda sônica em depósito de carvão obtida por perfilagem geofísica

A modelagem de um depósito mineral é realizada por meio da combinação de diversas fontes de informações. Dentre estas fontes pode-se citar a sísmica de reflexão. A sísmica de reflexão fornece dados de tempos de propagação de ondas sísmicas até essas serem refletidas pelas estruturas dos depósitos minerais. As profundidades dessas estruturas podem ser obtidas multiplicando-se os tempos pelas velocidades de propagação das ondas sísmicas. Normalmente, a velocidade de uma onda sísmica é determinada indiretamente por meio do processamento dos próprios dados sísmicos, o que pode gerar erros na interpretação de seções geológicas. A perfilagem geofísica é uma alternativa na determinação dessa velocidade, uma vez que a velocidade de onda acústica é obtida ao longo do furo perfilado e a velocidade de onda acústica pode ser relacionada com a velocidade de onda sísmica. As estimativas de valores de velocidade na região entre os furos perfilados permite as estimativas de valores de profundidade nessa região. Neste estudo, foram analisadas possibilidades de se estimar, em um grid, valores de velocidade e incertezas associadas a esses valores, utilizando-se ferramentas geoestatísticas. A simulação seqüencial Gaussiana, dentre as ferramentas analisadas, pareceu ser a mais adequada para obtenção de velocidades e incerteza dos valores estimados em cada nó do grid considerado. Para o caso abordado, alguns valores de profundidade da estrutura de interesse apresentaram variações significativas em função da incerteza de velocidade. Essas variações são muito importantes na execução de certos métodos de lavra subterrânea, o que enfatiza a importância da determinação da incerteza das velocidades estimadas. A metodologia é apresentada e ilustrada em um importante depósito de carvão em Queensland, Austrália.

Análise multivariada de dados espaciais na classificação interpretativa de solos

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Avaliação do impacto do reservatório de Três Irmãos sobre a superfície potenciométrica do aquífero livre na cidade de Pereira Barreto (SP): uma abordagem numérica e geoestatística

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

O emprego da geoestatística na determinação do tamanho "ótimo" de amostras aleatórias com vistas à obtenção de estimativas dos volumes dos fustes de espécies florestais em Paragominas, estado do Pará

O objetivo do trabalho foi determinar o tamanho adequado de amostra para estimar o volume de fustes de espécies florestais de uma população de árvores a serem cortadas no sistema de manejo florestal da empresa Cikel Brasil Verde Madeiras - Pará. Utilizaram-se as metodologias da amostragem sistemática e do estimador geoestatístico da krigagem ordinária com simulação sequencial, respectivamente para a escolha das amostras e estimação dos volumes dos fustes das árvores. Os resultados mostraram que os métodos podem ser utilizados no cálculo dos volumes de fustes de árvores. Entretanto, o método da krigagem apresenta um efeito de suavização, tendo como conseqüência uma subestimação dos volumes calculados. Neste caso, um fator de correção foi aplicado para minimizar o efeito da suavização. A simulação sequencial indicativa apresentou resultados mais precisos em relação à krigagem, uma vez que tal método apresentou algumas vantagens, tal como a não exigência de amostras com distribuições normais e ausência de efeito de suavização, característico dos métodos de interpolação.

Simulação numérica de esteiras em transição utilizando o método dos contornos virtuais

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

Obtenção e identificação de grânulos de polihidroxialcanoatos (PHA) presentes em biomassa de um reator descontínuo sequencial (SBR) com tratamento de águas residuais

As implicações ambientais dos plásticos convencionais, produzidos através de combustíveis fósseis, não permitem a sua utilização de forma sustentável pela sociedade devido aos seus problemas de degradabilidade e à sua dependência de matérias-primas não renováveis. Surgem então os polihidroxialcanoatos (PHA) como uma alternativa ambientalmente sustentável, já que conseguem reproduzir grande parte das aplicações dos plásticos ditos normais, tendo as grandes vantagens de serem biodegradáveis e produzidos a partir de uma fonte de carbono renovável. Mais recentemente, de forma a tornar a produção do PHA mais sustentável, tem-se desenvolvido o conceito de associar a produção destes polímeros ao tratamento de efluentes líquidos através do processo de lamas ativadas. Apesar destas vantagens todo o processo que envolve a produção do PHA é ainda dispendioso. Esta dissertação insere-se na estratégia de tentar minorar esses custos. Uma das formas de o fazer será na otimização de técnicas de quantificação do PHA em culturas mistas. Neste trabalho o principal objetivo reside na avaliação da capacidade da metodologia de análise de imagem, acoplada a microscopia de fluorescência, como método de quantificação de PHA. Pretende-se também explorar a produção de PHA a partir da simulação de um sistema de tratamento de águas residuais implementado num reator descontínuo sequencial (SBR) em regimes cíclicos de feast/famine, aliando assim a função de tratamento de efluentes líquidos com a produção de um produto de valor acrescentado como é o PHA. Deste modo, foi montada uma instalação com um sistema de lamas ativadas num SBR cujo principal objetivo era produzir PHA. Este objetivo de produzir PHA através do regime de alimentação dinâmica aeróbia com uma alimentação com acetato como fonte de carbono, e sem controlo do pH foi conseguido, atingindo-se um máximo de concentração mássica de 0,37 mg PHA/mg SST e de concentração volumétrica de 1785 mg PHA/L. Em relação à capacidade de tratamento, comprovou-se a adequabilidade do sistema utilizado para a remoção da carga orgânica devido a uma elevada captação de nutrientes do afluente pela biomassa. Foram encontradas algumas limitações na metodologia de análise de imagem ao nível da transformação da bidimensionalidade da área projetada dos grânulos de PHA nas imagens para a sua tridimensionalidade real, na gama de concentração de PHA na qual a metodologia é válida e na amostragem relativamente curta. Contudo, foi possível comprovar a potencialidade deste método obtendo-se um R2 de 0,80 para o modelo de previsão da concentração mássica de PHA. Deste modo, pode-se concluir que esta é uma metodologia promissora para quantificar PHA em amostras de licor misto.

Educação para cidadania na simulação parlamentar : o caso do programa “deputado por um dia” da ALRS

Analisa o projeto "Deputado por um Dia", mantido desde 1999 pela Assembleia Legislativa do Rio Grande do Sul, que promove para os jovens estudantes a vivência da rotina dos deputados estaduais. A partir da recuperação de documentos e de entrevistas, analisa efeitos do programa sobre os jovens. Conclui-se que a experiência da decisão coletiva ajuda na percepção da complexidade do processo deliberativo e na criação de uma visão menos ingênua da política.

Modelagem, simulação e otimização operacional de unidades de síntese de metanol

A indústria de processos químicos tem sofrido consideráveis transformações devido ao acirramento da competitividade. Importantes progressos tecnológicos têm sido atingidos através de técnicas de modelagem, simulação e otimização visando o aumento da lucratividade e melhoria contínua nos processos industriais. Neste contexto, as plantas de metanol, um dos mais importantes produtos petroquímicos, podem ser destacadas. Atualmente, a principal matéria-prima para obtenção de metanol é o gás natural. A produção do metanol é caracterizada por três etapas: geração de gás de síntese, conversão do gás de síntese em metanol (unidade de síntese ou loop de síntese) e purificação do produto na especificação requerida. Os custos fixos e variáveis da unidade de síntese são fortemente dependentes das variáveis operacionais, como temperatura, pressão, razão de reciclo e composição da carga. Desta forma, foi desenvolvido um conjunto de modelos e algoritmos computacionais para representar matematicamente unidades de síntese de metanol. O modelo apresenta operações unitárias associadas aos seguintes equipamentos: divisores de correntes, misturadores de correntes, compressores, trocadores de calor, vasos de flash e reatores. Inicialmente, foi proposto um simulador estacionário, que serviu como base para um pseudo-estacionário, o qual contempla a desativação do catalisador por sinterização térmica. Os simuladores foram criados segundo uma arquitetura seqüencial modular e empregou-se o método de substituição sucessiva para a convergência dos reciclos. O estudo envolveu dois fluxogramas típicos, um constituído por reatores adiabáticos em série, enquanto o outro constituído por um reator tipo quench. Uma análise do efeito das principais variáveis operacionais foi realizada para o reator e para o loop de síntese. Estudou-se também o efeito da desativação do catalisador ao longo do tempo. Uma ferramenta de otimização operacional foi empregada para alcançar a máxima produção manipulando as injeções de carga fria na entrada dos leitos catalíticos. Formulou-se também um problema de maximização do lucro em função da otimização de paradas da unidade para troca do catalisador. Os resultados obtidos apontam que a ferramenta desenvolvida é promissora para a compreensão e otimização da unidade objeto deste estudo

Sistemas complexos, séries temporais e previsibilidade

Para qualquer sistema observado, físico ou qualquer outro, geralmente se deseja fazer predições para sua evolução futura. Algumas vezes, muito pouco é conhecido sobre o sistema. Se uma série temporal é a única fonte de informação no sistema, predições de valores futuros da série requer uma modelagem da lei da dinâmica do sistema, talvez não linear. Um interesse em particular são as capacidades de previsão do modelo global para análises de séries temporais. Isso pode ser um procedimento muito complexo e computacionalmente muito alto. Nesta dissertação, nos concetraremos em um determinado caso: Em algumas situações, a única informação que se tem sobre o sistema é uma série sequencial de dados (ou série temporal). Supondo que, por detrás de tais dados, exista uma dinâmica de baixa dimensionalidade, existem técnicas para a reconstrução desta dinâmica.O que se busca é desenvolver novas técnicas para poder melhorar o poder de previsão das técnicas já existentes, através da programação computacional em Maple e C/C++.

Simulação biofisicoquímica em reservatório tropical polimítico

O reservatório do Lobo, localizado no estado de São Paulo, é um sistema dinâmico no qual se desenvolve um ciclo diurno de estratificação e mistura, de modo similar ao que tem sido observado em outros lagos tropicais. Utilizou-se simulação 3D computacional com os softwares ELCOM (Estuary and Lake Computer Model) acoplado ao CAEDYM (Computacional Aquatic Ecosystem Dynamics Model), ambos desenvolvidos pelo CWR (Center for Water Research) da Universidade da Austrália. Foram realizadas cinco simulações: Piloto Primavera baseada em dados reais da estação no ano primavera no reservatório para o ano de 2007; Primavera-P em que as concentrações de fósforo total, fosfato inorgânico e fosfato total dissolvido foram aumentadas em 100% no reservatório (coluna de água e sedimento) e nos rios tributários; Primavera-V na qual a intensidade dos ventos foi aumentada em 50%; Primavera-T onde a temperatura da água (reservatório e tributários) e do ar foram aumentadas em 10C e, Primavera-X, onde a temperatura da água (reservatório e tributários) e do ar sofreu aumento em 10C, as concentrações de fósforo total, fosfato inorgânico e fosfato total dissolvido foram aumentadas em 100% e a velocidade do vento aumentada em 50%. A concentração de clorofila a foi representada pelos grupos cianobactérias e clorofíceas. O espaço de tempo das simulações representou 90 dias. As clorofíceas apresentaram maior desenvolvimento populacional do que as cianobactérias em todas as simulações. No reservatório, a mistura vertical é ocasionada diariamente pelo vento ou por processos convectivos causados pela perda de calor no corpo de água. A oxigenação do reservatório é maior com a ocorrência de ventos e de grupos fotossintéticos. As concentrações totais de fósforo e nitrogênio apresentaram aumento em todas as simulações.

Uso de ambientes virtuais para a simulação de acidentes radiológicos

Em um ambiente virtual, construído com o uso de tecnologia computacional, encontram-se presentes entidades virtuais inseridas em um espaço tridimensional, que é utilizado para a simulação de processos críticos, como os acidentes radiológicos. A pronta detecção de um acidente radiológico e a determinação da sua possível extensão são fatores essenciais para o planejamento de respostas imediatas e de ações de emergência. A integração das representações georeferenciadas do espaço tridimensional, com modelos baseados em agentes autônomos, com o objetivo de construir ambientes virtuais que tenham a capacidade de simular acidentes radiológicos é a proposta dessa tese. As representações georeferenciadas do espaço tridimensional candidatas são: i)as representações espaciais usadas nos sistemas de informações geográficas (SIG) e ii) a representação adotada pelo Google MapsTM. Com o uso deste ambiente pode-se: quantificar as doses recebidas pelas pessoas; ter uma distribuição espacial das pessoas contaminadas; estimar o número de indivíduos contaminados; estimar o impacto na rede de saúde; estimar impactos ambientais; gerar zonas de exclusão; construir cenários alternativos; treinar pessoal técnico para lidar com acidentes radiológicos.

Simulação com o código GEANT4 de medida de espessura de revestimento metálico em metal por XRF

Nesta dissertação são apresentados resultados de simulações Monte Carlo de fluorescência de raios X (XRF), utilizando o programa GEANT4, para medidas de espessura de revestimento metálico (Ni e Zn) em base metálica (Fe). As simulações foram feitas para dois tamanhos de espessura para cada metal de revestimento, (5μm e 10μm), com passos de 0,1 μm e 0,001 μm e com 106 histórias. No cálculo da espessura do revestimento foram feitas as aproximações de feixe de raios X monoenegético, com a análise da transmissão apenas da energia do K-alfa e para uma geometria compatível com um sistema real de medição (ARTAX-200). Os resultados mostraram a eficiência da metodologia de simulação e do cálculo da espessura do revestimento, o que permitirá futuros cálculos, inclusive para multirevestimentos metálicos em base metálica.