Estudos de inclusões fluidas e de isótopos estáveis no depósito Moreira Gomes do campo mineralizado do Cuiú-Cuiú, Província Aurífera do Tapajós, Estado do Pará

Moreira Gomes é um dos depósitos do campo mineralizado do Cuiú-Cuiú, província Aurífera do Tapajós, com recursos de 21,7 t de ouro. A zona mineralizada, com 1200 metros de comprimento, 30-50 metros de largura e, pelo menos, 400 metros de profundidade é controlada por uma estrutura subvertical de orientação E-W, associada a um sistema de falhas transcorrentes sinistrais. As rochas hospedeiras nesse depósito são predominantemente tonalitos de 1997 ± 2 Ma (Suite Intrusiva Creporizão). O estilo da alteração hidrotermal relacionado à mineralização é predominantemente fissural e localmente pervasivo. Os tipos de alteração hidrotermal são sericitização, carbonatação, cloritização, sulfetação, silicificação e epidotização, além da formação de veios de quartzo de espessuras variadas. Pirita é principal sulfeto e contém inclusões de galena, esfalerita, calcopirita e, em menor quantidade, de hessita e bismutinita. O ouro ocorre mais comumente como inclusão em cristais de pirita e, secundariamente, na forma livre em veios de quartzo. Ag, Pb e Bi foram detectados por análise semi-quantitativa como componentes das partículas de ouro. Estudo de inclusões fluidas identificou fluidos compostos por CO₂ (Tipo 1), H₂O-C O₂-sal (Tipo 2) e H₂O-sal (Tipo 3). O volátil CO₂ é predominante na fase carbônica. O fluido do Tipo 2 apresenta densidade baixa a moderada, salinidade entre 1,6 e 11,8 % em peso equivalente de NaCl e foi aprisionado principalmente entre 280° e 350°C. No fluido do Tipo 3 o sistema químico pode conter aCl₂ e, talvez, MgCl₂, e a salinidade varia de zero a 10,1% em peso equivalente de NaCl. Apenas localmente a salinidade atingiu 25% em peso equivalente de NaCl. Esse fluido foi aprisionado principalmente entre 120° e 220°C e foi interpretado como resultado de mistura de fluido aquoso mais quente e levemente mais salino, com fluido mais frio e diluído. Globalmente, o estudo das inclusões fluidas indica estado heterogêneo durante o aprisionamento e ocorrência de separação de fases, mistura, flutuação de pressão e reequilíbrio das inclusões durante aprisionamento. A composição isotópica do fluido em equilíbrio com minerais hidrotermais (quartzo, clorita e calcita e pirita) e de inclusões fluidas apresenta valores de δ¹⁸O e δD entre +0,5 e +9,8 ‰, e -49 a -8 ‰, respectivamente. Os valores de ³⁴S de pirita (-0,29 ‰ a 3,95 ‰) são provavelmente indicativos da presença de enxofre magmático. Pares minerais forneceram temperaturas de equilíbrio isotópico em geral concordante com as temperaturas de homogeneização de inclusões fluidas e compatíveis com as relações texturais. Os resultados isotópicos, combinados com os dados mineralógicos e de inclusões fluidas são interpretados como produto da evolução de um sistema magmático hidrotermal em três estágios. (1) Exsolução de fluido magmático aquoso e portador de CO₂ entre 400°C e 320-350°C, seguido de separação de fases e precipitação principal da assembleia clorita-sericita-pirita-quartzo-ouro sob pressões menores que 2,1 kb e a 6-7 km de profundidade. (2) Resfriamento e continuação da exsolução do CO₂ do fluido magmático geraram fluido aquoso, mais pobre a desprovido de CO₂ e levemente mais salino, com aprisionamento dominantemente a 250°-280°C. A assembleia hidrotermal principal ainda precipitou, mas epidoto foi a principal fase nesse estágio. (3) Mistura do fluido aquoso do estágio 2, mais quente e mais salino, com um fluido aquoso mais frio e menos salino, de origem meteórica. Carbonatação está associada com esse estágio. A assembleia hidrotermal e os valores isotópicos indicam que fluido foi neutro a levemente alcalino e relativamente reduzido, que H₂S (ou HS-) pode ter sido a espécie de enxofre predominante, e que Au(HS) ^-2 deve ter sido o complexo transportador de ouro. A deposição do ouro em Moreira Gomes ocorreu em resposta a diversos mecanismos, envolvendo a separação de fases, mistura e reações fluido-rocha. O depósito Moreira Gomes é interpretado como o produto de um sistema magmático-hidrotermal, mas não possui feições clássicas de depósitos relacionados a intrusões graníticas, tanto oxidadas como reduzidas. A idade de deposição do minério (1,86 Ga) sugere que o sistema magmático-hidrotermal pode estar relacionado com a fase final do extenso magmatismo cálcio-alcalino da Suíte Intrusiva Parauari, embora o magmatismo transicional a alcalino da Suíte Intrusiva Maloquinha não possa ser descartado.

