Carta de sensibilidade ambiental para derramamentos de óleo em rodovias: proposta aplicada a Estrada dos Tamoios (SP-099)

Pós-graduação em Geociências e Meio Ambiente - IGCE

Avaliação da vulnerabilidade natural à contaminação do solo e aqüífero do reservatório Batatã - São Luís (MA)

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

Falhas de transferência da porção norte da Bacia de Santos interpretadas a partir de dados sísmicos: sua influência na evolução e deformação da bacia

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Wave equation depth migration using complex Padé approximation

Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domínio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido às ondas evanescentes e produz imagens com menos ruídos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos íngremes. Os resultados foram comparados àqueles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD).

Riscos potenciais de impacto ambiental negativo da suinocultura na bacia do córrego Bebedouro, Uberlândia-MG

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

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.

Uma nova abordagem para interpretação de anomalias gravimétricas regionais e residuais aplicada ao estudo da organização crustal: exemplo da Região Norte do Piauí e Noroeste do Ceará

A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas.

Análise de velocidade por otimização do semblance na reflexão sísmica

Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, v_int,n, as espessuras intervalares, z_n, e as velocidades média-quadrática, v_RMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática.

Estudo da contaminação por hidrocarbonetos em postos de serviços de Abaetetuba-PA utilizando georadar

O trabalho apresenta os resultados de um levantamento Geofísico utilizando o método Radar de Penetração de Solo (GPR) ou Georadar para detectar possíveis zonas de contaminação provocadas por vazamentos de derivados de hidrocarbonetos em postos de serviços da região urbana do município de Abaetetuba, no estado do Pará. A metodologia foi aplicada em postos de serviços porque eles constituem uma das principais fontes potenciais urbanas de contaminação do solo e de aqüíferos rasos por combustíveis. Conceitos básicos sobre a contaminação por derivados do petróleo e sua interação com o subsolo são apresentados, juntamente com os princípios básicos que permitem o entendimento do funcionamento do método GPR para o problema abordado. Durante o trabalho, foram realizadas medidas do nível de água em poços rasos, visando a elaboração de um mapa de fluxo subterrâneo. A interpretação geofísica foi auxiliada pelo conhecimento do comportamento do fluxo hídrico subterrâneo local, que mostra o sentido de movimentação da provável contaminação. A correlação entre os dados de GPR, os dados de fluxo e as informações sobre o histórico dos postos levaram a classificá-los em suspeitos de produzir contaminação e possivelmente contaminados.

Geofísica aplicada à Arqueologia: investigação no Sítio Histórico Engenho Murutucu, em Belém, Pará

Esta pesquisa visou à indicação de possíveis localizações, através da prospecção geofísica, de vestígios materiais de ocupação e de alicerces da edificação antiga da senzala no Engenho Murutucu, local de grande importância e relevância para a história do Estado do Pará. Tendo como base a orientação de um mapa magnético realizado num levantamento anterior, foram realizados 26 perfis de GPR, utilizando uma antena monoestática de 400 MHz. Os perfis alcançaram uma profundidade aproximada de 1,75 m. Após a análise dos dados magnéticos e registros de GPR, as informações obtidas foram correlacionadas, resultando nas indicações almejadas. Tais indicações servirão de guia para futuras intervenções arqueológicas no sítio.

Mapeamento de uma interface separando dois meios homogêneos através da reconstituição gravimétrica compacta: aplicação a dados do LEPLAC-IV

