Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.

Empilhamento sísmico por superfície de reflexão comum: um novo algoritmo usando otimização global e local

O método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC) produz a simulação de seções com afastamento nulo (NA) a partir dos dados de cobertura múltipla. Para meios 2D, o operador de empilhamento SRC depende de três parâmetros que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀), o raio de curvatura da onda ponto de incidência normal (R_NIP) e o raio de curvatura da onda normal (R_N). O problema crucial para a implementação do método de empilhamento SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros ótimos associados a cada ponto de amostragem da seção AN a ser simulada. No presente trabalho foi desenvolvido uma nova sequência de processamento para a simulação de seções AN por meio do método de empilhamento SRC. Neste novo algoritmo, a determinação dos três parâmetros ótimos que definem o operador de empilhamento SRC é realizada em três etapas: na primeira etapa são estimados dois parâmetros (β°₀ e R°_NIP) por meio de uma busca global bidimensional nos dados de cobertura múltipla. Na segunda etapa é usado o valor de β°₀ estimado para determinar-se o terceiro parâmetro (R°_N) através de uma busca global unidimensional na seção AN resultante da primeira etapa. Em ambas etapas as buscas globais são realizadas aplicando o método de otimização Simulated Annealing (SA). Na terceira etapa são determinados os três parâmetros finais (β₀, R_NIP e R_N) através uma busca local tridimensional aplicando o método de otimização Variable Metric (VM) nos dados de cobertura múltipla. Nesta última etapa é usado o trio de parâmetros (β°₀, R°_NIP, R°_N) estimado nas duas etapas anteriores como aproximação inicial. Com o propósito de simular corretamente os eventos com mergulhos conflitantes, este novo algoritmo prevê a determinação de dois trios de parâmetros associados a pontos de amostragem da seção AN onde há intersecção de eventos. Em outras palavras, nos pontos da seção AN onde dois eventos sísmicos se cruzam são determinados dois trios de parâmetros SRC, os quais serão usados conjuntamente na simulação dos eventos com mergulhos conflitantes. Para avaliar a precisão e eficiência do novo algoritmo, este foi aplicado em dados sintéticos de dois modelos: um com interfaces contínuas e outro com uma interface descontinua. As seções AN simuladas têm elevada razão sinal-ruído e mostram uma clara definição dos eventos refletidos e difratados. A comparação das seções AN simuladas com as suas similares obtidas por modelamento direto mostra uma correta simulação de reflexões e difrações. Além disso, a comparação dos valores dos três parâmetros otimizados com os seus correspondentes valores exatos calculados por modelamento direto revela também um alto grau de precisão. Usando a aproximação hiperbólica dos tempos de trânsito, porém sob a condição de R_NIP = R_N, foi desenvolvido um novo algoritmo para a simulação de seções AN contendo predominantemente campos de ondas difratados. De forma similar ao algoritmo de empilhamento SRC, este algoritmo denominado empilhamento por Superfícies de Difração Comum (SDC) também usa os métodos de otimização SA e VM para determinar a dupla de parâmetros ótimos (β₀, R_NIP) que definem o melhor operador de empilhamento SDC. Na primeira etapa utiliza-se o método de otimização SA para determinar os parâmetros iniciais β°₀ e R°_NIP usando o operador de empilhamento com grande abertura. Na segunda etapa, usando os valores estimados de β°₀ e R°_NIP, são melhorados as estimativas do parâmetro R_NIP por meio da aplicação do algoritmo VM na seção AN resultante da primeira etapa. Na terceira etapa são determinados os melhores valores de β°₀ e R°_NIP por meio da aplicação do algoritmo VM nos dados de cobertura múltipla. Vale salientar que a aparente repetição de processos tem como efeito a atenuação progressiva dos eventos refletidos. A aplicação do algoritmo de empilhamento SDC em dados sintéticos contendo campos de ondas refletidos e difratados, produz como resultado principal uma seção AN simulada contendo eventos difratados claramente definidos. Como uma aplicação direta deste resultado na interpretação de dados sísmicos, a migração pós-empilhamento em profundidade da seção AN simulada produz uma seção com a localização correta dos pontos difratores associados às descontinuidades do modelo.

Inversão de dados de ondas de superfície Rayleigh

