Atenuação de múltiplas e compressão do pulso fonte em dados de sísmica de reflexão utilizando o filtro Kalman-Bucy

O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função b_i(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas.

Estimativa de parâmetros elásticos em meios anisotrópicos

As medidas de amplitude, polarização e vagarosidade contem informações sobre o meio onde a propagação de onda ocorre. Esta tese investiga esses dados com objetivo de estimar as propriedades elásticas deste meio. Coeficientes de reflexão podem ser estimados das amplitudes dos dados e dependem de forma não linear dos contrastes dos parâmetros elásticos e do contraste de densidade entre os meios separados por uma interface. Quando o contraste de impedância é fraco, as aproximações lineares para a refletividade qP são mais convenientes para inversão da densidade e dos parâmetros elásticos usando as análises de amplitude versus ângulo de incidência (AVO) e amplitude versus a direção do plano de incidência (AVD). Escrevendo as equações de Zoepprittz de forma separada nos permite escrever uma solução destas equações em termos das matrizes de impedância e polarização. Usando esta solução são determinadas aproximações lineares para a refletividade da onda qP considerando fraco contraste de impedância, fraca anisotropia mas com classe de simetria de arbitrária. As linearizações são avaliadas para diferentes geometrias de aquisição e várias escolhas do meio de referência. Estas aproximações apresentam bom desempenho comparado com o valor exato do coeficiente de reflexão da onda qP e de suas ondas convertidas para incidências de até 30° e meios que obedecem à hipótese de fraca anisotropia. Um conjunto de fraturas orientado é representado efetivamente por um meio transversalmente isotrópico (TI), as aproximações lineares da refletividade da onda qP podem ser usadas para estimar a orientação de fratura. Partindo deste pressuposto este problema consiste em estimar a orientação do eixo de simetria a partir de dados de refletividade de onda qP. Este trabalho mostra que são necessários múltiplos azimutes e múltiplas incidências para se obter uma estimativa estável. Também é mostrado que apenas os coeficientes das ondas qS e qT são sensíveis ao mergulho da fratura. Foi investigada a estimativa da anisotropia local através de dados de VSP multiazimutal dos vetores de polarização e vagarosidade. Foram usadas medidas da componente vertical do vetor de vagarosidade e o vetor de polarização de ondas qP diretas e refletidas. O esquema de inversão é validado através de exemplos sintéticos considerando diferentes escolhas do vetor normal à frente de onda no meio de referência, meios de referências e geometria de aquisição. Esta análise mostra que somente um subgrupo dos parâmetros elástico pode ser estimado. Uma importante aplicação desta metodologia é o seu potencial para a determinação de classes de anisotropia. A aplicação desta metodologia aos dados do mar de Java mostra que os modelos isotrópicos e TIV são inadequados para o ajuste desses dados.

Migração por equação de onda em meios anisotrópicos com correção de amplitude

Este trabalho discute dois aspectos da migração em profundidade através da continuação para baixo dos campos de onda: o tratamento de modos evanescentes e a correção da amplitude dos eventos migrados. Estes dois aspectos são discutidos em meios isotrópicos e para uma classe de meios anisotrópicos. Migrações por diferenças finitas (FD) e por diferenças finitas e Fourier (FFD) podem ser instáveis em meios com forte variação lateral de velocidade. Estes métodos utilizam aproximações de Padé reais para representar o operador que descreve a propagação de ondas descendentes. Estas abordagens não são capazes de tratar corretamente os modos evanescentes, o que pode levar à instabilidades numéricas em meios com forte variação lateral de velocidade. Uma solução possível para esse problema é utilizar aproximação de Padé complexa, que consegue melhor representar os modos evanescentes associados às reflexões pós-críticas, e neste trabalho esta aproximação é utilizada para obter algoritmos FD e híbrido FD/FFD estáveis para migração em meios transversalmente isotrópicos com eixo de simetria vertical (VTI), mesmo na presença de forte variação nas propriedades elásticas do meio. A estabilidade dos algoritmos propostos para meios VTI foi validada através da resposta ao impulso do operador de migração e pela sua aplicação na migração de dados sintéticos, em meios fortemente heterogêneos. Métodos de migração por equação de onda em meios heterogêneos não tratam corretamente a amplitude dos eventos durante a propagação. As equações de onda unidirecionais tradicionais descrevem corretamente apenas a parte cinemática da propagação do campo de onda. Assim, para uma descrição correta das amplitudes deve-se usar as equações de onda unidirecionais de amplitude verdadeira. Em meios verticalmente heterogêneos, as equações de onda unidirecionais de amplitude verdadeira podem ser resolvidas analiticamente. Em meios lateralmente heterogêneos, essas equações não possuem uma solução analítica. Mesmo soluções numéricas tendem a ser instáveis. Para melhorar a compensação de amplitude na migração, em meios com variação lateral de velocidade, é proposto uma aproximação estável para solução da equação de onda unidirecional de amplitude verdadeira. Esta nova aproximação é implementada nas migrações split-step e diferenças finitas e Fourier (FFD). O algoritmo split-step com correção de amplitude foi estendido para meios VTI. A migração pré e pós-empilhamento de dados sintéticos, em meios isotrópicos e anisotrópicos, confirmam o melhor tratamento das amplitudes e estabilidade dos algoritmos propostos.

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Análise das aproximações RPP e RSP para meios isotrópicos

Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de incidência mostram-se mais precisas que as aproximações lineares em função da média angular. Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento).

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

