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.

Tomografia eletromagnética para caracterização de reservatórios de hidrocarbonetos

Na produção de petróleo é importante o monitoramento dos parâmetros do reservatório (permeabilidade, porosidade, saturação, pressão, etc) para o seu posterior gerenciamento. A variação dos parâmetros dinâmicos do reservatório induz variações na dinâmica do fluxo no reservatório, como por exemplo, perdas na pressão, dificultando o processo de extração do óleo. A injeção de fluidos aumenta a energia interna do reservatório e incrementa a pressão, estimulando o movimento do óleo em direção aos poços de extração. A tomografia eletromagnética poço-a-poço pode se tomar em uma técnica bastante eficaz no monitoramento dos processos de injeção, considerando-se o fato de ser altamente detectável a percolação de fluidos condutivos através das rochas. Esta tese apresenta o resultado de um algoritmo de tomografia eletromagnética bastante eficaz aplicado a dados sintéticos. O esquema de imageamento assume uma simetria cilíndrica em torno de uma fonte constituída por um dipolo magnético. Durante o processo de imageamento foram usados 21 transmissores e 21 receptores distribuídos em dois poços distanciados de 100 metros. O problema direto foi resolvido pelo método dos elementos finitos aplicado à equação de Helmhotz do campo elétrico secundário. O algoritmo resultante é válido para qualquer situação, não estando sujeito às restrições impostas aos algoritmos baseados nas aproximações de Born e Rytov. Por isso, pode ser aplicado eficientemente em qualquer situação, como em meios com contrastes de condutividade elétrica variando de 2 a 100, freqüências de 0.1 a 1000.0 kHz e heterogeneidades de qualquer dimensão. O problema inverso foi resolvido por intermédio do algoritmo de Marquardt estabilizado. A solução é obtida iterativamente. Os dados invertidos, com ruído Gaussiano aditivo, são as componentes em fase e em quadratura do campo magnético vertical. Sem o uso de vínculos o problema é totalmente instável, resultando em imagens completamente borradas. Duas categorias de vínculos foram usadas: vínculos relativos, do tipo suavidade, e vínculos absolutos. Os resultados obtidos mostram a eficiência desses dois tipos de vínculos através de imagens nítidas de alta resolução. Os tomogramas mostram que a resolução é melhor na direção vertical do que na horizontal e que é também função da freqüência. A posição e a atitude da heterogeneidade é bem recuperada. Ficou também demonstrado que a baixa resolução horizontal pode ser atenuada ou até mesmo eliminada por intermédio dos vínculos.

Estudo da contaminação por hidrocarbonetos em postos de serviços de Abaetetuba-PA utilizando georadar

O trabalho apresenta os resultados de um levantamento Geofísico utilizando o método Radar de Penetração de Solo (GPR) ou Georadar para detectar possíveis zonas de contaminação provocadas por vazamentos de derivados de hidrocarbonetos em postos de serviços da região urbana do município de Abaetetuba, no estado do Pará. A metodologia foi aplicada em postos de serviços porque eles constituem uma das principais fontes potenciais urbanas de contaminação do solo e de aqüíferos rasos por combustíveis. Conceitos básicos sobre a contaminação por derivados do petróleo e sua interação com o subsolo são apresentados, juntamente com os princípios básicos que permitem o entendimento do funcionamento do método GPR para o problema abordado. Durante o trabalho, foram realizadas medidas do nível de água em poços rasos, visando a elaboração de um mapa de fluxo subterrâneo. A interpretação geofísica foi auxiliada pelo conhecimento do comportamento do fluxo hídrico subterrâneo local, que mostra o sentido de movimentação da provável contaminação. A correlação entre os dados de GPR, os dados de fluxo e as informações sobre o histórico dos postos levaram a classificá-los em suspeitos de produzir contaminação e possivelmente contaminados.

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

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

Hidrocarbonetos policíclicos aromáticos em sedimentos e organismos bentônicos do Terminal de Miramar (Baía do Guajará - Belém - Pará - Amazônia)

