Desenvolvimento de um magnetometro a precessão nuclear

O objetivo desta tese foi desenvolver um magnetômetro à precessão nuclear para prospecção geofísica e estações-base magnéticas. O magnetômetro à precessão nuclear mede a intensidade total do campo magnético. Seu funcionamento é baseado na ressonância magnética nuclear. A medida de campo é feita pela de terminação da freqüência de precessão de núcleos de hidrogênio – prótons - de líquidos não viscosos no campo magnético terrestre. O magnetômetro é constituído de duas partes: o sensor e o instrumento de medida. O sensor é uma bobina solenoidal, cujo núcleo é preenchido com o líquido. Três líquidos diferentes foram testados; água, propanol e um querosene sintético. Optou-se pelo uso do querosene porque oferece maior amplitude no sinal de precessão, dando, conseqüentemente, maior relação sinal/ ruído. O sistema de medida contém os circuitos de sintonia e amplificação do sinal e, os circuitos lógicos para a programação da operação e contagem da freqüência de precessão. Cada ciclo de medida tem duração de 3 segundos, sendo 2,3s para a polarização e 0,7s para a recepção do sinal. São possíveis dois modos de operação: manual, reciclando automaticamente e por controle remoto. O sinal de precessão é amplificado seletivamente em uma das 14 faixas de sintonia, que cobrem medidas entre 22000 e 95000 gammas. A freqüência de precessão é multiplicada por um fator de 64 e contada durante um tempo igual a 0,36699s, determinado com base na razão giromagnética do próton. O número de pulsos contados é numericamente igual ao valor do campo magnético em gammas. A resposta pode ser lida em mostradores digitais ou na saída BCD paralela quando operando por controle remoto. A precisão da medida é de 1 gamma. O instrumento foi testado no campo para avaliar a relação sinal/ruído, gradiente suportável e consumo de potência. Nos testes de aplicação do protótipo, foram obtidos dados de variação diurna e realizaram-se levantamentos magnético de reconhecimento e detalhe em um sítio arqueológico na Ilha de Marajó, Pará. As respostas dos testes foram comparados com dois magnetômetros comerciais - o GP-70, McPhar e o G-816, Geometrics e, ainda, com dados do Observatório Magnético de Tatuoca-Pa. Em todos os casos, a comparação dos dados mostrou bom desempenho do magnetômetro em teste.

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.

Efeitos do manto de intemperismo sobre anomalias eletromagnéticas provocadas por corpos tabulares inclinados: um estudo através de modelamento reduzido

Electromagnetic methods have been extensively applied in the prospecting of sulphide bodies and other conducting materials. The interpretation of e. m. data is based on the results obtained either with reduced scale or analytical modelling. In most models, the host rocks, the overburden as well as the halo of disseminated sulphides are considered highly resistive although in nature they are often weakly conductive. Presence of a well-developed conducting overburden in tropical and sub-tropical regions and a saline rich crustal layer in semi-arid regions have been found to modify significantly the e. m. anomalies in practice. Therefore, the parameters of the target, determined on the basis of simple models where the presence of the conducting environment is neglected, are found to be in considerable error. The effects of the overburden on the e. m. anomalies of a dipping tabular ore body were studied with reduced scale models for varying response parameters of both the overburden and the orebody, and also for different depths and dips of the target. The overburden and the orebody were represented by metallic sheets of varying thickness in the scale model developed in accordance with the law of electromagnetic similitude. The results of these investigations show that the overburden affects the anomaly by causing: a) phase rotation; b) amplitude reduction; c) base level displacement; d) reversion of the quadrature; and e) appearance of an extra peak in the quadrature in the case of low dipping models. The last two effects complicate the quantification of the anomalies but, on the other hand, provide a qualitative indication of the response parameters of both the ore body and the overburden. The results were assembled in the form of Argand diagrams and, finally, an interpretation scheme is suggested for the e. m. field data on the basis of these diagrams.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Empilhamento sísmico pela composição de ondas planas

