Estudo de pluma de contaminação por combustível com o método GPR e desenvolvimento de simulador geofísico para fins educativos

Este trabalho consistiu na investigação geofísica da possível contaminação de um posto de combustíveis em Belém-PA, bem como no desenvolvimento de um simulador geofísico para fins educativos. A investigação consistiu no levantamento de 13 perfis com o método geofísico eletromagnético Ground Penetrating Radar (GPR). Os radargramas obtidos, conjuntamente com os dados pré-existentes de sondagens mecânicas, Volatile Organics Compounds (VOC) e Benzeno-Tolueno-Etil-benzeno-Xileno (BTEX), permitiram identificar zonas de baixa reflexão do sinal eletromagnético (zonas atenuadas) em áreas com predominância de areia, o que indicaria contaminação. Como as medidas de GPR foram repetidas no período chuvoso, foi possível observar a redução da atenuação em perfil rico em areia, indicando a lixiviação dos ácidos graxos oriundos da biodegradação dos hidrocarbonetos, bem como o aumento da atenuação provocada pela argila, devido a sua elevada porosidade, e também, possivelmente, pelo aumento da biodegradação dos hidrocarbonetos que ficam retidos na argila. Este resultado constitui o primeiro caso inserido no simulador geofísico, experimento que permite mostrar o uso de diferentes métodos geofísicos. O simulador possui duas telas computacionais, a debaixo para simulação do modelo de subsuperfície escolhido e, a tela de cima, para representação da resposta geofísica obtida com o auxílio de uma unidade ótica que percorre o perfil representado na tela debaixo. Trata-se de uma modelagem analógica com recursos computacionais ainda não reportada na literatura, que permitirá mostrar diversos casos investigados constantes em dissertações e teses. O simulador poderá ser deslocado para empresas, universidades e secretarias, que fazem ou podem fazer uso da Geofísica, além de praças, escolas, etc., contribuindo de forma incisiva para a difusão da Geofísica como ferramenta na gestão de áreas contaminadas bem como para chamar a atenção da população para diferentes tipos de problemas, de modo a instruí-la sobre os mesmos.

SwImax: um simulador em realidade virtual para auxiliar o ensino do padrão WiMAX

O padrão IEEE 802.16, também chamado de WiMAX, é uma tecnologia da rede banda larga sem fio para áreas metropolitanas, utilizado como alternativa para transmissão de sinal de Internet a regiões que não possuem infraestrutura de rede cabeada. Atualmente, o ensino desta tecnologia em sala de aula é meramente teórico, o que dificulta a compreensão dos alunos com relação a determinadas funcionalidades do WiMAX. Nesse sentido, a presente dissertação aborda o projeto de desenvolvimento de um simulador em Realidade Virtual, chamado SwImax, voltado para auxílio do ensino do padrão IEEE 802.16. Assim, o SwImax simula algumas características do funcionamento deste padrão, quais sejam: faixas de frequência de operação, área de cobertura, procedimento de handover, transmissão sem linha de visada, entre outros. A dissertação apresenta os trabalhos correlatos que influenciaram o desenvolvimento do projeto, além de um resumo acerca do padrão IEEE 802.16. O texto apresenta também as ferramentas utilizadas no desenvolvimento do SwImax e a implementação do simulador. Ao final do desenvolvimento, o software foi submetido a uma avaliação dos usuários, de modo que os resultados também são abordados nesta dissertação.

Estudo geofísico para a caracterização geológica e hidrogeológica da subsuperfície rasa em área de deposição de rejeitos sólidos no município de Barcarena-PA