ABSTRACT: In Guajará bay, mouth of the Amazon River, is located the capital of Pará state (Belém) and its metropolitan area. Exists in this area an intense boat traffic, as well as transportation and sale of fuels in floating docks and activities related to the storage and transportation of oil in Miramar Petrochemical Terminal (TEMIR). Small spills and discharges of oil in water can serve as sources of pollution by polycyclic aromatic hydrocarbons (PAHs). PAHs are organic compounds generated by incomplete combustion of organic matter (OM) and are among the contaminants of most interest in environmental studies due to their mutagenic and carcinogenic potential. One way to detect and evaluate the impact of PAHs in an environment is using biomonitors, however the qualitative and quantitative analysis in sediments are most widely used. This study aimed to evaluate, by Gas Chromatography/Mass Spectrometry (GC/MS), the 16 PAHs considered as priority by the Environmental Protection Agency of the United States in sediments and benthic organisms (Namalicastys abiuma) of TEMIR. Field expeditions were carried out in December 2012, March, May and June 2013 representing the dried, wet (two samples) and dried seasons successively. Fine sediments levels dominated in both sampling periods and the %Fines in Guajará bays has a direct relationship to the %OM. In relation to the 16 PAHs studied, 10 of them were detected in the wet season sediments samples and 8 in the dried season. Even with a low diversity of aromatic compounds, sediment samples of the dried season presented greater ΣHPAs (1.351,43 ng g^-1) than the sediments of the wet season (263,99 ng g^-1), which can be related to the increased hydrodynamic in Guajará bay this last period. Correlation analysis indicated that ΣHPAs not seem to be influenced by %Fine and %OM. Benzo(a)pyrene accounted for 87% of ΣHPAs during the wet season, the other PAHs showed percentage ≤ 3%. During the other period stood out: pyrene (18% of ΣHPAs) fluoranthene (16%), the chrysene, benzo(b)fluoranthene (15%) and benzo(a)pyrene (11%). The use of the geochemical reason to interpret likely sources of PAHs indicated that there is a predominance of aromatics from pyrolytic origin in TEMIR. The activities that may be responsible for the PAHs levels observed in sediments and benthic organisms in TEMIR are the discharge of untreated domestic sewage and the emissions of particulates from cars and small/medium boats. The ΣHPAs seems to influence the density of polychaetes because a reduction of about 50% in the number of organisms was observed during the dry season when there was the greater amount of ΣHPAs in their body structure. Among the 16 PAHs studied, 11 were detected in the polychaetes during the dry season (ΣHPAs_MED = 848,71 ng g^-1) and 10 in the wet season (ΣHPAs_MED = 141,85 ng g^-1). Were highlighted during the dried season: indeno(1,2,3-c,d)pyrene (47%) and pyrene (23%). While in the wet season: pyrene (23%), chrysene (17%), fluorene (17%) and fluoranthene (13%). The %Rec obtained was >>100% indicating a matrix effect and reducing the quantitative accuracy of results. The use of the polychaete N. abiuma as biomonitor of pollution by PAHs in estuaries not seem to be effective when the sample period is limited, because collecting it requires a great effort to obtain a small amount of mass and it still results in a matrix effect in the chromatographic analysis (method %Rec >> 100%) that may not be corrected due to the restricted availability of sample mass. Studies with a larger number of samples, enabling depth statistical analysis, and covering other collection points in Guajará bay are required to prove mathematically that has been stated in this dissertation.

Craqueamento termocatalítico da borra de neutralização do óleo de palma (Elaeis guineensis) em escala piloto

