Avaliação dinâmica da interligação Tucuruí-Manaus-Macapá frente a pequenas e grandes perturbações

Esta dissertação apresenta uma avaliação dinâmica da conexão Tucuruí-Manaus-Macapá ao Sistema Interligado Nacional frente a pequenas e grandes perturbações. Todos os dados utilizados foram obtidos da base de dados do Operador Nacional do Sistema Elétrico. No entanto, os parâmetros dos estabilizadores dos sistemas de potência das usinas de Macapá não estão sintonizados para a interconexão e as usinas de Manaus não apresentam dispositivos suplementares amortecedores. Sendo assim, o objetivo deste trabalho é projetar estabilizadores de sistemas de potência para as usinas hidrelétricas de Coaraci Nunes e Balbina, e também para as usinas térmicas de Mauá e Aparecida, através de técnicas de respostas em frequência baseada no critério de Nyquist. A utilização dos estabilizadores tem como finalidade o amortecimento dos modos eletromecânicos que surgem quando os sistemas isolados do Amapá e do Amazonas passam a ser integrados ao sistema interligado nacional. Os projetos são avaliados em diversos cenários operativos através de simulações no domínio do tempo.

Projeto e implementação em DSP de sistema de comunicação multiportadora baseado no padrão G.fast

As redes de acesso usando cabos de cobre atualmente utilizam bandas de frequência até 30 MHz, especificada no padrão VDSL2. À medida que arquiteturas híbridas de fibra e cobre se tornam mais proeminentes na indústria e academia, torna-se possível utilizar cabos metálicos mais curtos (i.e. até 250 metros) conectando o último ponto de distribuição aos usuários, de modo que frequências mais altas podem ser exploradas para se alcançar taxas de transmissão de dados de 500 Mbps ou mais, como é o caso do padrão G.fast atualmente em desenvolvimento no ITU-T. Nesse trabalho, um simulador no domínio do tempo foi desenvolvido para avaliar a capacidade do sistema G.fast com diferentes tamanhos de extensão cíclica e diferentes topologias de rede especificadas pelo ITU-T. Os resultados das simulações mostram que sistemas G.fast são robustos a bridged taps e capazes de atingir altas taxas de dados para todas as topologias simuladas, provendo suporte à próxima geração de serviços de banda larga. Além disso, esse trabalho descreve o progresso da implementação de um protótipo de modem baseado no padrão G.fast em um ambiente híbrido de DSP multicore e FPGA utilizando kits de avaliação adquiridos pela UFPA. Arquiteturas, protocolos de comunicação e benchmarks são apresentados e avaliados para se chegar à conclusão de que tal protótipo é factível e fornece suporte flexível a várias linhas de pesquisa em banda larga da próxima geração.

Análise bidimensional dos efeitos das descargas atmosféricas em cabos OPGW utilizando o método ADI-FDTD

Neste trabalho é analisado o comportamento da densidade de corrente em cabos OPGW quando submetidos a descargas atmosféricas. Para esta análise, foram aplicados dois modelos de corrente, uma vez que a descarga atmosférica tem duas componentes principais que causam mais danos aos cabos em geral: a componente impulsiva e a componente contínua. O método numérico utilizado para as simulações computacionais é o método ADI-FDTD (Diferenças Finitas no Domínio do Tempo de Direções Alternadas Implícitas), o qual é truncado pela técnica CPML (Convolutional Perfectly Matched Layers). Com base no comportamento da corrente no cabo, é feita uma avaliação para determinar qual das duas componentes de descarga atmosférica pode causar mais danos aos cabos OPGW.

Análise da influência do monitoramento OTDR em WDM-PON: um estudo baseado em simulação

Redes Ópticas Passivas estão cada vez mais difundidas como rede de acesso de banda larga. Devido à sua importância, tais redes necessitam constantemente de supervisão. Diversas técnicas vêm sendo empregadas para supervisão de redes ópticas, dentre as quais se destacam as que utilizam reflectometria, pois agregam vantagens, tais como: baixo custo, em comparação com outros métodos de supervisão, e não exigência de intervenção na casa do usuário. Esta dissertação apresenta uma estratégia capaz de quantificar, por meio de simulações, o impacto que a técnica de reflectometria no domínio do tempo gera como interferência em transmissões típicas de uma rede WDM-PON.