Este estudo foi realizado na Área de Depósito de Rejeitos Sólidos (ADRS) da ALBRAS (Alumínio Brasileiro S.A.), localizada no município de Barcarena, Estado do Pará. Os resultados alcançados permitiram estabelecer um padrão de caracterização da subsuperfície, a partir de medidas geofísicas, fornecendo informações sobre as camadas geológicas e sobre parâmetros hidrogeológicos, tais como direção e velocidade do fluxo subterrâneo. Na realização do trabalho foram utilizados três métodos geofísicos: os métodos elétricos (Potencial Espontâneo e Eletrorresitividade) e eletromagnético (Slingram). A distribuição dos potenciais medidos através do Método do Potencial Espontâneo indicaram o sentido do fluxo subterrâneo local. Durante o estudo, foi montado um experimento de infiltração de solução salina no solo e seu efeito foi monitorado através de imageamento elétrico. Esse experimento permitiu que se estimasse a velocidade e o sentido do fluxo subterrâneo local ao longo da linha de imageamento. Os modelos interpretativos obtidos através da inversão de dados de resistividade aparente obtidos em sondagens elétricas verticais apresentaram uma boa correlação com perfis de raios gama corridos em poços tubulares da área de estudo. Estas Sondagens Elétricas Verticais permitiram que se detalhasse as unidades litológicas rasas através da estimativa de seus valores de resistividade e espessura. A interpretação dos dados eletromagnéticos através da análise de perfis de medidas, sondagens eletromagnéticas e mapas de contornos permitiu identificar a presença de zonas mais condutivas (material argiloso) e zonas mais resistivas (material arenoso). Além disso, observou-se que o mapa de contornos para a freqüência de 14080 Hz apresentou uma boa correlação com o mapa de contornos do Potencial Espontâneo referente a mesma área.

Estudo geofísico sobre a erosão fluvial do Rio Guamá na orla da Universidade Federal do Pará, Belém (PA)

Este trabalho apresenta estudo geofísico de parte de orla fluvial sujeita à erosão, com o objetivo primordial de avaliar o seu uso na detecção de áreas degradadas pela erosão em estágio crítico, antes de seu desmoronamento. Esta área está situada ao longo de 600 m à margem do rio Guamá, entre a ponte do rio Tucunduba e o Porto de Canoagem, dentro do campus da Universidade Federal do Pará (UFPA), na cidade de Belém no Estado do Pará (Brasil). Conta, em alguns trechos, com diferentes tipos de contenção de erosão, entre os quais muro de concreto e sacos de cimento. A área é consequentemente, um laboratório para se testar a eficiência da Geofísica em detectar as zonas em que a ação da erosão, embora ainda não observável em superfície, atue, de modo a auxiliar obras de contenção da orla antes do seu colapso. Foram utilizados quatro métodos geofísicos: Potencial Espontâneo (SP), Eletrorresistividade, Slingram (LIN) e Radar de Penetração do Solo (GPR). As medidas foram levantadas durante a maré baixa e a maré alta na tentativa de mapear os caminhos preferenciais subsuperficiais para a entrada de água trazida pela maré alta e, por esse meio, para o trabalho erosivo. O efeito global da maré é aumentar as oscilações nas medidas de SP, Eletrorresistividade e LIN e aumentar a intensidade das reflexões obtidas com o GPR em relação às medidas obtidas com maré baixa. Na maré alta, as medidas mostram, em relação às medidas obtidas na maré baixa: I) Em zona com erosão, a redução da voltagem obtida com o método SP, a diminuição da resistividade obtida com imageamento realizado com o método da Eletrorresistividade, o aumento da condutividade LIN bem como tanto o aumento da intensidade das reflexões como a perturbação da continuidade dos refletores obtidos com o Método GPR II) Em zona com contenção, o aumento da voltagem SP, o aumento da resistividade obtida com Eletrorresistividade, a diminuição da condutividade LIN e, finalmente, a atenuação de reflexões obtidas com o GPR e redução do efeito global da maré. A despeito das dificuldades, os resultados demonstram que a Geofísica pode ser uma ferramenta auxiliar na previsão de locais onde a queda erosiva do terreno provocado por ela está prestes a ocorrer bem como na análise da eficácia da obra de contenção realizada.

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ó.

SIMPLARCOM: simulador para planejamento de redes de comunicação sem fio utilizando realidade virtual e modelos de propagação

Este simulador é formado pela junção de técnicas de Realidade Virtual com modelos de propagação, desenvolvidos através dos estudos de rádio enlace, que descrevem a perda que o sinal transmitido sofre ao longo do percurso no ambiente. O simulador possui dois módulos. O primeiro permite a criação do ambiente virtual com o posicionamento, sobre um terreno, de prédios, árvores, carros, antenas e outras primitivas que permitem a construção de um ambiente tridimensional customizável. O segundo módulo permite a configuração dos parâmetros relacionados a propagação de sinal de antenas como a potência, a frequência, o ganho, etc., e também selecionar o modelo de propagação para a execução da simulação. Dentro deste segundo módulo, existe um submódulo responsável pelo estudo do planejamento da área de cobertura composta pelas antenas, em outras palavras, este submódulo simula a distância que cada antena no cenário consegue atingir e gera a respectiva área de cobertura. Para demonstrar a eficiência do simulador foram criados dois ambientes virtuais para testes. Um cenário representando um ambiente urbano onde empregou-se um modelo de propagação clássico, Okumura-Hata para cidades pequenas e médias, e um ambiente tridimensional arborizado utilizando um modelo especifico para simulação de propagação para regiões densamente arborizadas, desenvolvido na Universidade Federal do Pará chamado de Lyra-Castro-UFPA.

