Marcação e atenuação de múltiplas de superfície livre, processamento e imageamento em dados sísmicos marinhos

O presente trabalho teve como objetivos a análise e atenuação de múltiplas de superfície livre, o processamento e o imageamento de dados sísmicos marinhos visando obter imagens migradas de utilidade à interpretação geológica voltada à exploração de petróleo. Foi dada uma atenção sistemática ao estudo de múltiplas de superfície livre do ponto de vista do filtro de predição baseado na teoria da comunicação com o objetivo de melhor aplicar o filtro WH de deconvolução preditiva, na etapa de processamento, posterior à correção NMO, embora outros métodos possam ser considerados mais competitivos. A identificação e a atenuação de reflexões múltiplas em dados sísmicos reais continua como um grande desafio no processamento de dados sísmicos, uma vez que elas são consideradas como ruído. No entanto, sendo este ruído classificado como coerente, várias técnicas foram desenvolvidas objetivando sua atenuação para evitar erros em cascata nas etapas posteriores como processamento, marcação de eventos, inversão tomográfica, imageamento, e finalmente na interpretação geológica das imagens obtidas. Outro aspecto do trabalho foi estabelecer um fluxograma de processamento-imageamento sendo a atenuação de múltiplas de superfície livre uma etapa central. Foram obtidas seções migradas em tempo e em profundidade onde se permite a interpretação que se deseja. O desenvolvimento deste trabalho foi realizado com os pacotes CWP/SU e MatLab.

Métodos sísmicos de refração e eletromagnético aplicados ao estudo de um paleocanal da Ilha do Marajó, Pará

Os métodos sísmico de refração e eletromagnético foram aplicados, em caráter experimental, à prospecção de água subterrânea. A eficiência destes métodos foi testada com cerca de 1.500m de perfis de refração rasa na delimitação geométrica em três dimensões de um paleocanal na região dos campos naturais da Ilha de Marajó e com um perfil eletromagnético (técnica Slingram) para a detecção de contatos laterais de água doce-água salobra. A interpretação dos dados sísmicos foi realizada com as equações de tempo de intersecção para n camadas inclinadas com o auxílio do controle geológico de 2 (dois) poços para contornar problemas de camadas não detectadas. Obteve-se uma baixa velocidade (0,81 km/s) para os sedimentos silte-arenosos saturados que preenchem o paleocanal. Em virtude da inversão de velocidades ter sido obtida como resultado para estes sedimentos, não se aconselha o mapeamento de paleocanais pela sísmica de refração, salvo com a utilização de técnicas especiais de campo (registro digital, 'up-hole survey', registro ondas S, reflexão simultânea, ...). Os dados eletromagnéticos foram analisados semi-quantitativamente com o auxílio de um modelo reduzido construído a partir de uma seção de resistividades aparentes. Um diagrama de Argand, traçado com os dados de laboratório, fornece subsídios para a determinação da condutância em ambientes de sedimentos não consolidados análogos ao de Marajó. Recomenda-se a utilização dos métodos eletromagnéticos para o reconhecimento lateral de zonas de água doce-água-salobra, uma vez que a componente quadratura da frequência alta (2640 Hz) define, com precisão suficiente, estas zonas. A comparação dos resultados dos dois métodos mostra que apenas parte do paleocanal se comporta como um aquífero.

Modelagem direta de dados de eletrorresistividade 3-D

Dentre os métodos geofísicos, o da Eletrorresistividade é um dos mais tradicionais, com o seu desenvolvimento ocorrido há mais de 80 anos. Durante esse tempo o seu uso acompanhou o avanço no poder de processamento numérico e mais recentemente, a modelagem e inversão tridimensional tornou-se uma possibilidade para o geofísico. Apresentamos, neste trabalho, a técnica de elementos finitos aplicada ao método da eletrorresistividade 3-D, através do cálculo do potencial secundário. Para o desenvolvimento da metodologia, simulamos o levantamento do método da eletrorresistividade 3-D com os arranjos Dipolo-Dipolo e Schlumberger, visando medir as variações laterais e verticais da resistividade aparente do solo. Estes arranjos consistem na injeção de corrente elétrica na superfície e de medidas de diferenças de potencial elétrico, resultante da interação da corrente elétrica com o solo. Sendo que, as fontes e receptores são localizados de acordo com os arranjos escolhidos para o levantamento. Neste trabalho, as curvas de sondagem e as pseudo-secções de resistividade aparente, são obtidas através da modelagem de eletrorresistividade 3-D, usando malha de elementos finitos regular. Para efeito de validação, os resultados são comparados com a resposta 3-D obtida a partir dos potenciais totais.

