Detecção de refletores sísmicos por rede neural discreta

As redes neurais artificiais têm provado serem uma poderosa técnica na resolução de uma grande variedade de problemas de otimização. Nesta dissertação é desenvolvida uma nova rede neural, tipo recorrente, sem realimentação (self-feedback loops) e sem neurônios ocultos, para o processamento do sinal sísmico, para fornecer a posição temporal, a polaridade e as amplitudes estimadas dos refletores sísmicos, representadas pelos seus coeficientes de reflexão. A principal característica dessa nova rede neural consiste no tipo de função de ativação utilizada, a qual permite três possíveis estados para o neurônio. Busca-se estimar a posição dos refletores sísmicos e reproduzir as verdadeiras polaridades desses refletores. A idéia básica desse novo tipo de rede, aqui denominada rede neural discreta (RND), é relacionar uma função objeto, que descreve o problema geofísico, com a função de Liapunov, que descreve a dinâmica da rede neural. Deste modo, a dinâmica da rede leva a uma minimização local da sua função de Liapunov e consequentemente leva a uma minimização da função objeto. Assim, com uma codificação conveniente do sinal de saída da rede tem-se uma solução do problema geofísico. A avaliação operacional da arquitetura desta rede neural artificial é realizada em dados sintéticos gerados através do modelo convolucional simples e da teoria do raio. A razão é para explicar o comportamento da rede com dados contaminados por ruído, e diante de pulsos fonte de fases mínima, máxima e misturada.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Refinamento da representação de raízes no modelo de biosfera SiB2 em área de floresta na Amazônia

O objetivo desta pesquisa foi estudar a sensibilidade dos fluxos de superfície e hidrologia do solo em relação à representação e distribuição de raízes no modelo de biosfera para uma floresta de terra firme na Amazônia. A finalidade foi avaliar o impacto na representatividade dos fluxos de energia considerando a sazonalidade da região amazônica, usando como suporte medidas intensivas realizadas em uma reserva biológica (Reserva Biológica do Cuieiras, em Manaus). Foram realizadas oito simulações com o modelo de biosfera SiB2 (“Simple Biosphere Model” – versão 2) , onde cada simulação representou um cenário diferente de distribuição de raízes em uma profundidade de 4 m de solo, dividido em três camadas: 0,5 m, 1,5 m e 2,0 m. As raízes foram distribuídas privilegiando a concentração de raízes na camada superficial, em seguida, na camada intermediária e, por fim, uma concentração de raízes abaixo de 2,0 m de profundidade. As simulações foram realizadas para o período de 2003 a 2006, enfatizando o ano de 2005 para avaliar o efeito da representação de raízes nos fluxos de energia (calor latente – LE e calor sensível – H) e de dióxido de carbono (CO₂). A partir da análise integrada dos fluxos simulados com dados observacionais medidos no sítio experimental foi possível perceber que uma redução na precipitação no ano de 2005, apesar de ter sido menor na parte central da Amazônia, implicou na diminuição da umidade do solo, mostrando que a floresta passou por um período de estresse hídrico maior do que os outros anos analisados. O modelo representou a energia disponível com valores muito próximos aos observados, variando sazonalmente em concordância com os dados medidos em 2005. No entanto, LE é superestimado durante a estação chuvosa, mas mostra juntamente com o fluxo de CO₂, a redução com a umidade do solo na estação seca, enquanto H é superestimado em até 20 W m^-2 durante todo o período simulado. Estes resultados mostram que, a consideração de raízes rasas é mais apropriada para regiões que possuem estação seca curta, conforme se caracteriza a área de estudo, e raízes profundas devem favorecer a modelagem dos processos para superfície de áreas com estação seca mais pronunciada. Com isto os resultados revelam que há necessidade de obter mais informações de propriedades físicas do solo, apropriadas às condições da região, para que outros refinamentos sejam efetivos na distinção do comportamento de florestas tropicais sob diferentes regimes de disponibilidade de água no solo.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Diet and niche breadth and overlap in fish communities within the area affected by an Amazonian reservoir (Amapá, Brazil)

Nós investigamos a amplitude e a sobreposição de nicho de espécies de peixes que ocorrem em quatro ambientes afetados pelo reservatório Coaracy Nunes, no Estado brasileiro do Amapá. Amostras sazonais de peixes foram coletadas usando uma configuração padrão com malhadeiras, tarrafas, linhas e armadilhas. Quinhentos e quarenta conteúdos estomacais, representando 47 espécies de peixes, foram analisados e quantificados. Amplitude e sobreposição de nichos foram estimadas usando os índices de Levin e Pianka, respectivamente, enquanto a competição interespecífica foi avaliada usando (RA3) modelo nulo. ANOVA e o teste de Kruskal¬Wallis foram usados para avaliar, respectivamente, as diferenças de amplitude e sobreposição de nichos entre áreas. Os dados indicaram que a maioria das espécies de peixes pertence às guildas de psicívoros, onívoros e detritívoros. Estas espécies provavelmente colonizaram os ambientes devido à disponibilidade de adequados recursos alimentares e às favoráveis condições físicas criadas pelo represamento do rio. De maneira geral, poucas espécies têm nichos amplos, mas muitas delas são altamente especializadas. Variação sazonal de recursos tem pouco efeito no comportamento alimentar da maioria das espécies nas áreas de estudo. Os modelos nulos indicaram que competição não foi um fator determinante na estrutura da comunidade.