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.

Modelagem sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

A sucessão siliciclástica-carbonática neocarbonífera da bacia do Amazonas, regiões de Monte Alegre e Itaituba (PA)

Registros sedimentares do Neocarbonífero, particularmente do Moscoviano, na Bacia do Amazonas, Norte do Brasil, caracterizam a zona de contato entre as formações Monte Alegre (rochas siliciclásticas) e Itaituba (rochas carbonáticas). A análise faciológica da sucessão Moscoviana de até 40 m de espessura, exposta na região de Monte Alegre e Itaituba, Estado do Pará, permitiu identificar 5 associações de fácies (AF), que correspondem a depósitos estabelecidos no ambiente costeiro, representados por dunas/intedunas eólicas (AF1), lençóis de areia/wadi (AF2), laguna/washover (AF3), praia/planície de maré (AF4) e laguna/delta de maré (AF5). A associação de campo de dunas/interdunas (AF1) é constituída por arenitos finos a médios, bimodais com estratificação cruzada de médio porte, laminação cavalgante transladante subcrítica e arenitos com gradação inversa. Arenitos finos com acamamento maciço, marcas de raízes e, subordinadamente, verrugas de aderência (adhesion warts), ocorrem nos limites dos sets de estratificação cruzada e indicam, respectivamente, paleossolos e migração de grãos por ação eólica sobre interduna úmida. Depósitos de lençóis de areia/wadi (AF2) são compostos de arenitos finos a médios com estratificação plano-paralela e laminação cavalgante transladante subcrítica, relacionados a superfícies de deflação (lençóis de areia), enquanto arenitos finos a médios com estratificações cruzadas tangencial e recumbente, e acamamento convoluto, caracterizam rios efêmeros com alta energia. Pelitos laminados e arenitos finos com laminação cruzada cavalgante, contendo o icnofóssil Palaeophycus, representam sedimentação em ambiente de baixa energia e foram agrupados na associação de laguna/washover (AF3). Os depósitos de praia/planície de maré (AF4) consistem em arenitos finos a médios, com estratificação plano-paralela a cruzada de baixo ângulo, intercalados com lentes de dolomito fino maciço, localmente truncados por arenitos finos a médios. Estas fácies foram formadas pelo fluxo-refluxo em ambiente de praia, localmente retrabalhadas por pequenos canais, enquanto o carbonato é interpretado como precipitado em poças (ponds). Na AF4 encontram-se também pelitos laminados com gretas de contração, lâminas curvadas de argila e arenitos com estratificação cruzada tabular de pequeno a médio porte, contendo filmes de argila sobre foresets e superfícies de reativação, sugerindo a migração de sandwaves na intermaré. A associação de laguna/delta de maré (AF5) é constituída por calcários dolomitizados (mudstones, wackestones, packstones e grainstones) com poros do tipo vug e móldicos e bioclastos de braquiópodes, equinodermas, foraminíferos, ostracodes, briozoários, trilobitas, moluscos e coral isolado não fragmentado, além do ichnofóssil Thalassinoides. Conglomerados com seixos de calcário dolomitizado, arenitos finos com estratificação cruzada de baixo ângulo e superfícies de reativação, localmente sobrepostos por arenitos finos com estratificação cruzada sigmoidal e laminação cruzada cavalgante, foram interpretados como depósitos de tidal inlet e delta de maré. As associações de fácies/microfácies e os dados paleontológicos descritos neste trabalho corroboram a predominância de ambientes lagunares, em parte, conectados a um ambiente desértico costeiro para o intervalo de transição entre as formações Monte Alegre e Itaituba. A abundância de grãos arredondados de areia fina nas fácies carbonáticas corrobora influxo siliciclástico advindo do ambiente desértico adjacente ao ambiente costeiro. Condições mais quentes e tropicais para a sucessão estudada são também indicadas pela presença de carbonatos e argilominerais como illita e, principalmente, esmectita, bem como uma fauna diversificada. Os litotipos siliciclásticos e carbonáticos intercalados que caracterizam a fase final da deposição Monte Alegre e o início da sedimentação Itaituba, justificam sua representação em um mesmo sistema deposicional costeiro.