Estudos isotópicos (Pb, O, H, S) em zonas alteradas e mineralizadas do depósito cupro-aurífero Visconde, Província Mineral de Carajás

O depósito cupro-aurífero Visconde está localizado na Província Mineral de Carajás, a cerca de 15 km a leste do depósito congênere de classe mundial Sossego. Encontra-se em uma zona de cisalhamento de direção WNW-ESE, que marca o contato das rochas metavulcanossedimentares da Bacia Carajás com o embasamento. Nessa zona ocorrem outros depósitos hidrotermais cupro-auríferos com características similares (Alvo 118, Cristalino, Jatobá, Bacaba, Bacuri, Castanha), que têm sido enquadrados na classe IOCG (Iron Oxide Copper-Gold), embora muitas dúvidas ainda existam quanto a sua gênese, principalmente no que diz respeito à idade da mineralização e fontes dos fluidos, ligantes e metais. O depósito Visconde está hospedado em rochas arqueanas variavelmente cisalhadas e alteradas hidrotermalmente, as principais sendo metavulcânicas félsicas (2968 ± 15 Ma), o Granito Serra Dourada (2860 ± 22 Ma) e gabros/dioritos. Elas registram diversos tipos de alteração hidrotermal com forte controle estrutural, destacando-se as alterações sódica (albita + escapolita) e sódico-cálcica (albita + actinolita ± turmalina ± quartzo ± magnetita ± escapolita), mais precoces, que promoveram a substituição ubíqua de minerais primários das rochas e a disseminação de calcopirita, pirita, molibdenita e pentlandita. Dados isotópicos de oxigênio e hidrogênio de minerais representativos desses tipos de alteração mostram que os fluidos hidrotermais foram quentes (410 – 355°C) e ricos em ¹⁸O (δ¹⁸O_H2O= +4,2 a 9,4‰). Sobreveio a alteração potássica, caracterizada pela intensa biotitização das rochas, a qual ocorreu concomitantemente ao desenvolvimento de foliação milonítica, notavelmente desenhada pela orientação de palhetas de biotita, que precipitaram de fluidos com assinatura isotópica de oxigênio similar à dos estágios anteriores (δ¹⁸O_H2O entre +4,8 e +7,2‰, a 355°C). Microclina e alanita são outras fases características desse estágio, além da calcopirita precipitada nos planos da foliação. A temperaturas mais baixas (230 ± 11°C), fluidos empobrecidos em ¹⁸O (δ¹⁸O_H2O = -1,3 a +3,7‰) geraram associações de minerais cálcico-magnesianos (albita + epidoto + clorita ± calcita ± actinolita) que são contemporâneas à mineralização. Valores de δ¹⁸D_H2O e δO_H2O indicam que os fluidos hidrotermais foram inicialmente formados por águas metamórficas e formacionais, a que se misturou alguma água de fonte magmática. Nos estágios tardios, houve considerável influxo de águas superficiais. Diluição e queda da temperatura provocaram a precipitação de abundantes sulfetos (calcopirita ± bornita ± calcocita ± digenita), os quais se concentraram principalmente em brechas tectônicas - os principais corpos de minério - que chegam a conter até cerca de 60% de sulfetos. Veios constituídos por minerais sódico-cálcicos também apresentam comumente sulfetos. A associação de minerais de minério e ganga indica uma assinatura de Cu-Au- Fe-Ni-ETRL-B-P para a mineralização. Os valores de δ³⁴S (-1,2 a +3,4‰) de sulfetos sugerem enxofre de origem magmática (proveniente da exsolução de magmas ou da dissolução de sulfetos das rochas ígneas pré-existentes) e precipitação em condições levemente oxidantes. Datação do minério por lixiviação e dissolução total de Pb em calcopirita forneceu idades de 2736 ± 100 Ma e 2729 ± 150 Ma, que indicam ser a mineralização neoarqueana e, a despeito dos altos erros, permite descartar um evento mineralizador paleoproterozoico. A idade de 2746 ± 7 Ma (MSDW=4,9; evaporação de Pb em zircão), obtida em um corpo granítico não mineralizado (correlacionado à Suíte Planalto) que ocorre na área do depósito, foi interpretada como a idade mínima da mineralização. Assim, a formação do depósito Visconde teria relação com o evento transpressivo ocorrido entre 2,76 e 2,74 Ga, reponsável pela inversão da Bacia Carajás e pela geração de magmatismo granítico nos domínios Carajás e de Transição. Esse evento teria desencadeado reações de devolatilização em rochas do Supergrupo Itacaiúnas, ou mesmo, provocado a expulsão de fluidos conatos salinos aprisionados em seus intertícios. Esses fluidos teriam migrado pelas zonas de cisalhamento e reagido com as rochas (da bacia e do embasamento) pelas quais se movimentaram durante a fase dúctil. As concentrações subeconômicas do depósito Visconde devem ser resultado da ausência de grandes estruturas que teriam favorecido maior influxo de fluidos superficiais, tal como ocorreu na formação dos depósitos Sossego e Alvo 118.

