Efeito da fragmentação florestal sobre a estrutura da comunidade de drosofilídeos (Diptera) no extremo sul da Amazônia

Perturbações ambientais freqüentemente provocam alterações na estrutura da paisagem e na diversidade faunistica. Dados referentes ao histórico de degradação, estrutura ambiental e coleta de espécies foram obtidas em 30 fragmentos florestais de vários tamanhos nos municípios de Alta Floresta e Apiacás, sul da Amazónia, com o intuito de examinar como a comunidade de Drosophilidae responde à alteração florestal. Os resultados mostraram grande heterogeneidade entre os fragmentos florestais, porém a diversidade de espécies foi semelhante entre os maiores e mais conservado e os menores com alto grau de degradação. Um total de 10.692 indivíduos de Drosophilidae, divididos em 62 táxons, foram coletados. Não foi observada relação entre tamanho dos fragmentos florestais e diversidade de espécies. Entretanto, a abertura do sub-bosque foi uma das variáveis ambientais que mais contribuiu para a diminuição da diversidade. A estimativa de riqueza de espécies para a região ficou em tomo de 100 espécies. A estrutura de abundância das comunidades foi mais similar entre os fragmentos florestais do que a composição de espécies. Cerca de 34% da riqueza total esteve restrita à borda florestal e 19% ao centro dos remanescentes, 47% das espécies foram comuns a ambos os habitats. Algumas espécies de drosofilideos responderam a um gradiente de distância matriz-centro florestal. O. maierkolliana, S. latifasciaefonnis e Z. indianus decresceram suas freqüências com a proximidade do centro, por outro lado, as espécies do grupo willistoni aumentaram em abundância. O principal efeito da fragmentação florestal sobre a comunidade de drosofilídeos foi a alteração drástica na composição de espécies, com a matriz do habitat sendo determinante na composição faunistica do fragmento.

Composição, riqueza e abundância de espécies de anfíbios na região do médio rio Xingú

Este estudo teve como objetivo inventariar a fauna de anfíbios nos diferentes hábitats da região do Médio Rio Xingu. Para tanto, foram selecionadas seis localidades, sendo as amostragens realizadas em três incursões, entre novembro de 2007 e março de 2008, totalizando 48 dias. Foram identificados seis tipos de hábitats na região amostrada: floresta de terra-firme, floresta sazonalmente alagável (várzea), margem do rio, “pedral”, lagoa e área alterada. Para a coleta de dados foram utilizadas duas metodologias: armadilhas de interceptação e queda e procura ativa. As armadilhas foram instaladas apenas nas áreas de floresta de terra-firme, enquanto a procura ativa foi empregada em todos os tipos de hábitats identificados. As amostragens resultaram no registro de 56 espécies de anfíbios e outras oito espécies foram identificadas e registradas em estudos anteriores na área. A diversidade e riqueza de espécies foram maiores na localidade Caracol, onde predominou a floresta de terra-firme, e menor na localidade Ilha Grande, dominada por floresta sazonalmente alagável. Foram registrados nove modos reprodutivos na área de estudo, sendo todos observados na floresta de terra-firme e apenas três nos “pedrais”, o que parece refletir a baixa heterogeneidade ambiental dessa área. Para comparar as localidades estudadas foram realizadas duas análises de similaridade, uma para cada metodologia de coleta utilizada. A análise de similaridade dos dados de procura ativa apontou maior semelhança na composição das espécies entre as áreas de floresta de terra-firme do que entre os outros hábitats. A análise de agrupamento realizada entre a composição de espécies obtida nesse estudo e outros levantamentos realizados na Amazônia agrupou esta área com outra também localizada no Médio Xingu.

Estoque de biomassa e carbono florestal em unidades de paisagem na Amazônia: uma análise a partir da abordagem metodológica ecologia da paisagem