Análise morfométrica do sistema auditivo periférico da preguiça (Bradypus variegatus)

A superordem Xenarthra é composta de 31 espécies viventes de tatus, tamanduás e preguiças. As arborícolas preguiças pertencem a dois gêneros, Choloepus e Bradypus, cuja divergência se deu a aproximadamente 40 milhões de anos atrás. As similaridades entre os dois taxa, tais como a presença de algas verdes nos pêlos e habilidade locomotora suspensória, são notáveis exemplos de evolução convergente. A exata posição da linhagem Xenarthra entre os mamíferos na árvore filogenética ainda não é completamente compreendida, com alguns rearranjos na árvore da família dos mamíferos placentários, considerando os Xenarthras mais relacionados entre Afrotheria (que inclui musaranhos, porcos-da-terra, peixes-boi e elefantes) ou a Boreoeutheria (que inclui primatas, roedores, carnívoros e ungulados). O objetivo deste trabalho é descrever pela primeira vez características morfológicas dos ouvidos médio e interno de Bradypus variegatus e compará-las a outros mamíferos placentários que possuam dados publicados na literatura. Nós usamos 13 espécimes adultas post-mortem (machos e fêmeas) e 15 crânios da coleção do Museu Paraense Emílio Goeldi. Além das medições, foram usadas técnicas de microscopia óptica, microscopia eletrônica de varredura e tomografia computadorizada. Através da árvore filogenética das preguiças, o gênero Bradypus é posicionado como táxon-irmão de todas as outras preguiças. Nossos resultados mostram que a morfologia do ouvido médio e interno de Bradypus variegatus é similar a de outros mamíferos com dados publicados na literatura e que apresentam escalonamento alométrico.

Qualidades intrínsecas para inovação de produtos em piso e revestimento de cerâmica vermelha do APL de São Miguel do Guamá

A presente dissertação representa uma contribuição na resolução de problema de inovação de produtos na indústria de cerâmica vermelha no Município de São Miguel do Guamá, no Estado do Pará. Apresenta-se uma metodologia, baseada na estratégia do QFD (Desdobramento da Função Qualidade), que associa os requisitos dos consumidores na elaboração dos projetos do produto e dos processos de produção. A falta de conhecimento dos produtos ou serviços a serem vendidos é uma das principais (e óbvias) razões do fracasso de pequenas empresas varejistas e de serviços. No atual nível de competitividade, a informação é o elemento fundamental para o sucesso empresarial e adaptação ao mercado, não só no sentido de obter e acumular estas informações, mas também de saber interpretá-las e utilizálas adequadamente e de forma criativa. A qualidade do produto é avaliada tanto no aspecto interno (especificações técnicas) quanto no aspecto externo (qualidades exigidas pelos clientes). A pesquisa mostra, além da realidade que vive o setor das indústrias de cerâmica, a importância e os atributos da qualidade definidos pelos consumidores relativo as cerâmicas vermelhas. O trabalho também mostra os resultados da aplicação da metodologia em uma das empresas do setor, mediante o desdobramentos da qualidade exigida, estabelece a qualidade planejada.

Pesca, alimentação, reprodução e crescimento do Amuré, Gobioides broussonnetii Lacepède, 1800 (Pisces : Gobiidae) no Estuário Amazônico, município de Vigia - Pará