Relações ordinais sob controle contextual em crianças surdas

Dificuldades apresentadas por crianças surdas na aprendizagem da matemática têm conduzido educadores ao desenvolvimento de procedimentos especiais de ensino. O paradigma de equivalência tem sido útil na explicação de comportamentos complexos, como comportamentos conceituais numéricos. Uma expansão desse paradigma envolve a formação de classes de estímulos equivalentes em seqüência. A emergência de novas relações através do responder ordinal já foi documentada em estudos com contingências de reforçamento de três termos. Há necessidade de verificar se esses resultados se mantém estáveis sob contingências de quatro e cinco termos. Três experimentos foram programados com o objetivo de investigar a emergência de relações ordinais com controle discriminativo simples, sob controle condicional (sem e com randomização das tentativas) e sob controle contextual em crianças surdas. No Experimento 1 participaram cinco crianças surdas, matriculadas numa Escola Pública Especializada. Um microcomputador com um software (REL 4.0, utilizado nos Experimentos 1, 2a e 2b e atualizado para a versão 5.0 no Estudo 3) foi utilizado. Nesse estudo foi ensinado aos participantes seqüências de pares de estímulos sobrepostos. Em seguida, foram realizados testes de transitividade e conectividade. Todos os participantes alcançaram o critério de acerto e responderam aos testes. Os resultados replicaram estudos da literatura confirmando a eficiência do procedimento de ensino por sobreposição de estímulos no estabelecimento de relações ordinais. No Experimento 2a, quatro novos participantes e um com história experimental, foram ensinados a selecionar estímulos, aos pares, na ordem crescente na presença da cor verde e na ordem decrescente na presença da cor vermelha. Foram aplicados testes de transitividade e conectividade sob controle condicional. Em seguida foi conduzido um teste de generalização com estímulos do ambiente escolar. Todos os participantes alcançaram o critério de acerto e responderam aos testes de transitividade e conectividade. Nos testes de generalização, três participantes responderam consistentemente aos novos estímulos, um respondeu parcialmente e um não respondeu ao teste. Os resultados corroboraram a eficiência do procedimento de ensino por sobreposição de estímulos sob controle condicional em crianças surdas. O Experimento 2b envolveu os mesmos participantes do Experimento 1 com história experimental e a randomização das tentativas com os estímulos condicionais. Todos os participantes alcançaram o critério de acerto. Nos testes demonstraram um responder consistente com a linha de base. No Experimento 3 participaram três crianças dos Experimento 1 e duas do Experimento 2a, que foram expostas ao procedimento de ensino por pares sobrepostos sob controle contextual de duas formas círculo e triângulo e sob controle condicional das cores verde e vermelha (ex. A1A2, na presença do círculo e da cor verde; ou A2A1, na presença do círculo e da cor vermelha). Todos os participantes alcançaram o critério de acerto e responderam aos testes de transitividade e conectividade. Os resultados indicaram a eficiência do procedimento de ensino por sobreposição de estímulos sob controle contextual, sugerindo que o ensino por contingências de reforçamento simples e sob controle condicional foram pré-requisitos para a emergência de classes ordinais sob controle contextual. Uma extensão deste estudo deve ampliar o número de membros na seqüência e investigar a emergência de novas relações ordinais com seqüências mais longas, e verificar se a ordem de treino em que a seqüência é ensinada interfere sobre o responder ordinal.

Jornal do Pará: o caminho literário entre espaços e diálogos na Belém oitocentista