Contribuição dos sedimentos e ictiólitos para a identificação do ambiente de formação do Sambaqui do Moa (Saquarema-RJ)

O Sambaqui do Moa é um sítio arqueológico localizado em Itaúna, Saquarema, litoral do Estado do Rio de Janeiro. Três momentos de ocupação foram reconhecidos neste sítio: o estrato 3 na base correspondente ao início da ocupação deste sítio; o 2, intermediário, aponta para uma ocupação mais intensa, com grande concentração de carapaças de moluscos, ossos de peixes e sepultamentos humanos; e o 1, o mais superficial, relacionado à última ocupação. Para identificar as condições ambientais de desenvolvimento do Sambaqui do Moa foram coletadas amostras de sedimentos e material zooarqueológico nestes três estratos. O material zooarqueológico está representado por restos microscópicos de peixes (ictiólitos), compostos por microdentes com diferentes morfologias: caninos, incisivos e molares. Os sedimentos segundo análises por DRX são constituídos além de quartzo e caulinita, calcita, aragonita e por fluorapatita. Este último é o principal mineral do material zooarqueológico, enquanto calcita e aragonita refletem os restos de conchas contidos neles, abundantes no sítio. As análises mineralógicas foram confirmadas pelas análises químicas, em que os teores elevados de P₂O₅, CaO e PF (H₂O, CO₂), respondem pela fluorapatita, calcita e aragonita. Modificações químicas representadas pelas variações nos conteúdos de C e P nos microdentes sugerem que estes experimentaram mineralização, num processo inicial de fossilização, pósdeposição. Os dados de isótopos estáveis ¹³C e ¹⁵N permitiram definir a fonte de matéria orgânica do sambaqui do Moa como marinha/salobra, em que a vegetação representava-se, predominantemente, por plantas do tipo C3 de floresta tropical. As razões ⁸⁷Sr/⁸⁶Sr nos ictiólitos confirmam que o ambiente no entorno do sambaqui tenha sido estuarino. A morfologia dentária permitiu reconhecer cinco famílias até então não registradas para o sítio, como Labridae, Serranidae, Ariidae, Erythrinidae e Characidae, que confirmam o ambiente estuarino. A idade do sambaqui do Moa por datação radiocarbono a partir dos sedimentos mostrou perturbação no estrato 2, ocasionando a inversão das idades entre os estratos o que pode ser justificado por processos de formação e/ou mudanças dos rios que modificaram as configurações geológico-geomorfológicas da área. Outra explicação poderia ser a interferência humana, devido ao grande número de enterramentos (mais de 30), tenha perturbado a ordem dos momentos de ocupação do sambaqui do Moa, além de possíveis processos erosivos. O sambaqui do Moa se instalou, portanto, em uma área de transição marinho-estuarina.

Evolução geológica das sequências do embasamento do Cinturão Araguaia na região de Paraíso do Tocantins (TO), Brasil

