Deconvolução de perfis de poços através do ajuste de energia

As medidas de resistividade são de fundamental importância para o cálculo da saturação de óleo em reservatórios potencialmente produtores. A combinação das medidas de resistividade rasa e profunda permite a obtenção dos parâmetros R_t, R_xo e d_i. Mas, em reservatórios complexos existem dificuldades em se obter leituras confiáveis de R_t, devido à baixa resolução vertical das ferramentas de investigação profunda. Em reservatórios laminados, por exemplo, as leituras obtidas pela ferramenta de indução profunda (ILD) podem levar a uma interpretação errônea das mesmas, levando a acreditar que as medidas obtidas do perfil referem-se a uma única camada. Este problema pode ser em parte resolvido através de uma metodologia que melhore a resolução vertical dos perfis de investigação profunda, valendo-se do uso de informações obtidas de um perfil de alta resolução vertical, i.e; a curva de resistividade rasa. Uma abordagem neste sentido seria usar um perfil de alta resolução que apresente uma boa correlação com o perfil de investigação profunda. Esta correlação pode ser melhor avaliada se aplicarmos um filtro no perfil de alta resolução, de tal forma que o perfil resultante tenha teoricamente a mesma resolução vertical do perfil de baixa resolução. A obtenção deste filtro, porém, recai na premissa de que as funções respostas verticais para as ferramentas de alta e baixa resolução são disponíveis, o que não ocorre na prática. Este trabalho se propõe mostrar uma nova abordagem onde o filtro pode ser obtido a partir de um tratamento no domínio da freqüência. Tal tratamento visa igualar a energia espectral do perfil de alta resolução à energia espectral do perfil de baixa resolução tendo como base o Teorema de Parseval. Será mostrado que a resolução vertical depende fundamentalmente da energia espectral do perfil em questão. A seguir, uma regressão linear é aplicada sobre os perfis de alta resolução filtrado e de baixa resolução. Para cada ponto amostrado dos perfis, uma rotina de minimização é aplicada visando escolher o melhor intervalo de correlação entre os perfis. Finalmente, um fator de correção é aplicado sobre cada ponto do perfil de baixa resolução. Os resultados obtidos com os perfis de indução são promissores, demonstrando a eficácia da abordagem e mesmo quando aplicada em perfis com diferentes propriedades petrofísicas, a metodologia funcionou satisfatoriamente, sem danificar os perfis originais.

Detecção de refletores sísmicos por rede neural discreta

As redes neurais artificiais têm provado serem uma poderosa técnica na resolução de uma grande variedade de problemas de otimização. Nesta dissertação é desenvolvida uma nova rede neural, tipo recorrente, sem realimentação (self-feedback loops) e sem neurônios ocultos, para o processamento do sinal sísmico, para fornecer a posição temporal, a polaridade e as amplitudes estimadas dos refletores sísmicos, representadas pelos seus coeficientes de reflexão. A principal característica dessa nova rede neural consiste no tipo de função de ativação utilizada, a qual permite três possíveis estados para o neurônio. Busca-se estimar a posição dos refletores sísmicos e reproduzir as verdadeiras polaridades desses refletores. A idéia básica desse novo tipo de rede, aqui denominada rede neural discreta (RND), é relacionar uma função objeto, que descreve o problema geofísico, com a função de Liapunov, que descreve a dinâmica da rede neural. Deste modo, a dinâmica da rede leva a uma minimização local da sua função de Liapunov e consequentemente leva a uma minimização da função objeto. Assim, com uma codificação conveniente do sinal de saída da rede tem-se uma solução do problema geofísico. A avaliação operacional da arquitetura desta rede neural artificial é realizada em dados sintéticos gerados através do modelo convolucional simples e da teoria do raio. A razão é para explicar o comportamento da rede com dados contaminados por ruído, e diante de pulsos fonte de fases mínima, máxima e misturada.

Espalhamento geométrico em modelos plano-estratificados

A medição de parâmetros físicos de reservatórios se constitui de grande importância para a detecção de hidrocarbonetos. A obtenção destes parâmetros é realizado através de análise de amplitude com a determinação dos coeficientes de reflexão. Para isto, faz-se necessário a aplicação de técnicas especiais de processamento capazes de corrigir efeitos de divergência esférica. Um problema pode ser estabelecido através da seguinte questão: Qual o efeito relativamente mais importante como responsável pela atenuação de amplitudes, o espalhamento geométrico ou a perda por transmissividade? A justificativa desta pergunta reside em que a correção dinâmica teórica aplicada a dados reais visa exclusivamente o espalhamento geométrico. No entanto, a análise física do problema por diferentes direções põe a resposta em condições de dúvida, o que é interessante e contraditório com a prática. Uma resposta embasada mais fisicamente pode dar melhor subsídio a outros trabalhos em andamento. O presente trabalho visa o cálculo da divergência esférica segundo a teoria Newman-Gutenberg e corrigir sismogramas sintéticos calculados pelo método da refletividade. O modelo-teste é crostal para que se possa ter eventos de refração crítica além das reflexões e para, com isto, melhor orientar quanto à janela de aplicação da correção de divergência esférica o que resulta em obter o então denominado “verdadeiras amplitudes”. O meio simulado é formado por camadas plano-horizontais, homogêneas e isotrópicas. O método da refletividade é uma forma de solução da equação de onda para o referido modelo, o que torna possível um entendimento do problema em estudo. Para se chegar aos resultados obtidos foram calculados sismogramas sintéticos através do programa P-SV-SH desenvolvido por Sandmeier (1998), e curvas do espalhamento geométrico em função do tempo para o modelo estudado como descrito por Newman (1973). Demonstramos como uma das conclusões que a partir dos dados do modelo (velocidades, espessuras, densidades e profundidades) uma equação para a correção de espalhamento geométrico visando às “verdadeiras amplitudes” não é de fácil obtenção. O objetivo maior então deveria ser obter um painel da função de divergência esférica para corrigir as verdadeiras amplitudes.