No presente trabalho de tese é apresentada uma nova técnica de empilhamento de dados sísmicos para a obtenção da seção de incidência normal ou afastamento fonte-receptor nulo, aplicável em meios bidimensionais com variações laterais de velocidade. Esta nova técnica denominada Empilhamento Sísmico pela Composição de Ondas Planas (empilhamento PWC) foi desenvolvida tomando como base os conceitos físicos e matemáticos da decomposição do campo de ondas em ondas planas. Este trabalho pode ser dividido em três partes: Uma primeira parte, onde se apresenta uma revisão da técnica de empilhamento sísmico convencional e do processo de decomposição do campo de ondas produzido a partir de fontes pontuais em suas correspondentes ondas planas. Na segunda parte, é apresentada a formulação matemática e o procedimento de aplicação do método de empilhamento sísmico pela composição de ondas planas. Na terceira parte se apresenta a aplicação desta nova técnica de empilhamento na serie de dados Marmousi e uma analise sobre a atenuação de ruído. A formulação matemática desta nova técnica de empilhamento sísmico foi desenvolvida com base na teoria do espalhamento aplicado a ondas sísmicas sob a restrição do modelo de aproximação de Born. Nesse sentido, inicialmente se apresenta a determinação da solução da equação de onda caustica para a configuração com afastamento fonte-receptor finito, que posteriormente é reduzido para a configuração de afastamento fonte-receptor nulo. Por outra parte, com base nessas soluções, a expressão matemática deste novo processo de empilhamento sísmico é resolvida dentro do contexto do modelo de aproximação de Born. Verificou-se que as soluções encontradas por ambos procedimentos, isto é, por meio da solução da equação da onda e pelo processo de empilhamento proposto, são iguais, mostrando-se assim que o processo de empilhamento pela composição de ondas planas produz uma seção com afastamento fonte-receptor nulo. Esta nova técnica de empilhamento basicamente consiste na aplicação de uma dupla decomposição do campo de ondas em onda planas por meio da aplicação de dois empilhamentos oblíquos (slant stack), isto é um ao longo do arranjo das fontes e outro ao longo do arranjo dos detectores; seguido pelo processo de composição das ondas planas por meio do empilhamento obliquo inverso. Portanto, com base nestas operações e com a ajuda de um exemplo de aplicação nos dados gerados a partir de um modelo simples, são descritos os fundamentos e o procedimento de aplicação (ou algoritmo) desta nova técnica de obtenção da seção de afastamento nulo. Como exemplo de aplicação do empilhamento PWC em dados correspondentes a um meio com variações laterais de velocidade, foi aplicado nos dados Marmousi gerados segundo a técnica de cobertura múltipla a partir de um modelo que representa uma situação geológica real. Por comparação da seção resultante com a similar produzida pelo método de empilhamento convencional, observa-se que a seção de afastamento nulo desta nova técnica apresenta melhor definição e continuidade dos reflectores, como também uma melhor caracterização da ocorrência de difrações. Por último, da atenuação de ruído aleatório realizada nos mesmos dados, observa-se que esta técnica de empilhamento também produz uma atenuação do ruído presente no sinal, a qual implica um aumento na relação sinal ruído.

Equipamento eletromagnético para prospecção geofísica e modelamento reduzido

Nas últimas décadas, a prospecção por métodos eletromagnéticos vem-se constituindo numa técnica eficiente para prospecção mineral. O objetivo deste trabalho foi desenvolver um equipamento para prospecção eletromagnética quantitativa de corpos condutores, através do método dipolo-dipolo, podendo ainda ser usado em modelos reduzidos. Eletricamente, o sistema mede grandezas relacionadas ao acoplamento indutivo entre duas bobinas: transmissora e receptora. Elas são dispostas na superfície da terra, afastadas entre si, e a terra, desse modo, constitui o núcleo acoplador. Quando existem corpos condutores nas proximidades, estes são denunciados por alterações no comportamento do sinal induzido na bobina receptora. O equipamento compreende dois conjuntos: o transmissor e o receptor, além de acessórios. O transmissor gera um campo eletromagnético nas freqüências de 520 e 3.090 Hz, e um sinal de referência para o receptor, o qual é enviado através de um cabo. O receptor, inicialmente, separa o sinal induzido pelos campos secundários gerados por condutores, do campo normalmente recebido, quando a condutividade da subsuperfície é relativamente uniforme (campo primário). Em seguida, decompõe esse sinal em duas componentes ortogonais, uma em fase, e outra em quadratura com o campo primário. Através de duas escalas de precisão, as amplitudes dessas componentes são mostradas como percentagens do campo primário, com precisão de 1%. A sensibilidade do receptor é de 0,5 μV. O circuito eletrônico foi rigorosamente testado com preciso instrumental de laboratório. Em seguida, testou-se sua aplicação no Laboratório de Modelo Reduzido Eletromagnético do NCGG, refazendo-se experiências clássicas, encontradas na literatura especializada. No campo, foi experimentado próximo da cidade de Araci no Estado da Bahia, em áreas prospectadas pela "Rio Doce Geologia e Mineração S/A-DOCEGEO". Em ambos os casos, verificaram-se bons resultados.

Espalhamento geométrico em modelos plano-estratificados