Estudos isotópicos baseados nas metodologias Pb-Pb em zircão e Sm/Nd (rocha total), permitiram avanços no entendimento do quadro geológico evolutivo e litoestratigráfico do embasamento do segmento sul do Cinturão Araguaia. Os processos geológicos identificados aconteceram a partir do Arqueano (2,6 Ga e T_DM 2,78 – 3,25 Ga) e se estenderam até o Neoproterozoico. Os ortognaisses do Complexo Rio_dos Mangues posicionam-se no Paleoproterozoico (2,05 – 2,08 Ga) e T_DM 2,35 – 2,21 Ga. Um forte encurtamento crustal e fusão parcial de compartimentos isolados e espessados, gerou corpos ígneos de 1,85 e 1,82 Ga e o Granito Serrote (1,86 Ga), que provêm de fontes entre 2,50 e 2,43 Ga. No final do Mesoproterozoico a região foi marcada por processos tafrogenéticos, evidenciados por magmatismo máfico e alcalino (1,05 Ga) e bacias deposicionais, como a que acolheu os sedimentos que originaram as supracrustais do Cinturão Araguaia. No Neoproterozoico, através da inversão nas condições geodinâmicas, ocorreu novo processo de encurtamento/espessamento crustal com fusões que originaram expressivas massas batolíticas (Granitos Matança e Santa Luzia). O Cinturão Araguaia foi edificado a partir dessa movimentação tectônica. O transporte de massas tectônicas no sentido do Cráton Amazônico teria ocorrido, resultando na atual arquitetura em que se encontram as várias unidades lito-estratigráficas, organizadas sob a forma de lascas imbricadas.

Basaltos almofadados da Suíte Ofiolítica Morro do Agostinho: registros de fundo oceânico na porção centro-oeste do Cinturão Araguaia

Ao longo do domínio de baixo grau metamórfico (porção centro-oeste) do Cinturão Araguaia, afloram dezenas de corpos máficos e/ou ultramáficos de natureza ofiolítica. Cita-se como exemplo a Suíte Ofiolítica Morro do Agostinho nos arredores da cidade de Araguacema (TO) que configura um pequeno corpo isolado que sustenta o Morro do Agostinho e encontra-se encaixado tectonicamente em metarenitos, ardósias e filitos da Formação Couto Magalhães (Grupo Tocantins). A Suíte Ofiolítica Morro do Agostinho é constituída por peridotitos serpentinizados, basaltos e cherts ferríferos todos afetados por incipiente metamorfismo. A associação de basaltos é caracterizada por um expressivo derrame submarino com estruturas em lavas almofadadas, sobrepostas aos peridotitos serpentinizados. Os basaltos foram classificados em tipos maciços e hipovítreos com esferulitos. Os basaltos maciços são homogêneos, com textura intersertal definida, essencialmente, por finas ripas de plagioclásio, clinopiroxênio e raramente olivina, calcocita e calcopirita. Os basaltos hipovítreos apresentam feições texturais formadas por ultrarresfriamento de lavas apresentando esferulitos de plagioclásio, feixes de cristais aciculares e esqueletais de clinopiroxênio e plagioclásio, e cristais com terminações tipo rabo-de-andorinha. Geoquimicamente, os basaltos revelaram natureza subalcalina toleítica, compatíveis com o tipo MORB. As razões (La/Yb)_n < 1 e (La/Sm)_n < 1 apontam, mais especificamente, para magmas do tipo N-MORB na evolução dessas rochas relacionadas ao ambiente de fundo oceânico. Estas rochas revelaram que nos estágios iniciais da evolução do Cinturão Araguaia houve uma fase importante de oceanização da Bacia Araguaia, com exposição de peridotitos do manto litosférico seguido de extravasamento de lavas e sedimentação de cherts e formações ferríferas bandadas em ambiente oceânico profundo. Após o preenchimento sedimentar da Formação Couto Magalhães (Grupo Tocantins), e o descolamento da litosfera oceânica, a fase tectônica principal propiciou a inversão tectônica que levou à exumação dos corpos ofiolíticos, principalmente ao longo de superfícies de cavalgamento, fragmentando-os e misturando-os tectonicamente às rochas supracrustais, acompanhado de metamorfismo regional em condições da fácies xisto verde baixo. A Suíte Ofiolítica Morro do Agostinho representa, assim, um pequeno fragmento alóctone de um segmento litosférico manto/crosta oceânica, bem preservado, do início da evolução da Bacia Araguaia, similar a outros no Cinturão Araguaia, que é um importante registro da fase de oceanização do Cinturão Araguaia, durante o Neoproterozoico.