Estudo do comportamento das feições lineares observado em dados aeromagnéticos da Bacia do Paraná

O processamento dos dados aeromagnéticos em imagens revelou-se como um bom método de interpretação, que aqui denominamos de método das imagens sombreadas. A aplicação deste método foi feita em dados aeromagnéticos de uma área da Bacia do Paraná, com o objetivo de realçar os lineamentos que não são evidentes nos mapas de contorno. Para subsidiar os nossos estudos, utilizamos imagens de LANDSAT porque muitos lineamentos presentes na bacia não possuem resposta ao levantamento aeromagnético. Com a integração destas duas imagens definimos duas direções dominantes para os lineamentos nos seguintes intervalos: N40° — 60°W e N40° — 65°E. Com a definição deste padrão, obtivemos informações do arcabouço estrutural e, consequentemente, uma idéia da evolução tectônica da bacia. A outra parte desenvolvida neste trabalho foi o cálculo do vetor de magnetização total nos principais lineamentos da bacia, e eles apresentaram um total de 5 inclinações de magnetização diferentes. Essas inclinações foram explicadas pela superposição do efeito da magnetização induzida, efeito de reversão do campo geomagnético, e o efeito da magnetização anisotrópica causada pela desmagnetização nos diques.

Estudo das feições lineares observadas em dados aeromagnéticos da Bacia do Solimões a partir de imagens digitais de sombreamento do relevo anômalo total

A Bacia do Solimões é caracterizada por um padrão magnético complexo, expresso ao longo de toda sua extensão, e produzido pela superposição de anomalias magnéticas que se apresentam sob a forma de feições lineares. Tais feições refletem as diferentes atividades tectônicas que atuaram na Região Amazônica através do Pré-Cambriano e Fanerozóico. Neste trabalho, foi aplicado nos dados aeromagnéticos da Bacia do Solimões um método de processamento de imagens digitais de sombreamento do relevo magnético anômalo total que, dadas suas características metodológicas, permitindo utilizá-lo como um filtro direcional, possibilitou a definição de aspectos relevantes no contexto das relações entre os lineamentos magnéticos E/W, NE/SW, e NW/SE. Nesse sentido, foram identificados padrões de lineamentos magnéticos que refletem a existência de zonas de cisalhamento transcorrentes dextrais, orientadas preferencialmente na direção E/W. A interação entre os diversos segmentos transcorrentes promoveu o desenvolvimento de um regime predominantemente transpressivo, representado por falhas reversas associadas aos lineamentos magnéticos E/W e NE/SW; e duplexes direcionais formando falhas em flor positivas e negativas associadas, respectivamente, aos lineamentos magnéticos orientados na direção N70-80E e N70-80W. A análise quantitativa permitiu explicar dois aspectos importantes em relação às feições lineares observadas nas imagens digitais. O primeiro mostra, através de modelamentos baseados na superposição de prismas bidimensionais, que estas feições lineares magnéticas podem ser explicadas pela superposição de fontes profundas intraembasamento altamente magnéticas, e fontes rasas de alta frequência, sendo estas associadas a falhas reversas ao longo dos níveis de diabásio, presentes em forma de intrusões nos sedimentos paleozóicos da Bacia do Solimões. O segundo aspecto, baseado na utilização dos conceitos da cross-covariância, constata a presença de 'offsets' dextrais, associados aos lineamentos magnéticos NE—SW, ao longo da direção E—W. Este fato mostra, quantitativamente, que o padrão magnético desta região pode ser explicado pela presença de zonas de cisalhamento transcorrestes dextrais, cujos processos tectônicos associados foram fortemente condicionados por zonas de fraquezas pré-existentes (Pré-Cambrianas, Paleozóicas) durante o Mesozóico e Cenozóico.

