Modelagem de sistemas de detecção de descargas atmosféricas na Amazônia

As descargas atmosféricas representam um dos maiores fatores de riscos para o setor elétrico, especialmente na Amazônia, uma região do mundo identificada com altas densidades de ocorrências e picos de corrente de raios. Nesse contexto, a rede STARNET (Sferics Timing And Ranging NETwork), único sistema de detecção de descargas elétricas atmosféricas, terrestre, gratuito, e contínuo, cobrindo toda a região Amazônica, foi escolhido para a geração de padrões de descargas atmosféricas necessários a otimização dos sistemas de proteção das linhas de transmissão por meio de sistema de monitoramento de descargas atmosféricas. Contudo, o funcionamento intermitente, observado nos diagramas operacionais das várias estações da rede STARNET, afeta o desempenho global do sistema, em especial a eficiência de detecção. Por isso, foi desenvolvido um modelo de uniformização dos dados de descargas atmosféricas em função da configuração da rede (número e localização dos sensores em operação), no objetivo final de editar mapas de índice ceráunico e densidades de raios confiáveis. Existem regiões da Amazônia que sempre apresentam anomalias positivas de densidade de raios como as regiões de Belém e Manaus, bem como o estado do Tocantins que afeitam as linhas de transmissão. Depois, a incidência de raios nas linhas de transmissão na Amazônia foi estimada a partir da distribuição de pico de correntes registrada pela RDR-SIPAM e da densidade de raios corrigida da rede STARNET. A avaliação do desempenho da rede STARNET tendo como referencia à rede de detecção de raios do Sistema de Proteção da Amazônia (RDR-SIPAM) mostrou também uma forte dependência da eficiência de detecção da rede STARNET em função do pico de corrente das descargas atmosféricas.

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

Resistividade e polarização induzida (IP): modelagem analógica

Foi desenvolvida neste trabalho uma nova técnica de modelagem analógica para IP-resistividade em laboratório. Para tanto, foram utilizados modelos geométricos cilíndricos e tabulares. Para a construção destes modelos usou-se cimento, areia quartzítica e grafita. Como a grafita apresenta um intenso fenômeno IP, variou-se o teor de grafita contido nos modelos. Assim, para diferentes formas e medidas dos modelos, experimentaram-se diferentes conteúdos de partículas polarizáveis. As medidas foram feitas no sentido de verificar o comportamento da resposta IP e resistividade dos modelos quanto à forma, conteúdo de grafita, profundidade e orientação do modelo estudado.

Modelagem da disponibilidade de uma unidade geradora da UHE-Tucuruí

Com o advindo do marco regulatório do Setor Elétrico Brasileiro, a partir de 2004, os agentes que atuam neste mercado têm experimentado um acirramento nas disputas por novos negócios, evidenciando um aumento de competitividade. A Disponibilidade dos Ativos Físicos e os Custos com Manutenção se apresentam como os pontos chave para a competitividade dos agentes. O presente trabalho tem por objetivo apresentar uma metodologia de Análise de Disponibilidade de Sistemas Reparáveis, durante as etapas de projeto ou de operação do sistema, contemplando a mensuração dos Custos com Manutenção versus o Desembolso com Aquisição para um nível esperado de desempenho. A metodologia para a Análise de Disponibilidade sugerida se utiliza da construção do Diagrama de Blocos do Sistema com respectivas descrições funcionais, exportação das informações para o formato de Árvore de Sucesso, composta de portas lógicas dos tipos "E" e "OU" as quais caracterizam um subsistema integrante do sistema principal. O analista pode reavaliar a topologia do sistema, agregando ou retirando redundâncias com a finalidade de ajustar o desempenho do projeto aos requisitos de Disponibilidade, Custo de Aquisição e Custos de Manutenção. Como resultados do trabalho foram identificadas lacunas normativas que definem a forma de controle do desempenho dos ativos, estabelecida uma sistemática de integração entre técnicas de modelagem de confiabilidade e disponibilidade, estabelecidos e incorporados indicadores de desempenho de Manutenção Programada em um agente do mercado, foram modelados e discutidos diferentes cenários para um Sistema de Circulação de Óleo de Mancal e foi aplicado o modelo a toda uma Unidade Geradora Hidráulica por meio da implementação computacional do modelo aos componentes críticos dos principais sistemas.