Os subprodutos das indústrias de refino de óleo e de biodiesel, como a borra de neutralização, são bastante extensivos e devido a sua baixa pureza e valor econômico constituem uma problemática para essas indústrias no descarte e destinação, sendo de fundamental interesse pesquisas para a utilização desse rejeito. Portanto, este trabalho visa investigar uma alternativa viável, econômica e ambiental para o destino deste subproduto usando o Processo de Craqueamento Termocatalítico da Borra de Neutralização do Óleo de Palma (Elaeis guineensis, Jaqc) em escala piloto utilizando-se 5, 10 e 15% de Carbonato de Sódio (Na²CO³) como catalisador e temperaturas finais de 440ºC e 420ºC. A borra foi obtida pelo processo de neutralização e submetida a uma desidratação e caracterizada assim como o catalisador foi desidratado em estufa e caracterizado em relação à Difração de Raio-X, Análise Térmica Gravimétrica (ATG) e à Análise Térmica Diferencial (TDA). Foram realizados cinco testes de craqueamento termocatalítico na Usina Piloto de Craqueamento (THERMTEK/FEQ/UFPA) o qual é constituído em um reator com agitação mecânica e capacidade de 125 litros, além da dinâmica do processo e destilação dos produtos do craqueamento. O produto líquido orgânico (PLO), amostras da dinâmica do processo e frações da destilação foram caracterizados de acordo com cada norma exigida pela ANP N°65 e analisadas por IV e RMN. Após as análises dos resultados verificou-se que a eficiência do processo aumenta com catalisador e temperatura e que a matéria-prima fornece produtos de baixa acidez e com boas características para uso como combustível. Pela análise da termodinâmica do processo percebeu que alguns parâmetros como viscosidade, densidade e ponto de fulgor diminuem com o tempo e aumento da temperatura, formando hidrocarbonetos mais leves. Com relação à destilação, as frações nas faixas mais pesadas se assemelham ao óleo diesel do petróleo na maioria dos parâmetros exigidos pela ANP N°65.

Craqueamento termocatalítico do óleo de fritura residual

Neste trabalho foi estudado o processo de craqueamento termocatalítico do óleo de fritura nas escalas de bancada e piloto, variando-se o percentual do catalisador carbonato de sódio de 5 e 10% m/m em relação a matéria prima utilizada e temperatura de 440 ºC. O objetivo foi obter misturas de hidrocarbonetos ricas na fração diesel. O óleo de fritura neutralizado e seco foi caracterizado em relação ao Índice de Acidez, Índice de saponificação, Viscosidade Cinemática, Densidade e Índice de Refração. Após o craqueamento, o produto líquido obtido foi purificado por decantação da fase aquosa e filtração simples em escala de bancada. Esse produto foi fracionado por destilação fracionada e os condensados foram coletados em um funil de decantação de acordo a faixa de destilação da gasolina (40ºC-175ºC), querosene (175ºC-235ºC), diesel leve (235°C-305ºC) e diesel pesado (305ºC-400 ºC). Foi realizada a caracterização tanto físico química quanto da composição dos produtos líquidos e suas respectivas frações. Também foi realizada a evolução do processo de craqueamento em escala piloto, acompanhando o comportamento das características físico químicas e de composição do produto formado no decorrer do processo de craqueamento. Os resultados mostraram que o catalisador carbonato de sódio forneceu produtos de baixa acidez e com boas características para uso como combustível. A variação do percentual de catalisador influencia significamente as propriedades físico químicas e composição tanto do produto quanto de suas frações. Verificou-se, ainda, que o craqueamento termocatalítico do óleo de fritura propicia a formação de hidrocarbonetos ricos na fração do diesel (19,16% diesel leve e 41,18% diesel pesado para o teste com 10% de Na₂CO₃ e de 13,53% leve e 52,73% diesel pesado para o teste com 5% de Na₂CO₃ ). Os intervalos de tempos finais do craqueamento geram um combustível com baixo teor de acidez e com propriedades físico químicas em conformidade a norma especificada para o diesel mineral.

Simulação numérica da gaseificação de biomassa em leito fixo