O presente trabalho tem por objetivo central demonstrar a variabilidade existente na floresta no que tange aos estoques de biomassa e carbono florestal acima do solo, a partir da identificação e caracterização, com base em técnicas de sensoriamento remoto, de unidades de paisagem em uma área situada no município de Belterra, região oeste do Estado do Pará, a partir da matriz teórico-conceitual da abordagem Ecologia da Paisagem. Para o alcance de tal proposição, a metodologia empregada partiu da revisão da literatura sobre o tema, aquisição de dados cartográficos e orbitais, uso de técnicas de sensoriamento remoto, coleta de dados em campo, tratamento e análise estatística. O trabalho está dividido em quatro capítulos, seguidos pelas considerações gerais da obra. Partindo da matriz teórico-metodológica da Ecologia da Paisagem, analisa-se a dinâmica socioambiental do município de Belterra, que atualmente experimenta a expansão das atividades agrícolas, com destaque para a agricultura mecanizada da soja. A partir da análise multitemporal de imagens Landsat do município pôde-se avaliar a distribuição da cobertura florestal existente no mesmo, bem como o padrão espacial de distribuição das principais unidades de paisagem identificadas. Considerando esse recorte, realizou-se a coleta de dados em campo via inventário florestal em quatro tipologias florestais (floresta de alto platô, floresta de baixo platô, vegetação secundária e tensão ecológica) para obtenção de parâmetros morfométricos da vegetação e posterior quantificação dos estoques de biomassa e carbono contidos em cada unidade, bem como observar o comportamento estrutural da floresta nas mesmas. A adoção da paisagem como escala espacial de análise mostrou-se bastante satisfatória na quantificação dos estoques de biomassa e carbono florestal ao permitir considerar a influência da dinâmica socioeconômica na redução desses estoques. Além disso, possibilitou constatar que o reconhecimento da heterogeneidade da cobertura florestal é um elemento fundamental para a obtenção de estimativas de carbono de acordo com as características estruturais da vegetação, que varia de acordo com a topografia do terreno, com as espécies existentes e com as características geográficas, o que envolve a tipologia climática, as características geomorfológicas, pedológicas e geológicas da área.

Estudo da interação genótipo x ambiente sobre a produção de leite em rebanhos da raça Parodo-Suiço no Brasil, utilizando inferência Bayesiana

Verifica a presença de heterogeneidade de variâncias sobre a produção de leite na primeira lactação de fêmeas da raça Pardo-Suíço e, seu impacto sobre a avaliação genética dos reprodutores, utilizando a inferência Bayesiana por meio de amostrador de Gibbs, foram utilizados 2981 registros referentes às produções de leite e idade da vaca ao parto, em primeiras lactações de vacas da raça Pardo-Suíço, distribuídos em 62 rebanhos. Os registros foram provenientes do serviço de controle leiteiro da Associação Brasileira de Criadores de Gado Pardo Suíço, com os partos ocorridos entre os anos de 1980 a 2002. Foram estabelecidas duas classes de desvio-padrão fenotípico para produção de leite. Posteriormente, os dados foram analisados desconsiderando e considerando as classes de desvio-padrão. As médias observadas e desvio-padrão para produção de leite nas classes de alto e baixo desviopadrão e em análise geral foram iguais a 5802,02 ± 1929,96, 4844,37 ± 1592,99, 5373,47 ± 1849,13, respectivamente. As médias posteriores para os componentes de variâncias foram maiores na classe de alto desvio-padrão. A herdabilidade obtida na classe de alto desviopadrão foi próxima do valor observado na análise geral e inferior ao valor encontrado na classe de baixo desvio-padrão fenotípico. A correlação genética para a produção de leite entre as classes de desvio-padrão foi igual a 0,48. As correlações de Pearson e Spearman entre os valores genéticos para a produção de leite obtidos na análise geral, com os valores obtidos em cada classe de desvio-padrão foram todas maiores que 0,80, quando se considerou todos os reprodutores. Porém, refinando a amostra de reprodutores verifica-se que, as correlações diminuem em magnitude. Existindo uma maior variabilidade nos rebanhos presentes na classe de alto desvio-padrão e, o impacto dessa heterogeneidade de variância sobre a avaliação genética de reprodutores, é pequeno, pois a fonte principal dessa heterogeneidade é decorrente de fatores genéticos confirmando a presença de heterogeneidade de variâncias.

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.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Informação a priori na inversão de dados magnetotelúricos

