Análise da sensibilidade do seqüenciamento de lavra em função da incerteza do modelo geológico

Quantificação de incerteza e risco tem aplicação direta no projeto de limites de cava e na programação de produção. Pequenas variações nas condições de contorno de um projeto podem ter impacto significativo sobre o retorno final do mesmo na ordem de milhões de dólares. Preço de mercado do bem mineral, o custo de capital, taxas de atratividade, etc, são fatores usualmente testados para analisar a viabilidade de um empreendimento, porém, raramente é considerada a incerteza relacionada ao atributo geológico em questão. O propósito de um planejamento de lavra tem sido prover subsidio para o engenheiro de minas decidir sobre a capacidade de minerar determinadas unidades de lavra do depósito a partir de teores estimados. Salienta-se porém que existe, a partir dos dados amostrais, incertezas a respeito do modelo geológico que devem ser consideradas ao se projetar uma cava e desenvolver a lavra. A simulação geoestatistica tem a capacidade de produzir múltiplos modelos equiprováveis, os quais podem ser avaliados independentementecomo possíveis cenários do depósito mineral. Simulação condicional, ao contrário de técnicas de interpolação e estimativa, provê respostas para questões associadas a risco devido à variações nos teores do modelo geológico. Ao gerar múltiplos cenários tem-se acesso à probabilidade e conseqüentemente às variações de retorno financeiro e rotas de extração de minério de um projeto. o presente trabalho tem como objetivo a investigação de novas metodologias que contribuam para a construção de cenários de planejamento de lavra e avaliação do impacto provocado nestes pela incerteza fornecida a partir de modelos simulados. As respostas buscadas dentro da abordagem sugerida por esse trabalho, compreendem a definição de uma metodologia para controle e planejamento de lavra a médio e longo prazo por incorporação da flutuabilidade local associada ao minério, avaliando a sensibilidade do retorno financeiro e o traçado de lavra em relação ao método de geração do modelo geológico de teores. Para solucionar o problema em questão, sugere-se a geração de modelos estocásticos e a alimentação de múltiplos modelos selecionados por critérios considerados relevantes para a delimitação do espaço de incerteza. A aplicação de funções de transferência tais como mudança de suporte e seqüenciamento de produção sobre os modelos simulados é a chave para a obtenção de respostas a respeito do impacto sobre a lucratividade de um projeto. Ao alimentar à essas funções modelos equiprováveis que preservam as características estatísticas e a conectividade espacial dos dados amostrais mas que invariavelmente possuem diferentes distribuições de teores tem-se a dimensão em termos econômicos do risco associado à incerteza geológica. Foi confrontado o retorno financeiro produzido por modelos simulados e verificou-se para o caso especifico que os métodos de simulação geoestatistica empregados não produziram diferenças significativas quando comparados casos selecionados segundo os mesmos critérios. Porém, entre cenários extremos gerados pelo mesmo algoritmo de simulação foram verificadas desigualdades relevantes. Verificou-se também a transferência de minério entre diferentes classes de qualidade ao proceder-se com variações nas dimensões dos blocos de lavra.

Integração semi-Lagrangeana num modelo oceano-atmosfera e ajuste geostrófico no domínio tempo

Este trabalho visa o uso da função de Green de valor inicial no ajuste geostrófico e do método Semi-Lagrangeano na integração de um modelo acoplado oceano-atmosfera descrito pelas equações de águas rasas. O ajuste geostrófico é considerado atravées de perturbações na pressão e do vento. No caso de sistemas sem rotação, é discutida a relação da equação hidrostática com ondas longas não-dispersivas. Com rotação, a conservação da vorticidade potencial permite escolher a elevação correspondente a um estado de equilíbrio geostrófico. O sistema de equações de águas rasas é desacoplado em equações de Klein-Gordon com valores iniciais e termos não-homogêneos acoplados. A resposta dinâmica formada pela resposta transiente e a resposta forçada é obtida para uma perturbação inicial da elevação. A ação do vento como forçante nas equações de momento 2D, através do transporte de Eckman, conduz a uma equação de águas rasas forçada. Uma decomposição da resposta forçada é realizada com uma resposta permanente, que satisfaz a equação de Helmholtz , e com o uso da base dinâmica gerada pela resposta impulso. Um modelo hidrodinâmico 3D introduzido por Casulli e governado por equações não-lineares de águas rasas é integrado na vertical para a obtenção de um modelo 2D. Com isto, as condições de contorno devido a tensão do vento e a fricção devido a topografia do fundo, transformam-se em forçantes do modelo. O modelo foi integrado com um método semi-implícito em diferenças finitas, utilizando-se o método Semi-Lagrangeano para a parte advectiva. Simulações simbólicas foram realizadas para o ajuste geostrófico devido a perturbações de duração infinita e finita para a elevação e para o efeito da tensão do vento. Foram realizadas simulações numéricas para variadas geometrias, em particular a Baia de Guanabara e a Lagoa do Patos.