Comparação do desempenho dos métodos eletroresistivo e transiente eletromagnético num caso específico de prospecção de águas subterrâneas

O método eletroresistivo é um dos métodos de prospecção geofísica mais empregados na exploração de água subterrânea. Todavia, ele apresenta a desvantagem de ser de difícil uso em locais habitados. De fato, para realizar sondagens elétricas é necessário, muitas vezes, expandir cabos elétricos até um quilômetro ou mais, dependendo da profundidade que se deseja investigar, o que nem sempre é possível em áreas urbanas. O método transiente eletromagnético no domínio do tempo, é usado também para a exploração de água subterrânea e não apresenta as mesmas desvantagens, porque este método indutivo emprega uma bobina quadrada cujo tamanho não é um impedimento para ser utilizado em locais habitados. Por outro lado, o método eletromagnético é mais sensível a ruídos devido a instalações elétricas. Para comparar as vantagens e desvantagens dos dois métodos em áreas urbanas foram utilizados programas para a modelagem sintética das sondagens elétricas verticais e das sondagens transientes eletromagnéticas, ambos desenvolvidos por RIJO (1994ab, 1996). Para tal propósito foram considerados dois casos: no primeiro, foram definidos quatro tipos de modelos (H, K, A e Q) com base nos horizontes geoelétricos em Palestina do Pará. No segundo caso, se considerou os modelos obtidos pelas 21 sondagens elétricas verticais tipo Schlumberger da mesma área. Com base nestes modelos, se utilizou um programa de transientes electromagnéticas para gerar as curvas sintéticas das resistividades aparentes. Como resultado, os valores de resistividades aparentes dados em função do espaçamento de eletrodos (método eletroresistivo) e em função do tempo (método eletromagnético), permitiu obter nos dois casos mencionados uma relação linear em escala bilogaritmica entre o espaçamento (metros) e o tempo (milisegundos). Deste modo é possível inferir a profundidade de investigação desejada em áreas urbanas usando o método eletromagnético, quando a aplicação do método eletroresistivo ficar limitada pela acessibilidade.

Comparação dos filtros de velocidade e do operador WHLP-CRS na atenuação de múltiplas

A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.

Imageamento multifoco de refletores sísmicos

A simulação de uma seção sísmica de afastamento-nulo (ZO) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhorar a relação sinal/ruído. Segundo Berkovitch et al. (1999) o método Multifoco está baseado na Teoria do Imageamento Homeomórfico e consiste em empilhar dados de cobertura múltipla com distribuição fonte-receptor arbitrária de acordo com uma nova correção de sobretempo, chamada Multifoco. Esta correção de sobretempo esta baseada numa aproximação esférica local da frente de onda focalizante na vizinhança da superfície da terra. Este método permite construir uma seção sísmica no domínio do tempo de afastamento nulo aumentando a relação sinal/ruído. A técnica Multifoco não necessita do conhecimento a priori de um macro-modelo de velocidades. Três parâmetros são usados para descrever a aproximação de tempo de trânsito, Multifoco, os quais são: 1) o ângulo de emergência do raio de afastamento nulo ou raio de reflexão normal (β₀), 2) a curvatura da frente de onda no Ponto de Incidência Normal (R_NIP) e 3) curvatura da frente de Onda Normal (R_N). Sendo também necessário a velocidade próximo a superfície da terra. Neste trabalho de tese aplico esta técnica de empilhamento Multifoco para dados de cobertura múltipla referidos a modelos de velocidade constante e modelo heterogêneos, com o objetivo de simular seções sísmicas afastamento-nulo. Neste caso, como se trata da solução de um problema direto, o macro-modelo de velocidades é considerado conhecido a priori. No contexto do problema inverso tem-se que os parâmetros R_NIP, R_N e β₀ podem ser determinados a partir da análise de coerência aplicada aos dados sísmicos de múltipla cobertura. Na solução deste problema a função objetivo, a ser otimizada, é definida pelo cálculo da máxima coerência existente entre os dados na superfície de empilhamento sísmico. Neste trabalho de tese nos discutimos a sensibilidade da aproximação do tempo de trânsito usado no empilhamento Multifoco, como uma função dos parâmetros R_NIP, R_N e β₀. Esta análise de sensibilidade é feita de três diferentes modos: 1) a primeira derivada da função objetivo, 2) a medida de coerência, denominada semelhança, e 3) a sensibilidade no Empilhamento Multifoco.