A divulgação de estudos em jornais do século XIX se torna cada vez mais necessária, pois contribui para um melhor entendimento do processo de produção, circulação e formação da História Literária Brasileira. Neste trabalho, analisa-se um dos principais jornais da imprensa paraense oitocentista - O Jornal do Pará (1862 -1878) - que destinou vários espaços para publicação literária e manteve estáveis diálogos com periódicos de outros estados, como o impresso carioca Jornal das Famílias - produzido pela Editora Garnier e dirigido sob a mesma linha conservadora do periódico paraense. Desse jornal familiar foram extraídos desde traduções a textos de autores consagrados, como Machado de Assis. Esses espaços reservados para a conversa entre os periódicos mostram que o Pará não ficou a parte das conjecturas de seu tempo, como também participou ativamente do cenário que ajudou a constituir a Literatura Brasileira.

Agricultura familiar na Amazônia Oriental: uma comparação dos resultados da pesquisa socioeconômica sobre fronteiras agrárias sob condições históricas e agroecológicas diversas

Na Amazônia Oriental, a agricultura familiar representa uma parte economicamente significativa do setor agrário. O enfoque dominante qualifica-a como uma agricultura nômade, ineficiente e insustentável, baseada no modelo do ciclo da fronteira (a agricultura familiar é expulsa pelas grandes fazendas de gado) e no modelo clássico de ecologia tropical da floresta amazônica. Resultados da pesquisa ecológica recente mostram que o modelo precisa ser revisado em vários aspectos relacionados diretamente com a sustentabilidade da agricultura familiar. Uma reconstrução histórica do desenvolvimento dos sistemas de produção de Igarapé-Açu e uma comparação com dados sobre outras fronteiras mais recentes no Pará mostram que, dependendo do desenvolvimento histórico específico da fronteira agrária, dos fatores condicionantes ecológicos e econômicos e das condições macroeconômicas, uma ampla diversidade de trajetórias relativamente estáveis do desenvolvimento agrário pode ser identificada.

Co-infecção vírus da imunodeficiência humana e vírus da hepatite C (HIV/HCV): aspectos epidemiológicos, clínicos e laboratoriais de uma população atendida em um serviço de hepatopatias na cidade de Belém-Pará

A infecção pelo Vírus da imunodeficiência humana 1 (HIV -1) em associação com a do Vírus da hepatite C (HCV) representa, atualmente, uma comorbidade que pode interferir principalmente na história natural da hepatite C. Este trabalho tem como objetivo descrever aspectos demográfico, clínico e laboratorial, incluindo exame histopatológico de pacientes coinfectados HIV/RCV. No período de agosto de 2004 a dezembro de 2006, 36 pacientes coinfectados, foram selecionados para o estudo. Noventa e dois por cento desses pacientes eram procedentes de Belém, com média de idade de 42 anos. Entre as principais informações demográficas da população estudada, foram identificados 72,52% solteiros, 83,5% do sexo masculino e 61,1% relataram ser heterossexuais. Entre os fatores de risco para o HCV o uso de drogas ilícitas injetáveis foi identificado em 41,7% dos casos, o uso de cocaína intranasal foi relatado por 38,9% dos pacientes, e o compartilhamento de seringas e material pessoal, em 38,9% dos casos. A história de etilismo em 77,8% e o uso de TARV foram os possíveis fatores agravantes mais frequentes para a doença hepática. Apenas um paciente apresentou sinais clínicos de insuficiência hepática crônica. Entre os testes bioquímicos hepáticos, a mediana de ALT e AST foi de 68UI/L e 61UI/L, respectivamente. Os níveis de linfócitos T CD4+ apresentaram mediana de 327 células/mm³, a carga viral do HIV com mediana de 2,53 logl0 cópias/mL (ep=0,34), carga viral do HCV com mediana de 5,9 log10UI/mL. O genótipo 1 do HCV foi o mais frequente (58,82%). Cinqüenta e sete por cento dos pacientes submetidos à biópsia hepática apresentavam fibrose de moderada a severa, e 11% não apresentaram fibrose pela classificação MET AVIR. Houve associação entre níveis de linfócitos T CD4+ e níveis de ALT e de AST (p=0,0009 e p=0,0002, respectivamente), assim como associação entre genótipo 1 do HCV e HCV-RNA maior ou igual a 6 log10 UI/mL (p=0.0039). Foi observada também associação entre HCV-RNA e HIV-RNA (p=0,039). Os pacientes apresentam estado geral bom, imunologicamente estáveis, sem sinais de descompensação hepática, mas com alterações estruturais hepáticas importantes, sendo portanto bons candidatos à terapia antiviral para o HCV. Futuros estudos, talvez de caso controle, com casuística maior são necessários para melhor entendimento da co-infecção HIV/HCV.

