Avaliação geofísica e geológica de um porção de quebra de talude da Bacia do Jequitinhonha

A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies.

Correlação de poços com múltiplos perfis através da rede neural multicamadas

A correlação estratigráfica busca a determinação da continuidade lateral das rochas, ou a equivalência espacial entre unidades litológicas em subsuperfície, a partir de informações geológico-geofísicas oriundas de poços tubulares, que atravessam estas rochas. Normalmente, mas não exclusivamente, a correlação estratigráfica é realizada a partir das propriedades físicas registradas nos perfis geofísicos de poço. Neste caso, busca-se a equivalência litológica a partir da equivalência entre as propriedades físicas, medidas nos vários poços de um campo petrolífero. A técnica da correlação estratigráfica com perfis geofísicos de poço não é uma atividade trivial e sim, sujeita a inúmeras possibilidades de uma errônea interpretação da disposição geométrica ou da continuidade lateral das rochas em subsuperfície, em função da variabilidade geológica e da ambigüidade das respostas das ferramentas. Logo, é recomendável a utilização de um grande número de perfis de um mesmo poço, para uma melhor interpretação. A correlação estratigráfica é fundamental para o engenheiro de reservatório ou o geólogo, pois a partir da mesma, é possível a definição de estratégias de explotação de um campo petrolífero e a interpretação das continuidades hidráulicas dos reservatórios, bem como auxílio para a construção do modelo geológico para os reservatórios, a partir da interpretação do comportamento estrutural das diversas camadas em subsuperfície. Este trabalho apresenta um método de automação das atividades manuais envolvidas na correlação estratigráfica, com a utilização de vários perfis geofísicos de poço, através de uma arquitetura de rede neural artificial multicamadas, treinada com o algoritmo de retropropagação do erro. A correlação estratigráfica, obtida a partir da rede neural artificial, possibilita o transporte da informação geológica do datum de correlação ao longo do campo, possibilitando ao intérprete, uma visão espacial do comportamento do reservatório e a simulação dos possíveis paleoambientes. Com a metodologia aqui apresentada foi possível a construção automática de um bloco diagrama, mostrando a disposição espacial de uma camada argilosa, utilizando-se os perfis de Raio Gama (RG), Volume de Argila (Vsh), Densidade (ρ_b) e de Porosidade Neutrônica (φ_n) selecionados em cinco poços da região do Lago Maracaibo, na Venezuela.

Deconvolução de perfis de poço

A maioria dos perfis de poço utilizados nas avaliações petrofísicas de reservatórios possuem uma resolução vertical na ordem de um metro. Isto cria um problema quando as espessuras típicas das camadas são inferiores a um metro, uma vez que não há correção das leituras. Os perfis de alta resolução vertical como da ferramenta de propagação eletromagnética (EPT, Schlumberger), o dipmeter (SHDT, Schlumberger) ou das ferramentas de varredura acústica ou elétrica possuem uma resolução vertical da ordem de centimetros, mas apresentam uma limitada aplicação para as avaliações petrofísicas. Nós apresentamos um método para a deconvolução de um perfil de baixa resolução vertical que utiliza informações de um perfil de alta resolução vertical para identificar uma nítida interface entre camadas que apresentam valores da propriedade petrofísica contrastante, mas localmente constante em ambos os lados. A partir desse intervalo de controle, nós determinamos a função resposta vertical da ferramenta sob as condições atuais do poço com base no teorema da convolução. Utilizamos várias interfaces de modo a obter valores mais representativos da resposta da ferramenta. O perfil de baixa resolução é então deconvoluido utilizando a transformada discreta de Fourier (FFT) sobre todo o intervalo de interesse. É importante destacar que a invasão do filtrado da lama e a presença do bolo de lama não produzem efeitos danosos sobre o método, que foi aplicado a perfis sintéticos e a dados de campo, onde a aplicação de filtros com um correto ajuste de profundidade, bem como a própria escolha do intervalo de controle, antes da deconvolução, são de extrema importância para o sucesso do método.

Deconvolução de perfis de poços através do ajuste de energia