Esta Tese apresenta dois estudos aplicados à inversão de dados magnetotelúricos. No primeiro deles, os parâmetros obtidos na inversão são as dimensões da malha de parametrização da subsuperfície, sendo conhecida, a priori, a resistividade de uma heterogeneidade e a do seu meio envolvente; no outro estudo, é abordado o uso de operadores de derivadas de ordem maior do que um com a finalidade de estabilizar o problema inverso. No primeiro estudo, os resultados podem ser considerados satisfatórios somente se a informação sobre as resistividades tem erro menor do que 20%. No segundo estudo, os resultados demonstram que o uso de operadores de ordem maior do que um podem ser mais eficazes do que o uso convencional do operador de primeira derivada, pois além de estabilizarem o problema inverso, esses operadores contribuem para melhorar a resolução das heterogeneidades de resistividade da subsuperfície. Ambos os estudos são inéditos, pois a prática de inversão de dados magnetotelúricos consiste de obter como resultado do problema inverso a resistividade dos prismas de uma malha de parametrização de dimensões fixas, usando como estabilizador o operador de primeira derivada. Os modelos usados nos estudos são bidimensionais e representam uma subsuperfície com uma e duas heterogeneidades de forma prismática envolvidas por ambiente homogêneo. O desempenho das técnicas foi testado com dados sintéticos com e sem ruído gaussiano, bem como dados reais do perfil COPROD2. Durante o trabalho, são, ainda, descritas as técnicas de inversão denominadas creeping e jumping e feita uma comparação e avaliação sobre elas. Mostra-se aqui que, ao contrário do que afirmam muitos pesquisadores, a inclusão de informação a priori sore os parâmetros pode ser feita na técnica do creeping com a mesma facilidade com que é feita na técnica do jumping.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.

Região do espaço que mais influencia em medidas eletromagnéticas no domínio da frequência: caso de uma linha de corrente sobre um semi-espaço condutor

