Critérios para locação amostral baseados em simulação estocástica

Em planejamento de lavra, freqüentemente as campanhas de sondagem são complementadas por falta de segurança nas estimativas. Na maioria dos casos, essas sondagens adicionais são locadas sem o uso de qualquer procedimento matemático que indique os locais mais carentes de informação. Essa dissertação propõe uma metodologia fundamentada em simulação estocástica para tratar esse problema. Com a utilização desse método, é possível que seja considerada a variabilidade entre as informações, além da configuração espacial apresentada. Essa metodologia permite gerar várias estimativas equiprováveis do atributo em questão, reproduzindo o histograma e o variograma das informações originais. A variabilidade apresentada pelos conjuntos de n valores equiprováveis, em cada bloco discretizado da jazida, é medida por meio do cálculo de índices de incerteza. Esses índices são calculados após o processo de simulação estocástica por meio da compilação das n estimativas realizadas para cada bloco. Os índices de incerteza considerados nessa dissertação são: (i) variância condicional; (ii) coeficiente de variação condicional e (iii) intervalo entre quartis. A partir da classificação desses índices de incerteza, são identificados os setores que apresentam as maiores variabilidades em suas estimativas. Com essa classificação é possível a confecção de mapas de incerteza que auxiliam na visualização dos setores/blocos mais variáveis da jazida Nesses setores, são adicionadas informações oriundas da escolha aleatória de um dos n valores equiprováveis disponíveis para cada bloco simulado. Por meio dos índices de incerteza e de acordo com a metodologia, são propostas informações adicionais na área considerada, que contribuem na diminuição das incertezas globais e locais. Essa metodologia usa técnicas extremamente robustas, com embasamento teórico muito consistente, atacando diretamente os setores/blocos que apresentam as maiores flutuabilidades nas estimativas. Um estudo de caso em uma jazida de carvão verifica a aplicabilidade do método. Esse estudo, considera a variabilidade nas estimativas de espessura de uma camada carvão.

Modelos dinâmicos e simulação estocástica

This paper presents new methodology for making Bayesian inference about dy~ o!s for exponential famiIy observations. The approach is simulation-based _~t> use of ~vlarkov chain Monte Carlo techniques. A yletropolis-Hastings i:U~UnLlllll 1::; combined with the Gibbs sampler in repeated use of an adjusted version of normal dynamic linear models. Different alternative schemes are derived and compared. The approach is fully Bayesian in obtaining posterior samples for state parameters and unknown hyperparameters. Illustrations to real data sets with sparse counts and missing values are presented. Extensions to accommodate for general distributions for observations and disturbances. intervention. non-linear models and rnultivariate time series are outlined.

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)

Simulação estocástica da transmissividade em formações graníticas alteradas e fracturadas

Considerando como caso de estudo o do maciço envolvente e subjacente à antiga mina de urânio de Quinta do Bispo (Concelho de Mangualde), composto essencialmente por granitos hercínicos e metasedimentos do Complexo Xisto-Grauváquico, estabeleceu-se uma metodologia que permitiu a construção de modelos 3D de transmissividade, equiprováveis, mediante o cruzamento de propriedades intrínsecas do maciço rochoso – neste caso particular, a litologia, o grau de alteração e a densidade de fracturação - passíveis de modelação estocástica a 3D, com propriedades hidrogeológicas do meio, medidas in situ e avaliadas deterministicamente, por via da realização e interpretação de ensaios de bombagem direcionados. A interpretação dos resultados dos ensaios de bombagem realizou-se por patamares, recorrendo-se aos métodos de “Porosidade Dupla” e de “Theis com correcção de Jacob”, tendo-se processado curvas cumulativas de valores de transmissividade em função da litologia e do grau de alteração, conforme os distintos sectores comportamentais do maciço. A metodologia adoptada levou a que os valores de transmissividade de cada célula dos modelos 3D tenham sido simulados com condicionamento aos valores experimentais dos ensaios de bombagem, ponderados de acordo com as probabilidades dessa mesma célula representar, numa dada localização no espaço, uma certa Litologia, sob um determinado Grau de Alteração que possui, quando aplicável, um certo Número de Fracturas não preenchidas. As variáveis representativas dos atributos geológicos Litologia, Grau de Alteração e Número de Fracturas, foram simuladas em cadeia pelo que os modelos de transmissividades integram a variabilidade e heterogeneidade locais destes atributos, os quais condicionam o fluxo tridimensional da água.

Simulação de estruturas meandriformes por objectos e estatísticas multiponto e avaliação da porosidade – aplicação a reservatórios siliciclásticos