Análise numérica de seções reforçadas de concreto armado considerando o estado inicial de tensões

O presente trabalho trata da formulação e da implementação computacional, em MATLAB^®, para a análise numérica de seções reforçadas de concreto armado, submetidas à flexão composta, considerando o estado de tensões anterior ao reforço. A referida análise se dá com a geração de diagramas de interação momento fletor esforço normal por dois métodos, quais sejam: a) varredura dos domínios de deformação da NBR6118/2003; b) determinação dos picos de diagramas momento fletor – esforço normal – curvatura. Em ambos os procedimentos utiliza-se o método numérico do ponto médio na integração do cálculo dos esforços resistentes, e o método de Newton Raphson para a obtenção de raízes é usado na determinação da deformação no eixo de referência da seção, durante a determinação dos diagramas momento fletor -esforço normal - curvatura. Preliminarmente, concluiu-se que a primeira das duas metodologias aplicadas é inviável. Posteriormente, com a confirmação da eficácia da segunda metodologia, foi possível expandir o escopo do trabalho de modo a permitir a análise de seções de formatos quaisquer executadas em várias etapas, considerando o estado de tensões inicial em cada uma das etapas. A implementação computacional referente a este trabalho se baseou no programa para análise numérica de seções SECLAB, desenvolvido pelo professor Remo Magalhães de Souza.

Modelagem direta bidimensional do método magnetotelúrico com o método dos elementos finitos de arestas

Fizemos a modelagem direta 2D do método magnetotelúrico (MT) com o método dos elementos finitos (MEF) de arestas em termos dos campos primários e secundários. Para usarmos modelos de maior complexidade e diminuirmos o custo computacional utilizamos malhas não estruturadas. Nas malhas utilizadas, introduzimos quatro nós em torno de cada estação MT, constituindo um quadrado alinhado nas direções dos eixos cartesianos

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.

Modelagem de eletrorresistividade 2-D a partir do potencial elétrico secundário

Um dos métodos clássicos da geofísica de exploração é o Método de Eletrorresistividade, estabelecido há um século pelos irmãos Schlumberger e desde então amplamente empregado em prospecção mineral, estudos ambientais e hidrogeologia e em pesquisa de fontes geotermais. Conceitualmente o método consiste de injeção de corrente elétrica na subsuperfície e de medida de diferença de potencial elétrico, resultante da interação da corrente com o meio. As localizações dessas fontes e receptores são determinadas pelo arranjo escolhido para o levantamento. Após o processamento, obtém-se pseudo-seções de resistividade aparente que indicam a distribuição de condutividade em subsuperfície. Devido à simplicidade dos fundamentos físicos de sua formulação, o método apresenta fácil implementação computacional quando comparado aos métodos eletromagnéticos de fonte controlada. Na literatura há inúmeros trabalhos de modelagem computacional, onde se calcula a resposta para problemas 2-D e 3-D. Nestes trabalhos, as pseudo-seções são obtidas a partir do cálculo do potencial elétrico total. Neste trabalho, apresentaremos a resposta da modelagem de eletroresistividade 2-D com o arranjo dipolo-dipolo, obtida a partir do potencial elétrico secundário. A solução é calculada através do método de elementos finitos usando malhas não estruturadas. Para efeito de validação, os resultados são comparados com a resposta 2-D obtida a partir dos potencias totais.

Modelagem do mCSEM no domínio do tempo usando transformada discreta de Fourier

