Wave equation depth migration using complex Padé approximation

Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domínio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido às ondas evanescentes e produz imagens com menos ruídos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos íngremes. Os resultados foram comparados àqueles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD).

Linearização dos coeficientes de reflexão de ondas qP em meios anisotrópicos

Os coeficientes de reflexão de onda qP em uma interface plana que separa dois meios anisotrópicos dependem dos parâmetros elásticos dos meios envolvidos de maneira não linear. Aproximações lineares para a refletividade de ondas qP pressupondo fraco contraste de impedância e fraca anisotropia levam a uma forma simples para se fazer à análise de AVO/AVD. Neste trabalho a solução das equações de Zoeppritz é rescrita explicitamente em função de matrizes de impedância e polarização. Além disso, utilizando-se esta abordagem, é apresentada uma metodologia geral mais simples para se obter as formas linearizadas. Estas formas linearizadas para os coeficientes de reflexão da onda qP e das ondas convertidas para uma onda qP incidente apresentam resultados muito próximos da formulação exata para ângulos de incidência menores que 30° considerando-se um contraste de impedância moderado e anisotropia em limites geologicamente aceitáveis.

Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

Conforto ambiental de bezerros bubalinos (Bubalus bubalis Linnaeus, 1758) em sistemas silvipastoris na Amazônia Oriental

Foi analisado o efeito de dois sistemas silvipastoris nos índices de conforto ambiental e alterações nos parâmetros fisiológicos de bezerros bubalinos criados na Embrapa Amazônia Oriental, Belém–PA (clima Afi), no Período 1 (abril a setembro/2007) e Período 2 (outubro/2007 a março/2008). Foram inseridos onze bezerros no Sistema Silvipastoril 1 (SSP1), que apresentava sombreamento útil nas pastagens de 18 a 21%, e oito no Sistema Silvipastoril 2 (SSP2), sem sombreamento, com lago para banho. Foram mensuradas as variáveis fisiológicas: temperatura retal (TR), frequência respiratória (FR), temperatura da pele (TP), e calculados o Índice de Temperatura e Umidade (ITU) e Índice de Conforto de Benezra (ICB), nos dois períodos experimentais, comparados pelo Teste Tukey (P < 0,05). O ITU apresentou diferença estatística entre horários (P < 0,05) e período do ano (P < 0,05), e oscilou de 73,5 ± 1,3 até 82,2 ± 0,8. A TR apresentou diferença estatística entre horários e períodos do ano (P < 0,05), com amplitude de 38,3 ± 0,26 a 39,3 ± 0,38 °C. A FR apresentou diferença significativa entre horários (P < 0,05), com amplitude de 32,2 ± 9,2 a 56,5 ± 19,0 mov min^-1, consideradas acima dos níveis normais, enquanto a TP foi diferente estatisticamente entre períodos e horários (P < 0,05) e variou de 23,6 ± 8,3 a 31,7 ± 5,4 °C. Nos Períodos 1 e 2 e nos dois SSP's, os ICBs estiveram acima do valor ideal, variando de 2,46 ± 0,33 a 3,31 ± 0,62 (SSP1) e 2,42 ± 0,30 a 3,45 ± 0,66 (SSP2).

Perfil hematológico de búfalas da raça Murrah, criadas ao sol e à sombra, em clima tropical da Amazônia Oriental

Foram avaliadas as respostas hematológicas de 20 búfalas, criadas ao sol (grupo - SS) e à sombra (grupo - CS), em Belém, Pará. Os animais do grupo CS (n = 10) estavam em sistema silvipastoril, com Racosperma mangium e os do SS (n = 10), em piquetes sem sombra, em pastagem de Brachiaria humidicola, água para beber e sal mineral. Foram mensuradas temperatura do ar (TA), umidade relativa do ar (UR) e temperatura de globo negro, em cada tratamento. A coleta de sangue para eritrograma e número total de leucócitos foi realizada a cada 14 dias, às 13 h, durante o ano de 2009. Através da análise de variância constatou-se que em todos os períodos do ano, a TA e índice de temperatura de globo e umidade (ITGU) foram diferentes (P < 0,05), com níveis mais elevados no grupo SS. No período mais chuvoso, o grupo CS apresentou valores elevados de leucócitos, enquanto nos períodos de transição e menos chuvoso, foram maiores no grupo SS. No período menos chuvoso do ano, o grupo SS apresentou maiores valores de hemácias. O teor de hemoglobina teve maiores níveis (P < 0,05) nos períodos de transição e menos chuvoso. Somente a hemoglobina teve correlação significativa e negativa (P < 0,05) com a UR. Conclui-se que as novilhas búfalas Murrah estão sujeitas a um ambiente hostil e o período menos chuvoso é o mais propício a causar estresse térmico.