Neste trabalho apresenta-se uma metodologia destinada a simular a morfologia e a porosidade de canais de areia em reservatórios fluviais. Como informação de partida utilizam-se um ou vários canais de areia de treino, em estrutura vectorial, que sejam representativos do reservatório, e leis de distribuição das dimensões largura e profundidade dos canais. Para a simulação da morfologia dos canais de areia, são calculados os ângulos azimutais de todos os segmentos de recta das linhas poligonais que constituem os canais de treino, a que se segue a codificação em classes de azimute e a determinação de estatísticas multiponto destas classes. Seguidamente, faz-se a simulação estocástica dos novos canais no volume do reservatório, ou seja, são geradas novas linhas poligonais condicionadas às estatísticas multiponto, a que se associam as dimensões largura e profundidade. Este algoritmo multiponto é muito eficiente, porque é a 1D, e contorna a problemática da variável azimute ser do tipo circular. A avaliação da porosidade é feita em simultâneo na conversão do modelo vectorial da morfologia dos canais de areia para o modelo matricial da malha de blocos do reservatório. Para cada bloco reservatório é contabilizada a proporção de não canal / canal e das fácies conforme um modelo conceptual de zonamento na secção do canal. A estimativa da porosidade média de cada bloco do reservatório é obtida pela soma do produto da proporção de cada fácies pela respectiva porosidade de referência. Para ilustrar as potencialidades da metodologia proposta apresenta-se um caso de estudo com dados sintéticos. De acordo com critérios geológicos e estratigráficos, consideraram-se quatro fácies na secção do canal que são discriminadas lateralmente e em profundidade por valores de referência de porosidade. O algoritmo de simulação morfológica capturou adequadamente a morfologia das imagens de treino e gerou várias imagens equiprováveis que depois foram convertidas em porosidade.

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.

Aumento da eficiência dos métodos següenciais de simulação condicional

O algoritmo de simulação seqüencial estocástica mais amplamente utilizado é o de simulação seqüencial Gaussiana (ssG). Teoricamente, os métodos estocásticos reproduzem tão bem o espaço de incerteza da VA Z(u) quanto maior for o número L de realizações executadas. Entretanto, às vezes, L precisa ser tão alto que o uso dessa técnica pode se tornar proibitivo. Essa Tese apresenta uma estratégia mais eficiente a ser adotada. O algoritmo de simulação seqüencial Gaussiana foi alterado para se obter um aumento em sua eficiência. A substituição do método de Monte Carlo pela técnica de Latin Hypercube Sampling (LHS), fez com que a caracterização do espaço de incerteza da VA Z(u), para uma dada precisão, fosse alcançado mais rapidamente. A técnica proposta também garante que todo o modelo de incerteza teórico seja amostrado, sobretudo em seus trechos extremos.

Simulação e análise econômica da produção in vivo e in vitro de embriões em bovinos

The objective of this work was to evaluate the effect of the variables number of recipients, synchronization protocol, reproductive efficiency indicators and pregnancy cost, in the economic effectiveness of in vivo and in vitro bovine embryo production. A simulation application was elaborated to allow the user to insert the input variable parameters. A basic scenario, from the efficiency traditional rates of in vivo (ET) and in vitro production (IVP) techniques of bovine embryos, was introduced in the software as a criterion to compare the results. This software was able to reproduce both ET and IVP scenarios. The embryo production was simulated through stochastic simulation. The optimal number of recipients using sensitivity analysis was determined. The net present value and cost per pregnancy were used as a decision parameter. The synchronization for fixed-time embryo transfer decreased the recipient idleness and, consequently, the final cost of pregnancy, in comparison to the traditional methodology. Foetal sexing must be associated to IVP of bovine embryos. In addition, the optimal recipient number per donor is variable and depends on data inserted in the system.

Uso combinado de sêmen sexado e acasalamento dirigido sobre uma população de bovinos de corte submetida a seleção: estudo de simulação

Desenvolveu-se um estudo de simulação estocástica com o objetivo de verificar as consequências do uso combinado de acasalamento dirigido e sêmen sexado em uma população de bovinos de corte sob seleção. Simularam-se seis gerações de seleção para três cenários de acasalamento e uso de sêmen sexado. O primeiro cenário foi caracterizado por acasalamento aleatório e uso exclusivo de sêmen convencional. O segundo cenário caracterizou-se pelo uso de acasalamento associativo positivo nas 40% melhores vacas e acasalamento associativo negativo nas demais, sem uso de sêmen sexado. O terceiro cenário seguiu o mesmo procedimento de acasalamento do segundo, combinando-o com uso de sêmen sexado nas vacas submetidas a acasalamento associativo positivo. O acasalamento associativo positivo teve maior impacto no progresso genético que o uso de sêmen sexado, apesar de ter aumentado a incidência de endogamia na população. O uso de acasalamento associativo negativo foi ineficiente em reduzir a variabilidade dos animais destinados ao abate. O uso combinado de acasalmento associativo positivo e sêmen sexado aumentou a produção de animais geneticamente superiores.

