Desenvolvimento de técnicas de filtragem e sua aplicação a dados aeromagnéticos da Bacia Amazônica

Dois métodos de filtragem são combinados neste trabalho com a finalidade de ressaltar determinadas feições, em um mapa aeromagnético, tais como lineamentos e anomalias de grande porte, provindos do embasamento. A filtragem visual consiste na observação estereoscópica de pares de mapas, especialmente elaborados para este fim. A filtragem numérica consiste na eliminação dos componentes espectrais indesejáveis (ruído) e no realce dos componentes espectrais de interesse (sinal). Ambos os métodos foram aplicados a mapas aeromagnéticos da Bacia do Médio Amazonas. As operações de filtragem numérica foram todas efetuadas no domínio dos espaços, através de convolução dos dados com filtros de dimensões 13x13. A análise preliminar do espectro é enfatizada como guia na elaboração de filtros. Os resultados da aplicação de cada filtro são discutidos em detalhe. Além das filtragens visual e numérica, vários modelos de falhas e de anomalias do intraembasamento e do supraembasemento são apresentados. A análise de modelos constitui um método de interpretação independente e complementar dos métodos de filtragens. Uma interpretação estrutural do embasamento da Bacia do Médio Amazonas é incluída. As direções principais dos lineamentos aeromagnéticos são E-W, NW60°, N-S, NF45° e NE6O°. Dois grabens sao bem evidenciados um com direção N-S, o mais recente próximo à borda do Alto do Purus. O outro, mais antigo, situa-se a leste do primeiro e tem direção NE45°. Duas anomalias de grande e médio porte (15-60km) possuem magnetização remanente com altos valores de inclinação negativa, o que corrobora a posição muito mais ao sul da Bacia Amazônica, no Paleozóico. Na integração dos dados aeromagnéticos com outros dados geofísicos, como sísmica e gravimetria, e resultados de fotointerpretação geológica, a consistência das informações independentes, em relação à interpretação aeromagnética, concorre para aumentar a viabilidade e a confiabilidade das técnicas aplica das neste trabalho.

Estudo comparativo entre estereotomografia e da tomografia da onda NIP: aplicação em dados sintéticos e reais

A determinação de um acurado modelo de velocidades é um requisito fundamental para a realização do imageamento sísmico. Métodos novos como a Estereotomografia préempilhamento e a Tomografia da onda NIP são ferramentas poderosas e bastante sugestivas para este propósito. Basicamente, a Estereotomografia pré-empilhamento se baseia no conceito de eventos localmente coerentes interpretados como reflexões primárias e associados com pares de segmentos de raios, que tem um mesmo ponto de reflexão em profundidade. Na Tomografia da onda NIP um evento sísmico é representado por uma onda hipotética NIP, que está relacionada a um ponto de reflexão em profundidade. Os atributos da onda NIP são determinados no decorrer do Empilhamento de Superfície de Reflexão Comum (empilhamento CRS). Este trabalho tem como objetivo, fazer um estudo comparativo de ambos os métodos de determinação do modelo de velocidades em profundidade. Assim, é realizada uma revisão dos fundamentos teóricos de ambos os métodos tomográficos, destacando as principais diferenças e aplicando estas aproximações em um dado sintético e um dado real marinho (linha sísmica 214-2660 da Bacia do Jequitinhonha). Para avaliar os modelos de velocidades encontrados pelas aproximações, foi utilizada a migração pré-empilhamento em profundidade do tipo Kirchhoff e também as famílias de imagem comum (CIG). Os resultados mostraram que ambos os métodos tomográficos fornecem modelos de velocidades representativos. Contudo, constatou-se que a estereotomografia tem melhor desempenho em meios com variações laterais de velocidades, porém, aplicável somente em dados pré-empilhados com alta razão sinal/ruído.

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.

Interpretação aeromagnética automática com uso da equação homogênea de Euler e sua aplicação na Bacia do Solimões

ABSTRACT: We present here a methodology for the rapid interpretation of aeromagnetic data in three dimensions. An estimation of the x, y and z coordinates of prismatic elements is obtained through the application of "Euler's Homogeneous equation" to the data. In this application, it is necessary to have only the total magnetic field and its derivatives. These components can be measured or calculated from the total field data. In the use of Euler's Homogeneous equation, the structural index, the coordinates of the corners of the prism and the depth to the top of the prism are unknown vectors. Inversion of the data by classical least-squares methods renders the problem ill-conditioned. However, the inverse problem can be stabilized by the introduction of both a priori information within the parameter vector together with a weighting matrix. The algorithm was tested with synthetic and real data in a low magnetic latitude region and the results were satisfactory. The applicability of the theorem and its ambiguity caused by the lack of information about the direction of total magnetization, inherent in all automatic methods, is also discussed. As an application, an area within the Solimões basin was chosen to test the method. Since 1977, the Solimões basin has become a center of exploration activity, motivated by the first discovery of gas bearing sandstones within the Monte Alegre formation. Since then, seismic investigations and drilling have been carried on in the region. A knowledge of basement structures is of great importance in the location of oil traps and understanding the tectonic history of this region. Through the application of this method a preliminary estimate of the areal distribution and depth of interbasement and sedimentary magnetic sources was obtained.

