Métodos estatísticos recursivos aplicados ao problema de estimação de vazão de gás em plantas de Plunger-lift

This work has as main objective to find mathematical models based on linear parametric estimation techniques applied to the problem of calculating the grow of gas in oil wells. In particular we focus on achieving grow models applied to the case of wells that produce by plunger-lift technique on oil rigs, in which case, there are high peaks in the grow values that hinder their direct measurement by instruments. For this, we have developed estimators based on recursive least squares and make an analysis of statistical measures such as autocorrelation, cross-correlation, variogram and the cumulative periodogram, which are calculated recursively as data are obtained in real time from the plant in operation; the values obtained for these measures tell us how accurate the used model is and how it can be changed to better fit the measured values. The models have been tested in a pilot plant which emulates the process gas production in oil wells

Previsão de inflação utilizando modelos de séries temporais

Este trabalho compara modelos de séries temporais para a projeção de curto prazo da inflação brasileira, medida pelo Índice de Preços ao Consumidor Amplo (IPCA). Foram considerados modelos SARIMA de Box e Jenkins e modelos estruturais em espaço de estados, estimados pelo filtro de Kalman. Para a estimação dos modelos, foi utilizada a série do IPCA na base mensal, de março de 2003 a março de 2012. Os modelos SARIMA foram estimados no EVIEWS e os modelos estruturais no STAMP. Para a validação dos modelos para fora da amostra, foram consideradas as previsões 1 passo à frente para o período de abril de 2012 a março de 2013, tomando como base os principais critérios de avaliação de capacidade preditiva propostos na literatura. A conclusão do trabalho é que, embora o modelo estrutural permita, decompor a série em componentes com interpretação direta e estudá-las separadamente, além de incorporar variáveis explicativas de forma simples, o desempenho do modelo SARIMA para prever a inflação brasileira foi superior, no período e horizonte considerados. Outro importante aspecto positivo é que a implementação de um modelo SARIMA é imediata, e previsões a partir dele são obtidas de forma simples e direta.

Impondo mais restrições ao modelo de apreçamento vetorial com séries temporais

Este Trabalho se Dedica ao exercício empírico de gerar mais restrições ao modelo de apreçamento de ativos com séries temporais desenvolvido por Hansen e Singleton JPE 1983. As restrições vão, desde um simples aumento qualitativo nos ativos estudados até uma extensão teórica proposta a partir de um estimador consistente do fator estocástico de desconto. As estimativas encontradas para a aversão relativa ao risco do agente representativo estão dentro do esperado, na maioria dos casos, já que atingem valores já encontrados na literatura além do fato destes valores serem economicamente plausíveis. A extensão teórica proposta não atingiu resultados esperados, parecendo melhorar a estimação do sistema marginalmente.

Análise de séries temporais de coordenadas estimadas com GPS: uma proposta metodológica para eliminação de efeitos sazonais

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

Estimativa do expoente de Hurst, por meio da transformada Wavelet, de séries temporais de precipitação de chuvas das regiões climáticas do estado de São Paulo no período de 1978 a 1997

Pós-graduação em Física - IGCE

Aplicação de análise de séries temporais em perfis de poço

Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções.

Análise da evolução de software com séries temporais

Trabalho apresentado no âmbito do Mestrado em Engenharia Informática, como requisito parcial para obtenção do grau de Mestre em Engenharia Informática

Detecção de periodicidades em séries temporais: o Periodograma

Neste texto, depois de uma breve introdução à análise de Fourier, com especial incidência na apresentação de Fourier de sucessões periódicas, definimos o periodograma e deduzimos algumas das suas propriedades estatísticas. Mostramos como este pode ser utilizado para a detecção de periodicidades numa série temporal, e damos alguns testes nas ordenadas do periodograma. Por fim, aplicamos a teoria desenvolvida a um conjunto de observações referentes à temperatura da água do mar.

Detecção de correlação e causalidade em séries temporais não categóricas

Dissertação de Mestrado em Engenharia Informática

A metodologia Bootstrap associada ao método de Holt-Winters na previsão de séries temporais

Os métodos de alisamento exponencial são muito utilizados na modelação e previsão de séries temporais, devido à sua versatilidade e opção de modelos que integram. Na estatística computacional, a metodologia Bootstrap é muito aplicada em inferência estatística no âmbito de séries temporais. Este estudo teve como principal objectivo analisar o desempenho do método de Holt-Winters associado à metodologia Bootstrap, como um processo alternativo na modelação e previsão de séries temporais.

Exploração de séries temporais em processos de previsão de vendas

Dissertação de mestrado em Engenharia de Sistemas

Análise de séries temporais da mortalidade por doenças isquêmicas do coração e cerebrovasculares, nas cinco regiões do Brasil, no período de 1981 a 2001