Neste trabalho compilamos informações sobre um grande número de medidas de velocidade de grupo para ondas Rayleigh do modo fundamental, com período até 100 segundos. Tais dados consistiram de informações retiradas da literatura geofísica e cobriram toda a Terra. Parte dos dados foi organizada em trabalhos anteriores e uma segunda parte foi apresentada aqui de forma inédita. Para a América do Sul, selecionamos os principais conjuntos de dados de tais ondas e elaboramos diversos perfis onde a distribuição de velocidade de ondas cisalhantes foi obtida a partir da inversão das curvas de dispersão de velocidade de grupo. Tais perfis serviram para termos uma ideia inicial da estrutura interna da Terra em nosso continente. Com o conjunto global de dados de velocidade de grupo foi possível obtermos os mapas de distribuição lateral de valores de velocidade para cada período referencial entre 20 e 100 segundos. Tais mapas foram produzidos da mesma forma que os mapas de velocidade de fase de ROSA (1986), onde a amostragem for para realizada para blocas medindo 10x10 graus, englobando toda a Terra, em projeção mercator. O valor de velocidade de grupo em cada bloco, para cada período, foi obtido a partir da inversão estocástica dos dados de anomalia de velocidade em relação aos modelos regionalizados de JORDAN (1981) com os valores de velocidade de grupo de ROSA et al. (1992). Os mapas de velocidade de grupo obtidos aqui foram então empregados, na América do Sul, com os valores de velocidade de fase dos mapas obtidos por ROSA (1986). Assim, foi possível determinarmos, em profundidade, os mapas de variação de velocidade de onda cisalhante e os mapas de distribuição de valores de densidade. Com isto, pudemos construir o primeiro mapa de profundidade do Moho (todo do Manto Superior) da América do Sul.

Avaliação do efeito da combinação de pectina, gelatina e alginato de sódio sobre as características de gel de fruta estruturada a partir de "mix" de polpa de cajá e mamão, por meio da metodologia de superfície de resposta

Este trabalho teve por objetivo estabelecer procedimento tecnológico para produção de fruta estruturada a partir de "mix" de polpa de cajá e mamão, procurando unir as propriedades sensoriais de cada uma das frutas e potencializar as características funcionais do produto final. Avaliou-se o efeito da combinação de pectina, gelatina e alginato de sódio, via metodologia de superfície de resposta, nas características do gel de fruta. As polpas de mamão e cajá e os estruturados obtidos foram caracterizados com relação aos compostos funcionais, avaliando-se o teor de taninos e carotenóides totais, além das análises de composição centesimal, pH, acidez titulável, sólidos solúveis, açúcares, atividade de água, carboidratos e valor energético total. Os resultados obtidos através do planejamento experimental indicam que para o estruturado misto de cajá e mamão, somente o aumento da concentração de gelatina afeta a firmeza do produto final. Os estruturados de frutas desenvolvidos apresentaram boa aceitação sensorial para todos os atributos avaliados. Com relação à intenção de compra, 70% dos provadores responderam que provavelmente ou certamente, comprariam o estruturado misto de cajá e mamão se o encontrassem à venda.

Marcação e atenuação de múltiplas de superfície livre, processamento e imageamento em dados sísmicos marinhos

O presente trabalho teve como objetivos a análise e atenuação de múltiplas de superfície livre, o processamento e o imageamento de dados sísmicos marinhos visando obter imagens migradas de utilidade à interpretação geológica voltada à exploração de petróleo. Foi dada uma atenção sistemática ao estudo de múltiplas de superfície livre do ponto de vista do filtro de predição baseado na teoria da comunicação com o objetivo de melhor aplicar o filtro WH de deconvolução preditiva, na etapa de processamento, posterior à correção NMO, embora outros métodos possam ser considerados mais competitivos. A identificação e a atenuação de reflexões múltiplas em dados sísmicos reais continua como um grande desafio no processamento de dados sísmicos, uma vez que elas são consideradas como ruído. No entanto, sendo este ruído classificado como coerente, várias técnicas foram desenvolvidas objetivando sua atenuação para evitar erros em cascata nas etapas posteriores como processamento, marcação de eventos, inversão tomográfica, imageamento, e finalmente na interpretação geológica das imagens obtidas. Outro aspecto do trabalho foi estabelecer um fluxograma de processamento-imageamento sendo a atenuação de múltiplas de superfície livre uma etapa central. Foram obtidas seções migradas em tempo e em profundidade onde se permite a interpretação que se deseja. O desenvolvimento deste trabalho foi realizado com os pacotes CWP/SU e MatLab.