Avaliação dos níveis de ruído e das soluções arquitetônicas adotadas em ambientes das unidades de urgência e emergência no município de Belém-PA

Nos grandes centros urbanos pode-se constatar que o nível de ruído é tão elevado que ameaça a integridade psicológica e física dos habitantes. Entretanto, observa-se que muito do ruído no ambiente hospitalar origina-se de dentro do hospital, tendo como uma das principais fontes geradoras de ruído nas unidades os equipamentos e a conversação entre a equipe hospitalar. A perda auditiva induzida pelo ruído (PAIR) é uma das consequências que o ruído elevado pode acarretar à saúde, entre outros danos. Assim, o ambiente hospitalar que deveria ser sereno e silencioso transforma-se em um espaço ruidoso e estressante, aumentando a ansiedade e a percepção dolorosa do paciente, reduzindo o sono e estendendo a convalescença. Objetivo: Avaliar as condições acústicas dos principais ambientes nas unidades de urgência e emergência no município de Belém, Estado do Pará, Brasil. Verificar a adequação às normas sanitárias e técnicas vigentes, com a simulação de um ambiente piloto validado com a utilização do software ODEON e apresentar soluções virtuais de condicionamento acústico. Material e Métodos: foi utilizado o equipamento Medidor de Nível Sonoro com certificado de calibração da Rede Brasileira de Calibração (RBC) para aferição dos NPS nos períodos diurno e noturno, de acordo com o estabelecido pela NBR 10.151 (2000), além da observação e identificação das fontes de ruído dos ambientes. Resultados: Os valores obtidos durante as aferições dos NPS nos ambientes do Hospital A foram de 58 dB(A) a 70 dB(A) e do Hospital B foram de 62 dB(A) a 70 dB(A). O ambiente piloto foi selecionado de acordo com os resultados e com o tempo de permanência da equipe de saúde e dos pacientes. Conclusão: foi desenvolvido o modelo computacional validado do ambiente piloto, gerando um modelo proposto com intervenções arquitetônicas visando o conforto acústico da equipe de saúde e dos pacientes. Os ambientes pesquisados não apresentaram níveis de ruído superiores a 85 dB(A), comprovando que não há risco ocupacional para as equipes de saúde. Entretanto, níveis de ruídos elevados, por mais que não causem PAIR nos trabalhadores, são responsáveis por estimular ou até mesmo por desencadear diversos sintomas que prejudicam a capacidade laborativa da equipe de saúde, além de aumentar a convalescença e prejudicar a recuperação dos pacientes.

Inversão em geofísica de poço: um estudo sobre ambiguidade

A ambiguidade na inversão de dados de geofísica de poço é estudada através da análise fatorial Q-modal. Este método é baseado na análise de um número finito de soluções aceitáveis, que são ordenadas, no espaço de soluções, segundo a direção de maior ambiguidade. A análise da variação dos parâmetros ao longo dessas soluções ordenadas permite caracterizar aqueles que são mais influentes na ambiguidade. Como a análise Q-modal é baseada na determinação de uma região de ambiguidade, obtida de modo empírico a partir de um número finito de soluções aceitáveis, é possível analisar a ambiguidade devida não só a erros nas observações, como também a pequenos erros no modelo interpretativo. Além disso, a análise pode ser aplicada mesmo quando os modelos interpretativos ou a relação entre os parâmetros não são lineares. A análise fatorial é feita utilizando-se dados sintéticos, e então comparada com a análise por decomposição em valores singulares, mostrando-se mais eficaz, uma vez que requer premissas menos restritivas, permitindo, desse modo, caracterizar a ambiguidade de modo mais realístico. A partir da determinação dos parâmetros com maior influência na ambiguidade do modelo é possível reparametrizá-lo, agrupando-os em um único parâmetro, redefinindo assim o modelo interpretativo. Apesar desta reparametrização incorrer na perda de resolução dos parâmetros agrupados, o novo modelo tem sua ambiguidade bastante reduzida.