Investigação com os métodos GPR e magnético de urnas funerárias do Sítio Nossa Senhora do Perpétuo Socorro (Pará)

Este trabalho apresenta os resultados da aplicação de métodos geofísicos na eleição de áreas potenciais para salvamento arqueológico no Sítio Nossa Senhora do Perpétuo Socorro no estado do Pará. O Sítio possui urnas funerárias indígenas aflorantes e sub-aflorantes de bastante interesse arqueológico. Os métodos utilizados foram Magnetometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais mais apropriados para escavação reduzindo bastante as áreas potenciais. O GPR foi usado para confirmar as anomalias detectadas pelas medidas magnéticas e proporcionar uma melhor avaliação espacial tanto horizontal quanto verticalmente, diminuindo os erros cometidos quando se identificam anomalias magnéticas que não são causadas por feições arqueológicas. A metodologia de estudo baseou-se na comparação de anomalias obtidas em terrenos desconhecidos com terrenos onde comprovadamente havia a existência de elementos mapeáveis como urnas aflorantes e raízes, caracterizando e agrupando essas anomalias.

Métodos geofísicos aplicados a prospecção de água subterrânea no município de Ponta de Pedras - Ilha de Marajó - Pará

Este trabalho foi realizado numa área de 25 km² em Ponta de Pedras, Ilha de Marajó. Tem por objetivo o dimensionamento geométrico e seleção de zonas favoráveis à água subterrânea. Foram utilizados dois métodos geofísicos: sísmica de refração e eletroresistividade. O método sísmico foi empregado em caráter experimental, visando verificar a viabilidade de sua aplicação na área. Os resultados da interpretação confirmaram a suspeita inicial, da existência de horizontes geológicos pouco profundos, não possíveis de serem distinguidos pelo método sísmico. Com o método de eletroresistividade foram feitas sondagens verticais e perfis horizontais. Os resultados da interpretação indicaram áreas promissoras na parte central e numa faixa ao norte da área prospectada. A espessura do aquífero raso varia entre 0 e 50 m.

Métodos geofísicos aplicados à arqueologia no Estado do Pará

Este trabalho constitui uma primeira tentativa de aplicação de métodos geofísicos à pesquisa de sítios arqueológicos no Brasil. Esses métodos foram usados com a finalidade de detetar feições arqueológicas no subsolo, em dois sítios arqueológicos no estado do Pará: o sítio PA-SA-4, da fase arqueológica Areão e o sítio Ilha dos Bichos da fase Marajoara. No PA-SA-4, foi empregado o método magnético, enquanto que na Ilha dos Bichos, além do método magnético, foi usado o eletro-resistivo. O método eletro-resistivo possibilitou não só a determinação da camada de ocupação, pois esta camada se comporta geoeletricamente diferente das camadas que compõem as formações geológicas da região, como também permitiu detetar camadas compactas de argila, planas e horizontais, de origem arqueológica, situadas no lado oriental da Ilha dos Bichos. O método magnético levou à deteção de antigos locais de queima, devido a magnetização termo-remanente adquirida pelas argilas no processo de aquecimento/resfriamento, permitindo a descoberta de antigos fornos, bem como de uma urna funerária sob espessa camada de argila queimada.

Métodos geofísicos aplicados ao mapeamento geológico numa área da Região Bragantina-Pará

Este trabalho trata da aplicação de métodos geofísicos, numa área piloto da região Bragantina, NE do Pará, com a finalidade de auxiliar no mapeamento geológico. A área acha-se coberta por Sedimentos Terciários espessos e apresenta ocorrências de rochas industriais como calcário e granito, representando uma superfície aflorante menor que 5% da área. Estudos geoquímicos, realizados por Goulart (1978), revelaram que no horizonte A, é impossível diferenciar os materiais de superfície quanto às suas origens. Isto torna difícil o mapeamento geológico de superfície, podendo levar a extrapolações exageradas dos contatos, o que ocorre no mapeamento pré-existente. Os métodos geofísicos aplicados foram: Eletroresistividade e Magnetometria. Também foram realizados alguns furos de sondagens mecânicas. Os resultados obtidos contribuíram para um melhor conhecimento das feições estruturais mapeadas, bem como tornou possível a descoberta de outras estruturas importantes. Definiram-se também duas áreas para subsequente estudo de detalhe, visando a separação do calcário, em melhor e pior qualidade para a fabricação do cimento. Quanto à região dos granitos, recomenda-se que seja feito o seu mapeamento sob forma de "ilhas de granito" em meio à espessas coberturas de solos e sedimentos.