Modelos acoplados do IPCC-AR4 e o gradiente meridional de temperatura da superfície do mar no Atlântico Tropical: relações com a precipitação no norte do nordeste do Brasil

Este artigo mostra como três modelos acoplados do Intergovernmental on Panel Climate Change - (IPCC-AR4), o FGOALS1.0G - LASG do Institute of Atmospheric Physics of China, o GISSER da National Aeronautics Space Admnistration (NASA) e o GFDL_CM2 da National Oceanic and Atmospheric Administration (NOAA), simularam a variabilidade do gradiente meridional de Temperatura da Superfície do Mar (TSM), entre os meses de fevereiro a maio, no Atlântico Tropical (1901-1999). A precipitação durante a estação chuvosa (fevereiro a maio) no setor norte do Nordeste do Brasil (NEB) foi também analisada pelos três modelos e comparada com as observações. Os modelos GISSER e FGOALS1.0G mostraram melhor desempenho na simulação do sinal do gradiente meridional de TSM no Atlântico Tropical para o período de 1901 a 1999. Destaca-se que os modelos apresentaram um melhor desempenho na simulação da tendência decadal, conseguindo explicar entre 50% a 80% da variabilidade do gradiente, com a TSM do setor sul sendo melhor simulada.

Obtenção e caracterização de carvão ativado de caroço de buriti (Mauritia flexuosa L. f.) para a avaliação do processo de adsorção de cobre (II)

Na região amazônica, algumas indústrias despejam cobre nos corpos receptores que, em elevadas concentrações, é tóxico para os seres vivos. A remoção de cobre de efluentes industriais é realizada por diversos processos como a adsorção. Neste trabalho mostrou-se o resultado da adsorção de cobre (II) em carvão ativado de caroço de buriti carbonizado a 400 °C e ativado a 900 °C. O carvão ativado foi caracterizado em termos de área específica, tamanho dos poros, densidades aparente e real, porosidade, microscopia eletrônica de varredura, conteúdo de cinzas, pH, umidade, carbono fixo e grupos funcionais de superfície. O estudo de equilíbrio de adsorção avaliou a influência do diâmetro das partículas do carvão, do tempo de contato adsorvente/adsorbato, do pH e da concentração inicial da solução de cobre (II) sobre a remoção de cobre (II). Com base nos resultados, concluiu-se que há uma maior eficiência de remoção de cobre (II) para diâmetro < 0,595 mm, tempo de contato de 300 minutos, pH 4,0 e concentrações iniciais de cobre (II) de 50 e 80 mg L^-1. O modelo matemático de Langmuir foi o que melhor se ajustou aos dados de equilíbrio de adsorção. A partir do tempo de contato de 15 minutos todas as concentrações de equilíbrio ficaram abaixo do máximo permitido de 1,0 mg L^-1 previsto pela legislação vigente para lançamento de efluentes em corpos receptores.

Otimização da desidratação osmótica de filés de mapará (Hypophthalmus edentatus) através da metodologia de superfície de resposta

A desidratação osmótica de filés de mapará tem sido utilizada como pré-tratamento antes da secagem e refrigeração de alimentos como um meio de reduzir a umidade dos produtos. O objetivo do trabalho foi otimizar a transferência de massa ocorrida durante a desidratação osmótica de mapará através da metodologia de superfície de resposta. A desidratação osmótica foi realizada de acordo com um planejamento fatorial completo 2³, com oito pontos fatoriais, três centrais e seis axiais, totalizando 17 ensaios, em que o valor das variáveis dependentes, a perda de peso (PP), a perda de água (PA), o ganho de sólidos (GS), GS/PA e ácido tiobarbitúrico (TBA), são funções das variáveis independentes, temperatura, concentração de NaCl e tempo de imersão. A análise estatística aplicada aos dados experimentais foi realizada através do erro puro e pelo SS residual. Aplicando a metodologia de superfície de resposta, a condição ótima de convergência de menor perda de peso, maior perda de água e menor ganho de sal correspondeu a uma concentração de cloreto de sódio de 22%, temperatura de 46 °C e tempo de 7 horas.

Análise dos impactos das diferentes formas de ocupação da superfície sobre as condições meteorológicas na região de Santarém, Pará