Abundância, composição e diversidade de Arctiinae (Lepidoptera, Arctiidae) em um fragmento de floresta na Amazônia Oriental em Altamira, PA, Brasil

Este estudo avaliou a fauna de Arctiinae em um fragmento de floresta primária em Altamira, Pará, na Amazônia Oriental brasileira. As mariposas foram amostradas durante dois anos (de agosto de 2007 a julho de 2009), com auxílio de armadilha luminosa. Foram medidos os seguintes parâmetros: riqueza, abundância, constância, índices de diversidade e uniformidade de Shannon (H' e E') e de Brillouin (H e E) e o índice de dominância de Berger-Parker (BP). As estimativas de riqueza, foram efetuadas através dos procedimentos não paramétricos, "Bootstrap", "Chao 1", "Chao 2", "Jackknife 1", "Jackknife2" e "Michaelis-Mentem". Foram capturados 466 exemplares pertencentes a 78 espécies de Arctiinae, das quais 12 são novos registros para o Estado. Os valores dos parâmetros analisados para todo o período foram: H'= 3,08, E'= 0,708, H= 2,86, E= 0,705 e BP= 0,294. As comunidades dos meses menos chuvosos foram mais diversas. Os estimadores previram o encontro de 17 a 253 espécies a mais.

Metodologia geofísica aplicada ao estudo arqueológico dos sítios Bittencourt e Jambuaçu, Estado do Pará

Medidas magnéticas, de radiação gama e eletromagnéticas com o Radar de Penetração no Solo (GPR) foram realizadas em dois sítios arqueológicos localizados nas proximidades de uma área destinada a construção de um mineroduto. O objetivo das medidas foi identificar locais para escavação arqueológica visando a descoberta de remanescentes de habitantes que povoaram a região em tempos pretéritos, auxiliando nos estudos para recompor a sua história. As medidas magnéticas evidenciaram anomalias que levaram a descoberta de uma grande quantidade de fragmentos cerâmicos e de outros artefatos. As medidas de radiação gama permitiram a discriminação entre solo arenoso e argiloso, o que representa um resultado importante, pois o solo arenoso geralmente é correlacionado a camada de ocupação arqueológica onde são encontrados os remanescentes. O Radar de Penetração no Solo confirmou as anomalias magnéticas e mostrou outras feições anômalas não relacionadas aos remanescentes arqueológicos.

Caracterização geoelétrica de área de curtume localizada no distrito industrial de Icoaraci, Bélem-Pará

Foi desenvolvido um estudo da subsuperfície rasa para investigar a possibilidade de contaminação do subsolo por líquidos usados no tratamento de couros. O estudo foi realizado na área da indústria Couro do Norte, Bairro de Maracacuera, no Distrito Industrial de Icoaraci em Belém-PA. As metodologias geofísicas utilizadas foram: eletrorresistividade (na forma de imageamento e através de sondagens verticais), potencial espontâneo, eletromagnético (slingram) e radar de penetração no solo (GPR). A interpretação conjunta dos dados geofísicos indicou a direção do fluxo subterrâneo e uma pluma de contaminação rasa na área do estudo causada provavelmente pelo despejo dos efluentes do curtume Couro do Norte.

Mecanismos de controle da variação sazonal da transpiração de uma floresta tropical no Nordeste da Amazônia

