Estudo metaloproteômico em amostras de tecido muscular de Tilápia do Nilo alimentadas com fontes de zinco orgânico e inorgânico

Pós-graduação em Zootecnia - FMVZ

Análise de algumas proteinas miofibrilares envolvidas na maciez da carne em bovinos de corte

Pós-graduação em Zootecnia - FMVZ

Análise comparativa dos efeitos alelopáticos das substâncias químicas titonina e titonina acetilada

Este trabalho teve por objetivo estabelecer as variações na atividade alelopática da substância química titonina, em função da acetilação de sua molécula. Bioensaios de germinação (25 ºC de temperatura constante e fotoperíodo de 12 horas) e de desenvolvimento da radícula e do hipocótilo (25 ºC de temperatura constante e fotoperíodo de 24 horas) foram desenvolvidos. Como planta receptora, utilizou-se a planta daninha Mimosa pudica (malícia). Análise de espectros RMN ¹H e ¹³C e técnicas de RMN bidimensionais foram realizadas na molécula acetilada. O processo de acetilação produziu a molécula 3'-acetil-7,4'-dimetoxiflavona, que diferiu da molécula original, identificada como 3'-hidroxi-7,4'-dimetoxiflavona. A estrutura da titonina-Ac foi confirmada pelos espectros de RMN ¹H, ¹³C, DEPT, COSY e HETCOR. A titonina foi acetilada com anidrido acético em piridina. A análise comparativa da atividade alelopática das duas substâncias revelou que titonina-Ac apresentou maior potencial para inibir tanto a germinação das sementes como o desenvolvimento da radícula e do hipocótilo da planta daninha malícia. A intensidade dos efeitos alelopáticos das duas substâncias esteve positivamente associada à concentração. O conjunto das informações obtidas permite sugerir a possibilidade de se aumentar a atividade alelopática de uma substância química sem comprometer suas peculiaridades biológicas desejáveis à natureza e aos interesses da sociedade.

Entropia conjunta de espaço e reqüência espacial estimada através da discriminação de estímulos espaciais com luminância e cromaticidade moduladas por funções de Gábor: implicações para o processamento paralelo de informação no sistema visual humano

