Aproveitamento dos resíduos gerados pelo processamento industrial/artesanal de pedra-sabão da região de Ouro Preto

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Desenvolvimento de cerâmicas de mulita a partir de alumina, ácido silícico e aerosil

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Impactos ambientais provocados pela exploração de granito na Serra da Meruoca - CE

Pós-graduação em Geografia - IGCE

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

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

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

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

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.

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

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

Aplicação dos métodos elétrico e eletricomagnético para o estudo hidrogeológico no município de Altamira-Pará

Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea.

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.

Comparação dos filtros de velocidade e do operador WHLP-CRS na atenuação de múltiplas

A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.

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

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

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

Geoquímica das águas do médio e baixo rio Madeira e seus principais tributários - Amazonas - Brasil

Este trabalho teve como objetivo estudar as águas do rio Madeira e seus principais tributários entre a cidade de Humaitá e sua foz no rio Amazonas. Foram analisados pH, condutividade, turbidez, íons maiores, elementos traço e isótopos de Sr nos períodos de seca, cheia e transição para a seca entre 2009 e 2010. As águas do Madeira, classificadas com brancas, são bicarbonatadas-cálcicas, têm pH entre 5 e 6 e são mais concentradas que as dos tributários. Estes têm águas de cor preta, mais ácidas e quimicamente heterogêneas, os da margem esquerda são quimicamente mais semelhantes as do Madeira, enquanto os da margem direita têm alta concentração em SiO2. Os cátions, Cl- e NO3- são mais concentrados na cheia o que sugere influência do solo, da vegetação e da composição da água da chuva (Cl-), enquanto HCO3-, SO42-, Al, Br e P, com maiores concentrações na seca, devem estar relacionados com a química das rochas. A SiO2 e os elementos terras raras (ETR) com concentrações elevadas na seca e na cheia, estão associados tanto a vegetação e ao solo como as rochas. A interação desses fatores é a causa da heterogeneidade química das águas. Contudo, a semelhança entre as águas dos tributários da margem esquerda e as do Madeira são consequência das rochas dos Andes serem a fonte dos sedimentos cenozóicos percolados por elas, enquanto a química das águas dos tributários da margem direita retrata a estabilidade tectônica, o intenso intemperismo e a baixa taxa de erosão das rochas do cráton Amazônico.

The ceramic artifacts in archaeological black earth (terra preta) from Lower Amazon Region, Brazil: chemistry and geochemical evolution

Neste trabalho realizou-se a caracterização química de fragmentos de artefatos cerâmicos encontrados em sítios arqueológicos com terra preta no Baixo Amazonas (Cachoeira-Porteira, Pará, Brasil), representativos da cultura Konduri (de 900 a 400 anos AP). Esses fragmentos são constituídos de SiO₂, Al₂O₃, Fe₂O₃, Na₂O e P₂O₅, sendo que SiO₂ e Al₂O₃, juntos, perfazem mais de 80 % em peso. Os teores de P₂O,₅ são relativamente elevados (2,37 % em média) sob a forma de (Al,Fe)-fosfatos, incomuns em cerâmicas vermelhas primitivas, mas encontrados em algumas cerâmicas arqueológicas egípcias e romanas. As concentrações dos elementos traços são comparáveis ou mesmo inferiores ao nível crustal, embora a composição total seja próxima a mesma. A composição química (exceto P₂O₅) em conjunto com os dados mineralógicos e texturais indicam material saprolítico derivado de rochas ígneas félsicas ou rochas sedimentares como matéria-prima das cerâmicas. Os teores de K, Ca e Na mostram que os feldspatos e fragmentos de rochas foram adicionados ao material argiloso, como sugerido pela mineralogia. Os altos teores de sílica respondem pela presença de cauixi, cariapé e/ou areias quartzosas. Fósforo deve ter sido incorporadoà matriz argilosa da cerâmica, quando do cozimento de alimentos nos vasos cerâmicos, e ainda, em parte, durante a formação do perfil de solo tipo ABE sobre Latossolos Amarelos. A matéria prima e os temperos (cauixi, cariapé, rochas trituradas e fragmentos de vasos cerâmicos descartados) encontram-se disponíveis próximos aos sítios até a atualidade, e, portanto foram a área fonte dos mesmos para a confecção dos artefatos cerâmicos.

The ceramic artifacts in archaeological black earth (terra preta) from lower Amazon region, Brazil: Mineralogy

Sítios arqueológicos com Terra Preta, denominados de Terra Preta de Índio ou ainda Terra Preta Arqueológica (TPA) são muito freqüentes na Amazônia. As TPA geralmente contém fragmentos de vasos cerâmicos, por vezes abundantes, além de líticos, que são materiais de grande importância para os estudos arqueológicos. Para consubstanciar esses estudos, realizou-se pesquisas mineralógicas e químicas em fragmentos cerâmicos provenientes de dois sítios arqueológicos da região de Cachoeira-Porteira, Estado do Pará. Os fragmentos foram classificados segundo seus principais temperos em: cauixi, cariapé, areia+feldspatos e caco de vaso cerâmico. Mineralogicamente são compostos de quartzo, minerais de argila calcinados (especialmente caulinita), feldspatos (albita e microclínio), hematita, goethita, maghemita, variscita-estrengita, fosfatos amorfos, anatásio, e raramente apatita, rhabdophana e óxidos de Mn e Ba. Cauixi e cariapé são componentes orgânicos silicosos e amorfos a DRX. A composição mineralógica e a morfologia dos seus grãos indicam saprólito (material argiloso rico em quartzo) derivado de rochas ígneas félsicas de granulação fina ou rochas sedimentares ricas em argilominerais como matéria-prima dos vasos cerâmicos. Neste material argiloso cauixi, cariapé e/ou areias, ricas em sílica, foram intencionalmente adicionados. O elevado conteúdo de fosfatos de Al-Fe, amorfos ou como de baixa cristalinidade, originou-se a partir do contato entre a matriz argilosa da parede do vaso cerâmico com a solução aquosa quente durante o cozimento diário de alimentos de origem animal (principal fonte de fósforo). A cristalização dos fosfatos deve ter prosseguida mesmo depois que os vasos foram descartados, e juntos com os restos de matéria orgânica vegetal e animal incorporaram-se aos solos residuais. Participaram desta forma na formação dos solos tipo TPA.