Análise estratigráfica dos sedimentos eo/mesodevonianos da porção ocidental da Bacia do Amazonas sob a ótica da estratigrafia de seqüências no interior cratônico

O trabalho foi desenvolvido com a finalidade de verificação da aplicabilidade dos conceitos da moderna Estratigrafia de Seqüências no interior cratônico e a possibilidade do refinamento cronoestratigráfico da seção eo/mesodevoniana da Bacia do Amazonas. O autor utilizou as seguintes ferramentas disponíveis para a interpretação e elaboração de um modelo geológico para a seção sedimentar estudada: a) a análise da Paleogeografia, do Paleomagnetismo e da Paleoecologia, através do estudo da Tectônica e do Clima atuantes no Eo/mesodevoniano, no Supercontinente Gondwana, retratadas em reconstituições do mundo devoniano e apoiada em extensiva consulta bibliográfica e em correlação com os conteúdos faunístico, icnológico e litológico da seção que compõe o intervalo pesquisado, com o auxílio de testemunhos, amostras de calha e afloramentos, b) o estabelecimento de superfícies-chave da Estratigrafia de Seqüências, definidas com o apoio de perfis elétrico-radioativos, notadamente o perfil de raios-gama, c) a Cicloestratigrafia química com a utilização dos teores de carbono orgânico e o índice de hidrogênio, d) a Cicloestratigrafia orbital e climática, mediante a análise espectral do perfil de raios-gama da seção estudada e a definição do controle da sedimentação, influenciada pela excentricidade curta da órbita terrestre, em ciclos dentro da banda de freqüências de Milankovitch. O encadeamento dessas análises levou o autor a montar um arcabouço cronoestratigráfico para o Eo/mesodevoniano da porção ocidental da Bacia do Amazonas. Para tanto, foram consideradas a hierarquização das unidades, a definição dos tratos de sistemas deposicionais, limites de seqüências e outras superfícies-chave estratigráficas, e a duração temporal dos eventos. Além disso, são discutidas as possíveis causas principais da evolução tectono-estratigráfica da seção estudada e sua associação com a curva de variação do nível do mar devoniano. Realizou-se ainda, tentativamente, a correlação com outras bacias intracratônicas gondwânicas através da comparação do conteúdo faunístico, da xvii paleoclimatologia, da tectônica, da análise da variação da curva do nível do mar devoniano, da posição geográfica e da relação com o Pólo Sul devoniano, que serviram de base para a compreensão do evento relacionado à passagem do Givetiano ao Frasniano, nessas bacias. Como contribuições principais o autor aponta a aplicabilidade da Estratigrafia de Seqüências no interior cratônico e a possibilidade de refinamento cronoestratigráfico em pelo menos uma ordem de grandeza (10 6 anos para 10 5 anos) com a utilização da metodologia da Cicloestratigrafia orbital.

Aperfeiçoamento da estratégia de homogeneização de minérios utilizando simulação geoestatística