A medição de parâmetros físicos de reservatórios se constitui de grande importância para a detecção de hidrocarbonetos. A obtenção destes parâmetros é realizado através de análise de amplitude com a determinação dos coeficientes de reflexão. Para isto, faz-se necessário a aplicação de técnicas especiais de processamento capazes de corrigir efeitos de divergência esférica. Um problema pode ser estabelecido através da seguinte questão: Qual o efeito relativamente mais importante como responsável pela atenuação de amplitudes, o espalhamento geométrico ou a perda por transmissividade? A justificativa desta pergunta reside em que a correção dinâmica teórica aplicada a dados reais visa exclusivamente o espalhamento geométrico. No entanto, a análise física do problema por diferentes direções põe a resposta em condições de dúvida, o que é interessante e contraditório com a prática. Uma resposta embasada mais fisicamente pode dar melhor subsídio a outros trabalhos em andamento. O presente trabalho visa o cálculo da divergência esférica segundo a teoria Newman-Gutenberg e corrigir sismogramas sintéticos calculados pelo método da refletividade. O modelo-teste é crostal para que se possa ter eventos de refração crítica além das reflexões e para, com isto, melhor orientar quanto à janela de aplicação da correção de divergência esférica o que resulta em obter o então denominado “verdadeiras amplitudes”. O meio simulado é formado por camadas plano-horizontais, homogêneas e isotrópicas. O método da refletividade é uma forma de solução da equação de onda para o referido modelo, o que torna possível um entendimento do problema em estudo. Para se chegar aos resultados obtidos foram calculados sismogramas sintéticos através do programa P-SV-SH desenvolvido por Sandmeier (1998), e curvas do espalhamento geométrico em função do tempo para o modelo estudado como descrito por Newman (1973). Demonstramos como uma das conclusões que a partir dos dados do modelo (velocidades, espessuras, densidades e profundidades) uma equação para a correção de espalhamento geométrico visando às “verdadeiras amplitudes” não é de fácil obtenção. O objetivo maior então deveria ser obter um painel da função de divergência esférica para corrigir as verdadeiras amplitudes.

Estudo de águas subterrâneas com métodos elétricos na região do semi-árido de Patos - Pb

O Nordeste do Brasil, em particular a região semi-árida do Estado da Paraíba, é muito castigada pelos períodos de grandes estiagens e escassez de recursos hídricos. A água de superfície não supre a demanda da região, portanto, a busca de água subterrânea e sua exploração é uma constante nos programas dos Governos Federal e Estadual. Nesta região, ocorrem geralmente apenas aquíferos aluviais rasos e aquíferos nas fraturas do embasamento cristalino. Tradicionalmente, as fraturas no cristalino são prospectados através de estudos aerofotogramétricos e geoestruturais, com média de 30% de furos secos. O dimensionamento dos aquíferos aluviais é feito através de perfurações sistemática do subsolo, que é uma metodologia muito demorosa e cara. O objetivo deste trabalho é mostrar que o emprego de métodos geofísicos, especialmente os métodos geoelétricos, são eficazes e econômicos para localizar e avaliar áreas promissoras de água subterrânea, tanto nos aquíferos aluviais como nos cristalinos. Neste trabalho, foram aplicados o convencional método da eletroresistividade e o método eletromagnético, ainda pouco usado na prospecção de água subterrânea, em três áreas selecionadas na região semi-árida de Patos - Pb, sendo duas representantes dos aquíferos aluviais e uma representante dos aquíferos cristalinos. Nos aquíferos aluviais, foi determinada a profundidade até o embasamento cristalino através de sondagens elétricas verticais (SEVs), e os respectivos limites laterais foram mapeados com perfilagens de eletroresistividade. Os resultados das SEVs são coerentes com as profundidades obtidas pelas sondagens mecânicas de confirmação mostrando, assim, a aplicabilidade desta metodologia na pesquisa deste tipo de aquífero. No aquífero cristalino, inicialmente determinou-se a orientação geral das zonas fraturadas do embasamento cristalino pela técnica modificada da sondagem elétrica vertical (SEV Radial) e, em seguida, a localização exata em planta das mesmas foram determinadas com perfilagens de eletroresistividade e eletromagnéticas no sistema horizontal coplanar nas frequências 880 Hz e 2640 Hz. Os resultados dos dois métodos utilizados foram satisfatórios, porém, a perfilagem eletromagnética é mais rápida, necessita de menos mão de obra, reduzindo assim os custos de prospecção de água subterrânea no embasamento cristalino.

Estudo de anomalias eletromagnéticas de um condutor tabular vertical sob camadas parcialmente condutiva em multifrequência e multiseparação através de modelamento analógico