Estudo de algoritmos para classificação de séries temporais: uma aplicação em qualidade de energia elétrica

Aborda a classificação automática de faltas do tipo curto-circuito em linhas de transmissão. A maioria dos sistemas de transmissão possuem três fases (A, B e C). Por exemplo, um curto-circuito entre as fases A e B pode ser identicado como uma falta\AB". Considerando a possibilidade de um curto-circuito com a fase terra (T), a tarefa ao longo desse trabalho de classificar uma série temporal em uma das 11 faltas possíveis: AT, BT, CT, AB, AC, BC, ABC, ABT, ACT, BCT, ABCT. Estas faltas são responsáveis pela maioria dos distúrbios no sistema elétrico. Cada curto-circuito é representado por uma seqüência (série temporal) e ambos os tipos de classificação, on-line (para cada curto segmento extraído do sinal) e off-line (leva em consideração toda a seqüência), são investigados. Para evitar a atual falta de dados rotulados, o simulador Alternative Transient Program (ATP) é usado para criar uma base de dados rotulada e disponibilizada em domínio público. Alguns trabalhos na literatura não fazem distinção entre as faltas ABC e ABCT. Assim, resultados distinguindo esse dois tipos de faltas adotando técnicas de pré-processamento, diferentes front ends (por exemplo wavelets) e algoritmos de aprendizado (árvores de decisão e redes neurais) são apresentados. O custo computacional estimado durante o estágio de teste de alguns classificadores é investigado e a escolha dos parâmetros dos classificadores é feita a partir de uma seleção automática de modelo. Os resultados obtidos indicam que as árvores de decisão e as redes neurais apresentam melhores resultados quando comparados aos outros classificadores.

DesG: uma nova metodologia para intervenção e coleta de dados magnetotelúricos

Esta tese apresenta duas contribuições para a interpretação geofísica, voltadas ao Método Magnetotelúrico. A primeira trata de uma nova abordagem para a interpretação MT, denominada MÉTODO DESCRITIVO-GEOLÓGICO (DesG), em alusão à incorporação explícita de informação a priori de correlação fácil com a descrição geológica tradicional. O intérprete define por meio de elementos geométricos (pontos e linhas) o arcabouço de feições geológicas bem como fornece valores de resistividade aos corpos geológicos presumidos. O método estima a distribuição de resistividade subsuperficial em termos de fontes anômalas próximas aos elementos geométricos, ajustando as respostas produzidas por estes corpos às medidas de campo. A solução obtida fornece então informações que podem auxiliar na modificação de algumas informações a priori imprecisas, permitindo que sucessivas inversões sejam realizadas até que a solução ajuste os dados e faça sentido geológico. Entre as características relevantes do método destacam-se: (i) os corpos podem apresentar resistividade maior ou menor do que a resistividade do meio encaixante, (ii) vários meios encaixantes contendo ou não corpos anômalos podem ser cortados pelo perfil e (iii) o contraste de resistividade entre corpo e encaixante pode ser abrupto ou gradativo. A aplicação do método a dados sintéticos evidencia, entre outras vantagens, a sua potencialidade para estimar o mergulho de falhas com inclinação variável, que merece especial interesse em Tectônica, e delinear soleiras de diabásio em bacias sedimentares, um sério problema para a prospecção de petróleo. O método permite ainda a interpretação conjunta das causas do efeito estático e das fontes de interesse. A aplicação a dados reais é ilustrada tomando-se como exemplo dados do COPROD2, cuja inversão produziu soluções compatíveis com o conhecimento geológico sobre a área. A segunda contribuição refere-se a desenho de experimento geofísico. Por meio de indicadores diversos, em especial a matriz densidade de informação, é mostrado que a resolução teórica dos dados pode ser estudada, que guia o planejamento da prospecção. A otimização do levantamento permite determinar os períodos e as posições das estações de medida mais adequados ao delineamento mais preciso de corpos cujas localizações são conhecidas aproximadamente.