As medidas de resistividade são de fundamental importância para o cálculo da saturação de óleo em reservatórios potencialmente produtores. A combinação das medidas de resistividade rasa e profunda permite a obtenção dos parâmetros R_t, R_xo e d_i. Mas, em reservatórios complexos existem dificuldades em se obter leituras confiáveis de R_t, devido à baixa resolução vertical das ferramentas de investigação profunda. Em reservatórios laminados, por exemplo, as leituras obtidas pela ferramenta de indução profunda (ILD) podem levar a uma interpretação errônea das mesmas, levando a acreditar que as medidas obtidas do perfil referem-se a uma única camada. Este problema pode ser em parte resolvido através de uma metodologia que melhore a resolução vertical dos perfis de investigação profunda, valendo-se do uso de informações obtidas de um perfil de alta resolução vertical, i.e; a curva de resistividade rasa. Uma abordagem neste sentido seria usar um perfil de alta resolução que apresente uma boa correlação com o perfil de investigação profunda. Esta correlação pode ser melhor avaliada se aplicarmos um filtro no perfil de alta resolução, de tal forma que o perfil resultante tenha teoricamente a mesma resolução vertical do perfil de baixa resolução. A obtenção deste filtro, porém, recai na premissa de que as funções respostas verticais para as ferramentas de alta e baixa resolução são disponíveis, o que não ocorre na prática. Este trabalho se propõe mostrar uma nova abordagem onde o filtro pode ser obtido a partir de um tratamento no domínio da freqüência. Tal tratamento visa igualar a energia espectral do perfil de alta resolução à energia espectral do perfil de baixa resolução tendo como base o Teorema de Parseval. Será mostrado que a resolução vertical depende fundamentalmente da energia espectral do perfil em questão. A seguir, uma regressão linear é aplicada sobre os perfis de alta resolução filtrado e de baixa resolução. Para cada ponto amostrado dos perfis, uma rotina de minimização é aplicada visando escolher o melhor intervalo de correlação entre os perfis. Finalmente, um fator de correção é aplicado sobre cada ponto do perfil de baixa resolução. Os resultados obtidos com os perfis de indução são promissores, demonstrando a eficácia da abordagem e mesmo quando aplicada em perfis com diferentes propriedades petrofísicas, a metodologia funcionou satisfatoriamente, sem danificar os perfis originais.

Detecção do contato entre camadas através do perfil de indução

Dois dos principais objetivos da interpretação petrofísica de perfis são a determinação dos limites entre as camadas geológicas e o contato entre fluidos. Para isto, o perfil de indução possui algumas importantes propriedades: É sensível ao tipo de fluido e a distribuição do mesmo no espaço poroso; e o seu registro pode ser modelado com precisão satisfatória como sendo uma convolução entre a condutividade da formação e a função resposta da ferramenta. A primeira propriedade assegura uma boa caracterização dos reservatórios e, ao mesmo tempo, evidencia os contatos entre fluidos, o que permite um zoneamento básico do perfil de poço. A segunda propriedade decorre da relação quasi-linear entre o perfil de indução e a condutividade da formação, o que torna possível o uso da teoria dos sistemas lineares e, particularmente, o desenho de filtros digitais adaptados à deconvolução do sinal original. A idéia neste trabalho é produzir um algoritmo capaz de identificar os contatos entre as camadas atravessadas pelo poço, a partir da condutividade aparente lida pelo perfil de indução. Para simplificar o problema, o modelo de formação assume uma distribuição plano-paralela de camadas homogêneas. Este modelo corresponde a um perfil retangular para condutividade da formação. Usando o perfil de entrada digitalizado, os pontos de inflexão são obtidos numericamente a partir dos extremos da primeira derivada. Isto gera uma primeira aproximação do perfil real da formação. Este perfil estimado é então convolvido com a função resposta da ferramenta gerando um perfil de condutividade aparente. Uma função custo de mínimos quadrados condicionada é definida em termos da diferença entre a condutividade aparente medida e a estimada. A minimização da função custo fornece a condutividade das camadas. O problema de otimização para encontrar o melhor perfil retangular para os dados de indução é linear nas amplitudes (condutividades das camadas), mas uma estimativa não linear para os contatos entre as camadas. Neste caso as amplitudes são estimadas de forma linear pelos mínimos quadrados mantendo-se fixos os contatos. Em um segundo passo mantem-se fixas as amplitudes e são calculadas pequenas mudanças nos limites entre as camadas usando uma aproximação linearizada. Este processo é interativo obtendo sucessivos refinamentos até que um critério de convergência seja satisfeito. O algoritmo é aplicado em dados sintéticos e reais demonstrando a robustez do método.