O objetivo desta dissertação foi obter informações referentes ao uso e ocupação do solo na região de Santarém, em diferentes anos das últimas décadas, para melhor representar os efeitos causados pelas modificações das propriedades da superfície sobre as condições atmosféricas simuladas por modelos numéricos de tempo e clima. As superfícies continentais caracterizam-se por causar efeitos substanciais sobre a atmosfera e, consequentemente, influir na qualidade das previsões de tempo e de clima. Por outro lado, o desmatamento contribui com as mudanças climáticas, por eliminar grandes quantidades de gases de efeito estufa para a atmosfera. Estas atividades também causam efeitos na saúde publica, na agricultura, nos recursos florestais, nos recursos faunísticos e nos recursos hídricos. Além disso, a substituição da superfície natural por pastagem ou agricultura altera as propriedades térmicas e radiativas da superfície, gerando modificações nas condições atmosféricas locais, regionais e globais. Neste trabalho foram analisados períodos representativos de possíveis mudanças climáticas na região, identificados a partir do tratamento e analise estatística de dados climatológicos de estações meteorológicas de superfície, bem como a evolução temporal e quantitativa do desmatamento na região de estudo com os dados do Projeto PRODES Monitoramento da Floresta Amazônica Brasileira por Satélite). Para avaliar os efeitos atmosféricos das mudanças no uso e ocupação do solo utilizou-se como base o mapa de vegetação do IBGE, e a inclusão da classe “desmatamento” ao mesmo em diferentes períodos analisados (anos de 1997e 2009) trabalhadas no software Arc. Gis. 9.2. Foram criados arquivos de dados de tipos de superfície compatíveis com a leitura do modelo BRAMS, que foi então utilizado para simular os diferentes efeitos desses mapas temáticos de uso e ocupação do solo na atmosfera local. Os resultados indicam uma tendência de aumento da precipitação média anual e da frequência média de dias com precipitação, diminuição da temperatura média das máximas e aumento da temperatura média das mínimas ao longo dos anos na região de Santarém. A área de estudo, até o ano de 1997, registrou um desmatamento de 19,44% e até o ano de 2009 passou para 25,54%. As simulações com os arquivos gerados de uso e ocupação do solo para 1997 e 2009 apresentaram poucas variações para os diferentes mapas temáticos em suas variáveis (temperatura, umidade e fluxos de calor sensível e latente), quando considerado os valores médios da área total simulada. Porém, quando se considera pequenas áreas localizadas somente sobre as regiões que sofreram maiores modificações, observam-se maiores influências com o aumento do desmatamento.

Análise dos efeitos térmicos de superfície na cidade de Belém-Pará-Brasil utilizando imagens de satélite

Utilizando dados observacionais de precipitação e temperatura do ar para o período de 1967 a 2008, do qual se elaborou a climatologia destas variáveis, e analisando as condições atmosféricas influenciadas pelo fenômeno ENOS, pretende-se identificar, caracterizar e analisar os efeitos térmicos espaciais na cidade de Belém-PA, a partir de dados meteorológicos de superfície e de imagens do sensor MODIS, o qual está disposto sobre a plataforma do satélite Aqua. Observou-se que, a média da precipitação anual foi de 2978,6 mm/ano, e que apresenta tendência de aumento ao longo dos anos, comportamento semelhante observou-se para a temperatura do ar. Em geral, os resultados mostram dois núcleos de maiores intensidade de temperaturas da superfície, um na cidade de Belém e outro na cidade vizinha, Ananindeua. Estes variam espacialmente e temporalmente de intensidade. Durante eventos de La Niña, o núcleo da ilha de calor fica localizado em bairros mais próximos a baía do Guajará, enquanto que durante eventos de El Niño estes bairros apresentam temperaturas mais amenas do que os bairros mais afastados dos corpos hídricos. Observou-se ainda que, a amplitude térmica da temperatura superficial entre áreas urbana e rural variam bastante, com a maior variação de 30,8°C e a menor de 16,8°C. Neste sentido, as maiores temperaturas da superfície foram observadas nos bairros com baixo NDVI, conseqüência de uma urbanização mais densa. As superfícies urbanas e as superfícies vegetadas apresentam relações de causa e efeito muito próximas, principalmente, durante o período menos chuvoso, isto pode ser percebido pela correlação que apresenta valor acima de 50%. Este estudo apresenta resultados que auxiliam no melhor entendimento do comportamento e dos efeitos térmicos espaciais e temporais na cidade de Belém, pois o uso de imagens do satélite é de fundamental importância para a identificação e caracterização das condições ambientais climáticas e ilhas de calor urbanas.