Gobioides broussonnetii é o maior membro da família Gobiidae, e está distribuído no Atlântico Ocidental dos Estados Unidos até o Rio Grande do Sul. É abundante na foz do Rio Amazonas onde habita as águas salobras dos manguezais. Além de ser um importante elo na cadeia trófica do estuário amazônico, é um recurso pesqueiro explorado por pescadores artesanais. Estudaram-se aspectos da pesca, alimentação, reprodução e crescimento de Gobioides broussonnetii na baía de Marajó, município de Vigia, Pará, Brasil, entre setembro de 2003 e agosto de 2004. Exemplares para estudo foram adquiridos de pescadores locais e informações sobre a pesca foram obtidas através da observação de campo e entrevistas com os pescadores. No estudo da alimentação foi analisada a composição dos itens alimentares encontrados nos estômagos. A análise reprodutiva abordou aspectos do período reprodutivo e tamanho da 1a maturação gonadal. Foram estimados os parâmetros de crescimento utilizando o método indireto baseados nas análises das medidas de comprimento total de 1155 exemplares. G. broussonnetii é conhecida como amurés no estuário amazônico e é capturada no fundo lodoso em águas rasas sendo utilizada como isca na pesca comercial. De maneira geral é uma espécie fitófaga. Análises de conteúdo estomacal revelaram que algas da espécie Coscinodiscus concinus são alimentos mais consumidos. Há maior conteúdo estomacal durante a estação seca (julho a dezembro). O período entre janeiro e junho corresponde a época de atividade reprodutiva para essa espécie na área de estudo. O comprimento médio da primeira maturação gonadal foi de 23,9 cm. A razão sexual observada foi de 1:1. Os resultados obtidos neste estudo indicam que G. Broussonnetii usa o estuário amazônico como local de alimentação e reprodução. A análise macroscópica das gônadas indicou uma época de reprodução durante o inverno (janeiro a junho), com apenas um período de desova por ano entre fevereiro e abril. Foram encontrados nove grupos etários de G. btoussonnetii no estuário amazônico. Estimaram-se os valores de L∞ = 67,36 ( K=0,205 ano^-1, Ø = 3,014). A relação peso-comprimento apresentou diferenças significativas entre machos e fêmeas, com as fêmeas que pesam mais que os machos.

CRS seismic data imaging system: a case study for basib reevaluation

Este trabalho apresenta resultados práticos de uma atenção sistemática dada ao processamento e à interpretação sísmica de algumas linhas terrestres do conjunto de dados do gráben do Tacutu (Brasil), sobre os quais foram aplicadas etapas fundamentais do sistema WIT de imageamento do empilhamento CRS (Superfície de Reflexão Comum) vinculado a dados. Como resultado, esperamos estabelecer um fluxograma para a reavaliação sísmica de bacias sedimentares. Fundamentado nos atributos de frente de onda resultantes do empilhamento CRS, um macro-modelo suave de velocidades foi obtido através de inversão tomográfica. Usando este macro-modelo, foi realizado uma migração à profundidade pré- e pós-empilhamento. Além disso, outras técnicas baseadas no empilhamento CRS foram realizadas em paralelo como correção estática residual e migração de abertura-limitada baseada na zona de Fresnel projetada. Uma interpretação geológica sobre as seções empilhadas e migradas foi esboçada. A partir dos detalhes visuais dos painéis é possível interpretar desconformidades, afinamentos, um anticlinal principal falhado com conjuntos de horstes e grábens. Também, uma parte da linha selecionada precisa de processamento mais detalhado para evidenciar melhor qualquer estrutura presente na subsuperfície.

Filtros ótimos: efeitos de janelas e de tipo de fase

O objetivo central deste trabalho é o estudo do desempenho do operador deconvolucional WHL na compressão do pulso-fonte sísmico, sob condições especiais de fase não-mínima e da densidade de eventos no traço, como casos advogados para dados reais e processamento em rotina. O método de ataque ao problema construído é centrado no conteúdo da informação da função autocorrelação submetida a diferentes condições: (a) de truncamento e tipo de janelas; (b) das características da fase do operador (se mínima ou não-mínima); (c) da medida de qualidade; (d) do nível de embranquecimento; (e) do ruído presente e da equalização; (f) do balanceamento do traço; (g) dos princípios físicos da propagação expressos e limitados pelo modelo convolutional. Os resultados obtidos são apenas na forma numérica, organizados na forma de álbuns com dificuldades crescentes, e demonstram como o uso de janelas na autocorrelação serve para diagnosticar e melhorar a performance dos operadores. Concluímos que muitas perguntas ainda surgem quando técnicas de deconvolução são aplicadas a seções sísmicas de bacias sedimentares, e que o modelo de Goupillaud é conveniente para simulações e análises devido a sua descrição matemática simples e completa.