Determinação automática da porosidade e zoneamento de perfis através da rede neural artificial competitiva

Duas das mais importantes atividades da interpretação de perfis para avaliação de reservatórios de hidrocarbonetos são o zoneamento do perfil (log zonation) e o cálculo da porosidade efetiva das rochas atravessadas pelo poço. O zoneamento é a interpretação visual do perfil para identificação das camadas reservatório e, consequentemente, dos seus limites verticais, ou seja, é a separação formal do perfil em rochas reservatório e rochas selante. Todo procedimento de zoneamento é realizado de forma manual, valendo-se do conhecimento geológico-geofísico e da experiência do intérprete, na avaliação visual dos padrões (características da curva do perfil representativa de um evento geológico) correspondentes a cada tipo litológico específico. O cálculo da porosidade efetiva combina tanto uma atividade visual, na identificação dos pontos representativos de uma particular rocha reservatório no perfil, como a escolha adequada da equação petrofísica que relaciona as propriedades físicas mensuradas da rocha com sua porosidade. A partir do conhecimento da porosidade, será estabelecido o volume eventualmente ocupado por hidrocarboneto. Esta atividade, essencial para a qualificação de reservatórios, requer muito do conhecimento e da experiência do intérprete de perfil para a efetiva avaliação da porosidade efetiva, ou seja, a porosidade da rocha reservatório, isenta do efeito da argila sobre a medida das propriedades físicas da mesma. Uma forma eficiente de automatizar estes procedimentos e auxiliar o geofísico de poço nestas atividades, que particularmente demandam grande dispêndio de tempo, é apresentado nesta dissertação, na forma de um novo perfil, derivado dos perfis tradicionais de porosidade, que apresenta diretamente o zoneamento. Pode-se destacar neste novo perfil as profundidades do topo e da base das rochas reservatório e das rochas selante, escalonado na forma de porosidade efetiva, denominado perfil de porosidade efetiva zoneado. A obtenção do perfil de porosidade efetiva zoneado é baseado no projeto e execução de várias arquiteturas de rede neural artificial, do tipo direta, com treinamento não supervisionado e contendo uma camada de neurônios artificiais, do tipo competitivo. Estas arquiteturas são projetadas de modo a simular o comportamento do intérprete de perfil, quando da utilização do gráfico densidade-neutrônico, para as situações de aplicabilidade do modelo arenito-folhelho. A aplicabilidade e limitações desta metodologia são avaliadas diretamente sobre dados reais, oriundos da bacia do Lago Maracaibo (Venezuela).

Determinação estatística dos contatos fluidos em perfis de poço

No interior dos reservatórios petrolíferos, o efeito da gravidade produz naturalmente a segregação dos fluidos presentes. Em razão basicamente das forças capilares não existe uma superfície definida de separação entre cada dois fluidos. No entanto, é comum na engenharia de petróleo admitir-se uma superfície convencional de separação dos fluidos. À essa superfície dá-se o nome de interface fluida ou contato fluido. A obtenção da localização dos contatos fluidos no interior dos reservatórios, como os contatos óleo-água, contatos gás-água e os contatos gás-óleo, sem a interpretação completa dos perfis geofísicos ou a realização dos procedimentos diretos dos testes de formação é um problema de difícil solução para a indústria do petróleo. O posicionamento em profundidade dos contatos pode ser utilizado para a definição dos procedimentos de completação do poço e para o posicionamento do canhoneio nos horizontes produtores dos poços verticais. Temos também a aplicação do mapeamento deste parâmetro, para os projetos e controle da perfuração no caso dos poços direcionais e horizontais, que utilizam tal informação para o posicionamento do eixo do poço no interior da zona de hidrocarboneto no reservatório, a fim de se evitar, principalmente, a produção de água. Este trabalho apresenta uma metodologia para a identificação e posicionamento dos contatos fluidos, através da aplicação da técnica da estatística multivariada, chamada de análise discriminante. Nesta aplicação, faz-se uso direto dos perfis de resistividade (R_T), de raios gama natural (RG), de densidade (ρ_b), de porosidade neutrônica (Φ_N), de caliper (Cal) e do perfil de volume de argila (Vcla), como dados de entrada no programa de análise discriminante. Para o caso de depósitos clásticos, com seqüências arenito-folhelho, a análise discriminante fornece ainda a indicação da litologia e as espessuras aparentes dos reservatórios. As curvas de perfis geofísicos aplicadas para a avaliação desta metodologia são provenientes de poços perfurados na bacia do lago Maracaibo, na Venezuela, onde, segundo estudos geológicos, existe a ocorrência de seqüências estratigráficas, com camadas de folhelhos e arenitos, sendo também confirmada a presença de hidrocarbonetos nos poços utilizados neste trabalho.

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.