No presente trabalho foram estudadas a variação sazonal da transpiração, de uma floresta tropical, e sua dependência com fatores bióticos e abióticos. Utilizaram-se dados do projeto CARBOPARÁ, parte integrante do Experimento de Grande Escala da Biosfera-Atmosfera na Amazônia (LBA), coletados na reserva florestal de Caxiuanã, região nordeste da Amazônia. A evapotranspiração total num intervalo de 39 dias para o período chuvoso foi 108,2 mm, com valor médio de 2,9 mm dia^-1, enquanto, durante o período menos chuvoso, a evapotranspiração total num intervalo de 29 dias foi 128,8 mm, com média de 4,3 mm dia^-1 para o período. Os valores máximos da condutividade de superfície (Cs), nos dois períodos, ocorreram às 08:00 hl, sendo estes valores de 0,060 m s^-1 e 0,045 m s^-1 para o período chuvoso e menos chuvoso, respectivamente. A condutância aerodinâmica média (Ca) foi 0,164 m s^-1 e 0,210 m s^-1, para os períodos chuvoso e menos chuvoso, respectivamente. Os valores máximos da Ca observados para os períodos chuvoso e menos chuvoso foram, respectivamente, 0,220 e 0,375 m s^-1. Verificou-se que Cs guarda uma relação exponencial inversa com o déficit de vapor de água atmosférico, para diferentes intervalos de irradiância solar global. A análise horária do fator de desacoplamento sugere que a evapotranspiração, durante a manhã, tem um maior controle realizado pela disponibilidade de energia, quando comparado ao período menos chuvoso. Durante a tarde verifica-se que o dossel da floresta progressivamente tende a estar mais acoplado à atmosfera, para ambos os períodos estudados, demonstrando maior controle superficial na transpiração.

Interações entre nuvens, chuvas e a biosfera na Amazônia

Apresenta-se uma revisão dos resultados já obtidos no Projeto LBA, na componente de Física do Clima, referente ao efeito da interação entre a biosfera e a formação de nuvens e de chuva em diferentes estações do ano e regiões da Bacia Amazônica. O efeito do desmatamento e das queimadas na atmosfera é abordado assim como mecanismos atmosféricos que foram caracterizados a partir das diversas campanhas do projeto LBA. Uma visão da complexidade das interações e do progresso que tem sido feito com base nos dados colhidos é apresentada.

CRS seismic data imaging system: a case study for basib reevaluation

Este trabalho apresenta resultados práticos de uma atenção sistemática dada ao processamento e à interpretação sísmica de algumas linhas terrestres do conjunto de dados do gráben do Tacutu (Brasil), sobre os quais foram aplicadas etapas fundamentais do sistema WIT de imageamento do empilhamento CRS (Superfície de Reflexão Comum) vinculado a dados. Como resultado, esperamos estabelecer um fluxograma para a reavaliação sísmica de bacias sedimentares. Fundamentado nos atributos de frente de onda resultantes do empilhamento CRS, um macro-modelo suave de velocidades foi obtido através de inversão tomográfica. Usando este macro-modelo, foi realizado uma migração à profundidade pré- e pós-empilhamento. Além disso, outras técnicas baseadas no empilhamento CRS foram realizadas em paralelo como correção estática residual e migração de abertura-limitada baseada na zona de Fresnel projetada. Uma interpretação geológica sobre as seções empilhadas e migradas foi esboçada. A partir dos detalhes visuais dos painéis é possível interpretar desconformidades, afinamentos, um anticlinal principal falhado com conjuntos de horstes e grábens. Também, uma parte da linha selecionada precisa de processamento mais detalhado para evidenciar melhor qualquer estrutura presente na subsuperfície.

Filtros ótimos: efeitos de janelas e de tipo de fase

O objetivo central deste trabalho é o estudo do desempenho do operador deconvolucional WHL na compressão do pulso-fonte sísmico, sob condições especiais de fase não-mínima e da densidade de eventos no traço, como casos advogados para dados reais e processamento em rotina. O método de ataque ao problema construído é centrado no conteúdo da informação da função autocorrelação submetida a diferentes condições: (a) de truncamento e tipo de janelas; (b) das características da fase do operador (se mínima ou não-mínima); (c) da medida de qualidade; (d) do nível de embranquecimento; (e) do ruído presente e da equalização; (f) do balanceamento do traço; (g) dos princípios físicos da propagação expressos e limitados pelo modelo convolutional. Os resultados obtidos são apenas na forma numérica, organizados na forma de álbuns com dificuldades crescentes, e demonstram como o uso de janelas na autocorrelação serve para diagnosticar e melhorar a performance dos operadores. Concluímos que muitas perguntas ainda surgem quando técnicas de deconvolução são aplicadas a seções sísmicas de bacias sedimentares, e que o modelo de Goupillaud é conveniente para simulações e análises devido a sua descrição matemática simples e completa.