Alterações clínicas e metabólicas em portadores de hanseníase multibacilares

A hanseníase é doença infecto-contagiosa crônica causada pelo Mycobacterium leprae. Caracteriza-se por acometimento dermatoneurológico, variando em espectro entre dois polos estáveis (tuberculoide e virchoviano), com formas intermediárias instáveis. Uma classificação operacional, para fins de tratamento, reúne os doentes em dois grupos: paucibacilares (PB) que correspondem a formas clínicas que possuem 1-5 lesões e baciloscopia negativa; multibacilares (MB) que correspondem a formas clínicas com mais de 5 lesões e com ou sem baciloscopia positiva. Apesar de curável, a hanseníase ainda representa relevante problema de saúde pública. Sua maior morbidade associa-se aos estados reacionais e ao acometimento neural que podem causar incapacidades físicas e deformidades permanentes, comprometendo significativamente a qualidade de vida dos pacientes. Consequências clínicas, no que diz respeito as alterações oftalmológicas, endócrinas e cardiovasculares podem advir da etiopatogenia do processo infeccioso e imunopatológico, assim como dos efeitos adversos medicamentosos, desse modo tais eventos necessitam de esclarecimento afim de que o planejamento em saúde possa minimizar tais agravos. Um pronto diagnóstico, possibilita um tratamento precoce e eficaz, evitando com isso alterações clínicas importantes e sequelas. Foi realizado um estudo descritivo do tipo série de casos com 68 pacientes com alta terapêutica da hanseníase maior ou igual a 2 anos e com tratamento dos estados reacionais, tendo como objetivo descrever os aspectos clínicos, epidemiológicos e laboratoriais em pacientes multibacilares após alta com enfoque no diagnóstico de hipertensão, diabetes, osteoporose e discutir possíveis relações com tratamento dos episódios reacionais hansênicos. Observamos que a maioria dos pacientes eram do sexo masculino, com faixa etária acima de 45 anos, procedente de Belém, baixa escolaridade e de baixa renda familiar, o tipo de hanseníase predominante foi a forma virchowiana, a reação reversa foi o estado reacional mais prevalente, o corticoide foi o medicamento mais utilizado para o tratamento nos estados reacionais, os pacientes que não usaram corticoide apresentaram maior percentagem de densitometria normal, as comorbidades diabetes, hipertensão e osteoporose foram mais frequente em pacientes que usaram corticoide e com idade acima de 45 anos. O grau 2 foi o grau de incapacidade mais prevalente.

Expressão gênica e viabilidade de folículos pré-antrais inclusos em fragmentos de córtex ovarianos cultivados in vitro de macaco-prego (Sapajus apella)