Novas centrais térmicas utilizam carvão ROM (Run of mine) na sua alimentação. Sem o processamento do minério as flutuações de qualidade do carvão serão transferidas da mina para o consumidor final. Freqüentemente, contratos entre produtor-consumidor estabelecem limites pré-estabelecidos de parâmetros geológicos e tecnológicos como enxofre, cinzas e poder calorífico. Lotes de minério com qualidade fora dos limites estabelecidos podem ser rejeitados ou penalizados. O custo elevado dessas penalizações recomenda a homogeneização do minério ROM. As principais estratégias de homogeneização são as baseadas em técnicas geoestatísticas, as pilhas de homogeneização e as usinas de beneficiamento ou metalúrgicas. A estratégia de homogeneização baseada em técnicas geoestatísticas apresenta os menores custos de implementação, no entanto o conhecimento do depósito deverá ser fundamental. Tradicionalmente, essa estratégia de homogeneização utiliza um modelo de blocos gerado por um estimador tradicional, geralmente krigagem ordinária. A estimativa dos blocos por krigagem não considera a incerteza associada ao modelo geológico, não sendo adequada para prever flutuações locais Diferente da krigagem, os métodos de simulação geoestatística têm como objetivo reproduzir a variabilidade in situ e a continuidade espacial dos dados originais. Uma vez que os teores e sua variabilidade são estimados, flutuações da qualidade do minério podem ser previstas para determinada rota de lavra e certo tamanho de lote entregue ao consumidor. Diferentes tamanhos de lote são testados a fim de obter-se o controle das flutuações da qualidade do minério. Este trabalho apresenta um procedimento para acessar a variabilidade dos teores e predizer a sua flutuabilidade para diferentes tamanhos de lote em um estudo de mineração subterrânea de carvão no sul do Brasil. Simulação geoestatística e planejamento de lavra apropriado proporcionam uma solução para o problema de homogeneização de minérios.

Mapeamento da Bacia do Camaquã com a utilização de dados geofísicos, eologia e sensoriamento remoto

A análise tectono-estratigráfica da Bacia do Camaquã, uma sequência vulcanosedimentar do Neoproterozóico ao Eoproterozóico com depósitos de Cu (Au, Ag), Zn e Pb, é aqui apresentada com a utilização de ferramentas de sensoriamento remoto, gravimetria e perfilagem de poço. Nas imagens LANDSAT TM demarcou-se as concentrações de lineamentos junto as principais estruturas regionais e delimitou-se quatro domínios estruturais de acordo com a orientação dos trends dos lineamentos estruturais. Os perfis de poços que abrangem as formações Guaritas e Bom Jardim evidenciam eventos tectônicos com deformação rúptil e dúctil-ruptil, estabelecendo-se diferentes fácies tectonoestratigráficas de seqüências deposicionais (ambiente deltáico) e sequência deformacionais. Nos poços observa-se a variação da densidade com a profundidade entre poços, indicando a presença de duas aloformações de compactação distintas. Com base nos dados gravimétricos locais e regionais pode-se delimitar anomalias gravimétricas do embasamento, pacotes sedimentares de espessuras distintas em subsuperfície, com espessamento para NE, como também valores diferenciados para as principais unidades sedimentares da região estudada, bem como uma compartimentação escalonada da bacia do Camaquã. Ferramentas computacionais complementam a análise gravimétrica e de perfilagem geofísica, possibilitando a integração das técnicas já relacionadas e a formatação de dois perfis esquemáticos EW e SW-NE da bacia. Estes perfis auxiliam na visualização dos limites estruturais e formato da bacia, trazendo importantes informações para o modelo geológico da Bacia do Camaquã.

Modelo matemático de um agitador eletromagnético

O objetivo deste trabalho é estudar os efeitos eletromagnéticos e fluido-dinâmicos induzidos no aço, decorrentes do uso de um agitador eletromagnético. Para tal, foi proposta a construção de um modelo numérico que resolva, de forma acoplada, os problemas de eletromagnetismo e fluido-dinâmica. O modelo numérico do problema eletromagnético, em elementos finitos, foi construído utilizando-se o software Opera-3d/Elektra da Vector Fields. O mesmo foi validado com medidas experimentais de densidade de fluxo magnético feitas na usina. O escoamento decorrente da agitação eletromagnética foi resolvido fazendo-se o acoplamento das forças de Lorentz com as equações de Navier-Stokes. Essas últimas foram resolvidas pelo método de volumes finitos, usando-se o software CFX-4 da AEA Technology. O modelo eletromagnético mostrou que existe um torque máximo dependente da freqüência do campo magnético. Também foi observado que a força magnética aumenta em quatro vezes seu valor, quando a corrente é duplicada. O perfil de escoamento produzido no molde, sob agitação eletromagnética, indica, que as situações de lingotamento testadas, não propiciam o arraste da escória. A velocidade crítica de arraste, determinada via modelo físico, não foi atingida para nenhum caso testado. O modelo fluido-dinâmico e térmico apresentou um aumento do fluxo de calor cedido pelo fluido para a casca solidificada com o uso do agitador eletromagnético. Como conseqüência, observou-se uma queda na temperatura do banho. Também foi observado, que o uso do agitador propicia a remoção de inclusões das camadas mais externas do tarugo. Ao mesmo tempo, notou-se que o uso do agitador aumenta o índice de remoção de inclusões para as duas seções de molde analisadas.