Desenvolvimento de Interface gráfica como suporte para soluções numéricas das equações de Maxwell em coordenadas gerais – 3D

Neste trabalho, é implementada uma interface gráfica de usuários (GUI) usando a ferramenta Qt da Nokia (versão 3.0). A interface visa simplificar a criação de cenários para a realização de simulações paralelas usando a técnica numérica Local Nonorthogonal Finite Difference Time-Domain (LN-FDTD), aplicada para solucionar as equações de Maxwell. O simulador foi desenvolvido usando a linguagem de programação C e paralelizado utilizando threads. Para isto, a biblioteca pthread foi empregada. A visualização 3D do cenário a ser simulado (e da malha) é realizada por um programa especialmente desenvolvido que utiliza a biblioteca OpenGL. Para melhorar o desenvolvimento e alcançar os objetivos do projeto computacional, foram utilizados conceitos da Engenharia de Software, tais como o modelo de processo de software por prototipagem. Ao privar o usuário de interagir diretamente com o código-fonte da simulação, a probabilidade de ocorrência de erros humanos durante o processo de construção de cenários é minimizada. Para demonstrar o funcionamento da ferramenta desenvolvida, foi realizado um estudo relativo ao efeito de flechas em linhas de baixa tensão nas tensões transitórias induzidas nas mesmas por descargas atmosféricas. As tensões induzidas nas tomadas da edificação também são estudadas.

Análise de um regulador de velocidade digital para a Usina Hidroelétrica de Curuá-Una

Analisa uma proposta de um controlador digital para o regulador de velocidade da Usina Hidroelétrica de Curuá-Una, utilizando duas estratégias de projeto, a estratégia indireta com o projeto do controlador sendo feito no plano contínuo e depois discretizado e a estratégia direta, em que o controlador é projetado inteiramente no plano discreto. O conhecimento gerado no simulador é de grande importância, pois é uma das grandes ferramentas para avaliar o comportamento dos controladores digitais propostos, em ambiente seguro. O estudo dos controladores também permite uma substituição de equipamentos antigos, com produção descontinuada, por equipamentos novos que permitem a confecção de controladores modernos, digitais e inteligentes, que proporcionam uma substituição de técnicas antigas de controle por estratégias avançadas de controle, maximizando o rendimento do sistema em condições adversas de operação. O trabalho levanta as diferenças, vantagens e desvantagens de cada controlador, com o objetivo de auxiliar na escolha do controlador mais adequado para projetos de controladores avançados: tipo um controle adaptativo, controle fuzzy ou controle neural. As dificuldades de projeto e os resultados das simulações foram os principais indicadores na avaliações dos dois RV digitais projetados. De uma forma geral o RV digital direto apresentou melhor desempenho, estabilidade e menor esforço computacional, entretanto, o RV digital indireto mostrou desempenho similar, menor degradação devido às não linearidades e ao menor esforço de projeto.

Projeto e implementação em DSP de sistema de comunicação multiportadora baseado no padrão G.fast

As redes de acesso usando cabos de cobre atualmente utilizam bandas de frequência até 30 MHz, especificada no padrão VDSL2. À medida que arquiteturas híbridas de fibra e cobre se tornam mais proeminentes na indústria e academia, torna-se possível utilizar cabos metálicos mais curtos (i.e. até 250 metros) conectando o último ponto de distribuição aos usuários, de modo que frequências mais altas podem ser exploradas para se alcançar taxas de transmissão de dados de 500 Mbps ou mais, como é o caso do padrão G.fast atualmente em desenvolvimento no ITU-T. Nesse trabalho, um simulador no domínio do tempo foi desenvolvido para avaliar a capacidade do sistema G.fast com diferentes tamanhos de extensão cíclica e diferentes topologias de rede especificadas pelo ITU-T. Os resultados das simulações mostram que sistemas G.fast são robustos a bridged taps e capazes de atingir altas taxas de dados para todas as topologias simuladas, provendo suporte à próxima geração de serviços de banda larga. Além disso, esse trabalho descreve o progresso da implementação de um protótipo de modem baseado no padrão G.fast em um ambiente híbrido de DSP multicore e FPGA utilizando kits de avaliação adquiridos pela UFPA. Arquiteturas, protocolos de comunicação e benchmarks são apresentados e avaliados para se chegar à conclusão de que tal protótipo é factível e fornece suporte flexível a várias linhas de pesquisa em banda larga da próxima geração.