Estimativa dos perfis de permeabilidade e de porosidade utilizando rede neural artificial

A permeabilidade e a porosidade são duas das mais importantes propriedades petrofísicas para a qualificação dos reservatórios de óleo e gás. A porosidade está relacionada à capacidade de armazenamento de fluidos e a permeabilidade, com a capacidade de produção destes fluidos. Suas medidas são, normalmente, realizadas em laboratório, através de testemunhos da rocha. Esses processos têm custos elevados e nem todos os poços são testemunhados. As estimativas da permeabilidade e da porosidade são de fundamental importância para os engenheiros de reservatório e geofísicos, uma vez que seus valores podem definir a completação ou não de um poço petrolífero. O perfil de porosidade e sua relação com o perfil de densidade, é bem conhecida na geofísica de poço. No entanto, existem poucas relações quantitativas e/ou qualitativas entre a porosidade e a permeabilidade, como por exemplo as relações de Kozeny. Sendo assim, este trabalho busca o estabelecimento do perfil de permeabilidade e do perfil de porosidade, a partir de informações do perfil de densidade. Para tanto, buscamos a relação entre a propriedade física da rocha (densidade) e as propriedades petrofísicas: permeabilidade e porosidade, utilizando como metodologia à técnica de redes neurais artificiais, como a rede neural artificial com função de base radial. A obtenção da permeabilidade e da porosidade a partir da rede neural artificial, que possui como entrada a informação da densidade possibilita um menor custo para a aquisição dessas importantes informações petrofísicas, permite ao intérprete de perfis de poço optar ou não pela exploração de uma unidade estudada, além de uma visão mais completa do reservatório. Os procedimentos para a estimativa da permeabilidade e da porosidade estão direcionados para uma única formação, mas os intérpretes de perfis poderão aplicar a diretriz apresentada no programa de rede neural artificial com função de base radial, utilizando a estimativa dessas propriedades petrofísicas para outras formações, inclusive de outros campos petrolíferos. Portanto, recomenda-se a utilização de um conjunto de dados completo, com quantidade de dados suficientes de um mesmo poço, a fim de viabilizar corretamente a melhor interpretação.

Estudo da resposta de diferentes arranjos de bobinas na perfilagem de indução de poço: modelamento analógico

A sonda de indução tem sido largamente utilizada na perfilagem de poço para detectar reservatórios através do contraste de condutividade existente entre as zonas contendo hidrocarbonetos e as zonas contendo água conata. A sonda de indução, além de auxiliar na localização de formações que contêm hidrocarbonetos, possibilita o cálculo de reservas para análise da viabilidade de exploração do reservatório. Porém, apesar do grande avanço da sonda de indução, a qual utiliza arranjo coaxial, existem várias restrições ao seu uso. Camadas delgadas, poços preenchidos com lama condutiva e invasão profunda são algumas das situações onde não é adequada a utilização da sonda de indução. Considerando as limitações da sonda de indução é que desenvolvemos em laboratório um estudo sobre as diferentes disposições de bobinas para a perfilagem de indução. O estudo consiste em comparar a resposta da sonda de arranjo de bobinas coaxiais, com a resposta da sonda de arranjo de bobinas coplanares e verificar se as respostas encontradas resolvem alguns dos problemas apresentados pela sonda de bobinas coaxiais. As respostas das sondas de diferentes disposições de bobinas foram obtidas através da montagem de um sistema de modelamento analógico simulando um poço. Para o sistema modelo foi utilizado um fator de escala linear igual a 20. Um tanque contendo soluções salinas de diferentes condutividades simulava algumas situações geológicas controladas. O sistema permitia uma variação na frequência de 50 kHz a 100 kHz. Os resultados obtidos foram satisfatórios, em virtude do arranjo de bobinas coplanares ter apresentado melhor resolução vertical que o arranjo de bobinas coaxiais.