Respostas dinâmicas em sistemas distribuídos e decomposição forçada da superfície livre para um modelo acoplado Oceano-Atmosfera

O objetivo deste trabalho é a introdução e desenvolvimento de uma metodologia analítico-simbólica para a obtenção de respostas dinâmicas e forçadas (soluções homogêneas e não homogêneas) de sistemas distribuídos, em domínios ilimitados e limitados, através do uso da base dinâmica gerada a partir da resposta impulso. Em domínios limitados, a resposta impulso foi formulada pelo método espectral. Foram considerados sistemas com condições de contorno homogêneas e não homogêneas. Para sistemas de natureza estável, a resposta forçada é decomposta na soma de uma resposta particular e de uma resposta livre induzida pelos valores iniciais da resposta particular. As respostas particulares, para entradas oscilatórias no tempo, foram calculadas com o uso da fun»c~ao de Green espacial. A teoria é desenvolvida de maneira geral permitindo que diferentes sis- temas evolutivos de ordem arbitrária possam ser tratados sistematicamente de uma forma compacta e simples. Realizou-se simulações simbólicas para a obtenção de respostas dinâmicas e respostas for»cadas com equações do tipo parabólico e hiperbólico em 1D,2D e 3D. O cálculo das respostas forçadas foi realizado com a determinação das respostas livres transientes em termos dos valores iniciais das respostas permanentes. Foi simulada a decomposição da resposta forçada da superfície livre de um modelo acoplado oceano-atmosfera bidimensional, através da resolução de uma equação de Klein-Gordon 2D com termo não-homogêneo de natureza dinâmica, devido a tensão de cisalhamento na superfície do oceano pela ação do vento.

Modelagem estrutural e tridimensional para a prospecção e avaliação dos depósitos de ágata do distrito mineiro de Salto do Jacuí (RS)

O Distrito Mineiro de Salto do Jacuí (DMSJ) abrange uma área aproximada de 250 km2 localizada no Município de Salto do Jacuí na região central do estado do Rio Grande do Sul (Brasil). A região é caracterizada pela ocorrência de derrames basálticos a dacíticos pertencentes a Fm. Serra Geral, uma das unidades de topo da Bacia do Paraná. O DMSJ é um dos maiores produtores mundiais de geodos de ágata que ocorrem em derrames basálticos. Em escala de mina, os depósitos de geodos de ágata do DMSJ ocorrem em uma seqüência vulcânica constituída por 3 unidades líticas principais: 1) o dacito semi-vítreo vesículoamigdaloidal inferior (DSVI), 2) o basalto vesículo-amigdaloidal mineralizado (BM, “tabatinga”), e 3) o dacito vesicular (DV, “cupim”). Essas unidades vulcânicas constituem a denominada Estrutura Jacuí. Essa tese propõe um modelo estrutural para o controle da distribuição dos geodos de ágata no DMSJ e um novo conjunto de procedimentos aplicados para a exploração dos geodos de ágata. A modelagem estrutural 3D e a distribuição espacial dos depósitos e dos tipos de geodos de ágata são avaliados para permitir futuras investigações detalhadas a respeito da estimativa de reservas remanescentes do DMSJ. Mapeamento geológico em escala de mina, análise estrutural e investigação da distribuição espacial dos geodos de ágata foram realizados para desenvolver um modelo 3D do controle estrutural dos geodos de ágata e definir parâmetros de prospecção de geodos de ágata. A modelagem 3D e o controle estrutural da distribuição dos geodos de ágata permitiram definir alguns critérios para investigações quanto à estimativa de recursos/reservas no DMSJ. Nessa linha de investigação, devem-se considerar principalmente os seguintes parâmetros: i) os novos procedimentos para exploração de geodos de ágata, ii) a diluição introduzida pela presença de diques de arenito Botucatu que cortam o Basalto Mineralizado, e iii) o controle estrutural proposto parta os depósitos de ágata do DMSJ. Esses parâmetros são úteis nas atividades de exploração e na avaliação de novas áreas de geodos de ágata. O controle estrutural e a modelagem 3D da distribuição dos geodos de ágata também permitem definir os limites das unidades vulcânicas mineralizadas e a extensão do corpo mineralizado.