Adsorção de cromo (VI) por carvão ativado granular de soluções diluídas utilizando um sistema batelada sob pH controlado

Na Amazônia o cromo é empregado principalmente na indústria de couro e de madeira, sendo responsável por vários problemas de saúde porque é tóxico para os seres vivos. A remoção de cromo de efluentes industriais é feita por meio de diversos processos como a adsorção. Este trabalho mostra os resultados da adsorção de Cr(VI) por carvão ativado granular comercial (CAG) como adsorvente de soluções diluídas empregando um sistema de adsorção batelada com controle de pH. Os grupos funcionais da superfície do CAG foram determinados pelo método de Boehm. Além disso, o efeito do pH na adsorção de Cr(VI), o equilíbrio e a cinética de adsorção foram estudados nas condições experimentais (pH = 6, MA = 6g, tempo de adsorção 90min.). Na superfície do CAG, os grupos carboxílicos foram determinados em maior concentração (MAS=0,43 mmol/gCAG), estes, presentes em concentrações elevadas aumentam a adsorção do metal, principalmente em valores de pH ácidos. A capacidade de adsorção é dependente do pH da solução, devido a sua influência nas propriedades de superfície do CAG e nas diferentes formas iônicas das soluções de Cr(VI). Os dados de equilíbrio da adsorção foram ajustados satisfatoriamente pela isoterma de Langmuir (R²=0,988), tipo favorável. A partir da cinética de adsorção a 5mg/L e 20mg/L, os resultados obtidos foram compatíveis com o valor limite preconizado na legislação nacional (Res. nº 357/05). Portanto, para o sistema experimental utilizando CAG foi eficiente na remoção de Cr(VI) a partir de correntes líquidas contendo baixas concentrações do metal.

Otimização multiobjetivo da compensação de potência reativa em redes de distribuição considerando restrições de distorção harmônica

A localização de bancos de capacitores nas redes de distribuição de energia elétrica, corretamente dimensionados, busca compensar eventuais excessos de circulação de potência reativa pelas linhas, o que implica a redução de custos operacionais pela redução das perdas de energia e um aumento da capacidade de transmissão de potência ativa assegurando os níveis estabelecidos de tensão e fator de potência simultaneamente. A proliferação das cargas não lineares provocou uma mudança nos cenários de estudo dos sistemas elétricos de potência devido aos efeitos nocivos que os harmônicos gerados por elas ocasionam sobre a qualidade da energia elétrica. Considerando este novo cenário, esta tese tem como objetivo geral desenvolver uma ferramenta computacional utilizando técnicas de inteligência computacional apoiada em algoritmos genéticos (AG), para a otimização multiobjetivo da compensação da potência reativa em redes elétricas de distribuição capaz de localizar e dimensionar de forma ótima as unidades de compensação necessárias para obter os melhores benefícios econômicos e a manutenção dos índices de qualidade da energia estabelecidos pelas normas brasileiras. Como Inovação Tecnológica do trabalho a ferramenta computacional desenvolvida permite otimizar a compensação da potência reativa para melhorar do fator de potência em redes de distribuição contaminadas com harmônicos que, diferentemente de métodos anteriores, não só emprega bancos de capacitores, mas também filtros de harmônicos com esse objetivo. Utiliza-se o algoritmo NSGA-II, que determina as soluções ótimas de Pareto para o problema e permite ao especialista determinar as soluções mais efetivas. A proposta para a solução do problema apresenta várias inovações podendo-se destacar que a solução obtida permite determinar a compensação de potência reativa com capacitores em sistemas com certa penetração harmônica, atendendo a normas de qualidade de energia pertinentes, com relação aos níveis de distorção harmônica tolerados.