Freqüência de anticorpos anti-Borrelia burgdorferi em eqüinos na mesorregião metropolitana de Belém, Estado do Pará

A Borreliose de Lyme é uma doença multissistêmica causada pela espiroqueta Borrelia burgdorferi latu sensu e transmitida por carrapatos ixodideos, acometendo seres humanos e animais domésticos, tendo como reservatórios naturais os animais silvestres, sendo ainda considerada uma zoonose de ampla distribuição geográfica. Foram coletadas 300 amostras sanguíneas de eqüinos aparentemente sadios, procedentes dos municípios de Ananideua, Belém, Benevides, Castanhal, Marituba e Santa Izabel do Pará da mesorregião metropolitana de Belém – Pará. O sangue foi coletado pela veia jugular e os soros foram analisados através do ensaio de imunoadsorção enzimática (ELISA) indireto no Laboratório de Doenças Parasitárias na Universidade Federal Rural do Rio de Janeiro. O objetivo do trabalho foi verificar a soroepidemiologia da Borrelia burgdorferi na mesorregião metropolitana de Belém. A freqüência de soropositividade foi de 26,67% (n=80), sendo 72 (24%) com título de 1:800, seis (2%) a 1:1600 e dois (0,6%) a 1:3200. A soropositividade entre os sexos foram valores aproximados sendo 13,67% nas fêmeas e 13% nos machos. A freqüência entre os animais de raças e mestiços foram 9% e 18%, respectivamente. Não houve diferença significativa de animais soropositivos quanto aos municípios, sexo, raça e faixa etária. A freqüência encontrada corrobora a hipótese da ocorrência de Borrelia sp. na região estudada.

O impacto da infecção por Chlamydia em populações indígenas da Amazônia brasileira

A disseminação das bactérias do gênero Chlamydia no Brasil, inclusive na região Amazônica, é pouco conhecida. Este estudo soroepidemiológico incluiu 2.086 amostras de soro de populações indígenas da Amazônia brasileira, empregando metodologia de triagem pela imunofluorescência indireta para pesquisa de anticorpos. Usou-se o sorotipo L2 da C. trachomatis como substrato; a seguir, para os quinze sorotipos de C. trachomatis e para a C. pneumoniae, discriminou-se a sororreatividade pela microimunofluorescência específica. A prevalência média de anticorpos para Chlamydia foi de 48,6%. Sua variação entre as comunidades indicou as que não tiveram contato com as bactérias e aquelas em que quase todos os testados tiveram. Por meio da titulação dos anticorpos IgG e a presença de IgM específica nas amostras com títulos altos viu-se que 6,1% dos infectados persistiam com a infecção, servindo de reservatórios à disseminação das espécies de Chlamydia. Pela resposta à C. trachomatis, evidenciou-se a circulação dos sorotipos A, B, Ba, D, E, G, H, I e L1. Ademais, constatou-se que há C. pneumoniae na região. As duas espécies causariam impacto significativo no hospedeiro humano.

Interpretação do gráfico de Hingle através de rede competitiva angular