O objetivo desse trabalho foi analisar a ativação e o crescimento de folículos préantrais de Sapajus apella, submetidos a um sistema de cultivo in vitro em curto-prazo. Nesse sentido, dois experimentos foram realizados neste trabalho. Experimento I: exposição à fresco e à crioprotetores de biopsias de fragmentos do córtex ovariano. Ambos os fragmentos ovarianos foram submetidos à extração total de RNA e síntese de cDNA. Após a amplificação do cDNA por PCR em tempo-real (RT-PCR), os software GeNorm, bestkeeper e Normfinder foram usados para avaliar a estabilidade dos genes GAPDH (glyceraldehyde-2-phosphate dehydrogenase), HPRT1 (hypoxanthine phosphoribosyltransferase 1) e TBP (TATA-binding protein). Experimento II: o tecido do córtex ovariano de quatro fêmeas foi oletado e dividido em nove pedaços de 1 mm³. Um fragmento ovariano (grupo controle) foi imediatamente dividido em dois pedaços, os quais foram destinados para análise da viabilidade ou por RT-PCR. Os 8 fragmentos restantes foram individualmente cultivados in vitro em um meio constituído de TCM suplementado com 100 ng/mL EGF (T1), com adição de 10 μM de BME (T2), 100 ng/mL de BMP4 (T3), 25 IU de PMSG (T4), 10 μM de BME e 100 ng/mL de BMP4 (T5), 10 μM de BME, 25 IU de PMSG (T6), 100 ng/mL de BMP4, 25 IU de PMSG (T7) ou 10 μM de BME, 100 ng/mL de BMP4, 25 IU de PMSG (T8). Os resultados demonstraram que no tecido do córtex ovariano de S. apella, os genes HPRT1 e TBP foram os mais apropriados como genes de referência, podendo ser usados como parâmetro para normalizar dados em estudos futuros. Ao contrário, o GAPDH se apresentou como menos estável dos genes de referencia testados. Após o cultivo in vitro, todos os tratamentos alcançaram percentuais similares da viabilidade de folículos pré-antrais. O tecido ovariano cultivado na presença de GF+BME/BMP4/PMSG resultou no aumento da taxa de ativação e crescimento folicular, assim como no aumento da expressão de AMH, BMP15 e GDF9, genes conhecidos como indicadores específicos de desenvolvimento folicular. Dessa forma, o HPRT1 e TBP são os genes de referência mais estáveis, na exposição à crioprotetores, a fresco e no o cultivo de tecidos de córtex ovariano de S. apella. Os folículos pré-antrais são capazes de desenvolverem-se in vitro quando cultivados em meio suplementado com PMSG, BME e BMP4. A viabilidade folicular, entretanto, permaneceu independentemente do meio de cultivo in vitro e o uso de fatores de crescimento, como marcadores de desenvolvimento folicular, foi crucial para identificar o melhor meio de cultivo in vitro.

Fonte escalar acoplada ao campo de Klein-Gordon orbitando um objeto estelar

Neste trabalho determinamos, utilizando Teoria Quântica de Campos em nível de árvore, a radiação escalar emitida por uma fonte em movimento circular uniforme no espaço-tempo plano de Minkowski, assumindo Gravitação Newtoniana, e no espaço-tempo curvo de um buraco negro sem carga e com momento angular nulo, assumindo Relatividade Geral. Efetuamos este cálculo analiticamente para o caso de Minkowski e numericamente no âmbito do espaço-tempo de Schwarzschild, sendo que neste espaço-tempo curvo obtivemos a forma analítica e a normalização dos modos nas regiões assintóticas. Verificamos que, para as órbitas circulares estáveis de acordo com a Relatividade Geral, a potência irradiada no caso de um buraco negro de Schwarzschild é menor do que a obtida no espaço-tempo de Minkowski assumindo a Gravitação Newtoniana. Obtemos também que apenas uma pequena parcela da radiação emitida é absorvida pelo buraco negro. Verificamos que a diferença entre as potências irradiadas em Schwarzschild e Minkowski diminui na medida em que aumentamos o valor da massa do campo. Em Schwarzschild, uma parcela cada vez maior da radiação emitida é absorvida pelo buraco negro na medida em que aumentamos o valor da massa do campo.

Planejamento de novos derivados da Edaravona