Imageamento homeomórfico de refletores sísmicos

Neste trabalho é apresentada uma nova técnica para a realização do empilhamento sísmico, aplicada ao problema do imageamento de refletores fixos em um meio bidimensional, suavemente heterogêneo, isotrópico, a partir de dados de reflexão. Esta nova técnica chamada de imageamento homeomórfico tem como base a aproximação geométrica do raio e propriedades topológicas dos refletores. São utilizados, portanto, os conceitos de frente de onda, ângulo de incidência, raio de curvatura da frente de onda, cáustica e definição da trajetória do raio; de tal modo que a imagem obtida mantém relações de homeomorfismo com o objeto que se deseja imagear. O empilhamento sísmico é feito, nesta nova técnica de imageamento, aplicando-se uma correção local do tempo, ∆ t, ao tempo de trânsito, t, do raio que parte da fonte sísmica localizada em x_o, reflete-se em um ponto de reflexão, C_o, sendo registrado como uma reflexão primária em um geofone localizado em x_g, em relação ao tempo de referência t_o no sismograma, correspondente ao tempo de trânsito de um raio central. A fórmula utilizada nesta correção temporal tem como parâmetros o raio de curvatura R_o, o ângulo de emergência β_o da frente de onda, no instante em que a mesma atinge a superfície de observação, e a velocidade v_o considerada constante nas proximidades da linha sísmica. Considerando-se uma aproximação geométrica seguido um círculo para a frente de onda, pode-se estabelecer diferentes métodos de imageamento homeomórfico dependendo da configuração de processamento. Sendo assim tem-se: 1) Método Elemento de Fonte (Receptor) Comum (EF(R)C). Utiliza-se uma configuração onde se tem um conjunto de sismogramas relacionado com uma única fonte (receptor), e considera-se uma frente de onda real (de reflexão); 2) Método Elemento de Reflexão Comum (ERC). Utiliza-se uma configuração onde um conjunto de sismogramas é relacionado com um único ponto de reflexão, e considera-se uma frente de onda hipoteticamente originada neste ponto; 3) Método Elemento de Evoluta Comum (EEC). Utiliza-se uma configuração onde cada sismograma está relacionado com um par de fonte e geofone coincidentemente posicionados na linha sísmica, e considera-se uma frente de onda hipoteticamente originada no centro de curvatura do refletor. Em cada um desses métodos tem-se como resultados uma seção sísmica empilhada, u(x_o, t_o); e outras duas seções denominadas de radiusgrama, R_o (x_o, t_o), e angulograma, β_o(x_o, t_o), onde estão os valores de raios de curvatura e ângulos de emergência da frente de onda considerada no instante em que a mesma atinge a superfície de observação, respectivamente. No caso do método denominado elemento refletor comum (ERC), a seção sísmica resultante do empilhamento corresponde a seção afastamento nulo. Pode-se mostrar que o sinal sísmico não sofre efeitos de alongamento como consequência da correção temporal, nem tão pouco apresenta problemas de dispersão de pontos de reflexão como consequência da inclinação do refletor, ao contrário do que acontece com as técnicas de empilhamento que tem por base a correção NMO. Além disto, por não necessitar de um macro modelo de velocidades a técnica de imageamento homeomórfico, de um modo geral, pode também ser aplicada a modelos heterogêneos, sem perder o rigor em sua formulação. Aqui também são apresentados exemplos de aplicação dos métodos elemento de fonte comum (EFC) (KEYDAR, 1993), e elemento refletor comum (ERC) (STEENTOFT, 1993), ambos os casos com dados sintéticos. No primeiro caso, (EFC), onde o empilhamento é feito tendo como referência um raio central arbitrário, pode-se observar um alto nível de exatidão no imageamento obtido, além do que é dada uma interpretação para as seções de radiusgrama e angulograma, de modo a se caracterizar aspectos geométricos do model geofísico em questão. No segundo caso, (ERC), o método é aplicado a série de dados Marmousi, gerados pelo método das diferenças finitas, e o resultado é comparado com aquele obtido por métodos convecionais (NMO/DMO) aplicados aos mesmos dados. Como consequência, observa-se que através do método ERC pode-se melhor detectar a continuidade de refletores, enquanto que através dos métodos convencionais caracterizam-se melhor a ocorrência de difrações. Por sua vez, as seções de radiusgrama e angulograma, no método (ERC), apresentam um baixo poder de resolução nas regiões do modelo onde se tem um alto grau de complexidade das estruturas. Finalmente, apresenta-se uma formulação unificada que abrange os diferentes métodos de imageamento homeomórfico citados anteriormente, e também situações mais gerais onde a frente de onda não se aproxima a um círculo, mas a uma curva quadrática qualquer.

