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 geofísico regional sobre águas subterrâneas na Ilha do Marajó-Pará-Brasil

Nos últimos dez anos foram realizadas na parte leste da Ilha de Marajó (região dos campos naturais) pelo IDESP e NCGG, mais de 800 SEVs para fins hidrogeológicos. Na época, grande parte dessas SEVs não foram totalmente interpretadas em forma quantitativa, devido à falta de recursos técnicos para fazê-lo de forma eficiente. Agora, usando meios mais modernos para interpretação automática de SEVs, voltou-se a interpretá-las com a finalidade de apresentar uma visão regional dos principais aquíferos da área, agrupar as SEVs em famílias características, testar até que ponto essa interpretação é confiável e propor o modelamento bidimensional como técnica alternativa para interpretar as SEVs realizadas em certos locais da área em questão. Como resultado dessa interpretação, com base na teoria convencional dos meios estratificados, foram definidos três tipos de sistemas de aquíferos. 1. O primeiro, denominado de aquífero profundo, situado a profundidades maiores que 50m, estende-se por toda a região prospectada, estando provavelmente associada às camadas superiores da Formação Marajó ou às litologias altamente resistivas das camadas mais profundas do Grupo Pará. 2. O segundo, denominado de aquífero raso e de média profundidade, localiza-se na parte sul e sudeste da região a profundidades compreendidas entre 10 a 50m, e está associado às lentes arenosas do Grupo Pará. 3. O terceiro, é constituído pelos paleocanais e estruturas similares, distribuídos aleatoriamente na região a pouca profundidade. A partir do estudo detalhado das SEVs, decidiu-se classificá-las em 3 famílias características com seus respectivos tipos e apresentar mapas de localização e da espessura dos aquíferos, bem como mapas de condutância longitudinal total e resistividade média da área. Estes últimos, permitem que se divida a região dos campos da Ilha de Marajó em três zonas principais: 1. Uma, altamente resistiva, situada ao sul e sudeste, a qual coincide com os terrenos aflorantes do Grupo Pará. 2. Outra, altamente condutiva, está localizada no centro e norte, onde se encontram aleatoriamente distribuídos os paleocanais e coincide com os terrenos topograficamente mais baixos, geralmente argilosos e embebidos de água salgada, que são procedentes da erosão dos terrenos circundantes topograficamente mais altos. 3. A última é medianamente resistiva e está relacionada com os terrenos vizinhos à cidade de Chaves (noroeste da região dos campos), os quais apresentam semelhanças com os do sul e sudeste da área. Usando-se a técnica de inversão na interpretação de uma SEV característica de cada família, testou-se, através do seu tratamento estatístico, até que ponto os modelos usados na interpretação dessas SEVs (teoria convencional dos meios estratificados) seriam confiáveis. Conclui-se, então, que a alta correlação existente entre os parâmetros dos modelos assumidos (camadas horizontais, isotrópicas e homogêneas) pode-se dever à utilização de modelos geofísicos muito simples para interpretar a complexa geologia de Marajó. Tendo-se verificado que nem sempre é possível aplicar a teoria das SEVs em meios horizontalmente estratificados para interpretar SEVs obtidas em certos locais de Marajó, os quais muitas vezes apresentam bruscas variações laterais de resistividade, passou-se a demonstrar que estas variações laterais afetam profundamente os dados das SEVs, utilizando-se para isto a técnica dos elementos finitos, a qual leva em conta essa variação bidimensional das propriedades físicas do meio. Foi também possível com esta técnica, modelar uma estrutura rasa, semelhante a um paleocanal, concluindo-se que estes resultados sugerem o emprego, duma forma mais profunda, deste tipo de tratamento para os dados obtidos na região dos campos da Ilha de Marajó.

Estudo da hidrogeologia no município de Palestina do Pará usando o método de eletroresistividade

Esse trabalho de geofísica aplicada à hidrogeologia, foi realizado numa área localizada na porção sudeste do estado do Pará, mas precisamente na cidade de Palestina à margem esquerda do rio Araguaia. Ironicamente existem problemas de abastecimento de água potável para a população da cidade. Esse problema deve-se ao fato que, no momento, é antieconômico o tratamento e transporte de água do rio para a cidade. Considerando que a extração de água subterrânea é economicamente mais viável para a solução desse problema, o Departamento de Geofísica e o curso de Pós-graduação de Geofísica da Universidade Federal do Pará (UFPa), por solicitação da prefeitura daquele município e da Fundação Nacional de Saúde (FNS), participaram no levantamento geofísico de eletroresistividade para investigar os possíveis locais de ocorrência de água subterrânea. Foram realizadas, na área, 21 Sondagens Elétricas Verticais (SEVs) medidas na superfície, utilizando o arranjo Schlumberger com abertura máxima dos eletrodos AB e MN de 500m e 50m, respectivamente. Para a sua interpretação foram usados os programas de modelagem direta (EGSLIB/SEV1D) e inversa (EGSLIB/SEV1DIN) de RIJO (1994), de modelos de n camadas horizontais, homogêneas e isotrópicas. Com base nos resultados da interpretação das 21 SEVs, e também das informações da geologia local foram confeccionadas três mapas de contorno e seis secções geoelétricas. Dos seis horizontes que aparecem nessas secções, o mais favorável para a exploração hidrogeológica é o quarto, com espessura variando de 10 a 79 m, estando a uma profundidade média de 50m e com valores de resistividade variando de 150 — 850 Ωm, sendo correlacionavel com a Formação Potí. A analise final dos mapas e das secções geoelétricas, indicou que o local mais apropriado para a perfuração de poços, é uma área próximo à atual área de captação.