Um estudo teórico detalhado da edaravona foi realizado usando o método TFD (Teoria do Funcional da Densidade) com o conjunto de base B3LYP/6-31G*, com o objetivo de esclarecer o mecanismo sequestraste de radical e a influência do tautomerismo da edaravona em condições ácidas, importantes para auxiliar na elucidação do mecanismo de sequestro de radical pela edaravona no processo isquêmico. Em estudos teóricos anteriores, a tautomerização, o efeito do solvente e a abstração do elétron não foram considerados. Em nosso estudo, a estabilidade e reatividade foram determinadas através de parâmetros geométricos e energéticos. Os cálculos foram realizados em fase gasosa e o método PCM em fase aquosa e metanólica. As condições ácidas ou básicas foram consideradas pelas energias de dissociação e protonação, com consequente formação das formas aniônicas e catiônicas, respectivamente. As propriedades antioxidantes foram calculadas com base nos valores de HOMO, potencial de ionização (PI) e energia de dissociação da ligação (EDL). Os valores de HOMO e PI mostraram que o tautomero N-H é melhor antioxidante pela saída do elétron e os valores de EDL mostraram que o tautomero O-H é melhor antioxidante pela saída de hidrogênio. A protonação é mais favorecida termodinamicamente que a desprotonação. Além disso, a energia de protonação explica, teoricamente, a baixa diferença entre a protonação do N e O. O efeito do solvente diminuiu as barreiras de energias para a isomerização nos tautomeros O-H ou N-H. Adicionalmente, foram avaliados três derivados de pirazolonas com suas capacidades antioxidantes comparadas com a edaravona, em um esforço para identificar o farmacóforo antioxidante. A capacidade antioxidante da antipirina, dipirona, fenilbutazona e edaravona foram determinadas pela inibição de dois radicais estáveis DPPH e ABTS. A edaravona e fenilbutazona foram as mais potente para a inibição de radical DPPH e ABTS que a dipirona, enquanto que a antipirina não mostrou atividade em todas as concentrações analisadas. Simultaneamente, o método TFD previu o valor do potencial antioxidante in silico para explicar a relação de estrutura-atividade (REA). Além disso, foram calculados seus valores de HOMO, EDL X-H, e energia de estabilização (DEiso). Todos os cálculos foram executados usando o programa Gaussian 03, Hyperchem 7.5, e ChemOffice 2005. Os resultados mostraram que o derivado com C-H na posição C-4 aumentou a abstração do elétron ou de hidrogênio. Finalmente, a estratégia geral empregada para planejar os novos derivados se baseou nos estudos de dezoito derivados da edaravona descritos na literatura. Os estudos de REA classificaram os compostos em três grupos, como mais ativo, ativo e menos ativo. Nove derivados foram planejados a partir do composto mais ativo selecionado pelo estudo teórico.

Dissociação entre observação e interação na modificação do display agressivo do Betta splendens

A observação de luta (ou seja, a presença de um observador anônimo de uma luta) é algo que altera os parâmetros do display de luta do Betta splendens. O objetivo desse trabalho é identificar se e como a observação isolada ou associada a interação modifica a emissão do display. Foram utilizados Betta splendens (N=28) machos em uma bateria de experimentos a fim de identificar se a observação de um co-específico em situação de confronto modifica o display agressivo (Experimento I) ou se observação de um co-específico em situação de confronto e posterior exposição ao mesmo modifica o display agressivo (Experimento II). Não foram encontradas diferenças estatísticas nos dados obtidos através do experimento I, entretanto, com relação aos dados obtidos com o experimento II pode se identificar uma diferença estatística significativa (F(1,6)= 6,002; p= 0,05) refletida em um aumento da latência do Display Horizontal. Concluímos com isso que no Experimento I os comportamentos se mantiveram estáveis, e que a observação associada a interação (Experimento II) é capaz de modificar a emissão do display agressivo do Betta splendens.

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

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

Mapeamento do relevo do embasamento de bacias sedimentares através da inversão gravimétrica vinculada

Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada.

Migração por equação de onda em meios anisotrópicos com correção de amplitude

Este trabalho discute dois aspectos da migração em profundidade através da continuação para baixo dos campos de onda: o tratamento de modos evanescentes e a correção da amplitude dos eventos migrados. Estes dois aspectos são discutidos em meios isotrópicos e para uma classe de meios anisotrópicos. Migrações por diferenças finitas (FD) e por diferenças finitas e Fourier (FFD) podem ser instáveis em meios com forte variação lateral de velocidade. Estes métodos utilizam aproximações de Padé reais para representar o operador que descreve a propagação de ondas descendentes. Estas abordagens não são capazes de tratar corretamente os modos evanescentes, o que pode levar à instabilidades numéricas em meios com forte variação lateral de velocidade. Uma solução possível para esse problema é utilizar aproximação de Padé complexa, que consegue melhor representar os modos evanescentes associados às reflexões pós-críticas, e neste trabalho esta aproximação é utilizada para obter algoritmos FD e híbrido FD/FFD estáveis para migração em meios transversalmente isotrópicos com eixo de simetria vertical (VTI), mesmo na presença de forte variação nas propriedades elásticas do meio. A estabilidade dos algoritmos propostos para meios VTI foi validada através da resposta ao impulso do operador de migração e pela sua aplicação na migração de dados sintéticos, em meios fortemente heterogêneos. Métodos de migração por equação de onda em meios heterogêneos não tratam corretamente a amplitude dos eventos durante a propagação. As equações de onda unidirecionais tradicionais descrevem corretamente apenas a parte cinemática da propagação do campo de onda. Assim, para uma descrição correta das amplitudes deve-se usar as equações de onda unidirecionais de amplitude verdadeira. Em meios verticalmente heterogêneos, as equações de onda unidirecionais de amplitude verdadeira podem ser resolvidas analiticamente. Em meios lateralmente heterogêneos, essas equações não possuem uma solução analítica. Mesmo soluções numéricas tendem a ser instáveis. Para melhorar a compensação de amplitude na migração, em meios com variação lateral de velocidade, é proposto uma aproximação estável para solução da equação de onda unidirecional de amplitude verdadeira. Esta nova aproximação é implementada nas migrações split-step e diferenças finitas e Fourier (FFD). O algoritmo split-step com correção de amplitude foi estendido para meios VTI. A migração pré e pós-empilhamento de dados sintéticos, em meios isotrópicos e anisotrópicos, confirmam o melhor tratamento das amplitudes e estabilidade dos algoritmos propostos.