Localizar em subsuperfície a região que mais influencia nas medidas obtidas na superfície da Terra é um problema de grande relevância em qualquer área da Geofísica. Neste trabalho, é feito um estudo sobre a localização dessa região, denominada aqui zona principal, para métodos eletromagnéticos no domínio da freqüência, utilizando-se como fonte uma linha de corrente na superfície de um semi-espaço condutor. No modelo estudado, tem-se, no interior desse semi-espaço, uma heterogeneidade na forma de camada infinita, ou de prisma com seção reta quadrada e comprimento infinito, na direção da linha de corrente. A diferença entre a medida obtida sobre o semi-espaço contendo a heterogeneidade e aquela obtida sobre o semi-espaço homogêneo, depende, entre outros parâmetros, da localização da heterogeneidade em relação ao sistema transmissor-receptor. Portanto, mantidos constantes os demais parâmetros, existirá uma posição da heterogeneidade em que sua influência é máxima nas medidas obtidas. Como esta posição é dependente do contraste de condutividade, das dimensões da heterogeneidade e da freqüência da corrente no transmissor, fica caracterizada uma região e não apenas uma única posição em que a heterogeneidade produzirá a máxima influência nas medidas. Esta região foi denominada zona principal. Identificada a zona principal, torna-se possível localizar com precisão os corpos que, em subsuperfície, provocam as anomalias observadas. Trata-se geralmente de corpos condutores de interesse para algum fim determinado. A localização desses corpos na prospecção, além de facilitar a exploração, reduz os custos de produção. Para localizar a zona principal, foi definida uma função Detetabilidade (∆), capaz de medir a influência da heterogeneidade nas medidas. A função ∆ foi calculada para amplitude e fase das componentes tangencial (H_x) e normal (H_z) à superfície terrestre do campo magnético medido no receptor. Estudando os extremos da função ∆ sob variações de condutividade, tamanho e profundidade da heterogeneidade, em modelos unidimensionais e bidimensionais, foram obtidas as dimensões da zona principal, tanto lateralmente como em profundidade. Os campos eletromagnéticos em modelos unidimensionais foram obtidos de uma forma híbrida, resolvendo numericamente as integrais obtidas da formulação analítica. Para modelos bidimensionais, a solução foi obtida através da técnica de elementos finitos. Os valores máximos da função ∆, calculada para amplitude de H_x, mostraram-se os mais indicados para localizar a zona principal. A localização feita através desta grandeza apresentou-se mais estável do que através das demais, sob variação das propriedades físicas e dimensões geométricas, tanto dos modelos unidimensionais como dos bidimensionais. No caso da heterogeneidade condutora ser uma camada horizontal infinita (caso 1D), a profundidade do plano central dessa camada vem dada pela relação p_o = 0,17 δ_o, onde p_o é essa profundidade e δ_o o "skin depth" da onda plana (em um meio homogêneo de condutividade igual à do meio encaixante (σ₁) e a freqüência dada pelo valor de w em que ocorre o máximo de ∆ calculada para a amplitude de H_x). No caso de uma heterogeneidade bidimensional (caso 2D), as coordenadas do eixo central da zona principal vem dadas por d_o = 0,77 r₀ (sendo d_o a distância horizontal do eixo à fonte transmissora) e p_o = 0,36 δ_o (sendo p_o a profundidade do eixo central da zona principal), onde r₀ é a distância transmissor-receptor e δ_o o "skin depth" da onda plana, nas mesmas condições já estipuladas no caso 1D. Conhecendo-se os valores de r₀ e δ_o para os quais ocorre o máximo de ∆, calculado para a amplitude de H_x, pode-se determinar (d_o, p_o). Para localizar a zona principal (ou, equivalentemente, uma zona condutora anômala em subsuperfície), sugere-se um método que consiste em associar cada valor da função ∆ da amplitude de H_x a um ponto (d, p), gerado através das relações d = 0,77 r e p = 0,36 δ, para cada w, em todo o espectro de freqüências das medidas, em um dado conjunto de configurações transmissor-receptor. São, então, traçadas curvas de contorno com os isovalores de ∆ que vão convergir, na medida em que o valor de ∆ se aproxima do máximo, sobre a localização e as dimensões geométricas aproximadas da heterogeneidade (zona principal).

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 um esquema linear para estimativa de anisotropia local a partir de dados de onda P em experimentos de VSP multiazimutal

Neste trabalho é apresentada uma análise do esquema de inversão linear para a estimativa de anisotropia na vizinhança de um receptor situado em um poço a partir de da componente vertical do vetor de vagarosidade e do vetor de polarização de ondas P medidops em experimentos de VSP walkaway multiazimutal. Independente do meio acima do geofone (homogêneo ou heterogêneo) e da forma do poço (pode ser direcional ou curvado, vertical e inclinado), a inversão é feita a partir de uma aproximação de primeira ordem em torno de um meio isotrópico de referência. O esquma da inversão é analisado considerando fatores como: o nível de ruído nos dados, o tipo de onda P, o grau de anisotropia do meio, a escolha dos parâmetros no meio isotrópico de referência e grau de heterogeneidade do meio. Os resultados são apresentados.

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.

Eletroresistividade na prospecção hidrogeológica no município de Ourilândia do Norte-PA