Aplicação dos métodos gravimétrico e magnetométrico para a definição do comportamento estrutural da faixa de dobramentos Araguaia

O presente trabalho consiste na interpretação de informações gravimétricas e aeromagnetométricas (Projeto Geofísico Brasil-Canadá - PGBC) da região setentrional da Faixa de Dobramentos Araguaia, envolvendo uma área de aproximadamente 129.000km², compreendida entre os meridianos 47°50’W e 50°30’W e paralelos 4°50'S e 9°00'S. Abrange porções sudeste do Estado do Pará, noroeste de Goiás e oeste do Maranhão. Os trabalhos de campo constaram de levantamentos gravimétrico e altimétrico, ao longo de rodovias que constituem a rede viária regional, e os resultados obtidos, após correções e reduções, foram então organizados de forma a constituírem um mapa de anomalias Bouguer, apresentado em escala 1:500.000. O padrão gravimétrico da faixa de dobramentos Araguaia se caracteriza por mostrar feições predominantemente longitudinais, com curvas isoanômalas de direções submeridianas, concordante com o comportamento litológico-estrutural conhecido para a área. Na porção centro-meridional desse mapa aparece uma zona de anomalia negativa, alongada e intensa (menos de -100mgal), a qual sofre duas importantes inflexões para NW, sendo uma correlacionável ao Lineamento Carajás, e a outra na altura do paralelo 5°30'S. Apresenta zonas de altos gravimétricos, notadamente na porção noroeste (domínio cratônico), relacionada à influência das metavulcânicas do Grupo Grão Pará, e na porção nordeste, devida a massas densas introduzidas na crosta. Destaca-se também o relativo alto gravimétrico acompanhando o flanco oriental do eixo de ocorrência das braquidobras, podendo ser devido a massas excedentes colocadas sob essa região, através de esforços tectônicos, produzindo elevações do nível de base da crosta e consequentemente do embasamento. De forma genérica o flanco oriental da faixa mostra um gradiente mais intenso que o ocidental, devido principalmente à configuração geométrica das estruturas nesse setor. Da mesma forma merece destaque o gradiente regional ascendente de sul para norte. A análise das cartas de intensidade magnética do PGBC fornecem condições à individualização de diversos domínios magnéticos, perfeitamente correlacionáveis litológico e estruturalmente com as unidades geotectônicas que compõem o arcabouço regional, assim como suas subunidades. A execução de modelamento simples, quer para as informações gravimétricas como magnetométricas, objetivaram fornecer subsídios semi-quantitativos à interpretação, de forma a auxiliar no estabelecimento do padrão geométrico do embasamento da faixa de dobramentos Araguaia, assim como das estruturas a ela impostas. Assim, por meio do processamento das informações gravimétricas, verificou-se um padrão geométrico na forma de fatias imbricadas, estabelecidas sobre o substrato. A grande anomalia negativa da porção centro-meridional da faixa pode estar relacionada a uma deficiência de massa profunda (depressão na base da crosta) ou ainda a uma associação entre contraste de densidade e espessamento do pacote de supracrustais, nessa porção. O baixo gravimétrico da porção sudeste da área, entre Guarai e Tupirama, é devido a uma depressão do embasamento, da mesma forma como a anomalia negativa associada ao Lineamento Carajás está relacionada a uma estrutura sinformal conhecida.

Delineamento do pé do talude na margem continental do Ceará através da integração de dados geológicos e geofísicos