Atenuação de múltiplas pelo método WHLP-CRS

Nas bacias sedimentares da região Amazônica, a geração e o acúmulo de hidrocarboneto estão relacionados com a presença das soleiras de diabásio. Estas rochas magmáticas intrusivas possuem grandes contrastes de impedância com as rochas sedimentares encaixantes, resultando em múltiplas externas e internas, com amplitudes semelhantes às das reflexões sísmicas primárias. Estas múltiplas podem predominar sobre as informações oriundas de interfaces mais profundas, dificultando o processamento, a interpretação e o imageamento da seção de sísmica. O objetivo da presente tese é realizar a atenuação de múltiplas em seções sintéticas fontecomum (CS), através da combinação dos métodos Wiener-Hopf-Levinson de predição (WHLP) e o do empilhamento superfície-de-reflexão-comum (CRS), aqui denominando pela sigla WHLPCRS. O operador de deconvolução é calculado com as amplitudes reais do sinal sísmico e traço-a-traço, o que consideramos como uma melhor eficiência para a operação de atenuação. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, utilizando o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, e usados para calcular o operador WHLP-CRS. No desenvolvimento do presente trabalho, visamos evitar a inconveniência da seção processada ZO; desenhar e aplicar operadores na configuração CS; e estender o método WHL para camadas curvas.

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

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

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

Adaptação do método de Bott para contraste de densidade variando com a profundidade

Apresentamos um novo método para estimar o relevo do embasamento de bacias sedimentares através da extensão analítica da expressão da placa Bouguer para contraste de densidade entre o pacote sedimentar e o embasamento decrescendo monotonicamente com a profundidade de acordo com uma lei hiperbólica. O método requer ruído contido nos dados gravimétricos tenha desvio padrão inferior a 0,01 mGal. As estimativas do relevo do embasamento são obtidas nas posições espaciais localizadas diretamente abaixo de cada observação. A metodologia foi aplicada a dados sintéticos de bacias sedimentares simuladas apresentando relevo do embasamento suave. O método produziu relevo do embasamento estimado próximo do relevo simulado. O método foi também aplicado a três conjuntos de dados reais: Bacia do Recôncavo (Brasil), Graben do Büyük Menderes (Turquia) e Graben de San Jacinto (Estados Unidos). As soluções produzidas pelo método proposto apresentaram estimativas de profundidades máximas em acordo com as informações geológicas disponíveis: 6 km para Bacia do Recôncavo, 1,6 km para Vale do Büyük Menderes e 2,2 km para Graben de San Jacinto.

Análise de feições geológicas na região de Tartarugalzinho, Amapá, através de métodos geofísicos

Dois métodos geofísicos foram aplicados para a detecção de estruturas geológicas na borda oriental do Escudo das Guianas. O método gravimétrico foi aplicado para encontrar falhas, discordâncias e descontinuidades no embasamento, associadas a cinturões de cisalhamento que possuem orientação NW-SE, e o posicionamento da profundidade com que o embasamento se encontra da superfície. Foram encontradas falhas no embasamento e também indicações de que a crosta apresenta um afinamento em certa região da borda. O método eletromagnético GPR foi utilizado para encontrar descontinuidades geológicas nas rochas sedimentares e averiguar se essas feições são o reflexo das falhas no embasamento e/ou reativação de algumas delas. O GPR localizou feições interpretadas como descontinuidades em rochas e paleoestruturas. O estudo teve o intuito de compreender como as descontinuidades do embasamento podem ter influenciado na formação de bacias sedimentares do Terciário e na formação de paleoestruturas rasas.

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

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