Estudo do método de empilhamento SRC e aplicação em dados sintéticos e reais

O empilhamento sísmico é um método de imageamento utilizado para simulação de seções sísmicas com afastamento nulo (AN) a partir de dados de cobertura múltipla, produzindo um aumento da razão sinal/ruído na seção simulada. Este método é realizado usando o processamento convencional denominado ponto médio comum (PMC), através das chamadas correções sobretempo normal e sobretempo oblíquo (STN/STO). Alternativamente, um novo método de empilhamento sísmico, denominado superfície de reflexão comum (SRC) pode ser usado para simular seções sísmicas em afastamento nulo (AN). O método SRC tem como principais características: 1) a construção de uma aproximação do tempo de trânsito de reflexões primárias na vizinhança de um raio de reflexão normal; 2) esta aproximação de empilhamento tem como parâmetros os ângulos de emergência do raio normal e os raios de curvaturas de duas ondas hipotéticas denominadas onda ponto de incidência normal (onda PIN) e onda normal (onda N). O método pressupõe ainda que a velocidade da onda longitudinal na parte superior do modelo próxima à superfície da Terra seja conhecida. Para determinar adequadamente os parâmetros supracitados, fundamentais para o método de imageamento SRC, são necessárias estratégias de busca que utilizem casos particulares da aproximação de segunda ordem hiperbólica dos tempos de trânsito. As estratégias de busca apresentadas são: SRC pragmático-extendito e SRC global-local. Neste trabalho é apresentada uma revisão extensa dos fundamentos teóricos de método de empilhamento SRC, como também das estratégias de busca dos parâmetros SRC, dos algoritmos de implementação, e da delimitação de aberturas do operador de empilhamento SRC usando conceitos de zonas de Fresnel e o fator de estiramento da correção STN. Como um exemplo de aplicação dos atributos ou parâmetros SRC é apresentado o resultado da correção STN, sem estiramento do sinal sísmico, obtido utilizando as velocidades estimadas dos parâmetros SRC. Também, são apresentados os resultados da aplicação dos três métodos de empilhamento (PMC (STN/STO), SRC pragmático-extendido, SRC global-local) nos dados sintéticos Marmousi e real terrestre da bacia paleozóica Amazônica. Finalmente apresentam-se as discussões sobre os resultados, quanto à resolução da imagem sísmica e custo computacional.

Identificação de reflexões múltiplas, utilizando os parâmetros das frentes de onda PIN e normal

As reflexões múltiplas presentes nos sismogramas ocultam informações importantes sobre os refletores em subsuperfície e, podem até tornar completamente invisíveis as reflexões primárias, como no caso dos sismogramas marinhos, que muitas das vezes, exibem uma aparência anelar com fortes superposições das reflexões múltiplas, sobre as reflexões primárias. Problema este que tem sido alvo de importantes pesquisas, com o intuito de identificar, atenuar e/ou eliminá-las, através de vários métodos populares. O objetivo principal deste trabalho é a identificação das reflexões múltiplas. Com essa finalidade foi gerada, por modelamento direto, uma seção sísmica, com afastamento nulo (AN), contendo reflexões primárias e múltiplas simétricas de primeira ordem. Posteriormente, foi aplicada a migração cinemática do tipo Kirchhoff para obter o modelo em profundidade, apresentando uma boa recuperação dos refletores, bem como a presença de um refletor fictício, quando comparado com a seção anteriormente especificada. Foi obtida uma seção sísmica AN, do modelo migrado, na qual não é observado o segundo refletor, devido à ausência de contraste de impedância, entre a segunda e terceira camada, sendo este o primeiro indício de que o refletor fictício deste modelo é uma múltipla. Outro indício sobre a existência da múltipla foi a simetria encontrada entre as curvaturas do primeiro e terceiro refletor. Finalmente, foram calculados os parâmetros das frentes de ondas Hipotéticas Ponto de Incidência Normal (PIN) e Normal (N), bem como a velocidade Normal Moveout (NMO), tanto para os eventos de reflexões primárias como os eventos de reflexões múltiplas, para o modelo direto e para o modelo migrado. Em seguida, foram realizadas as comparações destes parâmetros, o que permitiu confirmar a veracidade dos indícios anteriores para a identificação das reflexões múltiplas.

Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