Paleoambiente da Formação Prosperança, embasamento neoproterozóico da Bacia do Amazonas

Rochas de idade neoproterozóica da Formação Prosperança, cobertura sedimentar da porção sul do Escudo das Guianas, são pobremente expostas quando comparadas com o registro paleozóico das bacias do Amazonas e Solimões. A Formação Prosperança consiste em conglomerados, arenitos arcosianos e pelitos que preenchem grábens no embasamento. Esta unidade é sotoposta por rochas carbonáticas da Formação Acarí (Neoproterozóico), observada apenas em subsuperfície, que agrupadas, representam o embasamento sedimentar das bacias paleozóicas produtoras de óleo da Amazônia. A precisa caracterização e reconstrução paleoambiental da Formação Prosperança são essenciais para a distinção entre unidades do embasamento sedimentar e paleozóicas. A análise estratigráfica foi realizada na região do baixo rio Negro, Estado do Amazonas. A Formação Prosperança consiste em quatro associações de fácies que foram interpretadas como produto de um sistema flúvio-deltaico: prodelta/lacustre, frente deltaica, foreshore/shoreface e planície braided distal. Camadas tabulares de pelitos distribuídos por quilômetros sugerem uma bacia sedimentar de provável origem lacustre/mar restrito. Lobos deltaicos complexamente estruturados foram alimentados por distributários braided que migravam para SE. Arenitos gerados sob condições de fluxo oscilatório/combinado são compatíveis com depósitos de face litorânea. Arenitos com estratificação cruzada e planar estão relacionados com a migração de dunas subaquosas associadas a processos fluviais braided. Camadas lenticulares de conglomerados e arenitos com estratificação cruzada e planar, de possível idade paleozóica, sobrepõem a Formação Prosperança erosivamente. Essas camadas são produtos de um sistema fluvial braided proximal que migrava para NW, sentido inverso dos valores de paleocorrentes dos arenitos Prosperança.

Considerações filogenéticas e biogeografia histórica dos malacostráceos (decápodes e isópodes) cenozóicos do Brasil

Análises da Biogeografia Histórica da paleocarcinofauna cenozóica brasileira, encerrada nas formações Maria Farinha (Paleoceno), Tremembé (Oligoceno) e Pirabas (Mioceno), revelaram que as mesmas guardam afinidades com as cretáceas do mar de Tethys e possuem aspecto moderno. Sua distribuição paleobiogeográfica permite enquadrar o seu surgimento, dispersão e irradiação nos padrões fauna de origem tetiana, fauna de origem de alta latitude sul e fauna anfitropical, sendo que a maior parte tem afinidade com a fauna de origem tetiana, que teve basicamente dois caminhos de dispersão, migração no sentido oeste a partir do sul da Europa e Mar de Tethys, para leste atingindo a costa leste dos Estados Unidos e região caribeana, e em seguida, chegando na costa oeste dos Estados Unidos e América do Sul, pelo corredor americano central; e para leste da Europa e Mar de Tethys e Japão, a partir da costa leste dos Estados Unidos, em um caminho inverso ao primeiro. As relações filogenéticas definidas até o momento sugerem que a superfamília Thalassinoidea originou-se provavelmente de um grupo primitivo da infraordem Caridea, e que a família Retroplumidae (gênero Costacopluma) presente na Formação Maria Farinha, é ancestral direta dos ocipodídeos cenozóicos (representados na Formação Pirabas pelo gênero Uca), ocupando os mesmos nichos ecológicos e com tolerâncias ambientais similares. Também as similaridades morfológicas sugerem monofiletismo entre Xanthoidea e Portunoidea, que a partir do Paleoceno, os goneplacídeos, também monofilético, originou outro grupo-irmão, a família Hexapodidae, e o gênero Glyphithyreus possui caracteres afins da subfamília Eucratopsinae, membro do grande grupo Xanthoidea.

Sazonalidade no balanço de energia em áreas de cultivo de soja na Amazônia