A saturação de água é a principal propriedade petrofísica para a avaliação de reservatórios de hidrocarbonetos, pois através da análise dos seus valores é definida a destinação final do poço recém perfurado, como produtor ou poço seco. O cálculo da saturação de água para as formações limpas é, comumente, realizado a partir da equação de Archie, que envolve a determinação da resistividade da zona virgem, obtida a partir de um perfil de resistividade profunda e o cálculo de porosidade da rocha, obtida a partir dos perfis de porosidade. A equação de Archie envolve ainda, a determinação da resistividade da água de formação, que normalmente necessita de definição local e correção para a profundidade da formação e da adoção de valores convenientes para os coeficientes de Archie. Um dos métodos mais tradicionais da geofísica de poço para o cálculo da saturação de água é o método de Hingle, particularmente útil nas situações de desconhecimento da resistividade da água de formação. O método de Hingle estabelece uma forma linear para a equação de Archie, a partir dos perfis de resistividade e porosidade e a representa na forma gráfica, como a reta da água ou dos pontos, no gráfico de Hingle, com saturação de água unitária e o valor da resistividade da água de formação é obtido a partir da inclinação da reta da água. Independente do desenvolvimento tecnológico das ferramentas de perfilagem e dos computadores digitais, o geofísico, ainda hoje, se vê obrigado a realizar a interpretação de ábacos ou gráficos, sujeito a ocorrência de erros derivados da sua acuidade visual. Com o objetivo de mitigar a ocorrência deste tipo de erro e produzir uma primeira aproximação para a saturação de água em tempo real de perfilagem do poço, insere-se o trabalho apresentado nesta dissertação, com a utilização de uma conveniente arquitetura de rede neural artificial, a rede competitiva angular, capaz de identificar a localização da reta da água, a partir da identificação de padrões angulares presentes nos dados dos perfis de porosidade e resistividade representados no gráfico de Hingle. A avaliação desta metodologia é realizada sobre dados sintéticos, que satisfazem integralmente a equação de Archie, e sobre dados reais.

Modelamento de perfis de indução

Neste trabalho foi implementado um algoritmo que permite o modelamento das respostas das sondas de indução (2C40-6FF40-ILD) em meios estratificados. O algoritmo é baseado na solução exata do problema de Green para a obtenção do campo eletromagnético em meios estratificados com a fonte no interior dos estratos. O código computacional desenvolvido neste trabalho para meios horizontalmente estratificados tem o objetivo de auxiliar o desenvolvimento de códigos mais avançados, uma vez que é difícil o acesso para a pesquisa dos códigos citados na literatura. Os códigos disponíveis são comerciais, mas estes têm o objetivo único de produção, além de serem, bastante dispendiosos. O código também foi utilizado para a obtenção da resposta da sonda ILD em reservatórios do tipo laminado, onde a sonda tem baixa resolução. E verificamos, a exemplo de outros autores, que a condutividade obtida pela sonda pode ser expressa pela soma das condutividades dos constituintes litológicos do reservatório ponderada pelos seus percentuais em volume. Esta relação foi definida para reservatórios ainda não analisados na literatura, tipo arenito moderadamente silicificado/arenito/folhelho, onde a lâmina de interesse tem resistividade superior a resistividade de uma das lâminas, e inferior a resistividade da outra. Foi analisado também reservatórios tipo arenito moderadamente silicificado/arenito, onde a lâmina de interesse tem a menor resistividade. Com esta relação pode-se obter para estes reservatórios laminados estimativas mais precisas de Rt e consequentemente determinar as saturações água/óleo e produtibilidade, também de forma mais precisa.

Paleoambiente e proveniência da formação cabeças da bacia do Parnaíba: evidências da glaciação famenniana e implicações na potencialidade do reservatório

