Adequações tecnológicas no beneficiamento do caranguejo-uçá Ucides cordatus (Linnaeus, 1763) desenvolvido na comunidade de Caratateua, município de Bragança, estado do Pará, Brasil

No Estado do Pará, o beneficiamento do caranguejo-uçá ocorre de forma totalmente artesanal, sendo o seu produto frequentemente associado a eventos de intoxicação alimentar. Neste trabalho o autor tem por finalidade propor adequações tecnológicas no processamento deste crustáceo desenvolvido na comunidade de Caratateua, município de Bragança/Pará. Foram realizadas excursões mensais à comunidade de abril/2008 a maio/2009 com aplicação de 120 entrevistas direcionadas aos principais atores sociais envolvidos na atividade, constando de aspectos econômicos, sociais e tecnológicos, além de análises laboratoriais (composição química e microbiologia), observações de campo e revisão de literatura. Constatou-se que a quantidade de caranguejo desembarcada e o tamanho médio dos indivíduos vêm diminuindo nos últimos anos. A captura é realizada exclusivamente por indivíduos do sexo masculino, sendo o gancho o principal apetrecho utilizado. O beneficiamento é praticado principalmente por mulheres e crianças, porém homens eventualmente também o fazem. Não há uma padronização do processamento e as condições higiênico-sanitárias dos locais de beneficiamento são na grande maioria das vezes inadequadas, podendo este ser realizado no interior da residência, em “puxadinhas” ou nos quintais, geralmente na presença de animais domésticos, insetos e sem a adequada assepsia de manipuladores, equipamentos e superfícies de contato. Desta forma, conclui-se que tal realidade é ocasionada principalmente pelo desordenamento da pescaria e falta de preocupação dos catadores quanto à qualidade do produto, sendo de fundamental importância a criação de uma legislação específica para a carne de caranguejo e a adoção de boas práticas de fabricação levando em consideração a realidade local visando agregação de valor ao produto e maior segurança alimentar aos consumidores.

Estratégias de reprodução do capital e as novas espacialidades urbanas: o comércio de auto-serviço em Uberlândia (MG)

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Avaliação dos perigos microbiológicos em uma indústria de beneficiamento de pescados e sugestão de um sistema de gestão de qualidade

Entre os produtos de origem animal os peixes são um dos mais susceptíveis à deterioração e podem representar risco para a saúde, se estiverem contaminados por micro-organismos patogênicos. Desta forma, as indústrias de pescado precisam ter cuidados especiais durante todas as etapas do processamento. Este estudo avaliou os perigos microbiológicos na linha de beneficiamento e produtos em uma indústria de pescado e propõe a sua adequação aos pré-requisitos necessários para a implantação do sistema de Análise de Perigos e Pontos Críticos de Controle (APPCC). Foram realizadas análises microbiológicas e físico-químicas em 144 amostras (filés e postas) de peixes congelados de diferentes espécies coletadas diretamente na indústria; foram também submetidas a análises microbiológicas 12 amostras de gelo, 24 amostras de água e 340 amostras de superfícies (equipamentos, utensílios e luvas de manipuladores). As Boas Práticas de Fabricação (BPF) foram avaliadas por meio de uma lista de verificação aplicada durante cinco vistorias na indústria. Foram realizados treinamentos e avaliadas as temperaturas de matérias-primas e de produtos semi-elaborados. Foi então proposto um plano APPCC para a linha de processamento de filé de peixe congelado. Entre os patógenos avaliados somente Estafilococos coagulase positiva foi detectado em 2,7% das amostras de peixes, resultado que está dentro do limite estabelecido na legislação brasileira. entretanto valores acima do estabelecido na legislação nacional para pl-! e B VT foram detectados em 18,1% e em 3,5% das amostras, respectivamente. Nas amostras de gelo e água não foram detectadas a presença de coliformes a 35°C e a 45°C. Do total de amostras de equipamentos e utensílios analisadas, 50% e 81,82% estavam dentro da recomendação estabelecida pela OPAS, de 1,7 log UFC/cm2 para mesófilos aeróbios e ausência de coliformes a 45°C, no primeiro e segundo período de coleta, respectivamente. As luvas, no primeiro período de coleta apresentaram contagem acima de 4 log UFC/luva de mesófilos aeróbios em 76,67% das amostras e no segundo período em apenas 21,67%. a avaliação das BPF nas duas primeiras vistorias a indústria foi classificada como de Risco Alto, por ter apresentado percentual de conformidades abaixo de 50%, e como de Risco Médio nas vistorias subsequentes, por ter atendido de 51 a 75% dos itens preconizados na lista de verificação. Com os resultados obtidos foi possível constatar que as BPF não estão totalmente implantadas na indústria analisada, detectando-se diversas não conformidades. A elaboração do plano APPCC seguiu as etapas: formação da equipe, descrição do produto, elaboração do fluxograma de produção, análise dos perigos, identificação dos Pontos Críticos de Controle (PCCs), estabelecimento de limites críticos e ações corretivas, e procedimentos de monitoração, registros e de verificação. A implementação do plano APPCC proposto é aplicável e pode contribuir, significativamente, para a elaboração de produtos de pescado com melhor qualidade e segurança alimentar.