Métodos geofísicos aplicados ao ultramafito da Serra de Quatipuru - Estado do Pará

O objetivo deste trabalho foi a aplicação de métodos geofísicos para a avaliação econômica de ocorrência de níquel, cromo e amianto, em uma área de 6,25 km² - Malha 2 - Serra de Quatipuru, localizada a cerca de 50 km a SW de Conceição do Araguaia, no sul do Estado do Pará. Inicialmente foi aplicado o método magnético em caráter de reconhecimento, usando magnetômetros de campo total e componente vertical visando localizar zonas anômalas de amianto e cromita, uma vez que estes minerais se encontram associados com a magnetita. Em seguida foi usado o método eletromagnético de "dip angle" com o objetivo de definir zonas condutoras associadas a possíveis depósitos minerais ou características estruturais. Com base na interpretação destes dois métodos, definiu-se os locais de aplicação do método de polarização elétrica induzida-resistividade, que teve caráter de detalhe visando como alvo principal depósitos de sulfetos disseminados. As possíveis ocorrências minerais na área foram correlacionadas com anomalias geofísicas e confirmadas através de testemunhos de sondagem concluindo-se não serem de importância econômica. Contudo a aplicação destes métodos geofísicos foi útil, pois possibilitou o levantamento detalhado da área, que de outra forma só poderia ser realizado por meio dos métodos mais convencionais, tais como sondagem mecânica através de perfurações, os quais são muito caros e demorados.

Métodos de otimização global aplicados na busca dos parâmetros SRC

O método de empilhamento por Superfície de Reflexão Comum (SRC) produz seções simuladas de afastamento nulo (AN) por meio do somatório de eventos sísmicos dos dados de cobertura múltipla contidos nas superfícies de empilhamento. Este método não depende do modelo de velocidade do meio, apenas requer o conhecimento a priori da velocidade próxima a superfície. A simulação de seções AN por este método de empilhamento utiliza uma aproximação hiperbólica de segunda ordem do tempo de trânsito de raios paraxiais para definir a superfície de empilhamento ou operador de empilhamento SRC. Para meios 2D este operador depende de três atributos cinemáticos de duas ondas hipotéticas (ondas PIN e N), observados no ponto de emergência do raio central com incidência normal, 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_PIN) e o raio de curvatura da onda normal (R_N). Portanto, o problema de otimização no método SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros (β₀, R_PIN, R_N) ótimos associados a cada ponto de amostragem da seção AN a ser simulada. A determinação simultânea destes parâmetros pode ser realizada por meio de processos de busca global (ou otimização global) multidimensional, utilizando como função objetivo algum critério de coerência. O problema de otimização no método SRC é muito importante para o bom desempenho no que diz respeito a qualidade dos resultados e principalmente ao custo computacional, comparado com os métodos tradicionalmente utilizados na indústria sísmica. Existem várias estratégias de busca para determinar estes parâmetros baseados em buscas sistemáticas e usando algoritmos de otimização, podendo estimar apenas um parâmetro de cada vez, ou dois ou os três parâmetros simultaneamente. Levando em conta a estratégia de busca por meio da aplicação de otimização global, estes três parâmetros podem ser estimados através de dois procedimentos: no primeiro caso os três parâmetros podem ser estimados simultaneamente e no segundo caso inicialmente podem ser determinados simultaneamente dois parâmetros (β₀, R_PIN) e posteriormente o terceiro parâmetro (R_N) usando os valores dos dois parâmetros já conhecidos. Neste trabalho apresenta-se a aplicação e comparação de quatro algoritmos de otimização global para encontrar os parâmetros SRC ótimos, estes são: Simulated Annealing (SA), Very Fast Simulated Annealing (VFSA), Differential Evolution (DE) e Controlled Rando Search - 2 (CRS2). Como resultados importantes são apresentados a aplicação de cada método de otimização e a comparação entre os métodos quanto a eficácia, eficiência e confiabilidade para determinar os melhores parâmetros SRC. Posteriormente, aplicando as estratégias de busca global para a determinação destes parâmetros, por meio do método de otimização VFSA que teve o melhor desempenho foi realizado o empilhamento SRC a partir dos dados Marmousi, isto é, foi realizado um empilhamento SRC usando dois parâmetros (β₀, R_PIN) estimados por busca global e outro empilhamento SRC usando os três parâmetros (β₀, R_PIN, R_N) também estimados por busca global.