O histórico de prospecção de hidrocarbonetos da Bacia Paleozoica do Parnaíba, situada no norte-nordeste do Brasil, sempre foi considerado desfavorável quando comparado aos super-reservatórios estimados do Pré-Sal das bacias da Margem Atlântica e até mesmo interiores, como a Bacia do Solimões. No entanto, a descoberta de gás natural em depósitos da superseqüência mesodevoniana-eocarbonífera do Grupo Canindé, que incluem as formações Pimenteiras, Cabeças e Longá, impulsionou novas pesquisas no intuito de refinar a caracterização paleoambiental, paleogeográfica, bem como, entender o sistema petrolífero, os possíveis plays e a potencialidade do reservatório Cabeças. A avaliação faciológica e estratigráfica com ênfase no registro da tectônica glacial, em combinação com a geocronologia de zircão detrítico permitiu interpretar o paleoambiente e a proveniência do reservatório Cabeças. Seis associações de fácies agrupadas em sucessões aflorantes, com espessura máxima de até 60m registram a evolução de um sistema deltaico Devoniano influenciado por processos glaciais principalmente no topo da unidade. 1) frente deltaica distal, composta por argilito maciço, conglomerado maciço, arenito com acamamento maciço, laminação plana e estratificação cruzada sigmoidal 2) frente deltaica proximal, representada pelas fácies arenito maciço, arenito com laminação plana, arenito com estratificação cruzada sigmoidal e conglomerado maciço; 3) planície deltaica, representada pelas fácies argilito laminado, arenito maciço, arenito com estratificação cruzada acanalada e conglomerado maciço; 4) shoreface glacial, composta pelas fácies arenito com marcas onduladas e arenito com estratificação cruzada hummocky; 5) depósitos subglaciais, que englobam as fácies diamictito maciço, diamictito com pods de arenito e brecha intraformacional; e 6) frente deltaica de degelo, constituída pelas fácies arenito maciço, arenito deformado, arenito com laminação plana, arenito com laminação cruzada cavalgante e arenito com estratificação cruzada sigmoidal. Durante o Fammeniano (374-359 Ma) uma frente deltaica dominada por processos fluviais progradava para NW (borda leste) e para NE (borda oeste) sobre uma plataforma influenciada por ondas de tempestade (Formação Pimenteiras). Na borda leste da bacia, o padrão de paleocorrente e o espectro de idades U-Pb em zircão detrítico indicam que o delta Cabeças foi alimentado por áreas fonte situadas a sudeste da Bacia do Parnaíba, provavelmente da Província Borborema. Grãos de zircão com idade mesoproterozóica (~ 1.039 – 1.009 Ma) e neoproterozóica (~ 654 Ma) são os mais populosos ao contrário dos grãos com idade arqueana (~ 2.508 – 2.678 Ma) e paleoproterozóica (~ 2.054 – 1.992 Ma). O grão de zircão concordante mais novo forneceu idade ²⁰⁶Pb/²³⁸U de 501,20 ± 6,35 Ma (95% concordante) indicando idades de áreas-fonte cambrianas. As principais fontes de sedimentos do delta Cabeças na borda leste são produto de rochas do Domínio Zona Transversal e de plútons Brasilianos encontrados no embasamento a sudeste da Bacia do Parnaíba, com pequena contribuição de sedimentos oriundos de rochas do Domínio Ceará Central e da porção ocidental do Domínio Rio Grande do Norte. No Famenniano, a movimentação do supercontinente Gondwana para o polo sul culminou na implantação de condições glaciais concomitantemente com o rebaixamento do nível do mar e exposição da região costeira. O avanço das geleiras sobre o embasamento e depósitos deltaicos gerou erosão, deposição de diamictons com clastos exóticos e facetados, além de estruturas glaciotectônicas tais como plano de descolamento, foliação, boudins, dobras, duplex, falhas e fraturas que refletem um cisalhamento tangencial em regime rúptil-dúctil. O substrato apresentava-se inconsolidado e saturados em água com temperatura levemente abaixo do ponto de fusão do gelo (permafrost quente). Corpos podiformes de arenito imersos em corpos lenticulares de diamicton foram formados pela ruptura de camadas pelo cisalhamento subglacial. Lentes de conglomerados esporádicas (dump structures) nos depósitos de shoreface sugere queda de detritos ligados a icebergs em fases de recuo da geleira. A elevação da temperatura no final do Famenniano reflete a rotação destral do Gondwana e migração do polo sul da porção ocidental da América do Sul e para o oeste da África. Esta nova configuração paleogeográfica posicionou a Bacia do Parnaíba em regiões subtropicais iniciando o recuo de geleiras e a influência do rebound isostático. O alívio de pressão é indicado pela geração de sills e diques clásticos, estruturas ball-and-pillow, rompimento de camadas e brechas. Falhas de cavalgamento associadas à diamictitos com foliação na borda oeste da bacia sugerem que as geleiras migravam para NNE. O contínuo aumento do nível do mar relativo propiciou a instalação de sedimentação deltaica durante o degelo e posteriormente a implantação de uma plataforma transgressiva (Formação Longá). Diamictitos interdigitados com depósitos de frente deltaica na porção superior da Formação Cabeças correspondem a intervalos com baixo volume de poros e podem representar trapas estratigráficas secundárias no reservatório. As anisotropias primárias subglaciais do topo da sucessão Cabeças, em ambas as bordas da Bacia do Parnaíba, estende a influência glacial e abre uma nova perspectiva sobre a potencialidade efetiva do reservatório Cabeças do sistema petrolífero Mesodevoniano-Eocarbonífero da referida bacia.