A metodologia de interpretação integrada dos dados geológicos e geofísicos observados em um perfil da margem continental do Ceará possibilita a identificação e a integração de características peculiares a cada tipo de dado. Dessa forma, é possível se definir a localização mais provável de feições estruturais importantes, tais como a fronteira entre as crostas continental e oceânica e o pé do talude, objeto do presente estudo. Segundo o Artigo 76 (parágrafo 4, item b) da Convenção das Nações Unidas para o Direito do Mar, o pé do talude é definido como o ponto de variação máxima do gradiente do talude na sua base. Entretanto, essa definição, apesar de simples no contexto fisiográfico, não é suficiente para propiciar a localização do pé do talude como preconiza a Convenção, razão pela qual se aplicam os métodos geofísicos. Dentro do contexto geofísico-geológico está implícita a interpretação quantitativa das anomalias gravimétricas ar-livre, que possibilita o delineamento do modelo geofísico representando a subsuperfície, cuja finalidade é subsidiar geologicamente a interpretação integrada dos dados aludidos. Um procedimento automático de ajuste de curvas combinando as técnicas de inversão de busca sistemática e a que utiliza derivadas foi usado com o propósito de gerar o modelo geofísico. A aplicação rigorosa de vínculos preliminarmente e a constante reavaliação desses vínculos através de um processo interativo entre a sísmica e a gravimetria, gerado durante a interpretação quantitativa das anomalias ar-livre, possibilitaram que o modelo geofísico final estivesse dentro dos padrões geológicos para área, notadamente quanto ao equilíbrio isostático (Teoria de Airy). O objetivo do presente trabalho é se estudar as características geológicas e geofísicas observadas ao longo de um perfil da margem continental do Ceará (LEPLAC III), notadamente quanto ao pé do talude, buscando estabelecer a aplicabilidade de uma metodologia de interpretação integrada desses dados, cuja finalidade é se definir de forma sistemática a localização mais provável para esta feição fisiográfica. A metodologia de interpretação integrada dos dados geológicos e geofísicos empregada mostrou-se eficiente para este objetivo. Foi possível se integrar: (i) a localização fisiográfica (distância da costa e profundidade) do pé do talude; (ii) a zona de instabilidade tectônica evidenciada pelos falhamentos, comuns nesta região; (iii) ao fim de uma zona magnética perturbada, associado a um ponto de mínimo na curva de anomalia magnética, e que possivelmente delimita o início de uma zona magnética quieta, denominada de anomalia E e (iv) um ponto de inflexão na curva de anomalia ar-livre, associado ao efeito gravimétrico do contraste de densidades entre as crostas continental, os sedimentos e a água do mar, evidenciado pela geometria do talude. Foi possível ainda se definir a localização mais provável para a fronteira entre as crostas continental e oceânica. Dada a rigorosidade na aplicação das técnicas de inversão e dos vínculos é provável que as correlações das características intrínsecas a cada tipo de dado efetuadas na conclusão desse trabalho tenham fundamento e possam ser confirmadas. A condição para isto é a aplicação da metodologia aqui estabelecida em um número maior de perfis.

Detecção de refletores sísmicos por rede neural discreta

As redes neurais artificiais têm provado serem uma poderosa técnica na resolução de uma grande variedade de problemas de otimização. Nesta dissertação é desenvolvida uma nova rede neural, tipo recorrente, sem realimentação (self-feedback loops) e sem neurônios ocultos, para o processamento do sinal sísmico, para fornecer a posição temporal, a polaridade e as amplitudes estimadas dos refletores sísmicos, representadas pelos seus coeficientes de reflexão. A principal característica dessa nova rede neural consiste no tipo de função de ativação utilizada, a qual permite três possíveis estados para o neurônio. Busca-se estimar a posição dos refletores sísmicos e reproduzir as verdadeiras polaridades desses refletores. A idéia básica desse novo tipo de rede, aqui denominada rede neural discreta (RND), é relacionar uma função objeto, que descreve o problema geofísico, com a função de Liapunov, que descreve a dinâmica da rede neural. Deste modo, a dinâmica da rede leva a uma minimização local da sua função de Liapunov e consequentemente leva a uma minimização da função objeto. Assim, com uma codificação conveniente do sinal de saída da rede tem-se uma solução do problema geofísico. A avaliação operacional da arquitetura desta rede neural artificial é realizada em dados sintéticos gerados através do modelo convolucional simples e da teoria do raio. A razão é para explicar o comportamento da rede com dados contaminados por ruído, e diante de pulsos fonte de fases mínima, máxima e misturada.