O objetivo deste estudo foi estimar a entropia conjunta do sistema visual humano no domínio do espaço e no domínio das freqüências espaciais através de funções psicométricas. Estas foram obtidas com testes de discriminação de estímulos com luminância ou cromaticidade moduladas por funções de Gábor. A essência do método consistiu em avaliar a entropia no domínio do espaço, testando-se a capacidade do sujeito em discriminar estímulos que diferiam apenas em extensão espacial, e avaliar a entropia no domínio das freqüências espaciais, testando-se a capacidade do sujeito em discriminar estímulos que diferiam apenas em freqüência espacial. A entropia conjunta foi calculada, então, a partir desses dois valores individuais de entropia. Três condições visuais foram estudadas: acromática, cromática sem correção fina para eqüiluminância e cromática com correção para eqüiluminância através de fotometria com flicker heterocromático. Quatro sujeitos foram testados nas três condições, dois sujeitos adicionais foram testados na condição cromática sem eqüiluminância fina e um sétimo sujeito também fez o teste acromático. Todos os sujeitos foram examinados por oftalmologista e considerados normais do ponto de vista oftálmico, não apresentando relato, sintomas ou sinais de disfunções visuais ou de moléstias potencialmente capazes de afetar o sistema visual. Eles tinham acuidade visual normal ou corrigida de no mínimo 20/30. O trabalho foi aprovado pela Comissão de Ética em Pesquisa do Núcleo de Medicina Tropical da UFPA e obedeceu às recomendações da Declaração de Helsinki. As funções de Gábor usadas para modulação de luminância ou cromaticidade compreenderam redes senoidais unidimensionais horizontais, moduladas na direção vertical, dentro de envelopes gaussianos bidimensionais cuja extensão espacial era medida pelo desvio padrão da gaussiana. Os estímulos foram gerados usando-se uma rotina escrita em Pascal num ambiente Delphi 7 Enterprise. Foi utilizado um microcomputador Dell Precision 390 Workstation e um gerador de estímulos CRS VSG ViSaGe para exibir os estímulos num CRT de 20”, 800 x 600 pixels, 120 Hz, padrão RGB, Mitsubishi Diamond Pro 2070SB. Nos experimentos acromáticos, os estímulos foram gerados pela modulação de luminância de uma cor branca correspondente à cromaticidade CIE1931 (x = 0,270; y = 0,280) ou CIE1976 (u’ = 0,186; v’= 0,433) e tinha luminância média de 44,5 cd/m2. Nos experimentos cromáticos, a luminância média foi mantida em 15 cd/m2 e foram usadas duas series de estímulos verde-vermelhos. Os estímulos de uma série foram formados por duas cromaticidades definidas no eixo M-L do Espaço de Cores DKL (CIE1976: verde, u’=0,131, v’=0,380; vermelho, u’=0,216, v’=0,371). Os estímulos da outra série foram formados por duas cromaticidades definidas ao longo de um eixo horizontal verde-vermelho definido no Espaço de Cores CIE1976 (verde, u’=0,150, v’=0,480; vermelho, u’=0,255, v’=0,480). Os estímulos de referência eram compostos por redes de três freqüências espaciais diferentes (0,4, 2 e 10 ciclos por grau) e envelope gaussiano com desvio padrão de 1 grau. Os estímulos de testes eram compostos por uma entre 19 freqüências espaciais diferentes em torno da freqüência espacial de referência e um entre 21 envelopes gaussianos diferentes com desvio padrão em torno de 1 grau. Na condição acromática, foram estudados quatro níveis de contraste de Michelson: 2%, 5%, 10% e 100%. Nas duas condições cromáticas foi usado o nível mais alto de contraste agregado de cones permitidos pelo gamut do monitor, 17%. O experimento consistiu numa escolha forçada de dois intervalos, cujo procedimento de testagem compreendeu a seguinte seqüência: i) apresentação de um estímulo de referência por 1 s; ii) substituição do estímulo de referência por um fundo eqüiluminante de mesma cromaticidade por 1 s; iii) apresentação do estímulo de teste também por 1 s, diferindo em relação ao estímulo de referência seja em freqüência espacial, seja em extensão espacial, com um estímulo sonoro sinalizando ao sujeito que era necessário responder se o estímulo de teste era igual ou diferente do estímulo de referência; iv) substituição do estímulo de teste pelo fundo. A extensão espacial ou a freqüência espacial do estímulo de teste foi mudada aleatoriamente de tentativa para tentativa usando o método dos estímulos constantes. Numa série de 300 tentativas, a freqüencia espacial foi variada, noutra série também de 300 tentativas, a extensão espacial foi variada, sendo que cada estímulo de teste em cada série foi apresentado pelo menos 10 vezes. A resposta do indivíduo em cada tentativa era guardada como correta ou errada para posterior construção das curvas psicométricas. Os pontos experimentais das funções psicométricas para espaço e freqüência espacial em cada nível de contraste, correspondentes aos percentuais de acertos, foram ajustados com funções gaussianas usando-se o método dos mínimos quadrados. Para cada nível de contraste, as entropias para espaço e freqüência espacial foram estimadas pelos desvios padrões dessas funções gaussianas e a entropia conjunta foi obtida multiplicando-se a raiz quadrada da entropia para espaço pela entropia para freqüência espacial. Os valores de entropia conjunta foram comparados com o mínimo teórico para sistemas lineares, 1/4π ou 0,0796. Para freqüências espaciais baixas e intermediárias, a entropia conjunta atingiu níveis abaixo do mínimo teórico em contrastes altos, sugerindo interações não lineares entre dois ou mais mecanismos visuais. Este fenômeno occorreu em todas as condições (acromática, cromática e cromática eqüiluminante) e foi mais acentuado para a frequência espacial de 0,4 ciclos / grau. Uma possível explicação para este fenômeno é a interação não linear entre as vias visuais retino-genículo-estriadas, tais como as vias K, M e P, na área visual primária ou em níveis mais altos de processamento neural.

Ecobiometria e Triplex Doppler para a avaliação da vitalidade fetal em Cebus apella