A modelagem do mCSEM é feita normalmente no domínio da frequência, desde sua formulação teórica até a análise dos resultados, devido às simplificações nas equações de Maxwell, possibilitadas quando trabalhamos em um regime de baixa frequência. No entanto, a abordagem através do domínio do tempo pode em princípio fornecer informação equivalente sobre a geofísica da subsuperfície aos dados no domínio da frequência. Neste trabalho, modelamos o mCSEM no domínio da frequência em modelos unidimensionais, e usamos a transformada discreta de Fourier para obter os dados no domínio do tempo. Simulamos ambientes geológicos marinhos com e sem uma camada resistiva, que representa um reservatório de hidrocarbonetos. Verificamos que os dados no domínio do tempo apresentam diferenças quando calculados para os modelos com e sem hidrocarbonetos em praticamente todas as configurações de modelo. Calculamos os resultados considerando variações na profundidade do mar, na posição dos receptores e na resistividade da camada de hidrocarbonetos. Observamos a influência da airwave, presente mesmo em profundidades oceânicas com mais de 1000m, e apesar de não ser possível uma simples separação dessa influência nos dados, o domínio do tempo nos permitiu fazer uma análise de seus efeitos sobre o levantamento. Como parte da preparação para a modelagem em ambientes 2D e 3D, fazemos também um estudo sobre o ganho de desempenho pelo uso do paralelismo computacional em nossa tarefa.

Modelagem de perfis de indução no domínio do tempo por diferenças finitas

Apresentamos aqui uma metodologia alternativa para modelagem de ferramentas de indução diretamente no domínio do tempo. Este trabalho consiste na solução da equação de difusão do campo eletromagnético através do método de diferenças finitas. O nosso modelo consiste de um meio estratificado horizontalmente, através do qual simulamos um deslocamento da ferramenta na direção perpendicular às interfaces. A fonte consiste de uma bobina excitada por uma função degrau de corrente e o registro do campo induzido no meio é feito através de uma bobina receptora localizada acima da bobina transmissora. Na solução da equação de difusão determinamos o campo primário e o campo secundário separadamente. O campo primário é obtido analiticamente e o campo secundário é determinado utilizando-se o método de Direção Alternada Implícita, resultando num sistema tri-diagonal que é resolvido através do método recursivo proposto por Claerbout. Finalmente, determina-se o valor máximo do campo elétrico secundário em cada posição da ferramenta ao longo da formação, obtendo-se assim uma perfilagem no domínio do tempo. Os resultados obtidos mostram que este método é bastante eficiente na determinação do contato entre camadas, inclusive para camadas de pequena espessura.

Modelagem unidimensional de dados do método MCSEM com variação na resistividade da água do mar e na orientação dos sensores

O Marine Controlled Source ElectroMagnetic (MCSEM) é um novo método geofísico que tem sido empregado para a detecção de camadas resistivas no sub-solo oceânico e tem se mostrado uma ferramenta promissora para a indústria do petróleo. Neste trabalho, é modelado levantamentos do método MCSEM incluindo variações nos valores de resistividade da água do mar. Tais variações podem resultar da influência de correntes marinhas, gradientes de temperatura ou de qualquer fator que possa afetar a salinidade da água. Observamos que os dados do método são afetados de maneira significativa por tais variações. Tanto a variação na resistividade quanto na espessura da camada de água pode afetar os resultados do MCSEM de maneira semelhante. É modelado também o posicionamento aleatório dos receptores, mostrando um decaimento significativo da amplitude do sinal quando os receptores se posicionam a uma certa distância da linha de medidas.

Modelagem do tempo de execução de obras civis: estudo de caso na Universidade Federal do Pará