A prospecção geofísica para a exploração de água subterrânea em regiões cristalinas torna-se uma tarefa extremamente difícil pela presença de estruturas como falhas, fraturas e ocorrência de variação lateral de resistividade no manto intempérico. Para analisar a influência da variação bidimensional da propriedade física no comportamento das curvas de sondagem elétrica vertical Schlumberger, simulou-se modelos contendo tais estruturas encaixadas no manto de alteração com ângulo de mergulho de 30° e 90° e com a posição do centro de sondagem variando em relação a heterogeneidade. A modelagem numérica foi feita através do método dos elementos finitos com o programa EGSLIB/SEV2D. O resultado da simulação mostra que o emprego da técnica da sondagem elétrica vertical Schlumberger sobre meios bidimensionais não funciona eficientemente, como quando aplicada em regiões sedimentares e de aluvião, devido o fenômeno da ambigüidade que dificulta a interpretação das curvas. Dentre os métodos geofísicos para a investigação hidrogeológica, a aplicação dos métodos eletrorresistivos mostra-se mais simples, eficaz e econômica. Assim, o levantamento geofísico para orientar os trabalhos de perfurações de poços na cidade de Ourilândia do Norte contou de trinta e duas sondagens elétricas verticais tipo Schlumberger e um caminhamento elétrico dipolo-dipolo realizados em diversas ruas da cidade. As curvas de resistividade aparente, obtidas com as sondagens foram tratadas e processadas considerando um modelo de camadas planas, horizontais, homogêneas e isotropicas, inicialmente com o algoritmo EGSLIB/SEV1D e posteriormente com uso do programa EGSLIB/SEV1DINV, permitindo estabelecer modelos geoelétricos aproximados para a configuração da subsuperfície da cidade. Depois da interpretação quantitativa, foi possível separar duas zonas geoelétricas distintas. Esta diferenciação foi realizada com base nos valores de resistividade aparente e as informações geológicas de superfície. A primeira, com pouca representatividade na área da cidade, apresenta-se constituída por solo/aterro sobre saprolito de natureza argilo-arenosa recobrindo os granitóides tipo Rio Maria. A segunda composta pelo solo/aterro seguido de um horizonte preenchido por sedimentos arenosos superposto a camada argilo-arenosa e finalmente a rocha sã. Esta seqüência geoelétrica predomina na cidade. Os estratos arenoso e argilo-arenoso apresentam condições potenciais para armazenamento de água.

Comparação de peso e rendimento de carcaça entre vacas bovinas e bubalinas

Registros de peso vivo, peso de carcaça, peso de cortes primários e rendimento de carcaça, provenientes do Matadouro Frigorífico da Cooperativa da Indústria Pecuária do Pará Ltda. (SOCIPE), localizada no distrito de Icoaraci, município de Belém, Estado do Pará, foram utilizados para discriminar possíveis diferenças entre fêmeas bubalinas e bovinas, abatidas durante o ano de 2007. Os registros de peso vivo foram utilizados para separar a amostra em fêmeas leves e pesadas, separadamente em cada espécie, utilizando como limiar o valor correspondente à mediana do peso vivo em cada espécie. Assim, foram compostos quatro grupos de fêmeas, bubalinas e bovinas, leves e pesadas. Para comparação entre fêmeas, vetores de médias dos pesos vivo, pesos carcaça quente, carcaça quente esquerda e direita, peso da carcaça fria, peso dianteiro, peso traseiro, peso ponta-de-agulha e peso couro, entre vacas bubalinas e bovinas, leves e pesadas, foram utilizados procedimentos de inferência multivariada. Os resultados revelaram a heterogeneidade para o peso de abate dos animais, dentro e entre espécies, não havendo um padrão bem definido de peso ideal ao abate. Vacas bubalinas obtiveram maiores pesos para componente de carcaça, quando comparadas às fêmeas do grupo bovino. Excetuando o peso da carcaça fria e a ponta de agulha, bubalinas do grupo pesado, obtiveram maiores médias quando comparadas ao grupo leve. Dentro do grupo de fêmeas bovinas não houve diferença estatística para o peso da caraça fria (P>0,05). Verifica-se que, em média, fêmeas bovinas apresentaram maior rendimento, exceto para o dianteiro. Os grupos intra-raciais apresentaram-se significativamente diferentes entre si (P<0,01). Por outro lado, não foram verificadas diferenças intra-raciais para os rendimentos de carcaças quentes, carcaça fria, dianteiro e traseiro.