Os objetos do trabalho foram: 1) realizar a ecobiometria do concepto para acompanhar o crescimento fetal e determinar a idade gestacional em Cebus apella; 2) descrever o momento em que os órgãos do feto são observados; 3) realizar a sexagem fetal; 4) avaliar o fluxo sanguíneo das artérias uterinas (AU) e umbilical (Aum), mensurando os índices de resistividade (IR) e pulsatilidade (IP); 5) observar a presença ou ausência da incisura protodiastólica (IPD) e do componente diastólico nas ondas de fluxo sanguíneo das AU e Aum, respectivamente, durante a gestação em C. apella. Seis fêmeas adultas da espécie C. apella foram selecionadas e, posteriormente, condicionadas para os procedimentos de colpocitologia ou contenção química. Para o monitoramento do ciclo reprodutivo das fêmeas e crescimento folicular, foram realizados exames colpocitológicos e ultrassonografia (US), respectivamente, para a escolha do melhor dia da monta natural ou inseminação artificial. Vinte dias após a cópula ou inseminação, era realizado o diagnóstico da gestação por meio da US. O dia de cada exame em relação ao parto foi contado em retrospecto (nascimento = dia 0). Os exames ultrassonográficos foram feitos nos dias -133, -113, -83, -53, -21 e -1 antes do parto. A US bidimensional em modo B foi utilizada para mensurar o saco gestacional do dia - 133 até -113; os diâmetros biparietal, fronto-occiptal, circunferência da cabeça e circunferência abdominal do dia -113 a -1; e comprimento do fêmur do dia -53 até -23, mostrando o aumento desses parâmetros com o avanço da gestação. O coração e o estômago começaram a ser visualizados no dia -113 e, os demais órgãos e a sexagem no dia -83. O Triplex Doppler foi empregado para avaliar o fluxo sanguíneo durante o período gestacional, mostrando uma diminuição nos IR e IP das AU e Aum, do dia -133 a -1, bem como observar o desaparecimento da IPD (dia -1) e aparecimento do componente diastólico (dia -53) na onda de fluxo sanguíneo das AU e Aum, respectivamente. Entre os dias -113 a -1, a média da freqüência cardíaca fetal, obtida pelo Triplex Doppler, foi de 189 ± 2,43 bpm. O presente trabalho permitiu determinar a idade gestacional, avaliar o crescimento anatômico do feto, descrever o momento em que os órgãos são visualizados e realizar a sexagem em C. apella. Constatou-se, também, que o fluxo sanguíneo das AU e Aum são importantes parâmetros para avaliar a vitalidade fetal em C. apella.

Diversidade genética de isolados ambientais de Vibrio cholerae da Amazônia brasileira

Vibrio cholerae, agente etiológico da cólera, é uma bactéria nativa de ambientes aquáticos de regiões temperadas e tropicais em todo o mundo. A cólera é endemica e epidemica em países da África, Ásia e Americas Central e do Sul. Neste trabalho o objetivo foi estudar a diversidade genética de isolados desta espécie, de ambientes aquáticos da Amazônia brasileira. Um total de 148 isolados de V.cholerae não-O1 e não-O139 (NAGs) e O1 ambientais da Amazônia, obtidos entre 1977 e 2007, foram caracterizados e comparados a linhagens clínicas de V.cholerae O1 da sexta e sétima pandemias. Utilizou-se os perfis de macrorestrição definidos em eletroforese em gel de agarose em campo pulsado (PFGE), para determinar a relação clonal entre V.cholerae non-O1 e O1 ambientais e clínicos. A presença de genes de virulência (hlyA/hem, hlyB, hlyC, rtxA, rtxC, tcp, ctx, zot, ace, stn/sto) e integrons de classe 1, 2 e 3 (intI 1, 2 e 3), foi analisada utilizando-se a reação em cadeia da polimerase. A análise dos perfis de macrorestrição revelou que os NAGs apresentaram uma grande diversidade genética comparada aos V.cholerae O1. Isolados de NAGs e O1 segregaram em distintos grupos e a maioria dos O1 ambientais apresentou relação clonal com isolados clínicos da sétima pandemia de cólera. A distribuição dos genes de virulência entre os NAGs é diferente a dos O1, os quais, em geral, foram positivos para todos os genes de virulência estudados exceto stn/sto e integrons de classe 1, 2 e 3. Alguns V.cholerae O1 ambientais pertencentes a linhagem da sétima pandemia, apresentaram uma extensiva perda de genes. Diferentes NAGs foram stn/sto+ e intI 1+. Dois alelos do gene aadA foram encontrados: aadA2 e aadA7. De modo interessante os V.cholerae O1 ambientais pertencentes à linhagem pandêmica, só foram isolados durante o período da última epidemia de cólera na região Amazônica brasileira (1991-1996).