Sistemas de liberação controlada de quitosana contendo antigeno capsular Vi de Salmonella Typhi

A utilização do antígeno Vi em vacinas é bastante promissor devido a este proporcionar um alto nível de imunidade em vacinas parenterais. A necessidade pela busca por alternativas para administração de vacinas levou à aplicação da tecnologia da liberação controlada de fármacos no campo da imunização. Nestes sistemas de liberação controlada, as doses administradas são diminuídas, porém o período de imunidade aumenta, já que prolonga a quantidade liberada do antígeno ao longo do tempo. O presente estudo propôs desenvolver e caracterizar um sistema de liberação controlada contendo antígeno Vi, utilizando como polímero veiculador a quitosana. As técnicas de RMN H^-1 e espectroscopia de infravermelho mostraram que o método de extração do antígeno Vi empregado foi satisfatório qualitativamente. A caracterização da quitosana e das nanopartículas através de ensaios de análise térmica mostrou maior estabilidade das partículas em relação à quitosana, além do aumento da temperatura de degradação nas nanopartículas à medida que aumenta a concentração da quitosana. Em relação ao potencial zeta todas as nanopartículas tiveram cargas positivas em pH 7,2, enquanto que, no tamanho, as partículas foram menores à medida que aumentou a quantidade de quitosana no sistema. Na microscopia eletrônica de transmissão, as partículas mostraram-se morfologicamente homogêneas e com formato esférico. Na cinética de adsorção o antígeno, contido em solução, sofreu uma adsorção de 55% nas partículas de quitosana. Com isso, observou-se que é possível criar um sistema de liberação controlada envolvendo nanopartículas de quitosana e antígeno Vi de Salmonella enterica sorotipo Typhi.

Uma abordagem híbrida e semiautomática para estimativa de regiões cobertas por nuvens e sombras em imagens de satélite: análise e avaliação

Os principais objetivos deste trabalho são propor um algoritmo eficiente e o mais automático possível para estimar o que está coberto por regiões de nuvens e sombras em imagens de satélite; e um índice de confiabilidade, que seja aplicado previamente à imagem, visando medir a viabilidade da estimação das regiões cobertas pelos componentes atmosféricos usando tal algoritmo. A motivação vem dos problemas causados por esses elementos, entre eles: dificultam a identificação de objetos de imagem, prejudicam o monitoramento urbano e ambiental, e desfavorecem etapas cruciais do processamento digital de imagens para extrair informações ao usuário, como segmentação e classificação. Através de uma abordagem híbrida, é proposto um método para decompor regiões usando um filtro passa-baixas não-linear de mediana, a fim de mapear as regiões de estrutura (homogêneas), como vegetação, e de textura (heterogêneas), como áreas urbanas, na imagem. Nessas áreas, foram aplicados os métodos de restauração Inpainting por suavização baseado em Transformada Cosseno Discreta (DCT), e Síntese de Textura baseada em modelos, respectivamente. É importante salientar que as técnicas foram modificadas para serem capazes de trabalhar com imagens de características peculiares que são obtidas por meio de sensores de satélite, como por exemplo, as grandes dimensões e a alta variação espectral. Já o índice de confiabilidade, tem como objetivo analisar a imagem que contém as interferências atmosféricas e daí estimar o quão confiável será a redefinição com base no percentual de cobertura de nuvens sobre as regiões de textura e estrutura. Tal índice é composto pela combinação do resultado de algoritmos supervisionados e não-supervisionados envolvendo 3 métricas: Exatidão Global Média (EGM), Medida De Similaridade Estrutural (SSIM) e Confiança Média Dos Pixels (CM). Finalmente, verificou-se a eficácia destas metodologias através de uma avaliação quantitativa (proporcionada pelo índice) e qualitativa (pelas imagens resultantes do processamento), mostrando ser possível a aplicação das técnicas para solucionar os problemas que motivaram a realização deste trabalho.

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.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Green markov - uma abordagem híbrida de política markoviana e simulação discreta para planejamento de alocação de usuários em redes macro-femto

O uso da comunicação de voz e dados através de dispositivos móveis vem aumentando significativamente nos últimos anos. Tal expansão traz algumas dificuldades inerentes, tais como: ampliação constante de capacidade das redes e eficiência energética. Neste contexto, vem se consolidando o conceito de Green networks, que se concentra no esforço para economia de energia e redução de CO2. Neste sentido, este trabalho propõe validar um modelo de uma política baseado em processo markoviano de decisão, visando a otimizar o consumo de energia, QoS e QoE, na alocação de usuários em redes macrocell e femtocell. Para isso o modelo foi inserido no simulador NS-2, aliando a solução analítica markoviana à flexibilidade característica da simulação discreta. A partir dos resultados apresentados na simulação, a política obteve uma economia significativa no consumo energético, melhorando a eficiência energética em até 4%, além de melhorar a qualidade de serviço em relação às redes macrocell e femtocell, demonstrando-se eficaz, de modo a alterar diretamente as métricas de QoS e de QoE.

Evolução das disparidades da extração vegetal e da silvicultura na Amazônia Legal: uma aplicação da cadeia de Markov

O artigo analisa a convergência municipal da produtividade vegetal (extração vegetal e silvicultura) na região da Amazônia Legal entre os anos de 1996 e 2006. Para analisar a convergência, optou-se pela metodologia da matriz de transição de Markov (Processo Estacionário de Primeira Ordem de Markov). Os resultados mostram a existência de 13 classes de convergência da produtividade vegetal. No longo prazo, a hipótese de convergência absoluta não se mantém, visto que 68,23% dos municípios encontram-se numa classe inferior à média municipal, 33,54% em uma classe intermediária acima da média e 13,41% em uma classe superior acima da média.