Análise de resolução em tomografia interpoços em 3-D

A aplicação de métodos de inversão aos dados de múltiplos levantamentos sísmicos interpoços permite a reconstrução de modelos de vagarosidade em 3-D de alta resolução adequados para monitoramento de processos de recuperação avançada de petróleo e caracterização de reservatórios. Entretanto, a falta de cobertura volumétrica uniforme de raios de levantamentos interpoços exige informação adicional ao sistema tomográfico para obtenção de soluções estáveis. A discretização do modelo em uma malha 3-D com células prismáticas triangulares e a decomposição em valores singulares são utilizadas para avaliar a reconstrução tomográfica em 3-D. O ângulo da projeção de modelos-alvo no subespaço ortogonal ao espaço nulo efetivo da matriz tomográfica é um critério adequado para se otimizar a malha de discretização do modelo interpretativo e a geometria de aquisição dos dados de modo a melhorar o condicionamento da reconstrução tomográfica. Esta abordagem pode ser utilizada durante as iterações lineares para redefinir a malha ou avaliar a necessidade de informação a priori adicional ao sistema tomográfico.

Inversão gravimétrica estável do relevo do embasamento e da variação da densidade com a profundidade em bacias sedimentares

O presente método postula uma variação hiperbólica para o contraste de densidade de uma bacia sedimentar em função da profundidade, e tem dois objetivos: (1) delinear o relevo do embasamento de uma bacia, conhecendo-se a anomalia gravimétrica, o contraste de densidade na superfície da bacia e o fator de decaimento do contraste de densidade com a profundidade; (2) estimar, além do relevo, o valor do contraste de densidade na superfície de uma bacia sedimentar e o fator de decaimento do contraste de densidade com a profundidade, sendo fornecida a anomalia gravimétrica e a profundidade do embasamento em alguns pontos da bacia. Nos dois casos o modelo interpretativo é um conjunto de prismas retangulares verticais justapostos, cujas espessuras, que são parâmetros a serem estimados, representam a profundidade da interface de separação entre os sedimentos e o embasamento. As soluções obtidas nos dois problemas acima formulados são estáveis devido à incorporação de informações adicionais sobre a suavidade do relevo estimado, e o conhecimento da profundidade do relevo do embasamento em alguns pontos, fornecido por furos de sondagem. O método foi testado em anomalias gravimétricas sintéticas produzidas pela simulação de bacias sedimentares com relevos suaves. Os resultados mostraram relevos com boa resolução e valores estimados do contraste de densidade na superfície da bacia e do fator de decaimento do contraste de densidade com a profundidade, próximos aos verdadeiros, indicando dessa maneira o potencial do método em interpretações gravimétricas de bacias sedimentares. O método foi aplicado à anomalia Bouguer da Bacia do Recôncavo, Brasil, delineando um relevo com um valor para a profundidade máxima de cerca de 6 km, semelhante ao estimado em interpretações sísmicas. As estimativas para o contraste de densidade na superfície da Bacia e o fator de decaimento com a profundidade foram -0,30 g/cm³ e 30 km, respectivamente, produzindo uma estimativa para a compactação máxima dos sedimentos da ordem de 4%.