Caracterização molecular e fenotípica de Salmonella Typhi isolada de casos de febre tifóide no Estado do Pará, no período de 1970 a 2009

A Salmonella Typhi é o agente etiológico da febre tifóide, uma doença infecciosa, sistêmica, que constitui importante problema de saúde pública principalmente nos países em desenvolvimento da África, Ásia, América Central e do Sul. Amostras de Salmonella Typhi isoladas de casos clínicos no período de 1970 a 2009 no Estado do Pará foram caracterizadas por meio da técnica de Eletroforese em Gel de Campo Pulsado (PFGE). O perfil de suscetibilidade aos antimicrobianos foi realizado através dos métodos manual e automatizado. Foi empregada a técnica de Reação em Cadeia Mediada pela Polimerase (PCR) para a pesquisa das integrases de classe 1 e 2 e do fator de virulência Vi. Em 151 amostras pôde-se observar 68 diferentes pulsotipos, sendo 66 deles agrupados em 5 clusters. As amostras independente de sua procedência, fonte de isolamento ou ano, apresentaram elevada similaridade genética que variou de 80 a 100%. Verificou-se resistência (1,99%) e resistência intermediária (6,62%) somente à nitrofurontoína, em nenhuma amostra foi detectado integrase de classe 1 e 2, demonstrando, que, no Estado do Pará, as cepas circulantes não apresentam multi-resistência como observado em várias regiões do mundo. Foram encontradas 4 (2,65%) amostras Vi-negativo, que pode ser devido ao longo período de armazenamento, pois a ilha de patogenicidade SPI-7 é geneticamente instável e pode ser perdida após sucessivos repiques.

Efeito Casimir em teorias de Chern-Simons

Nesta dissertação obtemos a força de Casimir a temperatura finita entre duas linhas paralelas sujeitas a condição de fronteira do tipo linhas mistas, no contexto da teoria de Maxwell- Chern-Simons em (2+1) dimensões. Além disso, analisamos a simetria de inversão de temperatura apresentada pela energia livre de Helmholtz do modelo para diferentes condições de fronteira. Iniciamos estudando aspectos gerais do formalismo de Matsubara no intuito de introduzirmos efeitos térmicos na teoria; também analisamos aspectos gerais da teoria de MCS em (2 + 1) dimensões. Posteriormente, revisitamos o cálculo da força de Casimir para o caso de duas linhas paralelas infinitamente permeáveis magneticamente a temperatura nula e finita, bem como o caso de linhas mistas a temperatura nula, onde tomamos uma linha perfeitamente condutora eletricamente e outra infinitamente permeável magneticamente. Em seguida, apresentamos novos resultados envolvendo a força de Casimir a temperatura finita com condições de fronteira do tipo linhas mistas. Por último, analisamos a simetria de inversão de temperatura associada a energia livre de Helmholtz do modelo, mostrando que mesmo para condições mistas e possível obter uma espécie de simetria residual, em analogia a resultados existentes para a eletrodinâmica em (3+1) dimensões.

Perfil de proteínas séricas em cães com linfoma multicêntrico de imunofenótipo B ou T: correlação com fatores prognósticos

Pós-graduação em Medicina Veterinária - FCAV

Transporte eletrônico em nanofitas de grafeno sob a influência de fatores externos, via primeiros princípios