Migração com amplitude verdadeira em meios bidimensionais (2-D) e introdução ao caso 2,5-D

Nos últimos anos tem-se verificado através de várias publicações um interesse crescente em métodos de migração com amplitude verdadeira, com o objetivo de obter mais informações sobre as propriedades de refletividade da subsuperfície da terra. A maior parte desses trabalhos tem tratado deste tema baseando-se na aproximação de Born, como em Bleistein (1987) e Bleistein et al. (1987), ou na aproximação do campo de ondas pela teoria do raio como Hubral et al. (1991), Schleicher et al. (1993) e Martins et al. (1997). Considerando configurações arbitrárias de fontes e receptores, as reflexões primárias compressionais podem ser imageadas em reflexões migradas no domínio do tempo ou profundidade de tal modo que as amplitudes do campo de ondas migrado são uma medida do coeficiente de reflexão dependente do ângulo de incidência. Para realizar esta tarefa, vários algoritmos têm sido propostos nos últimos anos baseados nas aproximações de Kirchhoff e Born. Essas duas abordagens utilizam um operador integral de empilhamento de difrações ponderado que é aplicado aos dados da seção sísmica de entrada. Como resultado obtém-se uma seção migrada onde, em cada ponto refletor, tem-se o pulso da fonte com amplitude proporcional ao coeficiente de reflexão naquele ponto. Baseando-se na aproximação de Kirchhoff e na aproximação da teoria do raio do campo de ondas, neste trabalho é obtida a função peso para modelos bidimensionais (2-D) e dois e meio dimensionais (2,5-D) que é aplicada a dados sintéticos com e sem ruído. O resultado mostra a precisão e estabilidade do método de migração em 2-D e 2,5-D como uma ferramenta para a obtenção de informações importantes da subsuperfície da terra, que é de grande interesse para a análise da variação da amplitude com o afastamento (ângulo). Em suma, este trabalho apresenta expressões para as funções peso 2-D e 2,5-D em função de parâmetros ao longo de cada ramo do raio. São mostrados exemplos da aplicação do algoritmo de migração em profundidade a dados sintéticos 2-D e 2,5-D obtidos por modelamento sísmico através da teoria do raio usando o pacote Seis88 (Cervený e Psencík, 1988) e os resultados confirmaram a remoção do espalhamento geométrico dos dados migrados mesmo na presença de ruído. Testes adicionais foram realizados para a análise do efeito de alongamento do pulso na migração em profundidade (Tygel et al., 1994) e a aplicação do empilhamento múltiplo (Tygel et al., 1993) para a estimativa de atributos dos pontos de reflexão - no caso o ângulo de reflexão e a posição do receptor.

Modelos paramétricos simples e causais para cabos de cobre apropriados para caracterização de enlaces em sistemas G.fast

O G.fast é um novo padrão da União Internacional de Telecomunicações que almeja atingir 1 Gb/s sobre enlaces de cobre curtos, utilizando freqüências de até 212 MHz. Essa nova tecnologia requer modelos paramétricos precisos de cabo para fins de projeto, simulação e testes de avaliação de desempenho. A maioria dos modelos de cabo de cobre foram desenvolvidos focando o espectro VDSL - isto é, freqüências de até 30 MHz - e adotam suposições que são violadas quando a faixa de freqüência é estendida para freqüências G.fast. Esta tese apresenta novos modelos de cabo simples e causais capazes de caracterizar com precisão enlaces de cobre compostos por segmentos individuais ou múltiplos, tanto no domínio do tempo quanto no domínio da freqüência. Resultados utilizando as topologias de referência do padrão G.fast mostram que, além de serem precisos, os novos modelos são atrativos devido ao baixo custo computacional e fórmulas fechadas para ajuste de seus parâmetros junto a dados medidos.