Aplicação de métodos estocásticos na gestão de custos de construção de condutas de saneamento

O trabalho apresentado centra-se na determinação dos custos de construção de condutas de pequenos e médios diâmetros em Polietileno de Alta Densidade (PEAD) para saneamento básico, tendo como base a metodologia descrita no livro Custos de Construção e Exploração – Volume 9 da série Gestão de Sistemas de Saneamento Básico, de Lencastre et al. (1994). Esta metodologia descrita no livro já referenciado, nos procedimentos de gestão de obra, e para tal foram estimados custos unitários de diversos conjuntos de trabalhos. Conforme Lencastre et al (1994), “esses conjuntos são referentes a movimentos de terras, tubagens, acessórios e respetivos órgãos de manobra, pavimentações e estaleiro, estando englobado na parte do estaleiro trabalhos acessórios correspondentes à obra.” Os custos foram obtidos analisando vários orçamentos de obras de saneamento, resultantes de concursos públicos de empreitadas recentemente realizados. Com vista a tornar a utilização desta metodologia numa ferramenta eficaz, foram organizadas folhas de cálculo que possibilitam obter estimativas realistas dos custos de execução de determinada obra em fases anteriores ao desenvolvimento do projeto, designadamente numa fase de preparação do plano diretor de um sistema ou numa fase de elaboração de estudos de viabilidade económico-financeiros, isto é, mesmo antes de existir qualquer pré-dimensionamento dos elementos do sistema. Outra técnica implementada para avaliar os dados de entrada foi a “Análise Robusta de Dados”, Pestana (1992). Esta metodologia permitiu analisar os dados mais detalhadamente antes de se formularem hipóteses para desenvolverem a análise de risco. A ideia principal é o exame bastante flexível dos dados, frequentemente antes mesmo de os comparar a um modelo probabilístico. Assim, e para um largo conjunto de dados, esta técnica possibilitou analisar a disparidade dos valores encontrados para os diversos trabalhos referenciados anteriormente. Com os dados recolhidos, e após o seu tratamento, passou-se à aplicação de uma metodologia de Análise de Risco, através da Simulação de Monte Carlo. Esta análise de risco é feita com recurso a uma ferramenta informática da Palisade, o @Risk, disponível no Departamento de Engenharia Civil. Esta técnica de análise quantitativa de risco permite traduzir a incerteza dos dados de entrada, representada através de distribuições probabilísticas que o software disponibiliza. Assim, para por em prática esta metodologia, recorreu-se às folhas de cálculo que foram realizadas seguindo a abordagem proposta em Lencastre et al (1994). A elaboração e a análise dessas estimativas poderão conduzir à tomada de decisões sobre a viabilidade da ou das obras a realizar, nomeadamente no que diz respeito aos aspetos económicos, permitindo uma análise de decisão fundamentada quanto à realização dos investimentos.

Abordagem probabilística na avaliação de depósitos de argila cerâmica

A argila é uma matéria-prima básica nos processos de fabricação de pisos e azulejos. Especificações industriais devem ser cumpridas para sua utilização como componente das formulações cerâmicas. A homogeneidade é vital. Nos depósitos minerais fonte desse material, as características cerâmicas da argila, sua representatividade, variabilidade e interação com a geologia local precisam ser bem conhecidas, a fim de propriamente quantificá-la e classificá-la de forma georeferenciada. Isso consubstancia a tomada de decisões mineiras relacionadas. Essa dissertação apresenta uma alternativa metodológica de avaliação de recursos geológicos e reservas recuperáveis de argila cerâmica, em uma perspectiva probabilística. Métodos geoestatísticos como krigagem e simulação estocástica são aplicados em suporte de blocos de lavra. Recursos geológicos ocorrentes nos diversos setores de um depósito estudado são estimados acuradamente. Em sua porção mais favorável reservas recuperáveis são avaliadas, sob critérios cerâmicos e de seletividade mineira, tomados separadamente e em conjunto, em escala global e local Para o estudo de caso apresentado, duas propriedades cerâmicas pós-queima são consideradas, absorção d’água e retração linear, além da espessura do horizonte argiloso. Modelos de cenários extremos e medianos de reservas são construídos, flexibilizando a apreciação do quadro geral obtido, e permitindo assim maior segurança na condução de decisões e eventuais operações mineiras. Comparativos com métodos tradicionalmente utilizados na avaliação deste tipo de depósitos, em suas diversas fases, comprovam as vantagens da metodologia proposta, a qual permite, além da quantificação e classificação cerâmica de reservas, uma incorporação de margem de riscos nas predições realizadas.

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.

Eutrofização e capacidade de carga de fósforo de seis reservatórios a bacia do Rio Seridó, região semi-árida do Estado do RN

