Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-14T12:15:27Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoRespostaDiferentes.pdf: 6343655 bytes, checksum: 9d709e578bbcbd072226e615633d6d33 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:46:31Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-20T16:12:29Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoRespostaDiferentes.pdf: 6343655 bytes, checksum: 9d709e578bbcbd072226e615633d6d33 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T16:56:35Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoRespostaDiferentes.pdf: 6343655 bytes, checksum: 9d709e578bbcbd072226e615633d6d33 (MD5) / Made available in DSpace on 2014-08-22T16:56:35Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoRespostaDiferentes.pdf: 6343655 bytes, checksum: 9d709e578bbcbd072226e615633d6d33 (MD5) Previous issue date: 1991 / 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. / Induction logs are frequently run in the drill holes to distinguish the zones containing hydrocarbons from those containing conate water, based on the conductivity contrast. Besides localizing these zones, the induction log can be used to calculate the hydrocarbon reserves reservoir evaluation. Most commercial induction logs use a coaxial coil system, which have poor response in thin beds, deep invasion and conductive muds. A comparative study of E.M. responses both coaxial and a coplanar coil systems was carried out for different geological situations to test the two systems in conditions where the coaxial system responds poorly. For this study, an existing analog model was modified. The model source frequency is in the range of 50 to 100 kHz. A scale factor of 20 was chosen for this study. Solutions of varying electrical conductivities, obtained using different concentrations of amonium chloride contained in various compartments in a fiber glass tank, simulated the geological models. The measurements obtained demonstrate that the coplanar coil system has a better vertical resolution than the coaxial coil system normaly used in induction logging.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-07-09T12:29:50Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:34:20Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T15:00:46Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-19T13:25:48Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) / Made available in DSpace on 2014-09-19T13:25:48Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EstudoEstabilidadeMecanica.pdf: 11477602 bytes, checksum: 6af6cc1a1e5c0885abc416d0e5bdc521 (MD5) Previous issue date: 1992 / FADESP - Fundação de Amparo e Desenvolvimento da Pesquisa / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / FINEP - Financiadora de Estudos e Projetos / A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação. / The stability of boreholes can be evaluated from the rocks' mechanical properties, which are compared to the acting stress field using an appropriate failure criterion. The elastic parameters can be measured on rock samples in the laboratory (static parameters) or can be calculated from acoustic velocities and the formation bulk density (dynamic parameters). Dynamic elastic parameters can be obtained from wireline logs, but in poorly consolidated rocks the arrival of the shear wave is often masked by the Stoneley wave, which may travel faster and has a higher amplitude. The giant Campos basin offshore Brazil contains such rocks. We performed laboratory measurements of acoustic velocities on samples from this basin and we use them to correct the shear velocities obtained from waveform processing. Based on this we calculate elastic properties (e.g. Poisson's ratio, bulk modulus) of the rocks which we assume to be homogeneous and isotropic. Mechanical properties such as the shear, compressional and tensional strength are obtained from empirical relations with mineralogy and the elastic parameters. The stress field was evaluated using an uniaxial strain model and actual formation pressure data. From this, near-wellbore stresses were calculated using the modified Kirsch equations. The stress concentrations around the wellbore are illustrated for a number of representative cases. Wellbore failure occurs when these local stresses exceed the formation strength. Three different failure criteria were tested: The Mohr-Coulomb, the Griffith and the Fairhurst criterion. They respectively represent failure due to shear, tension and a combination of compression and tension. Actual wellbore instability was determined from the differential caliper measurement. Two sequences were analyzed, both consisting of intercalations of sands and shales. We found that the Fairhurst criterion offered the best match with the actually observed wellbore failure, and we propose it as a method to predict borehole instability for this basin.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-18T14:01:22Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeterminacaoEstatisticaContatos.pdf: 6587528 bytes, checksum: c394cb002b28d37e58355d6f8383f245 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:18:34Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T15:14:41Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeterminacaoEstatisticaContatos.pdf: 6587528 bytes, checksum: c394cb002b28d37e58355d6f8383f245 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-22T11:22:23Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeterminacaoEstatisticaContatos.pdf: 6587528 bytes, checksum: c394cb002b28d37e58355d6f8383f245 (MD5) / Made available in DSpace on 2014-09-22T11:22:23Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DeterminacaoEstatisticaContatos.pdf: 6587528 bytes, checksum: c394cb002b28d37e58355d6f8383f245 (MD5) Previous issue date: 2001 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / 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 (RT), 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. / In oil reservoirs, the effect of gravity naturally produces the fluid segregation. Due to capillary forces, there is no defined surface of separation between two differents fluids. However, it is common in petroleum engineering to admit a conventional fluid separation surface, called fluid interface or fluid contact. The depth location of the fluid contacts (oil-water, gas-water and gas-oil) inside the reservoir rocks, without the complete well log interpretation or the accomplishment of the direct procedures involved in formation tests, is a petroleum industry problem. The knowledge of this parameter can be used for well completation procedures and for positioning perforating services in vertical wells. Another application of fluid contact mapping can be to control the drilling bit in directional and horizontal wells, where such information is used to maintain the well axis inside the hydrocarbon zone in order to avoid water production. We present a methodology which can identify and locate fluid contacts, through an application of the multivariate statistical technique called discriminant analysis. For clastics deposits, with sand-shale sequences, discriminant analysis may provide the indication of lithology and the apparent thicknesses of the reservoirs. The well logs applied for the evaluation of this methodology (resistivity (RT), gamma ray (RG), density (ρb), neutronic porosity (ΦN), caliper (Cal) and shaleness (Vcla) are from wells in Lake Maracaibo, Venezuela.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-03T14:22:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Rejected by Albirene Aires (albireneufpa@gmail.com), reason: Favor incluir as palavras chave on 2014-12-05T11:45:19Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-10T00:14:49Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-12-10T17:20:36Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) / Made available in DSpace on 2014-12-10T17:20:36Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPermissividadeDieletrica.pdf: 9420207 bytes, checksum: 28f0ebee4414f2a902ff9e96c0b0b87e (MD5) Previous issue date: 1990 / A ferramenta de propagação eletromagnética (EPT) fornece o tempo de propagação (T_pl) e a atenuação (A) de uma onda eletromagnética que se propaga num meio com perdas. Estas respostas da EPT são funções da permissividade dielétrica do meio. Existem vários modelos e fórmulas de misturas sobre a permissividade dielétrica de rochas reservatório que podem ser utilizados na interpretação da ferramenta de alta frequência. No entanto, as fórmulas de mistura não consideram a distribuição e a geometria do espaço poroso, e estes parâmetros são essenciais para que sejam obtidas respostas dielétricas mais próximas de uma rocha real. Foi selecionado um modelo baseado nos parâmetros descritos acima e este foi aplicado à dados dielétricos disponíveis na literatura. Foi obtida uma boa concordância entre as curvas teóricas e os dados experimentais, comprovando assim que a distribuição e a geometria dos poros têm que ser levadas em conta no desenvolvimento de um modelo realista. Foram conseguidas também funções de distribuição de razão de aspecto de poros, através das quais geramos várias curvas relacionando as respostas da EPT com diversas saturações de óleo/gás. Estas curvas foram aplicadas na análise de perfis. Como o modelo selecionado ajusta-se bem aos dados dielétricos disponíveis na literatura, torna-se atraente aplicá-lo à dados experimentais obtidos em rochas de campos brasileiros produtores de hidrocarbonetos para interpretação da EPT corrida em poços destes campos petrolíferos. / The electromagnetic propagation tool (EPT) provides the propagation time (T_pl) and the attenuation (A) of an electromagnetic wave propagating in a lossy medium. These EPT responses depend on the dielectric permittivity of the medium. There are several models and mixing equations concerning the dielectric permittivity of reservoir rocks that can be used in the interpretation of the high frequency tool. However, the mixing equations do not take into account the distribution and the geometry of the pore space, and these parameters are essential to obtaining dielectric responses approximating a true rock. A model based on the parameters described above was selected and this was applied to dielectric data available in the literature. A good agreement was reached between the theoretical curves and experimental data, confirming that the distribution and geometry of the pore space must be considered in the development of a realistic model. Aspect ratio distribution functions of the pores were also obtained, which were used for generating several curves relating the EPT responses to various oil/gas saturations. These curves were applied to the log analysis. The selected model fit the dielectric data available in the literature reasonably well, thus, making it suitable for application to experimental data of rock from Brazilian producing fields for the interpretation of the EPT in these fields.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-03T14:21:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Rejected by Albirene Aires (albireneufpa@gmail.com), reason: Favor incluir as palavras chave on 2014-12-05T11:44:55Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-12-10T00:15:50Z No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-12-10T17:23:04Z (GMT) No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) / Made available in DSpace on 2014-12-10T17:23:04Z (GMT). No. of bitstreams: 2 license_rdf: 22974 bytes, checksum: 99c771d9f0b9c46790009b9874d49253 (MD5) Dissertacao_ModelamentoPerfisInducao.pdf: 7044960 bytes, checksum: d5bab32aace9336f20b8ad536792da82 (MD5) Previous issue date: 1990 / 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. / This thesis applies an algorithm that allows modeling induction sondes (2C40-6FF40-ILD) responses in multi layered configuration. The algorithm is based on the exact solution of the Green's problem to obtain the electromagnetic field in layered configuration with the source located in the layers. The code, developed for horizontal layers, can be used to develop advanced codes, because there are not codes available for research. The codes available are commercial, but these are very expensive and primarily useful for production processing. The code was applied to obtain ILD sonde response in laminated reservoirs, where the sonde resolution is reduced. Using this code, the conductivity from the sonde response, can be expressed by the sum of the conductivities of the lithologic constituints of the reservoirs, weighed by their volume contribution. This relation was defined for reservoir type not yet analyzed. The models are for reservoirs of silicified sandstone/sandstone/shale where the resistivity of the layer of interest is superior to the resistivity of one layer and inferior to the resistivity of another layer, and for reservoirs of silidfied sandstone/sandstone where the layer of interest has the lower resistivity. These models allow the establishment of better estimates of resistivity (Rt) in layered reservoirs, and consequently, to the hydrocarbon saturation and producibility.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-02T15:50:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T13:51:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Made available in DSpace on 2014-06-30T13:51:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) Previous issue date: 1993 / PETROBRAS - Petróleo Brasileiro S.A. / A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P1, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm3 e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas. / The identification and description of lithological parameters of a formation are essential in the evaluation of complex formations. Based on this, the combination of the nuclear tool response in uncased wells has been used systematically. The resultant logs can be considered as the interaction between two distinct phases: • The radiation transport phase from a source to one or more detectors through the formation. • The detection phase that consists of the collection of radiation, its transformation into current pulses, and the spectral distribuition of these pulses. As the presence of the detector does not strongly affect the radiation transport result, each phase can be simulated independent of the other, which allows us to introduce a new type of model in which the transport phase and the detection phase are uncoupled. In this work, the final response is simulated combining transport numerical solutions with a library of the detector responses to different incident energies and for each specific source - detector array. The radiation is calculated by the Finite Elements Method (FEM), as a 2½-D scalar flux derived from the numerical solution of the multigroup diffusion approximation of the Boltzmann transport equation in phase space. This is known as P1 approximation, where the variable direction is expanded in terms of the Legendre orthogonal polymonials, leading to the dimensionality reduction of the problem in such a way as to let it be more consistent with the FEM, where the flux depends only on the spatial variable and the physical properties of the formation. The NaI(Tl) response function is obtained separately by the Monte Carlo method (MC) where the life of a particle within the scintillator crystal is reconstructed simulating interaction by interaction the position, direction and energy of the different particles using a random number technique with associated appropriate probabilities laws. Each type of interaction (e.g., Rayleigh, Photo-electric, Compton and Pair production) is determined similarly and the simulation is concluded when the detector response functions are convolved with the scalar flux. The final response is the pulse-height spectrum of the simulated system. From this spectrum, a set of channels called detection windows are then selected. The count rates in each window show different dependencies on density and lithology. This fact allows one to use a combination of these windows in determining the density and photoelectric absorption factor of the formation. According to the method developed in this work, the logs in both thin and thick layers can be simulated. The performance of the method has been tested in complex formations, mainly where the presence of clay minerals, feldspars and micas have produced considerable effects sufficient to perturb the final response of the sonde. The results show that it is possible to identify physical and lithological parameters in formations having densities between 1.8 and 4.0 g/cm3 and photoelectric absorption factors in the interval of 1.5 to 5.0 barns/e-. The concentrations of Potassium, Uranium and Thorium can be obtained through the introduction of a new system of calibration which corrects the effects due to high variances and negative correlations observed on the mass concentration of Uranium and Potassium. In the simulation of the CNL response using the Tittle polynomial regression algorithm, it is verified that due to the limited vertical resolution of this sonde, the porosity value is poorly measured for most layers of thickness less than the source - far detector spacing, thus it has application only in thick layers. A new method was developed to solve this problem; the contribution of the relative area of each layer within the maximum information zone is determined. Thus, this neutron porosity makes possible an in-depth evaluation of expected CNL porosity-lithology response, convolving that area factor with the local formation porosity index, considering only thick layers. The presence of perturbating minerals is solved by considering the formation as formed by a predominant base matrix mineral, totally saturated by fresh water; the rest of the components are then considered as a perturbation of this base case. These results enable the calculation of synthetic well logs that can be used in inversion schemes in order to get a more detailed quantitative evaluation of complex formations.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T12:28:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:13:21Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-14T14:23:37Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:04:51Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Made available in DSpace on 2014-08-22T17:04:51Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) Previous issue date: 1991 / Utilizando-se dados do campo de Camorim (Bacia de Sergipe-Alagoas), foi testado e aplicado um conjunto de técnicas estatísticas multivariantes (análises de agrupamentos, de componentes principais e discriminante) no intuito de identificar as fácies, previamente definidas em poços testemunhados, por meio dos perfis, viabilizando o reconhecimento das mesmas nos demais poços não testemunhados da área. A segunda etapa do processo de determinação das fácies consistiu no emprego de métodos auxiliares (análises composicional e de seqüência de fácies), que, combinados com as técnicas multivariantes, propiciaram melhores resultados na calibração rocha-perfil. A determinação das fácies, uma vez estabelecida, possibilitou o refinamento do processo de avaliação de formações ao viabilizar o exame de cada fácies-reservatório, isoladamente. Assim, esse procedimento tornou possível a escolha, para cada litologia, dos parâmetros utilizados na interpretação dos perfis ao mesmo tempo em que permitiu a totalização em separado dos valores de espessura, porosidade e saturação dos fluidos, bem como a adoção de diferentes valores de corte (cut-offs) para cada grupo considerado. Outras aplicações incluíram a melhoria na estimativa da porosidade e da permeabilidade, a adaptação de algoritmos para o cálculo preliminar de porosidade, a confecção de mapas de fácies e a geração automática de seções estratigráficas. Finalmente, foram destacadas a perspectiva de integração desse estudo com sistemas estatísticos de descrição de reservatórios, outras técnicas de determinação de fácies em desenvolvimento e a retomada da utilização de métodos estatísticos multivariantes em dados de perfis, como ferramenta de exploração. / A series of multivariate statistical techniques (cluster, principal component and discriminant analysis) was tested and applied to well log data from the Camorim field (offshore Sergipe State, Brazil) in order to identify facies previously defined through core description. The second step in the process of facies determination was supported by auxiliary methods (compositional and sequence facies analysis), which produced better results in the calibration of rock-log data, when combined with the multivariate techniques. The facies determination, once established, permits the refinement of the process of formation evaluation, enabling the examination of each facies separately. This procedure made it possible to choose, for each lithology, the parameters used in log interpretation. In parallel, this process allowed the summation of thickness, porosity, fluid saturation and the adoption of different cut-offs for each group, separately. Other applications included: enhancement in the estimation of porosity and permeability, the adaptation of algorithms designed for fast porosity estimation, the mapping of variables useful in the characterization of the vertical variability of the reservoir rocks and the automatic generation of stratigraphic sections. Finally, the possibility of integration of the work from this study with statistical systems of reservoir description, other facies determination techniques currently being developed and the utilization of multivariate statistical methods in well log data, as an auxiliary exploratory tool, were illustrated.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-21T15:39:26Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoPorosidadePerfis.pdf: 1768673 bytes, checksum: ba48f09908920b586e020f6650adeb26 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:15:33Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-12T15:24:10Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoPorosidadePerfis.pdf: 1768673 bytes, checksum: ba48f09908920b586e020f6650adeb26 (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-18T12:47:37Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoPorosidadePerfis.pdf: 1768673 bytes, checksum: ba48f09908920b586e020f6650adeb26 (MD5) / Made available in DSpace on 2014-09-18T12:47:37Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ImageamentoPorosidadePerfis.pdf: 1768673 bytes, checksum: ba48f09908920b586e020f6650adeb26 (MD5) Previous issue date: 2004 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / O imageamento da porosidade é uma representação gráfica da distribuição lateral da porosidade da rocha, estimada a partir de dados de perfis geofísicos de poço. Apresenta-se aqui uma metodologia para produzir esta imagem geológica, totalmente independente da intervenção do intérprete, através de um algoritmo, dito, interpretativo baseado em dois tipos de redes neurais artificiais. A primeira parte do algoritmo baseia-se em uma rede neural com camada competitiva e é construído para realizar uma interpretação automática do clássico gráfico o Pb - ΦN, produzindo um zoneamento do perfil e a estimativa da porosidade. A segunda parte baseia-se em uma rede neural com função de base radial, projetado para realizar uma integração espacial dos dados, a qual pode ser dividida em duas etapas. A primeira etapa refere-se à correlação de perfis de poço e a segunda à produção de uma estimativa da distribuição lateral da porosidade. Esta metodologia ajudará o intérprete na definição do modelo geológico do reservatório e, talvez o mais importante, o ajudará a desenvolver de um modo mais eficiente as estratégias para o desenvolvimento dos campos de óleo e gás. Os resultados ou as imagens da porosidade são bastante similares às seções geológicas convencionais, especialmente em um ambiente deposicional simples dominado por clásticos, onde um mapa de cores, escalonado em unidades de porosidade aparente para as argilas e efetiva para os arenitos, mostra a variação da porosidade e a disposição geométrica das camadas geológicas ao longo da seção. Esta metodologia é aplicada em dados reais da Formação Lagunillas, na Bacia do Lago Maracaibo, Venezuela. / Porosity images are graphical representations of the lateral distribution of rock porosity estimated from well log data. We present a methodology to produce this geological image entirely independent of interpreter intervention, with an interpretative algorithm approach, which is based on two types of artificial neural networks. The first is based on neural competitive layer and is constructed to perform an automatic interpretation of the classical Pb - ΦN cross-plot, which produces the log zonation and porosity estimation. The second is a feed-forward neural network with radial basis function designed to perform a spatial data integration, which can be divided in two steps. The first refers to well log correlation and the second produces the estimation of lateral porosity distribution. This methodology should aid the interpreter in defining the reservoir geological model, and, perhaps more importantly, it should help him to efficiently develop strategies for oil or gas field development. The results or porosity images are very similar to conventional geological cross-sections, especially in a depositional setting dominated by clastics, where a color map scaled in porosity units illustrates the porosity distribution and the geometric disposition of geological layers along the section. The methodology is applied over actual well log data from the Lagunillas Formation, in the Lake Maracaibo basin, located in western Venezuela.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T13:37:50Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoReservatoriosCampo.pdf: 10234818 bytes, checksum: 370a45261bcd9980b1496a9abf14d80c (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:22:28Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T17:12:04Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoReservatoriosCampo.pdf: 10234818 bytes, checksum: 370a45261bcd9980b1496a9abf14d80c (MD5) / Approved for entry into archive by Ana Rosa Silva (arosa@ufpa.br) on 2014-09-18T17:26:23Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoReservatoriosCampo.pdf: 10234818 bytes, checksum: 370a45261bcd9980b1496a9abf14d80c (MD5) / Made available in DSpace on 2014-09-18T17:26:23Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoReservatoriosCampo.pdf: 10234818 bytes, checksum: 370a45261bcd9980b1496a9abf14d80c (MD5) Previous issue date: 1990 / 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. / A new model for quantitative evaluation of shaly sand is proposed based on bibliography and clay mineral distribution in the pores of reservoir rocks. This model takes in account for two extra conductivity types besides the normal electrolytic conductivity one present in porous rocks. One conductivity type is due to surface conductance of clay minerals, which increases or decreases as result of their capability of cationic exchange. The second conductivity type behaves independently as a parallel system due to the microporous network formed by the clay minerals lining or bridging the grains of a shaly rock. The new equation was tested against several other equations and showed the best correlation with the Waxman & Smits equation. Both equations, the proposed one and the Waxman & Smits, have a similar approach to the clay conduction system. The new model was used in three different zones of the Ubarana oilfield (Potiguar basin, Brazil) It was verified that: (1) in two low resistivity zones (there was the possibility of) oil production was possible because the resistivities were supposedly affected by highly conductive minerals; and (2) one effective oil tested and producing tone of highly resistivity. The water saturations calculated with the new equation for the two low resistivity zones were above 50%, suggesting subsaturation and a continous resistivity phase through the salt water (about 85,000 ppm of NaCl) present in the pores. The practical results shown herein suggest that the proposed equation is theoretically consistent and must be exhaustively used in future evaluation projects of the Ubarana and similar oilfields in Brazilian sedimentary basins.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-20T12:04:23Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitosDelawareGroningen.pdf: 16863178 bytes, checksum: f72043e0b822489939d664df6c8bdf28 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:17:35Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-22T14:10:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitosDelawareGroningen.pdf: 16863178 bytes, checksum: f72043e0b822489939d664df6c8bdf28 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T15:19:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitosDelawareGroningen.pdf: 16863178 bytes, checksum: f72043e0b822489939d664df6c8bdf28 (MD5) / Made available in DSpace on 2014-08-22T15:19:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_EfeitosDelawareGroningen.pdf: 16863178 bytes, checksum: f72043e0b822489939d664df6c8bdf28 (MD5) Previous issue date: 1993 / Os efeitos Delaware e Groningen são dois tipos de anomalia que afetam ferramentas de eletrodos para perfilagem de resistividade. Ambos os efeitos ocorrem quando há uma camada muito resistiva, como anidrita ou halita, acima do(s) reservatório(s), produzindo um gradiente de resistividade muito similar ao produzido por um contato óleo-água. Os erros de interpretação produzidos têm ocasionado prejuízos consideráveis à indústria de petróleo. A PETROBRÁS, em particular, tem enfrentado problemas ocasionados pelo efeito Groningen sobre perfis obtidos em bacias paleozóicas da região norte do Brasil. Neste trabalho adaptamos, com avanços, uma metodologia desenvolvida por LOVELL (1990), baseada na equação de Helmholtz para HΦ, para modelagem dos efeitos Delaware e Groningen. Solucionamos esta equação por elementos finitos triangulares e retangulares. O sistema linear gerado pelo método de elementos finitos é resolvido por gradiente bi-conjugado pré-condicionado, sendo este pré-condicionador obtido por decomposição LU (Low Up) da matriz de stiffness. As voltagens são calculadas por um algoritmo, mais preciso, recentemente desenvolvido. Os perfis são gerados por um novo algoritmo envolvendo uma sucessiva troca de resistividade de subdomínios. Este procedimento permite obter cada nova matriz de stiffness a partir da anterior pelo cálculo, muito mais rápido, da variação dessa matriz. Este método permite ainda, acelerar a solução iterativa pelo uso da solução na posição anterior da ferramenta. Finalmente geramos perfis sintéticos afetados por cada um dos efeitos para um modelo da ferramenta Dual Laterolog. / The Delaware and Groningen effects are two different kinds of anomaly affecting electrode type resistivity logging tools. Both effects take place when there is a high-resistivity bed, like halite or anhydrite, above the reservoir(s), producing a gradient on the resistivity log, which resembles an oil-water contact. The misinterpretations therefore produced have caused considerable losses to oil industry. PETROBRAS, in particular, has faced problems caused by Groningen effect on logs obtained in paleozoic basins of northern Brazil. In this work we have adapted and improved upon a methodology developed by LOVELL (1990) based on Helmholtz's equation for HΦ for modeling Delaware and Groningen effects. We solve this equation by triangular and rectangular finite elements. The finite element linear system is solved by preconditioned bi-conjugate gradient, the preconditioner being obtained by incomplete LU (Low Up) decomposition. Voltages are calculated by a recently developed, and more precise, algorithm. Logs are generated by a new algorithm involving successive switching of subdomain resistivity. This procedure allows the computation of each new stiffness matrix from the previous one by the very fast computation of matrix variation. Moreover, this method allows rapid iterative solution by using the solution from the previous tool position. Finally we compute resistivity logs subject to each effect for a model of the Dual Laterolog Tool.