O objetivo deste trabalho foi avaliar a potencialidade das sementes de açaí rejeitadas após beneficiamento em Belém – Pará, a fim de ser usada como fonte de biomassa para a produção de energia por gaseificação. Resultados da análise elementar obtidas no LABCAT – UNB em base seca foram obtidas com percentuais de Carbono, 46,06 %, Hidrogênio 6,01 % Oxigênio 43,37 %, Nitrogênio 4,33 % e Enxofre 0,08%. Tais valores entre outros parâmetros obtidos em estudo experimental realizados por Santos 2011 e Cruz, 2010, foram utilizados para realizar uma simulação numérica do processo de gaseificação por meio de um código FORTRAN zero dimensional desenvolvido especificamente para esta finalidade. Variações da influência da razão de equivalência e teor de umidade nas concentrações de CO, H2 e CH4 foram comparados a resultados experimentais. Um modelo transiente unidimensional também foi analisado. Utilizando a aplicação de um mecanismo de pirólise de dois passos onde o alcatrão primário (ou oxigenados) são inicialmente formados e então craqueados em alcatrão secundário (ou hidrocarbonetos) e outros gases combustíveis. Assim se revela a complexidade dos processos físicos e químicos que ocorrem no reator pelo uso do balanço de energia e massa e com informações das taxas de reação das reções químicas e os processos de transporte físicos. O modelo computacional é capaz de prever o perfil da composição de gases, temperatura, alcatrão primário e secundário ao longo do eixo longitudinal do gaseificador. Foram obtidos resultados para simulações realizadas nos dois modelos e posteriormente analisado tal comportamento em comparação ao resultado experimental obtidos por Santos (2011).

Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Aplicação do gravímetro de poço em problemas de avaliação de formação no Brasil

No período de outubro a novembro de 1990, foi realizado na Bacia Sergipe-Alagoas o primeiro levantamento gravimétrico de poço no Brasil, através de um programa com a participação da Universidade Federal do Pará, a PETROBRÁS e o U.S.G.S. (United States Geological Survey). Este levantamento teve como objetivos o teste do equipamento, a resposta dos perfis gravimétricos em relação aos problemas encontrados nas bacias sedimentares brasileiras e a comparação com os dados de densidade obtidos com o perfil de densidade compensada (CDL). Os levantamentos foram realizados em três poços e os dados obtidos passaram por um processo de redução, onde os valores foram transformados para miligals e corrigidos dos efeitos de maré, deriva e terreno, para obter as densidades preliminares. Porém, a região onde foram realizados os levantamentos apresenta além do embasamento muito raso, uma seqüência evaporítica com contraste de densidade grande em relação ao resto do pacote sedimentar, gerando gradientes verticais anômalos. Estes efeitos foram então corrigidos e obtidos os valores finais de densidade. Com base em uma caracterização litológica prévia, os valores de densidade foram então comparados com o perfil CDL. Verificou-se a boa qualidade dos dados, o que permitiu o cálculo de valores de porosidade, mostrando assim possíveis intervalos de interesse à exploração de hidrocarbonetos, tanto no embasamento, reservatório mais importante, quanto no intervalo sedimentar.

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.

Avaliação dos reservatórios do Campo de Ubarana-RN, com baixas resistividades

Desenvolveu-se neste trabalho um novo modelo de avaliação quantitativa para arenitos argilosos, baseado na bibliografia e no estudo da disposição dos argilominerais dentro dos poros das rochas reservatório. Este novo modelo leva em consideração a contribuição de duas condutividades extras, além da condutividade eletrolítica das rochas: uma se deve à contribuição efetiva da condutância superficial dos argilominerais, aumentada ou diminuida em função da capacidade de troca catiônica dos mesmos, e outra independente e em paralelo, devido a condutância que se desenvolve em uma rede microporosa contínua formada por argilominerais sobre o arcabouço das rochas. Para o estudo da validade da equação proposta, foram feitos vários gráficos comparativos entre as equações e/ou modelos existentes n” bibliografia especializada, deles resultando ótimas correlações, principalmente com a tradicional equação de Waxman & Smits. Foi feita uma comparação do novo modelo em três seções distintas do Campo de Ubarana, Rio Grande do Norte, duas com baixas resistividades supostamente com hidrocarbonetos e uma outra efetivamente produtora. Verificou-se que as baixas resistividades resultam do fato das seções estudadas serem subsaturadas, com menos de 50% de óleo e aproximadamente igual ou maior proporção de água relativamente salgada (85.000 ppm de NaCl equivalente). Esta mesma água, sem dúvida, muito contribui para as baixas resistividades, por formar a fase eletricamente condutiva das referidas seções. Sendo a equação de Waxman & Smits mundialmente reconhecida, pode-se afirmar que a equação proposta neste trabalho tem consistência teórica e prática e para o caso particular do campo de Ubarana, mostrou-se mais coerente com o histórico de produção dos poços estudados, do que as demais equações existentes e testadas.