Este artigo propõe a adoção de um modelo numérico prognóstico para estimar a variável "tempo de execução" para empreendimentos públicos de forma objetiva. O trabalho de campo consistiu na aplicação de métodos estatísticos para analisar dados de obras licitadas e executadas durante o período de 2006 a 2009 na Universidade Federal do Pará (UFPA). A análise de dados envolveu cálculos de regressões lineares e transformadas das funções. Após estratificação e tratamento inicial dos dados, os elementos adotados para construção do modelo final se restringiram a 102 obras de um total de 225 originariamente pesquisadas, resultando nos seguintes parâmetros estatísticos: coeficiente de correlação (R) de 0,899; coeficiente de determinação (R²) de 0,808; coeficiente de determinação ajustado (R² ajustado) de 0,796 e erro padrão (Se) de 0,41. Estes parâmetros indicam forte correlação linear entre as variáveis, indicando que 79,60% da variação do tempo para executar uma obra pública podem ser causadas pela variação, em conjunto, das variáveis área construída, custo orçado, capacidade técnica operacional do contratante, capacidade operacional da empresa, tipologia de serviço e estação do ano.

Observações hidrológicas e resultados de modelagem no espalhamento sazonal e espacial da pluma de água Amazônica

Neste estudo são apresentadas as análises espaciais e temporais dos dados sinóticos coletados durante quatro campanhas oceanográficas do programa “REVIZEE/SCORE-N” realizadas em março-abril de 1995, outubro-novembro de 1997, abril-maio de 1999 e agosto-setembro de 2001. Em complemento às análises dos dados hidrológicos, foram utilizados resultados do modelo numérico CLIPPER - Experimento ATL6 (1/6 deg.). A ênfase deste artigo é avaliar o deslocamento sazonal e espacial da pluma de água Amazônica nas regiões localizadas à sudoeste, noroeste e em frente à foz do rio Amazonas. Dados hidrológicos e resultados do modelo demonstram que durante os meses de março-maio a pluma do rio Amazonas move-se para sudoeste, ao longo da costa norte brasileiro, tendendo a diminuir o deslocamento da pluma para sudoeste durante os meses de agosto-setembro. De outubro a novembro, não foi observado ou registrado o deslocamento da pluma para sudoeste. Resultados do modelo revelam anéis gerados pela retroflexão da Corrente Norte do Brasil (CNB) em torno de 7º N, e também a influência de água proveniente do rio Amazonas em volta desses anéis durante o mês de agosto.

Desenvolvimento de algoritmo para modelagem e simulação de sistemas por grafos de ligação

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

Modelagem de dinâmica do uso da terra e cobertura vegetal na região de Santarém, oeste do Pará

A região de Santarém, na última década, apresentou um aumento da área de agricultura de grãos, em especial, arroz, milho e soja. Na base da estrutura fundiária da região, essa dinâmica tem concorrido para a concentração fundiária devido à substituição da pequena propriedade pela grande propriedade capitalizada. Políticas de ordenamento territorial criaram um mosaico de unidades com regras específicas de uso da terra, são unidades de conservação e diferentes modalidades de projetos de assentamentos. Este trabalho tem dois objetivos: estudar o processo de transformação da paisagem após a introdução da agricultura capitalizada de grãos e construir cenários de futuro que analisem alternativas para conter o desflorestamento e a concentração fundiária em curso. O estudo foi realizado utilizando técnicas de sensoriamento remoto e geoprocessamento, a partir de imagens Landsat 5 TM dos anos de 1999, 2004 e 2007. Técnicas de modelagem dinâmica foram empregadas para explorar cenários de futuro (2015) considerando regras de uso do território. Os resultados obtidos mostram que, até 2004, a maior parte da agricultura mecanizada foi implantada em áreas onde anteriormente eram ocupadas pela agropecuária familiar, pastagem e capoeira. Após 2004, a sua expansão se deu, principalmente, sobre áreas de floresta, em especial, dentro de projetos de assentamento. A análise das transições de uso, em diferentes modalidades de assentamento, demonstra que as regras de uso do território estabelecidas por medidas de ordenamento territorial não têm sido seguidas em muitos casos. Este trabalho apresentou como principal contribuição metodológica a incorporação de questões institucionais relacionadas à estrutura fundiária na análise de transformação da paisagem e construção de cenários. Os resultados mostram que tal abordagem é essencial para entender os processos de transformação correntes à região em questão.