Aplicação da cromatografia líquida de alta eficiência na investigação de hidrocarbonetos policíclicos aromáticos em testemunhos sedimentares

Hidrocarbonetos policíclicos aromáticos (HPAs) são compostos orgânicos originários de fontes naturais e antrópicas e, por apresentarem potencial carcinogênico e mutagênico, são considerados poluentes prioritários por agências ambientais. Desta forma, métodos analíticos para investigação de tais compostos que sejam rápidos e de baixo custo são de relevância considerável para o monitoramento ambiental. O presente trabalho teve como objetivos otimizar um método analítico para HPAs utilizando cromatografia líquida com detecção por arranjo de diodos (HPLC-DAD) e aplicar em testemunho sedimentar de região estuarina. Para otimização e avaliação do método, uma coluna sedimentar de 46 cm de comprimento foi coletada na foz do Rio Tucunduba (Belém, Pará) e seccionada em porções de 2 cm (subamostras). Após secagem, 30 g de cada porção foram extraídos com mistura de diclorometano em acetona (1:1) em ultrassom por 40 min. Os extratos obtidos foram centrifugados, purificados em sílica-gel, adaptação em funil necessária principalmente para reter partículas finas, e em seguida concentrados em rotaevaporador à vácuo e, por fim, filtrados com membrana de nylon 0,22 μm antes da injeção no HPLC. Amostras fortificadas com padrões analíticos de 16 HPAs e brancos também foram processados da mesma maneira. Um conjunto de parâmetros para validação do método foi investigado e observou-se: (1) boa linearidade: as curvas de calibração (analíticas) apresentaram coeficientes de correlação elevadas; (2) precisão adequada: obteve-se desvio padrão relativo dentro do aceitável, sendo o mínimo de 2,1% para acenaftileno e máximo de 19,7% para o fluoranteno; (3) limites de detecção baixos: entre 0,004 a 1,085 ng g g^-1, viabilizando análises em concentrações reais in situ; (4) recuperação adequada para traços: sendo a mínima de 40,0% para o acenaftileno e máxima de 103,1% para o benzo(k)fluoranteno. As concentrações de HPAs totais variaram nas seções do testemunho sedimentar entre 60,77 - 783,3 ng g^-1 de sedimento seco. O método otimizado mostrou-se vantajoso com relação aos tradicionais que utilizam extrator soxhlet e colunas de adsorventes para purificação de extratos por minimizar o tempo de extração e reduzir custos com uso de volumes menores de solventes para purificação do extrato. A limitação do método, porém, foi a coeluição do criseno e do benzo(a)antraceno e a sobreposição do fluoreno e acenafteno, além da quantificação benzo(g,h,i)perileno. Essa limitação provavelmente está associada à eficiência da coluna cromatográfica disponível para a análise, que é para aplicação geral. O método mostrou-se aplicável a amostras estuarinas complexas e ricas em silte e argila. Razões diagnósticas de HPAs parentais indicam fontes petrogênicas a profundidades de 24 – 26 cm, 28 – 30 cm; e fontes pirolíticas a profundidades de 6 – 8 cm, 10 – 12 cm e 14 – 16 cm respectivamente.

Análise da incidência de casos de dengue na área urbana de Belém – PA: uma aplicação de modelos espaciais e temporais

Neste trabalho, foi realizado um estudo de mapeamento de áreas de incidência e previsões para os casos de dengue na área urbana de Belém. Para as previsões foi utilizada à incidência de dengue com a precipitação pluviométrica a partir de modelos estatísticos, baseados na metodologia de Box e Jenkins de series temporais. O período do estudo foi de 05 anos (2007-2011). Na pesquisa temos métodos multivariados de series temporais, com uso de função de transferência e modelos espaciais, em que se analisou a existência de autocorrelações espaciais na variável em estudo. Os resultados das análises dos dados de incidência de casos de dengue e precipitação mostraram que, o aumento no número de casos de dengue acompanha o aumento na precipitação, demonstrando a relação direta entre o número de casos de dengue e a precipitação nos anos em estudo. O modelo de previsão construído para a incidência de casos de dengue apresentou um bom ajuste com resultados satisfatórios podendo, neste caso, ser utilizado na previsão da dengue. Em relação à análise espacial, foi possível uma visualização da incidência de casos na área urbana de Belém, com as respectivas áreas de incidência, mostrando os níveis de significância em porcentagem. Para o período estudado observou-se o comportamento e as variações dos casos de dengue, com destaque para quatro bairros: Marco, Guamá, Pedreira e Tapanã, com possíveis influências destes bairros nas áreas (bairros) vizinhas. Portanto, o presente estudo evidencia a contribuição para o planejamento das ações de controle da dengue, ao servir de instrumento no apoio às decisões na área de saúde pública.