O problema do mapeamento de uma interface separando dois meios homogêneos tem sido tradicionalmente resolvido através da continuação para baixo da anomalia gravimétrica observada. Este procedimento requer a utilização de um filtro passa-baixa ou de um parâmetro de amortecimento, que tendem a diminuir a resolução do relevo estimado da interface. Além disso, o uso da continuação para baixo implica supor (desnecessariamente) que a interface é uma superfície harmônica. Desta forma, a utilização da continuação para baixo é restrita a interfaces que apresentam um relevo suave relativamente à sua profundidade média. As restrições impostas pela continuação para baixo do campo gravimétrico são reduzidas na técnica de inversão proposta para este estudo, cujo objetivo é mapear a interface crosta-manto baseada na incorporação de informações a priori sobre o mínimo momento de inércia da fonte causadora, em relação a um eixo de concentração de massa pré-estabelecido e coincidente com o topo da fonte anômala. Além da introdução de informações sobre o mínimo momento, foi incluída também a possibilidade de haver variação lateral e em profundidade de densidade em setores previamente estabelecidos, o que permite estender a aplicação desta técnica a situações geológicas mais complexas que envolvam a reconstituição de corpos justapostos, que sejam de diferentes naturezas, tais como aqueles presentes em zonas de transição crustal. O método de inversão apresentado foi avaliado em testes usando dados sintéticos e os resultados obtidos mostraram a eficiência desta técnica em recuperar corpos alongados horizontalmente, que apresentem como topo ou base, interfaces arbitrárias a serem mapeadas. Foi também constatada a eficácia do método em reconstituir corpos que comportem variação lateral de densidade que pudessem retratar situações em que se deseja investigar a natureza do embasamento em regiões de transição crustal. Com base nos resultados destes testes, a técnica foi então empregada na inversão gravimétrica de cinco perfis do LEPLAC IV que cruzam a região oceânica adjacente à margem continental sudeste brasileira. Os resultados desta aplicação permitiram inferir a natureza dos corpos que constituem o arcabouço estrutural da região onde encontram-se inseridas a Dorsal de São Paulo e a Zona de Fratura de Porto Alegre, assim como mapear a Moho subjacente.

Mapeamento de jazidas de bauxitas no Distrito Paragominas – Açailândia (PA-MA) pelo método da sísmica de refração

Nos últimos 30 anos no Brasil, o mapeamento de camada bauxítica tem sido sempre efetuado através de perfilagens de poços. Neste estudo, aplicou-se o método da sísmica de refração rasa pela primeira vez, em cerca de 3.000m de perfis, em caráter experimental em algumas localidades ao longo da rodovia BR-10, assim como no Platô Esperança (região de Jabuti-PA) e perto da cidade de Açailândia (MA); com objetivo de delimitar a camada de bauxita. Embora se tenha utilizado, em algumas ocasiões, durante a interpretação, o método de distância crítica no cálculo de profundidade, foi possível, mediante o emprego do Método Recíproco Generalizado (GRM), o processamento e interpretação dos dados, para se obter a delimitação contínua da camada investigada. O GRM apresenta vantagens sobre os métodos tradicionais, tais como poder delimitar camadas onduladas, detetar camadas escondidas, além da sensibilidade para camadas com mergulho até 20°. A partir dós gráficos da análise de velocidade, foi calculada para a camada de bauxita uma velocidade média de 840 m/s, e mediante o emprego dos tempos de profundidade chegou-se a sua delimitação: ela possui um comportamento irregular, tendo uma espessura média de 5m. A pesquisa geofísica realizada nas três áreas-piloto tem grande importância porque o método utilizado poderá ser aplicado a outras regiões, em condições geológicas semelhantes.

Variabilidade espacial dos atributos do solo, erosão e suscetibilidade magnética de uma vertente em Gilbués, Pi

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

Análise espacial do conteúdo de água da bacia hidrográfica do Rio Pardo

O presente trabalho tem por objetivo avaliar a distribuição espacial, por geoestatística, da produção de água gerada pelo modelo hidrológico matemático SWAT 2009 (Soil and Water Assessment Tool, versão 2009), da parte inicial da bacia hidrográfica do Rio Pardo – SP. Foi utilizado um Sistema de Informação Geográfica (SIG) associado a uma interface com o modelo SWAT para a confecção do banco de dados. Para isto, as informações de entrada necessárias para avaliar a produção de lâmina de água (mm), que infiltrou e armazenou em cada sub-bacia gerada pelo SWAT, referem-se a dados tabulares climáticos e de parâmetros físicos e químicos de solo e a planos de informações como: o Modelo Numérico do Terreno (MNT), Mapa de Uso do Solo e Mapa de Solos. A amostragem geoestatística foi representada por uma malha irregular georreferenciada com 43 pontos localizados na parte central de cada sub-bacia representando a quantidade de água produzida. A análise geoestatística foi realizada pela construção dos variogramas e posteriormente a confecção dos mapas interpolados por krigagem. Do resultado obtido observou-se que a produção de água apresentou dependência espacial e que esta ocorreu de forma homogênea, tanto para os maiores como para os menores valores de produção de água encontrados na bacia.

Bi-Lipschitz A-equivalence of K-equivalent map-germs

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