O grafeno é a primeira estrutura bidimensional que se obteve experimentalmente. Sua rede cristalina é uma rede hexagonal, conhecida como "Favo de Mel", possui apenas um átomo de espessura. Cortes em folhas de grafeno, privilegiando determinada direção, geram as chamadas nanofitas de grafeno. Embora o grafeno se comporte como um metal, é sabido que as nanofitas podem apresentar comportamentos semicondutor, metálico ou semimetálico, dependendo da direção de corte e/ou largura da fita. No caso de nanofitas semicondutoras, a largura da banda proibida (band gap), entre outros fatores, depende da largura da nanofita. Neste trabalho adotou-se métodos de primeiros princípios como o DFT (Density Functional Theory), afim de se obter as características tais como curvas de dispersão para nanofitas. Neste trabalho, primeiramente, são apresentados diagramas de bandas de energia e curvas de densidade de estados para nanofitas de grafeno semicondutoras, de diferentes larguras, e na ausência de influências externas. Utilizou-se métodos de primeiros princípios para a obtenção destas curvas e o método das funções de Green do Não Equilíbrio para o transporte eletrônico. Posteriormente foi investigado a influência da hidrogenização, temperatura e tensão mecânica sobre sistema, isso além, de se estudar o comportamento de transporte eletrônico com e sem influência destes fatores externos. Vale ressaltar que as nanofitas de grafeno apresentam possibilidades reais de aplicação em nanodispositivos eletrônicos, a exemplo de nanodiodos e nanotransistores. Por esse motivo, é importante se ter o entendimento de como os fatores externos alteram as propriedades de tal material, pois assim, espera-se que as propriedades de dispositivos eletrônicos também sejam influenciadas da mesma maneira que as nanofitas.

Avaliação de métodos eletromagnéticos aplicando campos polarizados e focalizados

Em investigações geofísicas rasas que empregam os métodos eletromagnéticos indutivos mais avançados, alvos com baixo número de indução (Low Induction Number – LIN) produzem anomalias eletromagnéticas muito baixas e de difícil interpretação. Para suprir esta deficiência, neste trabalho são estudados a aplicabilidade de campos eletromagnéticos polarizados e focalizados – POLFOCEM como fonte primária de indução. Os campos E.M. focalizados e polarizados, vertical e horizontalmente, são obtidos pelas combinações vetoriais de pares de dipolos transmissores e, ocorrem na região central entre eles. A focalização é observada nesta região na profundidade de 0,25 do espaçamento entre esses transmissores – L. Portanto, máximos acoplamentos podem ser obtidos através da seleção da polarização de acordo com a geometria do alvo, ocorrendo um aumento na densidade de fluxo magnético sobre ele e, máximas anomalias produzidas. É utilizada uma metodologia numérica para o cômputo dessas anomalias por meio da técnica dos elementos finitos para solução do problema 2,5-D. Em todos os experimentos numéricos são realizadas comparações qualitativas e quantitativas entre as respostas obtidas pelos sistemas POLFOCEM e convencional, o qual emprega um único dipolo como transmissor (dipolo-dipolo). As anomalias produzidas pelo sistema POLFOCEM, em que os dipolos transmissores são acionados simultaneamente, correspondem à soma das anomalias produzidas por cada um desses dipolos independentes, caracterizando, desta forma, a linearidade dos campos eletromagnéticos. Os experimentos numéricos são realizados para alvos prismáticos bidimensionais com três diferentes inclinações, inseridos num semi-espaço resistivo, e para as freqüências das fontes na faixa das ondas de rádio. As anomalias assimétricas no sistema convencional, que se tornam simétricas no sistema POLFOCEM, apresentam valores menores em amplitude. Contudo, aquelas anomalias tanto assimétricas quanto simétricas que se tornam anti-simétricas apresentam valores maiores. Em decorrência dessas diminuições e aumentos nas amplitudes ocorrem rotações nos diagramas de Argand, no sentido horário e anti-horário para alvos com baixos valores de condutividade, respectivamente. Em experimentos de identificação de presença de dois alvos próximos, o sistema convencional é capaz de identificá-los primeiramente, prevalecendo o seu uso.

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.

Imageamento homeomórfico de refletores sísmicos