A resposta eletromagnética (EM) de um corpo condutivo envolvido por uma zona parcialmente condutiva, torna-se bastante diferente daquela de um corpo condutivo em um meio altamente resistivo. As zonas parcialmente condutivas, como por exemplo, rocha encaixante, halo de sulfetos disseminados ou manto de intemperismo, que envolvem o corpo condutivo, afetam a resposta EM de diferentes maneiras, dependendo de suas características físicas e geométricas e, em particular, do sistema de prospecção EM utilizado. Neste trabalho em modelamento analógico, foi feita uma análise de anomalias EM provocadas por corpos condutivos tabulares verticais sob manto de intemperismo, em levantamentos terrestres para diferentes sistemas de bobinas - horizontal coplanar, vertical coplanar e vertical coaxial - em oito frequências na faixa de 250 Hz a 35 kHz e separações entre as bobinas de 0,15; 0,20 e 0,25m. O manto de intemperismo foi simulado por folhas de aço finas dispostas horizontalmente e o corpo condutor principal por folhas de alumínio finas colocadas verticalmente. As dimensões das folhas foram determinadas de acordo com as condições de modelamento para o plano e o semi-plano. Foram utilizados três corpos e três mantos com diferentes espessuras e condutividades, simulando, deste modo, diversas situações geológicas. Os resultados mostraram que cada sistema de bobinas é afetado diferentemente pela presença do manto de intemperismo. Para a análise dos resultados foi plotado um conjunto de diagramas considerando os valores pico-a-pico das anomalias em fase e em quadratura. Um outro conjunto de diagramas mostra as amplitudes máximas em fase, que ocorrem quando a componente em quadratura se anula em uma frequência relativamente baixa para um conjunto de corpo-manto, e as amplitudes máximas em quadratura, que ocorrem quando a resposta em fase atinge um mínimo próximo de zero, em frequências relativamente altas. Com isto foi possível conhecer a faixa de frequências para cada sistema de bobinas, onde a resposta EM se encontra o mínimo afetada pela presença do manto de intemperismo. A maior amplitude na resposta é obtida no sistema horizontal coplanar e a menor no sistema vertical coplanar. Um aumento na separação entre as bobinas é acompanhado por um deslocamento da anomalia para baixas frequências. A faixa de frequências, onde a presença do manto tem pouca influência na resposta do corpo condutivo, e maior para o sistema vertical coaxial e menor para o sistema horizontal coplanar. Esses resultados dão uma luz para o conhecimento da posição e da largura da banda de frequências utilizável, assim como as melhores separações entre transmissor-receptor, para auxiliar no planejamento de sistemas de prospecção EM, de modo que a resposta fique o mais livre possível de sinais indesejáveis, tais como os causados pela presença do manto de intemperismo.

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 da estrutura geo-elétrica da Região do Juruá, AM, pelo método magnetotelúrico

Utilizando-se dados magnetotelúricos (MT), foi obtida uma imagem geo-elétrica nítida da região do Juruá, Bacia do Solimões, na forma de seções geo-elétricas. Os dados de campo foram registrados ao longo de três linhas de 15 km, espaçadas de 3.5 km, recobrindo uma área de 100 km². O espaçamento entre as 35 estações é irregular, variando de 400 m a 3500 m. A faixa de freqüências utilizada cobriu de 0.001 Hz até 300 Hz, o que permitiu investigar de 100 m até 60 km de profundidade. Os dados apresentam-se afetados pelo efeito de distorção estática. Para corrigir este efeito foi utilizada a mediana da resistividade do primeiro condutor, correspondente à Formação Solimões. Foi utilizado o invariante do tensor MT para interpretar a estrutura geo-elétrica do Juruá. As seções geo-elétricas foram obtidas a partir do agrupamento dos dados resultantes da transformação de Bostick e da inversão 1D de Occam, para cada estação. Foi identificada uma seqüência de camadas condutivas e resistivas, correspondentes ao pacote sedimentar, uma zona de falhas e o topo do embasamento geo-elétrico, caracterizando a Bacia do Solimões. Abaixo do embasamento geo-elétrico foram também identificados uma zona condutora, seguida por uma camada de baixa condutividade, a profundidades iguais ou superiores a 20 km. Esta camada é interpretada como sendo de composição de gabro, estando associada a processos de acreção vertical, intimamente ligados à estabilização crustal e espessamento da litosfera. Os resultados apresentam uma boa concordância com os perfis de resistividade de poços e dados sísmicos de superfície.

Estudo de geração de campo uniforme para modelagem eletromagnética