Construção de um modelo experimental simulando condições de perfilagem de indução

A ferramenta de indução tem sido, desde 1949, largamente utilizada na prospecção de petróleo pois auxilia na localização de formações que contém hidrocarbonetos. Ainda existem, porém, várias restrições quanto ao uso deste equipamento. Camadas excessivamente finas, poços preenchidos com lamas muito condutivas e grandes diâmetros de invasão são exemplos de situações onde seu desempenho não é bom. Para que se possa minimizar estas limitações é necessário o conhecimento da resposta da sonda diante de condições conhecidas. Isto pode ser feito via modelamento analógico e/ou através de modelamento numérico. O presente trabalho consiste na montagem de um sistema de modelamento analógico capaz de simular o perfil de duas bobinas coaxiais frente a meios infinitos e camadas espessas invadidas. O meio geológico é representado através de solução salina de diferentes condutividades. O fator de escala linear é igual a 20 e o equipamento executa medidas na faixa dos 50 kHz aos 100 kHz. Este sistema será utilizado em futuros estudos envolvendo disposições de bobinas diferentes das normalmente utilizadas e em estudos de outros problemas encontrados na perfilagem de indução.

Deconvolução de perfis de poço através de rede neural recorrente

Para a indústria do petróleo, a interpretação dos perfis de poço é a principal fonte de informação sobre a presença e quantificação de hidrocarbonetos em subsuperfície. Entretanto, em duas situações as novas tecnologias, tanto em termos do processo construtivo das ferramentas, quanto da transmissão dos dados não têm justificativa econômica, ensejando a utilização de um conjunto de perfis convencionais: reavaliações de campos maduros e avaliações de campos marginais. Os procedimentos de aquisição dos perfis convencionais podem alterar o valor da propriedade física bem como a localização dos limites verticais de uma camada rochosa. Este é um antigo problema na geofísica de poço – o paradoxo entre a resolução vertical e a profundidade de investigação de uma ferramenta de perfilagem. Hoje em dia, isto é contornado através da alta tecnologia na construção das novas ferramentas, entretanto, este problema ainda persiste no caso das ferramentas convencionais como, a ferramenta de raio gama natural (GR). Apresenta-se, neste trabalho, um novo método para atenuar as alterações induzidas no perfil pela ferramenta, através da integração do clássico modelo convolucional do perfil com as redes neurais recorrentes. Assume-se que um perfil de poço pode ser representado através da operação de convolução em profundidade entre a variação da propriedade física da rocha (perfil ideal) e uma função que representa a alteração produzida sobre a propriedade física, chamada como resposta vertical da ferramenta. Assim, desenvolve-se um processamento iterativo dos perfis, o qual atua na forma da operação de deconvolução, composto por três redes neurais recorrentes. A primeira visa estimar a resposta vertical da ferramenta; a segunda procura definir os limites verticais de cada camada rochosa e a última é construída para estimar o valor real da propriedade física. Este processamento é iniciado com uma estimativa externa tanto para o perfil ideal, quanto para a resposta vertical da ferramenta. Finalmente, mostram-se as melhorias na resolução vertical e na avaliação da propriedade física produzida por esta metodologia em perfis sintéticos e em perfis reais da formação Lagunillas, bacia do Lago Maracaibo, Venezuela.

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