The eutrofization is a natural process of accumulation of nutrients in aquatic´s body that it has been accelerated for the human´s actives, mainly the related with the activities of camp, industrial and the inadequate disposition of the domestic sewage. The enrichment of the aquatic´s body with nutrients, mainly the nitrogen and the phosphorus, and the consequent proliferation of algae and Cyanobacteria can commit the quality of the water for the public provisioning, for the fish farming and for other ends. The problem becomes more critical when there is a shortage of water naturally as in the semi-arid area of the Brazilian northeast. Before that problem this work had as objective evaluates the trophic state of six reservoirs of the basin of River Seridó of Rio Grande of Norte and also estimate the capacity of load of match of the reservoir and risk probabilities based on the established limits by the resolution Conama 357/05. The results demonstrate that the six reservoirs are eutrofization, with concentration of total phosphorus and cloro a in the water upster to 50 e 12 μg l-1. The results show that space homogeneity exists in the state trophic of the reservoirs, but a significant variation interanual in function of the increase of the concentrations of nutrients and decrease of the transparency of the water with the reduction of the body of water accumulated in the reservoirs.The results of the simulation risk estocastic show that the reservoirs could receive annually from 72 to 216 Kg of P, assuming a risk of 10% of increasing in more than 30 μg l-1 the annual medium concentrations of total match in the water of these reservoirs. This load could be high in until 360 kg of P a year in case the managers assume a risk of 10% of increasing in more than 50 μg l-1 the annual medium concentrations of total phosphorus in the waters of these reservoirs

Caracterização estatística de processos sísmicos via Distribuição Generalizada de Pareto. Estudo de caso: João Câmara-RN

The work is to make a brief discussion of methods to estimate the parameters of the Generalized Pareto distribution (GPD). Being addressed the following techniques: Moments (moments), Maximum Likelihood (MLE), Biased Probability Weighted Moments (PWMB), Unbiased Probability Weighted Moments (PWMU), Mean Power Density Divergence (MDPD), Median (MED), Pickands (PICKANDS), Maximum Penalized Likelihood (MPLE), Maximum Goodness-of-fit (MGF) and the Maximum Entropy (POME) technique, the focus of this manuscript. By way of illustration adjustments were made for the Generalized Pareto distribution, for a sequence of earthquakes intraplacas which occurred in the city of João Câmara in the northeastern region of Brazil, which was monitored continuously for two years (1987 and 1988). It was found that the MLE and POME were the most efficient methods, giving them basically mean squared errors. Based on the threshold of 1.5 degrees was estimated the seismic risk for the city, and estimated the level of return to earthquakes of intensity 1.5°, 2.0°, 2.5°, 3.0° and the most intense earthquake never registered in the city, which occurred in November 1986 with magnitude of about 5.2º

Proposta metodológica para o imageamento, caracterização, parametrização e geração de modelos virtuais de afloramentos

The aim of this work was to describe the methodological procedures that were mandatory to develop a 3D digital imaging of the external and internal geometry of the analogue outcrops from reservoirs and to build a Virtual Outcrop Model (VOM). The imaging process of the external geometry was acquired by using the Laser Scanner, the Geodesic GPS and the Total Station procedures. On the other hand, the imaging of the internal geometry was evaluated by GPR (Ground Penetrating Radar).The produced VOMs were adapted with much more detailed data with addition of the geological data and the gamma ray and permeability profiles. As a model for the use of the methodological procedures used on this work, the adapted VOM, two outcrops, located at the east part of the Parnaiba Basin, were selected. On the first one, rocks from the aeolian deposit of the Piaui Formation (Neo-carboniferous) and tidal flat deposits from the Pedra de Fogo Formation (Permian), which arises in a large outcrops located between Floriano and Teresina (Piauí), are present. The second area, located at the National Park of Sete Cidades, also at the Piauí, presents rocks from the Cabeças Formation deposited in fluvial-deltaic systems during the Late Devonian. From the data of the adapted VOMs it was possible to identify lines, surfaces and 3D geometry, and therefore, quantify the geometry of interest. Among the found parameterization values, a table containing the thickness and width, obtained in canal and lobes deposits at the outcrop Paredão and Biblioteca were the more relevant ones. In fact, this table can be used as an input for stochastic simulation of reservoirs. An example of the direct use of such table and their predicted radargrams was the identification of the bounding surface at the aeolian sites from the Piauí Formation. In spite of such radargrams supply only bi-dimensional data, the acquired lines followed of a mesh profile were used to add a third dimension to the imaging of the internal geometry. This phenomenon appears to be valid for all studied outcrops. As a conclusion, the tool here presented can became a new methodology in which the advantages of the digital imaging acquired from the Laser Scanner (precision, accuracy and speed of acquisition) were combined with the Total Station procedure (precision) using the classical digital photomosaic technique