Para estudar os problemas de prospecção geofísica eletromagnética através de modelagem analógica, as condições em escala natural são representadas, no laboratório, em escala reduzida de acordo com a teoria da similitude. Portanto, para investigar os problemas de técnicas VLF, AFMAG e MT, frequentemente é necessário criar campo uniforme no arranjo experimental. Vários sistemas físicos para geração de campos uniformes são analisados teoricamente nesta tese. Os sistemas estudados aqui são a bobina circular, bobina de Helmholtz, solenóide, um plano de corrente e dois planos paralelos de correntes. As equações analíticas foram obtidas para campo magnético num ponto do espaço e subsequentemente as condições de campo uniforme. Nos casos em que as condições para o campo uniforme não puderam ser obtidas analiticamente, a porcentagem de desvio do campo em relação a um ponto pré-selecionado foi calculada. Contudo, os mapas de campo magnético, assim como o mapa de porcentagem de desvio, estão presentes para todos os sistemas estudados aqui. Também, foram calculados as áreas e os volumes espaciais de vários desvios de porcentagem do campo uniforme. Um estudo comparativo desses sistemas mostra que o solenóide é a maneira mais eficiente para criar um campo uniforme, seguido pelo sistema de bobinas de Helmholtz. Porém, o campo criado em um solenóide está em um espaço fechado onde é difícil colocar modelos e substituí-los para executar experimentos. Portanto, recomenda-se o uso de bobinas de Helmholtz para criar um campo uniforme. Este último sistema fornece campo uniforme com espaço aberto suficiente, o que facilita o experimento.

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.

Estudo do método de inversão tomográfica 2D baseada nos atributos cinemáticos da onda NIP e construção de interface gráfica para o processo de picking

Neste trabalho, estuda-se um novo método de inversão tomográfica de reflexão para a determinação de um modelo isotrópico e suave de velocidade por meio da aplicação, em dados sintéticos e reais, do programa Niptomo que é uma implementação do método de inversão tomográfica dos atributos cinemáticos da onda hipotética do ponto de incidência normal (PIN). Os dados de entrada para a inversão tomográfica, isto é, o tempo de trânsito e os atributos da onda PIN (raio de curvatura da frente de onda emergente e ângulo de emergência), são retirados de uma série de pontos escolhidos na seção afastamento nulo (AN) simulada, obtida pelo método de empilhamento por superfícies de reflexão comum (SRC). Normalmente, a escolha destes pontos na seção AN é realizada utilizando-se programas de picking automático, que identificam eventos localmente coerentes na seção sísmica com base nos parâmetros fornecidos pelo usuário. O picking é um dos processos mais críticos dos métodos de inversão tomográfica, pois a inclusão de dados de eventos que não sejam de reflexões primárias podem ser incluídos neste processo, prejudicando assim o modelo de velocidades a ser obtido pela inversão tomográfica. Este trabalho tem por objetivo de construir um programa de picking interativo para fornecer ao usuário o controle da escolha dos pontos de reflexões sísmicas primárias, cujos dados serão utilizados na inversão tomográfica. Os processos de picking e inversão tomográfica são aplicados nos dados sintéticos Marmousi e nos dados da linha sísmica 50-RL-90 da Bacia do Tacutu. Os resultados obtidos mostraram que o picking interativo para a escolha de pontos sobre eventos de reflexões primárias favorece na obtenção de um modelo de velocidade mais preciso.

Gradiometria magnética e GPR aplicados à arqueologia em sítio estuarino em Penalva - MA

Esta pesquisa apresenta os resultados da aplicação de métodos geofísicos não invasivos, visando a indicação de possíveis locais de deposição de materiais, no levantamento arqueológicos em um sítio estuarino localizado no município de Penalva, no estado do Maranhão. Os métodos utilizados foram Magnetometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais para escavação. O GPR foi utilizado para confirmar as anomalias geradas pelo método magnético, proporcionando a diminuição de erros quando as anomalias magnéticas não são provenientes de feições arqueológicas. Foram realizados 14 perfis de GPR, utilizando uma antena de 400 MHz e 1664 medidas magnéticas, utilizando o magnetômetro. Durante o levantamento arqueológico realizado, baseando-se nas anomalias fornecidas pela geofísica, a equipe liderada pelos arqueólogos Deusdédit Carneiro Leite Filho e Fernando Luiz Tavares Marques realizou cinco escavações no sítio, o que lhes permitiu o resgate de várias peças cerâmicas. Os resultados apresentados nesta pesquisa permitiram verificar e avaliar a potencialidade da prospecção geofísica como método auxiliar em estudos arqueológicos.