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

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

Caracterización petrográfica y petrofísica de la roca encajante de la Cueva del Rull (Vall d'Ebo, Alicante)

La Cueva del Rull se encuentra en el sector nororiental de la Cordillera Bética, en el denominado Prebético Externo de Alicante (Azema 1977). Regionalmente, la zona de estudio está dominada por la dinámica compresiva de los materiales calizos existentes (Cretácico Superior) afectados, desde el Mioceno Medio y durante el Mioceno Superior, por diversos movimientos tectónicos a partir de los cuales se origina la Depresión de la Vall d'Ebo. Esta fosa tectónica, cuyos bordes norte y sur quedan delimitados por fallas normales con dirección aproximada E-O, está rellena por materiales rudíticos de edad Mioceno Superior, predominantemente conglomeráticos, de espesor variable (decenas a más de 100 metros), localmente plegados y depositados sobre margas de facies “tap” (margas mal estratificadas de carácter arcillo-limoso, desagregadas y de color blanquecino en superficie, cuya edad se atribuye al Mioceno Medio).

Caracterização de reservatórios não convencionais Shale Gas na Formação Barreirinha Bacia do Amazonas

Os folhelhos pretos devonianos da Formação Barreirinha caracterizamse pela alta radioatividade na porção basal, grande extensão areal, espessura e profundidade de soterramento variável que vão de exposição na superfície até mais de 3000 m. Eles são as principais rochas geradoras do sistema petrolífero convencional da Bacia do Amazonas, e recentemente foram consideradas como promissores plays de gás não convencional. Folhelhos são geralmente caracterizados por uma matriz fechada, que faz com que sejam relativamente impermeáveis em relação ao fluxo de gás, a menos que ocorram fraturas, e dependendo das suas características geológicas e geoquímicas podem funcionar com um Sistema Petrolífero autossuficiente, atuando tanto como rocha fonte, quanto como reservatório de gás (reservatório Shale Gas). Assim, o gás natural termogênico ou biogênico gerado pode ser armazenado em folhelhos ricos em matéria orgânica na forma livre, adsorvida, ou em estado dissolvido. Em contraste com os sistemas petrolíferos convencionais, reservatórios Shale Gas, possuem mecanismos de aprisionamento e armazenamento únicos, sendo necessária a utilização de técnicas de avaliações específicas. No entanto, folhelhos prolíficos geralmente podem ser reconhecidos a partir de alguns parâmetros básicos: arquitetura geológica e sedimentar, propriedades geoquímicas e petrofísicas e composição mineralógica. Tendo em vista a carência de pesquisas de caráter descritivo, com cunho exploratório dos folhelhos geradores da Formação Barreirinha, esta dissertação tem como objetivo introduzir uma metodologia de identificação de intervalos de folhelho gerador com potencial para reservatório Shale Gas. Começando com uma investigação regional sobre o contexto geológico e sedimentar, seguido de uma avaliação abrangente enfocando as características geoquímicas, petrofísicas e litofácies dos folhelhos a partir da integração de parâmetros obtidos de perfis geofísicos de poço, análises geoquímicas e aplicação dos conceitos de Estratigrafia de Sequencia.

Proposta metodológica para avaliação preliminar de jazigos de shale gas

O potencial de um reservatório de shale gas e influenciado por um grande número de fatores, tais como a sua mineralogia e textura, o seu tipo e maturação de querogénio, a saturação de fluidos, os mecanismos de armazenamento de gás, a profundidade do reservatório e a temperatura e pressão de poros. Nesse sentido, o principal objetivo desta tese foi estabelecer uma metodologia de avaliação preliminar de potenciais jazigos de shale gas (estudo de afloramentos com base numa litoestratigrafia de alta resolução), que foi posteriormente aplicada na Formação de Vale das Fontes (Bacia Lusitânica, Portugal). Esta tese tem a particularidade de contribuir, não só para o aprofundamento da informação a nível geoquímico do local, mas também na abordagem inovadora que permitiu a caracterização petrofísica da Formação de Vale das Fontes. Para a aplicação da metodologia estabelecida, foi necessária a realização dos seguintes ensaios laboratoriais: Rock-Eval 6, picnometria de gás hélio, ensaio de resistência a compressão simples, Darcypress e a difracção de raios-X, aplicando o método de Rietveld. Os resultados obtidos na análise petrofísica mostram uma formação rochosa de baixa porosidade que segundo a classificação ISRM, e classificada como ”Resistente”, para alem de revelar comportamento dúctil e elevado índice de fragilidade. A permeabilidade média obtida situa a Formação no intervalo correspondente as permeabilidades atribuídas aos jazigos de tigh gas, indicando a necessidade de fracturação hidráulica, no caso de uma eventual exploração de hidrocarbonetos, enquanto a difracção de raios-X destaca a calcite, o quartzo e os filossilicatos como os minerais mais presentes na Formação. Do ponto de vista geoquímico, os resultados obtidos mostram que apesar do considerável teor médio de carbono orgânico total, a natureza da matéria orgânica analisada e maioritariamente imatura, composta, principalmente, por querogénio do tipo IV, o que indica a incapacidade de a formação gerar hidrocarbonetos em quantidades economicamente exploráveis.

Análise de agrupamentos dos dados de DFA oriundos de perfis elétricos de indução de poços de petróleo

The main objective of this study is to apply recently developed methods of physical-statistic to time series analysis, particularly in electrical induction s profiles of oil wells data, to study the petrophysical similarity of those wells in a spatial distribution. For this, we used the DFA method in order to know if we can or not use this technique to characterize spatially the fields. After obtain the DFA values for all wells, we applied clustering analysis. To do these tests we used the non-hierarchical method called K-means. Usually based on the Euclidean distance, the K-means consists in dividing the elements of a data matrix N in k groups, so that the similarities among elements belonging to different groups are the smallest possible. In order to test if a dataset generated by the K-means method or randomly generated datasets form spatial patterns, we created the parameter Ω (index of neighborhood). High values of Ω reveals more aggregated data and low values of Ω show scattered data or data without spatial correlation. Thus we concluded that data from the DFA of 54 wells are grouped and can be used to characterize spatial fields. Applying contour level technique we confirm the results obtained by the K-means, confirming that DFA is effective to perform spatial analysis

Determinação da distribuição de um nível de turbidito na formação Macaé no campo de Namorado através de uso de atributos sísmicos

This study was developed on the basis of petrophysics characterization, 3D seismic interpretation, and rock-wire-log integration. The results were integrated with amplitude attribute through cross plot (rock attribute x amplitude value). 3D seismic data also allowed inteipreting the top of Lagoa Feia, Macaé, and Carapebus formations, and Outeiro Member. The analyses of the maps of structures and attributes, integrated with wire-log data, suggest that Namorado turbidite level is candidate to be mapped with stratigraphic attributes like AverageTroughAmplitude, RMS velocity and Total AbsoluteAmplitude. Lagoa Feia Formation attribute maps (Dip-Azimuth and TotalEnergy) allows to interpret a low fault density at the carbonate coquina level in the Namorado field. This level is also considered one at the best seismic reflector and an important reservoir of the Campos Basin.

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.

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

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.

Caracterização dos reservatórios do grupo Macaé no Campo de Garoupa - Bacia de Campos - por meio de interpretação sísmica e perfis de poço

The study area is the Garoupa Field, located 75 km from the coast of Rio de Janeiro in bathymetric elevation of 125 m and it is the first important offshore discovery of the Campos Basin with hydrocarbon production in the Albian carbonates. The high significance of carbonate rock reservoirs for the global oil industry and its degree of complexity, combined with the fact that the Garoupa Field is a mature field, make this research project interesting to develop a re-exploratory vision of the study area. The objective of this paper is to characterize the physical properties of the Macae Group reservoirs in the Garoupa Field through well correlation techniques, petrophysical analysis and seismic interpretation (3D), in order to evaluate qualitatively and quantitatively the response of geophysical data to the presence of carbonate rock reservoirs. From the analysis of well profiles (gamma ray, resistivity, density and sonic) it was possible to identify six reservoir levels in the Quissamã Formation. As a result, the reservoirs have good correlation between the wells in the Garoupa Field and show lateral facies variation. In the seismic interpretation, structural contour maps were generated in time (ms) from the top of the Macae Group and two horizons that correspond to the reservoir levels. These maps show that the oil producing area of this field is a structural high located at the west of the interpreted seismic cube, and it was classified as a rollover limited by lateral faults forming a horst. The seismic attribute maps show a good response to the reservoirs of the Garoupa Field because they are associated with amplitude anomalies, verified by the correlation between the physical properties of the reservoir rock and seismic data through cross plots, with emphasis in a linear correlation between the resistivity profile and the values of Maximum Absolute Amplitude and RMS Amplitude attributes

Estudio de propiedades físicas de las rocas

El objetivo de este taller es mostrar a los alumnos mediante experiencias de laboratorio sencillas y de bajo coste algunas metodologías de estudio de las propiedades físicas de las rocas como material de construcción. Las propiedades físicas de las rocas determinan su uso y comportamiento tanto como materiales de construcción como soportes de obra civil. La determinación de las propiedades físicas complementan el estudio mineralógico y textural (petrográfico) de los materiales pétreos y naturales. Las propiedades petrofísicas más importantes que se abordan en este taller son el sistema poroso (porosidad); transporte de fluidos (permeabilidad, capilaridad); propiedades mecánicas (estáticas y dinámicas); la durabilidad de las rocas frente a las sales, hielo, ataque ácido, etc.

Desarrollo de electroimanes reconfigurables para uso en resonancia magnética nuclear y metrología magnética

Una técnica experimental de creciente interés en la comunidad internacional es la relaxometría magnética nuclear, la cual posee un amplio campo de aplicabilidad en la industria farmacéutica, alimentaria, petrofísica, caucho, cosmética y plástico, entre otras. Actualmente se desarrolla en el LaRTE (Laboratorio de Relaxometría y Técnicas Especiales) un prototipo alfa de un instrumento unico de resonancia magnética nuclear con campo magnético ciclado, que permitirá no solamente realizar estudios relaxométricos en muestras de hasta 35cm3, sino que posibilitará la obtención de imágenes y mediciones de difusión a diferentes valores de campo. El corazón de este aparato es un electroimán particular logrado con tecnología propia, el cual posee excelentes características al compararlo con el estado del arte. Este logro da sustento al desarrollo en marcha, el cual será posiblemente transferido a una empresa espín-off del laboratorio (Trovintek Advanced Magnetic Systems) en el futuro para un desarrollo beta. En esta dirección, este proyecto propone profundizar el desarrollo de los imanes logrados, tratando de alcanzar sistemas reconfigurables adaptivos según las características del campo requerido y según demande la aplicación (homogeneidad, velocidad de conmutación e intensidad). A su vez, estos mismos sistemas de imanes prometen excelente aplicabilidad en el campo de la metrología magnética, dirección en la cual el LaRTE ya ha comenzado a trabajar conjuntamente con la Unidad Técnica Electrónica del INTI (Córdoba). Dichos electroimanes pueden devenir, en si mismos, en dispositivos especialmente diseñados para tal aplicación.

Pirometamorfismo provocado por intrusões básicas cenozoicas em rochas da Bacia Potiguar, RN: integração de dados geológicos e petrofísicos

The Potiguar Basin is located in the Brazilian Equatorial Margin and presents sedimentary rocks affected by Cenozoic basic igneous intrusions, known as Macau Magmatism. The most prominent effect related to these intrusions is the formation of buchites, pyrometamorphic rocks that occur at very high temperatures and very low pressures in the sanidinite metamorphic facies. Through literature review, field observations, petrographic and petrophysical data, accessing the database of previous studies and results from this research, it was possible to characterize and estimate the effects produced in the thermal aureole of some hypabyssal bodies in the basin. The most relevant features associated with the intrusions are: compactation, hydraulic fracturing, partial melting and recrystallization of country rocks. According to the observed mineral occurrences, temperature of 800 to 1200 °C and pressure below 0,5 kbar were estimated at the contacts of the igneous bodies. The thermal modeling of the São João plug indicates thermal effects extending up to 150 m away from the contact and cooling time of approximately 265,000 years. After the peak of temperature, followed a cooling phase registered by remobilization and precipitation of minerals at low-temperature in faults, fractures and geodes, interpreted as derived from reactions with sedimentary rocks and metasomatic / hydrothermal fluids with abundant carbonatization and silicification.

Pirometamorfismo provocado por intrusões básicas cenozoicas em rochas da Bacia Potiguar, RN: integração de dados geológicos e petrofísicos

The Potiguar Basin is located in the Brazilian Equatorial Margin and presents sedimentary rocks affected by Cenozoic basic igneous intrusions, known as Macau Magmatism. The most prominent effect related to these intrusions is the formation of buchites, pyrometamorphic rocks that occur at very high temperatures and very low pressures in the sanidinite metamorphic facies. Through literature review, field observations, petrographic and petrophysical data, accessing the database of previous studies and results from this research, it was possible to characterize and estimate the effects produced in the thermal aureole of some hypabyssal bodies in the basin. The most relevant features associated with the intrusions are: compactation, hydraulic fracturing, partial melting and recrystallization of country rocks. According to the observed mineral occurrences, temperature of 800 to 1200 °C and pressure below 0,5 kbar were estimated at the contacts of the igneous bodies. The thermal modeling of the São João plug indicates thermal effects extending up to 150 m away from the contact and cooling time of approximately 265,000 years. After the peak of temperature, followed a cooling phase registered by remobilization and precipitation of minerals at low-temperature in faults, fractures and geodes, interpreted as derived from reactions with sedimentary rocks and metasomatic / hydrothermal fluids with abundant carbonatization and silicification.