Projeto de antenas multibandas integradas para terminais 4G

O presente trabalho apresenta um projeto e análise de uma antena monopolo planar multibanda operante nas frequências de 2000 MHz, 2500 MHz e 3500 MHz para aplicação em dispositivos móveis 3G e 4G. Dentre essas faixas estão contidos os espectros 3G (1895- 2165 MHz), 4G (2500-2690 MHz) e 4G (3400-3600 MHz) utilizadas pelo LTE e WIMAX 3,5. Para o desenvolvimento desse estudo foi tomado como referência uma antena dual-band de Hua- Ming Chen projetada para operar em sistemas ISM 1800 MHz, redes Wireless e Bluetooth 2400 MHz. O projeto proposto neste trabalho é uma modificação na estrutura da antena dual-band adaptada para tripla ressonância simultânea em frequências utilizadas em redes móveis 4G. Para alcançar o objetivo proposto utilizou-se para simulação o software CST® MICROWAVE STUDIO baseado no método Integração Finita no Domínio, as simulações foram feitas no domínio do tempo. Após obtenção satisfatória dos resultados simulados, um protótipo da antena foi confeccionado para medição e teste, os resultados apresentados são coerentes aos obtidos na simulação, o que valida o método proposto.

Detecção de potenciais corticais antecipatórios em sinais de eletroencefalografia (EEG) durante a condução de carros

O reconhecimento da intenção do condutor a partir de sinais de eletroencefalografia (EEG) pode ser útil no desenvolvimento de interfaces cérebro computador (BCI) para serem usadas em sinergia com veículos inteligentes. Isso pode ser benéfico para melhorar a qualidade de interação entre o motorista e o carro, por exemplo, fornecendo uma resposta do carro inteligente alinhada com a intenção do motorista. Neste estudo, considera-se a antecipação como sendo o estado cognitivo que leva a ações especificas durante a condução de um automóvel. Portanto, propomos investigar a presença de padrões antecipatórios em sinais EEG durante a condução de veículos para determinar duas ações especifícas (1) virar à esquerda e (2) virar à direita, alguns milissegundos antes que tais ações aconteçam. Um protocolo experimental foi proposto para gravar sinais EEG de 5 indivíduos enquanto eles operam um simulador de realidade virtual não invasiva - que foi projetado para tal experimento - que simula a condução de um carro virtual. O protocolo experimental é uma variante do paradigma da variação negativa contingente (CNV) com condições Go e No-go no sistema de condução de realidade virtual. Os resultados apresentados neste estudo indicam a presença de padrões antecipatórios em potenciais corticais lentos observados no domínio do tempo (medias dos sinais EEG) e da frequência (Power Spectra e coerência de fase). Isso abre um leque de possibilidades no desenvolvimento de sistemas BCI - baseados em sinais antecipatórios - que conectem o motorista ao veiculo inteligente favorecendo uma tomada de decisão que analise as intenções dos condutores podendo eventualmente evitar acidentes durante a condução.

Atenuação de múltiplas pelo método WHLP-CRS

Nas bacias sedimentares da região Amazônica, a geração e o acúmulo de hidrocarboneto estão relacionados com a presença das soleiras de diabásio. Estas rochas magmáticas intrusivas possuem grandes contrastes de impedância com as rochas sedimentares encaixantes, resultando em múltiplas externas e internas, com amplitudes semelhantes às das reflexões sísmicas primárias. Estas múltiplas podem predominar sobre as informações oriundas de interfaces mais profundas, dificultando o processamento, a interpretação e o imageamento da seção de sísmica. O objetivo da presente tese é realizar a atenuação de múltiplas em seções sintéticas fontecomum (CS), através da combinação dos métodos Wiener-Hopf-Levinson de predição (WHLP) e o do empilhamento superfície-de-reflexão-comum (CRS), aqui denominando pela sigla WHLPCRS. O operador de deconvolução é calculado com as amplitudes reais do sinal sísmico e traço-a-traço, o que consideramos como uma melhor eficiência para a operação de atenuação. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, utilizando o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, e usados para calcular o operador WHLP-CRS. No desenvolvimento do presente trabalho, visamos evitar a inconveniência da seção processada ZO; desenhar e aplicar operadores na configuração CS; e estender o método WHL para camadas curvas.

Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente.