Neste trabalho, avaliou-se o balanço de energia na cultura soja (Glycine max (L.) Merrill), variedade Tracajá, durante a safra e entressafra, em uma área de avanço de fronteira agrícola na Amazônia, por meio do método da razão de Bowen. Durante quase todo o ciclo da cultura, maior parte da energia foi consumida na forma de calor latente, principalmente durante as fases de florescimento e frutificação. Tal característica esteve relacionada à elevada condutância estomática foliar da soja e à disponibilidade de água na região. Próximo ao fim do ciclo ocorreu uma inversão na partição de energia entre os componentes H e LE, quando maior parte da energia foi usada no aquecimento do ar (79% do saldo de radiação). Durante a entressafra observou-se a redução de 75% no fluxo de calor latente e aumento considerável de 180% no fluxo de calor sensível comparado ao encontrada durante o ciclo.

Efeito do lodo de esgoto sobre a atividade microbiana e fitopatógenos habitantes do solo

Pós-graduação em Agronomia (Proteção de Plantas) - FCA

Influências das condições termodinâmicas nos eventos extremos de precipitação na cidade de Belém e região metropolitana

O objetivo desta dissertação foi verificar a importância dos fatores termodinâmicos na ocorrência de eventos extremos de precipitação na Cidade de Belém (PA) e região metropolitana, no período de agosto de 2008 a dezembro de 2009. Para tal, foram utilizados dados de precipitação e radiossondagens. Para análise das condições termodinâmicas e dinâmicas foram utilizadas imagens de satélite, cartas de superfície e análise do diagrama SKEW T LOG P. O estudo da termodinâmica da atmosfera foi realizado a luz das teorias da Energia Potencial Disponível para a Convecção (CAPE) e, também, da Energia de Inibição da Convecção (CINE). Foi utilizado o método dos decis para classificar os eventos extremos de precipitação a fim de associá-los aos valores da CAPE e da CINE com o objetivo de verificar o valor destes índices quando da ocorrência dos eventos extremos. Verificou-se que a região estudada possui forte atividade convectiva durante todo o ano, haja vista que seus valores médios mensais variam entre 900 J/kg e 1900 J/kg. Foi visto, também, que nem sempre CAPE alta e CINE baixa determinam precipitação. Esta situação determina o potencial para a convecção profunda, mas para converter este potencial em precipitação existe a necessidade da forçante dinâmica. Os resultados mostraram que quando o processo de precipitação dependeu, exclusivamente, da CAPE, foi necessário haver um valor alto para poder gerar convecção profunda e por consequência precipitação, enquanto, que no processo de precipitação com contribuição dinâmica não foi necessário um valor tão significativo da CAPE, neste caso, não ultrapassou a 1000 J/kg. A CINE esteve, sempre, menor no período chuvoso apresentando valores médios mensais menores que 300 J/kg. Isto não quer dizer que quanto menor a CINE maior será a precipitação. Quando a inibição está presente a instabilidade vai crescendo ao longo do dia determinando, com isso, nuvens com um desenvolvimento vertical mais acentuado, assim os pontos onde os inibidores enfraquecem primeiro, serão os pontos preferenciais para o disparo da tempestade. Logo, quando a instabilidade estiver alta e existir o mecanismo inibidor (CINE), em uma grande área, os locais mais propícios aos disparos das tempestades são os pontos onde a CINE e o NCE começam a diminuir e a inversão térmica, que por vezes acontece, começa a ser quebrada. Durante a execução desta pesquisa ficou claro que para a ocorrência de eventos extremos de precipitação, no período chuvoso, existe necessidade da influência da ITCZ e no período seco, conforme se observa no estudo de caso realizado para o mês de outubro o fator dinâmico que mais influencia é a Linha de Instabilidade (LI).

Inversion method for obtaining electrical parameters for soil in the rural region in Brazil

This paper presents an inversion methodology through weighted least squares to obtain the electrical parameters for the soil of a typical mid-western region in Brazil using the model based in the formalism of the parabolic equations to calculate the electric field intensity received. To validate this methodology, the results of the radio signal measurement campaign conducted in six radial routes leaving the city of Brasilia, Federal District, where the transmitter was located, were used. The measurements were compared to computer simulations and, thus, the optimal values for the electric conductivity and relative permittivity for the soil of the region could be estimated. Finally, a quantitative analysis of these parameters was performed with the values found in the literature, which demonstrated the effectiveness of the proposed methodology.

Complexos de lantanídeos contendo ligantes organofosforados como matrizes para materiais luminescentes

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

Satisfação com a vida em idosos

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