Neste trabalho é apresentada uma nova técnica para a realização do empilhamento sísmico, aplicada ao problema do imageamento de refletores fixos em um meio bidimensional, suavemente heterogêneo, isotrópico, a partir de dados de reflexão. Esta nova técnica chamada de imageamento homeomórfico tem como base a aproximação geométrica do raio e propriedades topológicas dos refletores. São utilizados, portanto, os conceitos de frente de onda, ângulo de incidência, raio de curvatura da frente de onda, cáustica e definição da trajetória do raio; de tal modo que a imagem obtida mantém relações de homeomorfismo com o objeto que se deseja imagear. O empilhamento sísmico é feito, nesta nova técnica de imageamento, aplicando-se uma correção local do tempo, ∆ t, ao tempo de trânsito, t, do raio que parte da fonte sísmica localizada em x_o, reflete-se em um ponto de reflexão, C_o, sendo registrado como uma reflexão primária em um geofone localizado em x_g, em relação ao tempo de referência t_o no sismograma, correspondente ao tempo de trânsito de um raio central. A fórmula utilizada nesta correção temporal tem como parâmetros o raio de curvatura R_o, o ângulo de emergência β_o da frente de onda, no instante em que a mesma atinge a superfície de observação, e a velocidade v_o considerada constante nas proximidades da linha sísmica. Considerando-se uma aproximação geométrica seguido um círculo para a frente de onda, pode-se estabelecer diferentes métodos de imageamento homeomórfico dependendo da configuração de processamento. Sendo assim tem-se: 1) Método Elemento de Fonte (Receptor) Comum (EF(R)C). Utiliza-se uma configuração onde se tem um conjunto de sismogramas relacionado com uma única fonte (receptor), e considera-se uma frente de onda real (de reflexão); 2) Método Elemento de Reflexão Comum (ERC). Utiliza-se uma configuração onde um conjunto de sismogramas é relacionado com um único ponto de reflexão, e considera-se uma frente de onda hipoteticamente originada neste ponto; 3) Método Elemento de Evoluta Comum (EEC). Utiliza-se uma configuração onde cada sismograma está relacionado com um par de fonte e geofone coincidentemente posicionados na linha sísmica, e considera-se uma frente de onda hipoteticamente originada no centro de curvatura do refletor. Em cada um desses métodos tem-se como resultados uma seção sísmica empilhada, u(x_o, t_o); e outras duas seções denominadas de radiusgrama, R_o (x_o, t_o), e angulograma, β_o(x_o, t_o), onde estão os valores de raios de curvatura e ângulos de emergência da frente de onda considerada no instante em que a mesma atinge a superfície de observação, respectivamente. No caso do método denominado elemento refletor comum (ERC), a seção sísmica resultante do empilhamento corresponde a seção afastamento nulo. Pode-se mostrar que o sinal sísmico não sofre efeitos de alongamento como consequência da correção temporal, nem tão pouco apresenta problemas de dispersão de pontos de reflexão como consequência da inclinação do refletor, ao contrário do que acontece com as técnicas de empilhamento que tem por base a correção NMO. Além disto, por não necessitar de um macro modelo de velocidades a técnica de imageamento homeomórfico, de um modo geral, pode também ser aplicada a modelos heterogêneos, sem perder o rigor em sua formulação. Aqui também são apresentados exemplos de aplicação dos métodos elemento de fonte comum (EFC) (KEYDAR, 1993), e elemento refletor comum (ERC) (STEENTOFT, 1993), ambos os casos com dados sintéticos. No primeiro caso, (EFC), onde o empilhamento é feito tendo como referência um raio central arbitrário, pode-se observar um alto nível de exatidão no imageamento obtido, além do que é dada uma interpretação para as seções de radiusgrama e angulograma, de modo a se caracterizar aspectos geométricos do model geofísico em questão. No segundo caso, (ERC), o método é aplicado a série de dados Marmousi, gerados pelo método das diferenças finitas, e o resultado é comparado com aquele obtido por métodos convecionais (NMO/DMO) aplicados aos mesmos dados. Como consequência, observa-se que através do método ERC pode-se melhor detectar a continuidade de refletores, enquanto que através dos métodos convencionais caracterizam-se melhor a ocorrência de difrações. Por sua vez, as seções de radiusgrama e angulograma, no método (ERC), apresentam um baixo poder de resolução nas regiões do modelo onde se tem um alto grau de complexidade das estruturas. Finalmente, apresenta-se uma formulação unificada que abrange os diferentes métodos de imageamento homeomórfico citados anteriormente, e também situações mais gerais onde a frente de onda não se aproxima a um círculo, mas a uma curva quadrática qualquer.

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.

Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.