Múltiplas trajetórias juvenis em Guarapuava: territórios e redes de sociabilidade

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

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

Caracterização e processamento de leite bubalino em pó em secador por nebulização

Avaliou-se o processo de secagem por nebulização em "spray dryer" para a obtenção de leite bubalino em pó. Foram realizadas análises físico-químicas no leite in natura e no leite desidratado, para averiguar perdas de macronutrientes ocorridas durante a secagem; análises microbiológicas para assegurar a qualidade higiênico-sanitária do processo e produto e análises sensoriais para estabelecer a aceitabilidade do produto. Foram obtidas isotermas de adsorção e dessorção de umidade a 25°C, para o produto em pó, visando observar o seu comportamento higroscópico. O leite bubalino em pó obtido apresentou como características físicoquímicas: 3,4% de umidade, 23,2% de proteínas, 46,1% de gordura, 4,2% de cinzas, 0,5% de lecitina de soja e 22,6% de carboidratos totais. De acordo com os exames microbiológicos está apto para o consumo humano. O produto foi caracterizado como integral por não ter sofrido qualquer padronização. As perdas ocorridas para a maioria dos macronutrientes durante o processamento, foram inferiores a 5% e apenas o teor de proteínas apresentou variação superior a 15%, devido à desnaturação da caseína em temperatura superior a 80°C. O produto apresentou isotermas de adsorção do tipo III, mesmo comportamento apresentado pelo leite bovino. Através da isoterma de dessorção determinou-se o valor de umidade equivalente à monocamada, obtendo-se o valor médio de 4,4g H₂O/100g b.s., o que validou o processo utilizado. O comportamento de histerese apresentado pelo produto foi similar ao do leite em pó integral, sendo a sua aceitabilidade de 74%, muito próxima a 80% obtido para um leite em pó integral comercial.

Processamento mínimo de pêssegos ‘Aurora-1’: estádio de maturação, embalagens, temperaturas de conservação e aditivos naturais

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

Injúrias mecânicas e processamento mínimo de goiabas: fisiologia e qualidade pós-colheita

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

Quantificação de cianeto total nas etapas de processamento das farinhas de mandioca dos grupos seca e d'água

O objetivo deste trabalho foi quantificar a concentração de cianeto total durante as etapas de produção da farinha de mandioca dos grupos seca e d'água. Em relação à farinha seca, a concentração de cianeto total na raiz de mandioca diminuiu de 160±11,8 mg HCN/kg para 149±12,3 mg HCN/kg após a trituração, 68±2,5 mg HCN/kg após a prensagem e chegando a 5±0,2 mg HCN/kg no produto final, após o processo de torração. Na produção da farinha d'água, a raiz de mandioca apresentava teor de cianeto total de 321±21,6 mg HCN/kg e durante o processo de fermentação da raiz, o teor de cianeto total nas primeiras 24 horas de fermentação era de 297±2,7 mg HCN/kg chegando a 64±2,3 mg HCN/kg após 96 horas em repouso no tanque. Após trituração e prensagem da massa fermentada, os valores diminuíram para 50±0,6 e 36±0,4 mg HCN/kg, respectivamente, obtendo-se no produto final a concentração de 9±0,1 mg HCN/kg, sendo evidenciado a eficiência do processo de destoxificação em ambos os processamentos.

Geração e otimização de trajetórias de um manipulador robótico utilizando algoritmos genéticos

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

Monitoramento por microscopia óptica e processamento digital de imagens do processo de conformação cerâmica por conformação com amidos comerciais

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

Processamento e caracterização de snack extrudado a partir de farinhas de quirera de arroz e de bandinha de feijão

Este trabalho teve por objetivo desenvolver nova formulação de snack por extrusão termoplástica a partir de mistura de farinhas de quirera de arroz e de bandinha de feijão, bem como avaliar o potencial nutricional, tecnológico e sensorial do novo produto. A farinha de bandinha de feijão carioca foi incorporada à farinha de quirera de arroz na proporção de 30%. O snack foi produzido em extrusora monorrosca, escala piloto. Os parâmetros de extrusão foram fixos, utilizando-se três zonas de extrusão com temperaturas de 40, 60 e 85 °C; velocidade do parafuso de 177 rpm; taxa de alimentação de 292 g.min-1, e matriz circular de 3,85 mm de diâmetro. A amostra de snack foi submetida a caracterizações fisicoquímica, tecnológica e sensorial. Observou-se efeito significativo da farinha de bandinha de feijão no aumento dos teores proteico e de fibras no snack obtido, quando comparada à farinha de quirera de arroz. Em relação às características tecnológicas do produto, obteve-se 0,17 g.cm-3 para densidade aparente, 7,75 para o índice de expansão e 435, g.f para a dureza instrumental. A formulação estudada foi aceita sensorialmente, com índice de aceitação para impressão global de 76%. Conclui-se que é possível produzir snacks por extrusão a partir da incorporação de 30% de farinha de bandinha de feijão à farinha de quirera de arroz, resultando em produto aceito sensorialmente e com adequado valor nutricional.