Estudo evolutivo do sistema morfoclimático e morfotectônico da bacia do Rio Verde (MG), sudeste do Brasil

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

Avaliação da rugosidade superficial e resistência à flexão de três tipos de resina composta indireta após variação de tratamento da superfície

Neste estudo avaliou-se a rugosidade superficial em corpos-de-prova confeccionados com três tipos de resinas indicadas para confecção de restaurações indiretas, SR Adoro (Ivoclar-Vivadent) (GA), Signum Matrix (Heraeus-Kulzer) (GSM) e Signum⁺ (Heraeus-Kulzer) (GSP), assim como a resistência à flexão após serem utilizados diferentes tratamentos de superfície com o propósito de se verificar a influência dos mesmos na resistência mecânica dos materiais. Obteve-se 120 corpos-de-prova de cada tipo de resina de acordo com as recomendações dos fabricantes, os quais foram divididos e tratados com: A- Jateamento com óxido de alumínio a 50μm por 6 s a 60-80 libras (Grupo Controle); B- Jateamento e silanização com Monobond (Ivoclar Vivadent) e Excite DSC (Ivoclar Vivadent); C- Jateamento e condicionamento com ácido ortofosfórico a 37% por 3 minutos e D- Jateamento, condicionamento e silanização. Após a execução do tratamento superficial, analisou-se os corpos-de-prova em rugosímetro Mytutoyo no parâmetro Ra e teste de flexão por compressão em três pontos em máquina de ensaio universal Kratos. Submeteu-se os dados observados à análise estatística, tendo como nível de significância 5% para a construção dos resultados. Os resultados mostraram um comportamento bastante coerente com o que é visto na literatura com relação à composição das resinas estudadas, os quais apresentaram inicialmente lisura de superfície similar (p>0.05). Com a aplicação dos tratamentos de superfície, houve um comportamento bastante diversificado entre os materiais. Os grupos que receberam a aplicação de silano e adesivo apresentaram superfícies mais lisas e maiores valores de resistência flexural em todas as resinas. No entanto não houve correlação dos valores de rugosidade com a resistência flexural, a qual demonstrou ser maior nas resinas SR Adoro e Signum Matrix, sendo que nesta última, apenas após a silanização é que se observou melhoria deste aspecto flexural.

O Afogar das veredas: uma análise comparativa espacial e temporal das veredas do Chapadão de Catalão (GO)

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

Avaliação de protocolos para colagem de bráquetes cerâmicos e remoção do adesivo remanescente com vistas á preservação da integridade de esmalte dentário

Pós-graduação em Ciências Odontológicas - FOAR

Osteogênese in vitro em amostras de titânio com diferente porosidades

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

Caracterização das propriedades de fluxo da lama vermelha para fins de dimensionamento de silos: aplicação para tremonhas cônicas