Deconvolução de perfis de poços através do ajuste de energia

As medidas de resistividade são de fundamental importância para o cálculo da saturação de óleo em reservatórios potencialmente produtores. A combinação das medidas de resistividade rasa e profunda permite a obtenção dos parâmetros R_t, R_xo e d_i. Mas, em reservatórios complexos existem dificuldades em se obter leituras confiáveis de R_t, devido à baixa resolução vertical das ferramentas de investigação profunda. Em reservatórios laminados, por exemplo, as leituras obtidas pela ferramenta de indução profunda (ILD) podem levar a uma interpretação errônea das mesmas, levando a acreditar que as medidas obtidas do perfil referem-se a uma única camada. Este problema pode ser em parte resolvido através de uma metodologia que melhore a resolução vertical dos perfis de investigação profunda, valendo-se do uso de informações obtidas de um perfil de alta resolução vertical, i.e; a curva de resistividade rasa. Uma abordagem neste sentido seria usar um perfil de alta resolução que apresente uma boa correlação com o perfil de investigação profunda. Esta correlação pode ser melhor avaliada se aplicarmos um filtro no perfil de alta resolução, de tal forma que o perfil resultante tenha teoricamente a mesma resolução vertical do perfil de baixa resolução. A obtenção deste filtro, porém, recai na premissa de que as funções respostas verticais para as ferramentas de alta e baixa resolução são disponíveis, o que não ocorre na prática. Este trabalho se propõe mostrar uma nova abordagem onde o filtro pode ser obtido a partir de um tratamento no domínio da freqüência. Tal tratamento visa igualar a energia espectral do perfil de alta resolução à energia espectral do perfil de baixa resolução tendo como base o Teorema de Parseval. Será mostrado que a resolução vertical depende fundamentalmente da energia espectral do perfil em questão. A seguir, uma regressão linear é aplicada sobre os perfis de alta resolução filtrado e de baixa resolução. Para cada ponto amostrado dos perfis, uma rotina de minimização é aplicada visando escolher o melhor intervalo de correlação entre os perfis. Finalmente, um fator de correção é aplicado sobre cada ponto do perfil de baixa resolução. Os resultados obtidos com os perfis de indução são promissores, demonstrando a eficácia da abordagem e mesmo quando aplicada em perfis com diferentes propriedades petrofísicas, a metodologia funcionou satisfatoriamente, sem danificar os perfis originais.

Efeito da topografia em sondagens eletromagnéticas bidimensionais nos domínios da frequência e do tempo

O efeito da topografia de um vale, uma colina, um declive e um aclive, em sondagens eletromagnéticas bidimensionais no domínio da freqüência, causa uma variação nos valores da amplitude e da fase da componente H_z devido a um meio homogêneo. A amplitude é menos afetada que a fase. A parede do vale mais próxima da linha de corrente causa uma forte diminuição dos valores da fase, enquanto que a parede do vale mais distante da linha causa um forte aumento. Os efeitos de uma colina são opostos aos do vale. Os efeitos do declive e do aclive num meio homogêneo, são similares, respectivamente, aos observados pelas paredes do vale e da colina mais próximas da linha de corrente. A resposta de um corpo condutivo retangular num meio homogêneo próximo a uma linha de corrente sofre pequenas variações devido a presença de um vale ou de uma colina situada longe da linha de corrente. Porém, se essas feições topográficas estiverem sobre o corpo, elas afetam fortemente a fase e a amplitude da componente H_z e apenas a amplitude, no caso da componente H_x. A resposta transiente não é afetada pela topografia para tempos muito baixos, pois nesse caso se investiga uma finíssima camada da superfície, assim como para tempos muito altos porque a profundidade de investigação é muito grande comparada com a dimensão da topografia. Para os modelos aqui estudados, a maior influência se dá para tempos intermediários, ao redor de 7 ms, ocasionando um retardo do ponto de "cross over" nas curvas de sondagens.