Deconvolução de processo sísmico não-estacionário

O presente trabalho trata da aplicação do filtro Kalman-Bucy (FKB), organizado como uma deconvolução (FKBD), para extração da função refletividade a partir de dados sísmicos. Isto significa que o processo é descrito como estocástico não-estacionário, e corresponde a uma generalização da teoria de Wiener-Kolmogorov. A descrição matemática do FKB conserva a relação com a do filtro Wiener-Hopf (FWH) que trata da contra-parte com um processo estocástico estacionário. A estratégia de ataque ao problema é estruturada em partes: (a) Critério de otimização; (b) Conhecimento a priori; (c) Algoritmo; e (d) Qualidade. O conhecimento a priori inclui o modelo convolucional, e estabelece estatísticas para as suas componentes do modelo (pulso-fonte efetivo, função refletividade, ruídos geológico e local). Para demostrar a versatilidade, a aplicabilidade e limitações do método, elaboramos experimentos sistemáticos de deconvolução sob várias situações de nível de ruídos aditivos e de pulso-fonte efetivo. Demonstramos, em primeiro lugar, a necessidade de filtros equalizadores e, em segundo lugar, que o fator de coerência espectral é uma boa medida numérica da qualidade do processo. Justificamos também o presente estudo para a aplicação em dados reais, como exemplificado.

Autonomia funcional de idosas praticantes de Pilates

O objetivo foi avaliar o efeito do método Pilates na autonomia funcional de idosas. Foram selecionadas 52 voluntárias, submetidas à avaliação geral e ao protocolo de avaliação da autonomia funcional do Grupo Latino-Americano de Desenvolvimento para Maturidade (GDLAM), que consiste em cinco testes (caminhada de 10 m, levantar-se da posição sentada, levantar-se da cadeira e locomover-se pela casa, levantar-se da posição de decúbito ventral e vestir e tirar uma camiseta); em seguida, foram divididas em dois grupos: grupo Pilates (GP, n=27; idade 66,9±5,3 anos) e grupo controle (GC, n=25; idade 65,2±3,9 anos). O GP foi submetido a uma série de dez exercícios de Pilates, por oito semanas, duas vezes por semana. Os dois grupos foram reavaliados após esse período. O nível de significância considerado foi de p?0,05. O GP obteve resultados significativamente melhores em todos os testes e no índice geral do GDLAM (p=0,035) após a intervenção. O GC obteve escore significativamente melhor (p=0,042) apenas no teste de caminhada de 10 m, tendo mantido sua classificação inicial de funcionalidade regular. Comparando-se os escores dos grupos após a intervenção, encontraram-se diferenças significativas em favor do GP, inclusive no índice GDLAM (p<0,05), tendo este passado da classificação de regular a bom. A prática do método Pilates gerou pois melhora significativa do desempenho funcional das idosas estudadas.

Chemical composition of phytoplankton from the estuaries of Eastern Amazonia

O fitoplancton é de grande importância para os estudos das interações solo-água, como indicadores de modificações químicas e biológicas nos ecossistemas naturais. O objetivo desse estudo foi determinar a composição química total do fitoplâncton nos estuários do rio Pará e Mocajuba (Pará, Brasil). As análises químicas foram realizadas na água superficial, sedimento de fundo (amostra total e fração biodisponível) e no fitoplancton, por espectrometria óptica e de massa com plasma indutivamente acoplado. A composição química elementar do fitoplancton é composta por elevadas concentrações de Ca, P, Mn, Fe, Zn, Al, Ba e Pb. O fitoplancton do estuário do rio Mocajuba é rico em Fe (2.967 a 84.750 µg g^-1) e do rio Pará rico em Al (1.216 a 15.389 µg g^-1), provavelmente com contribuição antropogênica. O material fitoplanctonico apresentou elevado fator de bioconcentração proveniente tanto da água quanto da fração biodisponível, e reflete a eficiência desses organismos em concentrar metais.