OBJETIVO: Analisar as tendências do risco de morte por doenças isquêmicas do coração e cerebrovasculares, nas regiões Norte, Nordeste, Centro-Oeste, Sudeste e Sul do Brasil, no período de 1981 a 2001. MÉTODOS: Dados de mortalidade por doenças isquêmicas do coração e cerebrovasculares nas cinco regiões brasileiras foram obtidos através do Ministério da Saúde. A fonte de dados foi o Sistema de Informações sobre Mortalidade (SIM), do Departamento de Análise da Informação de Saúde - Secretaria de Vigilância em Saúde do Ministério da Saúde. A população das regiões e estados teve por fonte o IBGE, censos 1991 e 2000, contagem populacional de 1996 e estimativas populacionais, disponibilizado pelo Datasus. Os dados do SIM foram relativos às seguintes causas de morte: doença cerebrovascular (Código Internacional de Doenças CID-9 430-438, CID-10 I60-I69) e a doença isquêmica do coração (CID-9 410-414, CID-10 I21-I25). O estudo estatístico utilizou-se para as análises inferenciais de modelos lineares generalizados ajustados. RESULTADOS: A tendência da mortalidade por doença cerebrovascular mostrou declínio nas regiões Sudeste, Sul e Centro-Oeste em todas as faixas etárias e sexo. Também a mortalidade por doença isquêmica do coração declinou nas regiões Sudeste e Sul, com estabilização do risco na região Centro-Oeste e aumento na região Nordeste. CONCLUSÃO: O risco de morte para as doenças circulatórias, cerebrovasculares e isquêmicas do coração diminuiu no Sul e no Sudeste, regiões mais desenvolvidas do país, e aumentou nas menos desenvolvidas, principalmente no Nordeste.

Séries temporais de NDVI do sensor SPOT Vegetation e algoritmo SAM aplicados ao mapeamento de cana‑de‑açúcar

O objetivo deste trabalho foi avaliar o mapeamento de área de cana‑de‑açúcar por meio de série temporal, de seis anos de dados do índice de vegetação por diferença normalizada (NDVI), oriundos do sensor Vegetation, a bordo do satélite "système pour l'observation de la Terre" (SPOT). Três classes de cobertura do solo (cana‑de‑açúcar, pasto e floresta), do Estado de São Paulo, foram selecionadas como assinaturas espectro‑temporais de referência, que serviram como membros extremos ("endmembers") para classificação com o algoritmo "spectral angle mapper" (SAM). A partir desta classificação, o mapeamento da área de cana‑de‑açúcar foi realizado com uso de limiares na imagem-regra do SAM, gerados a partir dos valores dos espectros de referência. Os resultados mostram que o algoritmo SAM pode ser aplicado a séries de dados multitemporais de resolução moderada, o que permite eficiente mapeamento de alvo agrícola em escala mesorregional. Dados oficiais de áreas de cana‑de‑açúcar, para as microrregiões paulistas, apresentam boa correlação (r² = 0,8) com os dados obtidos pelo método avaliado. A aplicação do algoritmo SAM mostrou ser útil em análises temporais. As séries temporais de NDVI do sensor SPOT Vegetation podem ser utilizadas para mapeamento da área de cana‑de‑açúcar em baixa resolução.

Tratamento de ruídos e caracterização de fisionomias do Cerrado utilizando séries temporais do sensor MODIS

O Cerrado é formado por um mosaico de fisionomias campestres, savânicas e florestais que possuem um típico ciclo fenológico. Nesse contexto, os dados do MODIS fornecem medidas diárias que permitem monitorar a sazonal fenologia da vegetação. O objetivo deste trabalho foi caracterizar formações savânicas, formações florestais e áreas de cerrado convertido pela ação antrópica, utilizando séries temporais de NDVI e EVI do sensor MODIS, após a suavização de ruídos. A metodologia adotada pode ser subdividida nos seguintes passos: (a) confecção do cubo temporal com NDVI e EVI, onde o perfil em z corresponde à assinatura temporal, (b) tratamento do ruído e (c) detecção das assinaturas temporais. A Transformação pela Fração Mínima de Ruído (MNF) foi aplicada para suavizar ruídos contidos nas séries temporais. Os resultados indicaram que os valores de NDVI foram maiores que os valores de EVI e existiu relação com as estações do ano. As formações florestais foram as que apresentaram maiores valores de NDVI e EVI, possuindo ainda as menores variações entre as estações. As áreas convertidas apresentaram os menores valores dos dois índices, tendo grande queda em seus valores no início da estação seca, provavelmente por se tratar do período de colheita. O estudo concluiu que as séries temporais NDVI e EVI são úteis na diferenciação dos tipos de vegetação.

Uso da função de transferência em problemas de condução do calor com a lei de Fourier modificada

Este trabalho visa o uso da função de transferência, a qual relaciona distribuição de temperatura e fluxo de calor, na comparação no domínio freqüência, entre o modelo de difusão usual (parabólico) e um modelo ondulatório (hiperbólico) que inclue o efeito de propagação do calor não instantâneo, sendo avaliados os casos de meio semi-infinito e finito. Para o caso de meio semi-infinito, são determinadas as expressões para as características de amplitude e de fase, considerando tanto a abordagem parabólica quanto hiperbólica. É observada a relação entre estas duas abordagens, mostrando que a abordagem parabólica é uma caso particular da abordagem hiperbólica, podendo ser obtida através de um processo de limite envolvendo o tempo de relaxação r . Para o caso de meio finito, são determinadas as expressões para as caracteríısticas de amplitude de ambas as faces da placa unidimensional, considerando tanto a abordagem parabólica quanto hiperbólica. Estas expressões são transformadas para a forma adimensional quando então são deduzidas as expressões correspondentes das características de amplitude. Mais uma vez, todos os resultados para o caso parabólico podem ser determinados a partir dos resultados do caso hiperbólico, através de um processo de limite envolvendo o tempo de relaxação r São apresentados resultados numéricos referentes às características de amplitude, onde é apontada a existência de uma freqüência limite, acima da qual a diferença entre os dois modelos, do tipo parabólico ou hiperbólico, aumenta rapidamente. Também é apresentada uma forma alternativa de cálculo da distribuição de temperatura transiente que faz uso da função de transferência do sistema.