Sustentabilidade em assentamentos rurais do estado do pará: uma aplicação do barômetro da sustentabilidade

Os assentamentos na Amazônia são fundamentais para a distribuição de terras, porém sua implementação ainda traz muitos problemas ambientais e sociais, o que tem levado a críticas sobre as condições de sustentabilidade dos mesmos. Nos últimos anos, têm-se buscado desenvolver métodos de monitoramento dos assentamentos rurais através de indicadores de sustentabilidades, porém os desafios ainda são enormes. Independente do método adotado, os estudos apontam que os principais passos na análise da sustentabilidade dos assentamentos são a identificação dos temas, o estabelecimento dos indicadores, a definição dos limites dos indicadores e a avaliação dos mesmos na construção do índice de sustentabilidade. Desta forma, o presente trabalho tem por objetivo analisar o nível de sustentabilidade dos assentamentos do Pará, por meio da ferramenta do Barômetro da Sustentabilidade (BS), comparando duas categorias de assentamentos: Projetos de Assentamentos Convencionais (PAC) e Projetos de Assentamentos Diferenciados (PAD). Foram selecionados 28 indicadores extraídos de um diagnóstico socioambiental efetuado pelo Instituo Nacional de Colonização e Reforma Agrária -INCRA e testou-se a ferramenta BS para 12 assentamentos. Os resultados mostram que há ligeira diferença entre os assentamentos, principalmente com relação à dimensão ambiental, mas as duas modalidades de assentamentos ocupam a mesma posição na escala do Barômetro da Sustentabilidade é Potencialmente Insustentável.

Influência de sistemas meteorológicos de grande escala sobre a distribuição de raios na Amazônia Oriental, com aplicação no setor elétrico

Foi feito um estudo da associação entre a densidade de ocorrência de raios com condições meteorológicas no leste da região Amazônica entre os anos de 2006 a 2008. A região estudada foi limitada por uma grade geográfica entre as latitudes 0º a -10º e longitudes -53º a -43º. Foram desenvolvidos métodos computacionais de processamento estatístico e geração de mapas de ocorrências de raios para diferentes intervalos de tempo. Esses métodos foram aplicados para determinar os pontos de altas ocorrências de raios ao longo de linhas de transmissão de energia elétrica, a fim de oferecer subsídios para a proteção e operação dos sistemas elétricos da região. Foram utilizados dados da rede de detecção de raios do SIPAM para redimensionar a detecção do sistema de alcance intercontinental STARNET, e subsequentemente, foram obtidos mapas de densidade de raios mais uniformes e realistas. Esses mapas foram interpretados em períodos semanais e sazonais, considerando as observações simultâneas de chuva, vento, radiação de onda longa e a presença de sistemas meteorológicos de grande escala. Zona de Convergência Intertropical e Zona de Convergência do Atlântico Sul (ZCIT, ZCAS) definidos pela convergência de ventos e anomalias negativas de Radiação de Ondas Longas ( ROL). Os resultados mostraram boa correspondência entre as áreas de convecção e a intensa ocorrência de raios. Nas regiões de mais alta atividade elétrica atmosférica, foi observada ocorrência de convergência dos ventos e anomalias negativas de ROL, para situações de presença de ZCAS e ZCIT sobre a região. A atividade de raios também coincidiu, com algumas exceções, com áreas de maior precipitação em períodos semanais. Foi observado também que, durante um trimestre seco, os segmentos de linhas de transmissão de Mojú-Tailândia e Jacundá-Marabá apresentaram maior incidência de raios. O pico da densidade de raios chegou a 18 raios/km²/trimestre em alguns locais. Durante o trimestre chuvoso, a densidade geral de raios foi 86% maior e apresentou distribuição espacial mais uniforme quando comparada com um trimestre seco. Este estudo mostrou não somente como as características meteorológicas influenciaram na distribuição e quantidade raios, mas também que a elevada atividade de descargas elétricas na região deve ser uma preocupação importante para os sistemas de distribuição de energia elétrica e outras atividades humanas na Amazônia Oriental.