Um dos fatores que afetam o fluxo de produtos sólidos em silos é o teor de umidade, o qual se destaca por favorecer a formação de arcos coesivos. Nesse sentido, o objetivo deste trabalho focou a determinação das propriedades de fluxo do produto lama vermelha e sua classificação quanto à sua escoabilidade em função da variação do teor de umidade. Tendo os teores de umidade influência no escoamento em tremonhas cônicas revestidas de chapas metálicas lisas e UHMW (Polietileno de Ultra Alto Peso Molecular). Nesse intuito se utilizou o aparelho de cisalhamento direto “Jenike Shear Cell”, juntamente com três superfícies (acrílico, metálica e de UHMW) de parede. Para determinar o fator de fluxo da tremonha e o ângulo de inclinação com a horizontal, recorreu-se a metodologia gráfica proposta por Jenike para fluxo mássico. E na determinação do diâmetro mínimo da moega, aplicaram-se as equações propostas por Jenike e Enstad. Verificou-se que a lama vermelha é coesiva e possui fluxo do tipo difícil, com pouca variação para os teores de umidades estudados. A rugosidade das paredes testadas teve uma suave influência na inclinação da tremonha. Enquanto os resultados do diâmetro da moega sofreram variação em relação aos níveis de umidade analisados para a lama vermelha. A partir das observações feitas para o escoamento no silo modelo de acrílico se verificou que a medida que se eleva o teor de umidade mais difícil se torna o fluxo do produto lama vermelha.

Comportamento mecânico de material compósito de matriz poliéster reforçado por sistema híbrido fibras naturais e resíduos da indústria madeireira

Neste trabalho, materiais compósitos de matriz poliéster reforçados por fibras curtas de sisal, por resíduo de madeira e por sistema híbrido sisal/resíduo de madeira, dispostos aleatoriamente foram produzidos, utilizando-se o menor nível possível de processamento tecnológico nas etapas produtivas, com vistas a se produzir um compósito tecnicamente viável a pequenos produtores. A matriz de poliéster utilizada foi a tereftálica pré-acelerada com naftenato de cobalto e curada a temperatura ambiente com peróxido de metil-etil-cetona (MEK) em diferentes proporções em relação à resina, 0,33%, 1,66%, 3,33% e 5,00% em volume, de forma a se avaliar a influência deste nas propriedades mecânicas. As fibras de sisal foram cortadas manualmente nos comprimentos de 5, 10 e 15mm e utilizadas da maneira como adquiridas, sem tratamento superficial. O resíduo de madeira utilizado foi o pó de lixadeira da madeira maçaranduba. Os compósitos foram fabricados por moldagem manual, sem pressão e a temperatura ambiente. Foram fabricados corpos de prova de matriz pura, compósitos reforçados por sisal, variando-se o comprimento das fibras, compósitos reforçados por pó de maçaranduba e compósitos de reforço híbrido, sisal/pó de madeira, em diferentes proporções entre os constituintes. As propriedades mecânicas foram avaliadas por ensaios de tração e impacto charpy e as superfícies de fratura geradas foram avaliadas por microscopia eletrônica de varredura de modo a se correlacionar os aspectos de fratura com as propriedades mecânicas. Foi determinada a massa específica de cada série de corpos de prova fabricada, bem como a fração volumétrica dos reforços nos compósitos. Os resultados demonstraram que com o aumento do comprimento da fibra de sisal a resistência à tração e ao impacto dos compósitos foi incrementada, alcançando, o compósito com fibras de sisal de 15 mm, o melhor desempenho mecânico dentre as séries testadas. Por outro lado, a heterogeneidade granulométrica do pó de maçaranduba teve efeito negativo sobre as propriedades mecânicas dos compósitos. Os compósitos híbridos sisal/pó de madeira com maior teor de fibras, alcançaram 80% do desempenho obtido para os compósitos de fibras de sisal.

Determinação das velocidades intervalares usando a teoria paraxial do raio: aproximação de segunda ordem dos tempos de trânsito

Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método.

Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.

Empilhamento sísmico por superfície de reflexão comum: um novo algoritmo usando otimização global e local

O método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC) produz a simulação de seções com afastamento nulo (NA) a partir dos dados de cobertura múltipla. Para meios 2D, o operador de empilhamento SRC depende de três parâmetros que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀), o raio de curvatura da onda ponto de incidência normal (R_NIP) e o raio de curvatura da onda normal (R_N). O problema crucial para a implementação do método de empilhamento SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros ótimos associados a cada ponto de amostragem da seção AN a ser simulada. No presente trabalho foi desenvolvido uma nova sequência de processamento para a simulação de seções AN por meio do método de empilhamento SRC. Neste novo algoritmo, a determinação dos três parâmetros ótimos que definem o operador de empilhamento SRC é realizada em três etapas: na primeira etapa são estimados dois parâmetros (β°₀ e R°_NIP) por meio de uma busca global bidimensional nos dados de cobertura múltipla. Na segunda etapa é usado o valor de β°₀ estimado para determinar-se o terceiro parâmetro (R°_N) através de uma busca global unidimensional na seção AN resultante da primeira etapa. Em ambas etapas as buscas globais são realizadas aplicando o método de otimização Simulated Annealing (SA). Na terceira etapa são determinados os três parâmetros finais (β₀, R_NIP e R_N) através uma busca local tridimensional aplicando o método de otimização Variable Metric (VM) nos dados de cobertura múltipla. Nesta última etapa é usado o trio de parâmetros (β°₀, R°_NIP, R°_N) estimado nas duas etapas anteriores como aproximação inicial. Com o propósito de simular corretamente os eventos com mergulhos conflitantes, este novo algoritmo prevê a determinação de dois trios de parâmetros associados a pontos de amostragem da seção AN onde há intersecção de eventos. Em outras palavras, nos pontos da seção AN onde dois eventos sísmicos se cruzam são determinados dois trios de parâmetros SRC, os quais serão usados conjuntamente na simulação dos eventos com mergulhos conflitantes. Para avaliar a precisão e eficiência do novo algoritmo, este foi aplicado em dados sintéticos de dois modelos: um com interfaces contínuas e outro com uma interface descontinua. As seções AN simuladas têm elevada razão sinal-ruído e mostram uma clara definição dos eventos refletidos e difratados. A comparação das seções AN simuladas com as suas similares obtidas por modelamento direto mostra uma correta simulação de reflexões e difrações. Além disso, a comparação dos valores dos três parâmetros otimizados com os seus correspondentes valores exatos calculados por modelamento direto revela também um alto grau de precisão. Usando a aproximação hiperbólica dos tempos de trânsito, porém sob a condição de R_NIP = R_N, foi desenvolvido um novo algoritmo para a simulação de seções AN contendo predominantemente campos de ondas difratados. De forma similar ao algoritmo de empilhamento SRC, este algoritmo denominado empilhamento por Superfícies de Difração Comum (SDC) também usa os métodos de otimização SA e VM para determinar a dupla de parâmetros ótimos (β₀, R_NIP) que definem o melhor operador de empilhamento SDC. Na primeira etapa utiliza-se o método de otimização SA para determinar os parâmetros iniciais β°₀ e R°_NIP usando o operador de empilhamento com grande abertura. Na segunda etapa, usando os valores estimados de β°₀ e R°_NIP, são melhorados as estimativas do parâmetro R_NIP por meio da aplicação do algoritmo VM na seção AN resultante da primeira etapa. Na terceira etapa são determinados os melhores valores de β°₀ e R°_NIP por meio da aplicação do algoritmo VM nos dados de cobertura múltipla. Vale salientar que a aparente repetição de processos tem como efeito a atenuação progressiva dos eventos refletidos. A aplicação do algoritmo de empilhamento SDC em dados sintéticos contendo campos de ondas refletidos e difratados, produz como resultado principal uma seção AN simulada contendo eventos difratados claramente definidos. Como uma aplicação direta deste resultado na interpretação de dados sísmicos, a migração pós-empilhamento em profundidade da seção AN simulada produz uma seção com a localização correta dos pontos difratores associados às descontinuidades do modelo.

Inversão de momentos de fonte em métodos potenciais

A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino.

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.