Investigação de estratégias de controle digital para regulação de velocidade e emulação da dinâmica de turbinas hidráulicas, com implementação e testes experimentais em uma micromáquina de 10 KVA

Devido não ser comum o estudo de sistemas de potência em plantas reais como usinas hidrelétricas por causa dos riscos e custos que envolvem sua operação, dá-se preferência pela realização computacional de modelos matemáticos desse tipo de planta na resolução de problemas, desenvolvimento de novas tecnologias e formação de recursos humanos. No entanto, modelos realizados computacionalmente não proveem a experiência, visual, auditiva e tátil que um modelo físico real pode oferecer. Portanto, neste trabalho, apresenta-se a descrição e a modelagem de um sistema de geração em escala reduzida de 10kVA, que é um sistema físico real, composto por um motor CC, um gerador síncrono e transformadores, chamado também de sistema micromáquina, o qual faz parte da infraestrutura do Laboratório de Engenharia Elétrica da UFPA. Para este sistema, por intermédio deste trabalho de mestrado e do trabalho de mestrado de Moraes (2011), foram desenvolvido subsistemas eletrônicos e computacionais de acionamento, automação e controle para operá-lo de forma segura resultando em uma excelente plataforma didática para dar suporte às pesquisas em dinâmica e controle de sistemas de potência, bem como o desenvolvimento de trabalhos acadêmicos e de ensino. Nesse mesmo contexto, é apresentada uma proposta de técnica de emulação de turbina hidráulica, tendo como base o controle de potência aplicado ao motor CC do sistema micromáquina. Tal técnica foi desenvolvida principalmente com o propósito de dar suporte ao estudo e desenvolvimento de técnicas de regulação de velocidade de turbinas hidráulicas. Consequentemente, também é apresenta uma proposta de um regulador de velocidade digital para turbinas hidráulicas baseado na estrutura canônica RST de controle digital, cujos parâmetros são projetados por duas técnicas de projeto estudadas neste trabalho: o método de alocação polinomial de polos e o projeto de compensadores por atraso de fase pelo método de resposta em frequência para sistemas discretos. Logo para comprovar a eficácia das ferramentas de hardware, software e teóricas desenvolvidas neste trabalho, resultados de experimentos realizados no sistema micromáquina são apresentados e analisados.

Em busca de carvão vegetal barato: o deslocamento de siderúrgicas para a Amazônia

O artigo analisa o processo de instalação, na Amazônia Oriental brasileira, de diversos produtores independentes de ferro-gusa, nas últimas duas décadas. Mesmo em número e velocidade menores do que indicava o planejamento estatal, foram implantadas no Corredor da Estrada de Ferro Carajás doze siderúrgicas. São empresas dedicadas somente à produção desta commoditie cuja produção requer elevada quantidade de energia suprida por carvão vegetal. Neste artigo, mostramos que houve deslocamento para a Amazônia de indústrias que até os anos 90 estavam concentradas exclusivamente no Sudeste brasileiro e que a demanda de carvão vegetal consolidou-se como principal elo de articulação da siderurgia com a socioeconomia da região de Carajás. A produção de carvão vegetal tem acarretado muitos impactos sociais e ambientais na região Amazônica. Impactos que se materializam pela ampliação da pressão exercida sobre a floresta amazônica, por práticas ambientalmente imprudentes e pela produção do carvão vegetal sustentada por trabalho precário, mal remunerado e insalubre. Exemplo de um cenário diante do qual se desenha como alternativa é a manutenção da fabricação de ferro-gusa nos moldes atuais e a incorporação de pequenos fornos elétricos a arco para a produção de aço. Contudo, a possível instalação de mini-aciarias na região, por si só, não resolve o principal problema socioambiental que envolve o beneficiamento do minério de ferro na região, vinculado aos efeitos deletérios da produção carvoeira, podendo inclusive agravá-los.

Estudo da transformação da austenita retida em martensita induzida por deformação plástica em aços multifásicos

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

Migração (2,5-D) com amplitudes verdadeiras em meios com gradiente constante de velocidade

A migração com amplitudes verdadeiras de dados de reflexão sísmica, em profundidade ou em tempo, possibilita que seja obtida uma medida dos coeficientes de reflexão dos chamados eventos de reflexão primária. Estes eventos são constituídos, por exemplo, pelas reflexões de ondas longitudinais P-P em refletores de curvaturas arbitrárias e suaves. Um dos métodos mais conhecido é o chamado migração de Kirchhoff, através do qual a imagem sísmica é produzida pela integração do campo de ondas sísmicas, utilizando-se superfícies de difrações, denominadas de Superfícies de Huygens. A fim de se obter uma estimativa dos coeficientes de reflexão durante a migração, isto é a correção do efeito do espalhamento geométrico, utiliza-se uma função peso no operador integral de migração. A obtenção desta função peso é feita pela solução assintótica da integral em pontos estacionários. Tanto no cálculo dos tempos de trânsito como na determinação da função peso, necessita-se do traçamento de raios, o que torna a migração em situações de forte heterogeneidade da propriedade física um processo com alto custo computacional. Neste trabalho é apresentado um algoritmo de migração em profundidade com amplitudes verdadeiras, para o caso em que se tem uma fonte sísmica pontual, sendo o modelo de velocidades em subsuperfície representado por uma função que varia em duas dimensões, e constante na terceira dimensão. Esta situação, conhecida como modelo dois-e-meio dimensional (2,5-D), possui características típicas de muitas situações de interesse na exploração do petróleo, como é o caso da aquisição de dados sísmicos 2-D com receptores ao longo de uma linha sísmica e fonte sísmica 3-D. Em particular, é dada ênfase ao caso em que a velocidade de propagação da onda sísmica varia linearmente com a profundidade. Outro tópico de grande importância abordado nesse trabalho diz respeito ao método de inversão sísmica denominado empilhamento duplo de difrações. Através do quociente de dois empilhamentos com pesos apropriados, pode-se determinar propriedades físicas e parâmetros geométricos relacionados com a trajetória do raio refletido, os quais podem ser utilizados a posteriori no processamento dos dados sísmicos, visando por exemplo, a análise de amplitudes.

Análise computacional dos campos de velocidade e temperatura do ambiente interno da usina termelétrica Santana - Amapá

A simulação numérica do escoamento de ar em ambientes internos é na atualidade o método mais apropriado para análise de conforto térmico em ambientes internos. O escoamento de ar nesses ambientes configura-se como um escoamento complexo, pois, em regra geral, é uma combinação de escoamentos cisalhantes livres (jatos) e cisalhantes de parede, além disso, esses escoamentos são governados por forças de inércia e forças de empuxo, caracterizando-o como de convecção mista. A combinação desses mecanismos cria um escoamento com características complexas, como zonas de recirculação, vórtices, descolamento e recolamento de camada-limite dentre outras. Portanto, a precisão da solução estará diretamente ligada, principalmente, na habilidade do modelo de turbulência adotado de reproduzir as características turbulentas do escoamento de ar e da transferência térmica. O objetivo principal do presente trabalho foi a simulação computacional do ambiente térmico interno do galpão que abriga os geradores e motores Wärtzilä da Usina Termelétrica Santana no estado do Amapá. A formulação matemática baseada na solução das equações gerais de conservação inclui uma análise dos principais modelos de turbulência aplicados ao escoamento de ar em ambientes internos, assim como os processos de transferência de calor associados. Na modelagem numérica o método de volumes finitos é usado na discretização das equações de conservação, através do código comercial Fluent-Airpak, que foi usado nas simulações computacionais para a análise dos campos de velocidade e temperatura do ar. A utilização correta do programa computacional foi testada e validada para o problema através da simulação precisa de casos retirados da literatura. Os resultados numéricos foram comparados a dados obtidos de medições experimentais realizados no galpão e apresentou boa concordância, considerando a complexidade do problema simulado, o objetivo da simulação em face da diminuição da temperatura no interior do galpão e, também, em função das limitações encontradas quando da tomada das medições experimentais. Além disso, foram feitas simulações de estratégias de melhoria do ambiente térmico da Usina, baseadas na realidade levantada e nos resultados da simulação numérica. Finalmente, foram realizadas simulações do protótipo de solução proposto para a diminuição da temperatura interna do galpão o que possibilitará um aumento, na faixa de 20 a 30%, do tempo de permanência no interior do galpão.

Análise experimental sobre velocidade de captura em sistemas de transporte pneumático

A velocidade de transporte é um dos mais importantes parâmetros no transporte pneumático de sólidos. O êxito no transporte de materiais particulados depende da predição ou da determinação da velocidade mínima de transporte. Uma velocidade acima daquela necessária ao transporte estável das partículas sólidas conduz a um grande consumo de energia devido ao aumento na perda de pressão do sistema, degradação dos sólidos e abrasão dos sólidos e da tubulação. Já uma velocidade abaixo desse valor limite certamente resultará na deposição das partículas sólidas para o fundo da tubulação e conseqüentemente o entupimento desta. Neste trabalho, uma técnica para medir a velocidade mínima de captura de partículas sólidas em uma tubulação na direção horizontal em um sistema de Transporte Pneumático é desenvolvida. Ela baseia-se em observações visuais do comportamento das partículas ao ocorrer a captura, em medidas da velocidade de operação do gás e da massa das partículas capturadas. É realizada ainda a análise qualitativa de alguns parâmetros que influenciam na velocidade de captura das partículas, permitindo uma maior compreensão dos fenômenos envolvidos.

Análise da velocidade crítica de escoamento e estudo reológico de polpas de caulim da bacia do rio capim em dutos circulares

A indústria mineral representa um grande estimulador da economia do estado, e como destaque nessa indústria podemos mencionar o setor do beneficiamento de Caulim. As soluções dos problemas relacionados a essa atividade mostram-se como uma oportunidade rara de interação universidade-indústria, além de ser um campo fértil para o desenvolvimento científico. Um ponto importante a ser estudado trata do escoamento da suspensão de caulim em dutos e equipamentos relacionados; problemas como incrustação e calculo da perda de carga não foram ainda resolvidos de forma definitiva pela indústria. Nesse contexto o presente trabalho se dispõe a desenvolver uma metodologia para o cálculo dos parâmetros de escoamento principalmente da velocidade mínima necessária para que as partículas em suspensão não depositem no interior dos dutos. Assim o objetivo macro desta dissertação reside em analisar através da implementação de metodologia, o cálculo da velocidade de deposição durante o escoamento de polpas de caulim com diferentes frações mássicas (% em peso de sólidos) no interior de dutos circulares; a partir do estudo reológico de amostras obtidas no seu processo de beneficiamento; propiciando assim, subsídio teórico aos técnicos e Engenheiros que desejam dimensionar um sistema de bombeamento aplicado a polpas de caulim; minimizando seus problemas de sub e superdimensionamento, que acarreta a deposição do material, causando danos à produção e meio ambiente, visto que se fará necessário à abertura da tubulação para remoção dos depósitos oriundos da baixa velocidade.

Avaliação das propriedades mecânicas e elétricas da liga 6201 não refinada, modificada com os teores de cobre (0,05 e 0,3)% para transporte de energia elétrica

No presente trabalho, foram feitos experimentos com a liga de Al – 0,6%Mg - 0,8%Si não refinada, com variações no teor de cobre (0,05% e 0,3%), que objetivaram avaliar dois aspectos do comportamento da liga. O primeiro aspecto diz respeito a afinidade entre esta e o molde, no qual busca se observar o efeito das variáveis térmicas de solidificação sobre a estrutura do fundido e sua correlação com propriedades mecânicas e elétricas, o segundo aspecto está relacionado com o comportamento intrínseco da variação da composição química quando submetida a baixas velocidades e taxas de resfriamento. Para a avaliação das variáveis térmicas de solidificação (velocidades das isotermas liquidus, taxas de resfriamento) utilizouse um dispositivo de solidificação unidirecional horizontal. Os lingotes produzidos a partir da solidificação da liga passaram pelos processos mecânicos de usinagem, laminação e trefilação, examinados em diferentes posições e diâmetros, através do ensaio de tração e condutividade elétrica. Foram analisadas as fraturas em função das microcavidades e diferentes concentrações de teores de cobre constatando-se que em relação ao material deformado a frio, o LRT tem uma tendência de crescimento para a liga de maior concentração de cobre. Os valores avaliados também possibilitaram concluir que a condutividade elétrica cresce quanto maior for o diâmetro das microcavidades e menor o teor de cobre. Com o objetivo de se avaliar somente o comportamento intrínseco da variação da composição química da liga, foi utilizado um segundo dispositivo de solidificação, molde em “U”, pintado internamente com solução de caulim que atribui ao molde baixas velocidades de resfriamento. Os resultados mostram que para os dois teores de cobre há uma tendência de crescimento do LRT com o aumento da redução sofrida pelas amostras, sendo o LRT maior para a liga com maior teor de cobre. Atribuímos este comportamento ao maior teor de soluto que por sua vez apresentou melhor resposta a deformação plástica, encruando mais o material. Entretanto, a caracterização elétrica ocorre de modo inverso, a exemplo do constatado na solidificação unidirecional, onde o menor teor de cobre foi mais eficiente.

Caracterização da liga 6101 refinada com a adição de diferentes teores de cobre e solidificadas em molde unidirecional horizontal e em molde "U"

A liga 6101 (série 6xxx) foi modificada através da adição de cobre nos percentuais de 0,05%Cu e 0,3%Cu, e também de refinadores de grão à base de uma pré-liga Al-Ti-B. As ligas foram caracterizadas termicamente (coeficiente de transferência de calor metal/molde, Velocidade de deslocamento da isoterma liquidus e taxa de resfriamento), mecanicamente (limite de resistência à tração e módulo de tenacidade) e eletricamente (resistência, resistividade e condutividade elétrica) utilizando duas formas de vazamentos: uma sendo através de solidificação em molde unidirecional horizontal e a outra através do molde em “U”. Para a solidificação em molde unidirecional horizontal as caracterizações mostraram-se bastante influenciadas pelos refinadores de grão e pela viscosidade que o cobre exerce na eficiência do contato metal/molde das ligas, alterando significativamente suas propriedades térmicas, mecânicas e elétricas. Por outro lado, na solidificação em molde “U” as ligas foram avaliadas apenas através de caracterização mecânica e elétrica, tendo em vista a preocupação sobre as variáveis do processo de conformação da liga metálica, tais como, a avaliação da taxa de deformação, do limite de resistência a tração, da tenacidade e da resistência, resistividade e condutividade elétrica, do perfil metálico produzido, para corpos de prova de diferentes diâmetros. Para o molde em “U” houve ainda a adição e avaliação de mais duas ligas sendo a Al-0,6Mg-0,4Si-0,1Cu e Al-0,6Mg-0,4Si-0,2Cu. As propriedades mecânicas e elétricas sofreram forte influência do encruamento, do teor de cobre e das microcavidades existentes nas estruturas das ligas.

Estudo da viabilidade operacional do processo de soldagem MAG com alimentação adicional de arame frio

Este trabalho apresenta um estudo experimental da viabilidade operacional, econômica, geométrica e da qualidade da solda com o uso do processo intitulado MAG-CW (Metal Active Gas–Cold Wire); entendido como: processo de soldagem MAG com adição de arame frio ou não energizado. Neste contexto, para o desenvolvimento deste trabalho utilizou-se os processos de soldagem MAG e o MAG-CW, sendo utilizados arames consumíveis sólidos de 1,2 mm de diâmetro, da classe AWS E-70S6 e o dióxido de carbono puro, CO², como gás de proteção com vazão de 17 l/min, para os dois processos empregados. A soldagem automatizada foi realizada com o auxílio de uma fonte eletrônica ajustada em CC⁺, no modo tensão constante e interfaceada a um microcomputador. Para ambos os processos, os cordões de solda foram depositados sobre corpos de prova de aço carbono SAE 1020 medindo 76x155x12,7 mm, com passes em simples deposição, na posição plana puxando e empurrando o arco voltaico. Para o processo MAG-CW adicionouse o arame frio, na região do arco elétrico, por um sistema de alimentação auxiliar de arame de forma seqüencial Tandem e Twin, também, na posição plana. As variáveis de influência foram a Velocidade de Alimentação do Arame-Eletrodo, VAE, que variou em três níveis: 6, 7,5 e 9 m/min, respectivamente, a Velocidade de Alimentação de Arame Frio, VAF, que também variou em três níveis, sendo: 3, 4 e 5 m/min combinadas com VAE. As variáveis de resposta para a avaliação experimental foram a estabilidade do arco elétrico, as características econômicas, as características geométricas e a qualidade superficial do cordão de solda. Os resultados mostraram satisfatória estabilidade do arco, bem como, elevada produção em relação ao processo MAG e bom aspecto superficial das soldas. Neste sentido, o processo MAG-CW, quando comparado ao processo MAG, foi superior na taxa de fusão (TF) em 57% para a VAE igual a 6 m/min, 59 % para VAE igual a 7,5 m/min e de 53 % para a VAE igual a 9 m/min. Observou-se, ainda, que o mesmo ocorreu para a taxa de deposição (TD) aumentando em todos os níveis de velocidade: 62% para a VAE na condição igual a 6 m/min e 7,5 m/min e de 57 % para a VAE na condição igual a 9 m/min. A boa performance no aumento das taxas obtidas pelo processo MAG-CW demonstra vantagem, com maior evidência no desempenho destas características, em função da adição de arame frio ao arco.

Análise do efeito da discretização do modelo de velocidades nas migrações Kirchhoff e Kirchhoff-Gaussian- Beam 2D pré-empilhamento em 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 no imageamento sísmicos. Nesta dissertação, 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).Afim de comparar os métodos Kirchhoff e KGB com respeito à sensibilidade em relação ao comprimento da discretização, aplicamos no conjunto de dados conhecido como Marmousi 2-D quatro grids de velocidade, ou seja, 60m, 80m 100m e 150m. Como resultado, temos que ambos os métodos apresentam uma imagem muito melhor para o menor intervalo de discretização da malha de velocidade. O espectro de amplitude das seções migradas nos fornece o conteúdo de frequência espacial das seções das imagens obtidas.

Análise de velocidade por otimização do semblance na reflexão sísmica

Este trabalho teve como objetivo geral desenvolver uma metodologia sistemática para a inversão de dados de reflexão sísmica em arranjo ponto-médio-comum (PMC), partindo do caso 1D de variação vertical de velocidade e espessura que permite a obtenção de modelos de velocidades intervalares, v_int,n, as espessuras intervalares, z_n, e as velocidades média-quadrática, v_RMS,n, em seções PMC individualizadas. Uma consequência disso é a transformação direta destes valores do tempo para profundidade. Como contribuição a análise de velocidade, foram desenvolvidos dois métodos para atacar o problema baseado na estimativa de velocidade intervalar. O primeiro método foi baseado na marcação manual em seções PMC, e inversão por ajuste de curvas no sentido dos quadrados-mínimos. O segundo método foi baseado na otimização da função semblance para se obter uma marcação automática. A metodologia combinou dois tipos de otimização: um Método Global (Método Price ou Simplex), e um Método Local (Gradiente de Segunda Ordem ou Conjugado), submetidos a informação à priori e vínculos. A marcação de eventos na seção tempo-distância faz parte dos processos de inversão, e os pontos marcados constituem os dados de entrada juntamente com as informações à priori do modelo a ser ajustado. A marcação deve, por princípio, evitar eventos que representem múltiplas, difrações e interseções, e numa seção pode ser feita mais de 50 marcações de eventos, enquanto que num mapa semblance não se consegue marcar mais de 10 eventos de reflexão. A aplicação deste trabalho é voltada a dados sísmicos de bacias sedimentares em ambientes marinhos para se obter uma distribuição de velocidades para a subsuperfície, onde o modelo plano-horizontal é aplicado em seções PMC individualizadas, e cuja solução pode ser usada como um modelo inicial em processos posteriores. Os dados reais da Bacia Marinha usados neste trabalho foram levantados pela PETROBRAS em 1985, e a linha sísmica selecionada foi a de número L5519 da Bacia do Camamu, e o PMC apresentado é a de número 237. A linha é composta de 1098 pontos de tiro, com arranjo unilateraldireito. O intervalo de amostragem é 4 ms. O espaçamento entre os geofones é 13,34 m com o primeiro geofone localizado a 300 m da fonte. O espaçamento entre as fontes é de 26,68 m. Como conclusão geral, o método de estimativa de velocidade intervalar apresentada neste trabalho fica como suporte alternativo ao processo de análise de velocidades, onde se faz necessário um controle sobre a sequência de inversão dos PMCs ao longo da linha sísmica para que a solução possa ser usada como modelo inicial ao imageamento, e posterior inversão tomográfica. Como etapas futuras, podemos propor trabalhos voltados direto e especificamente a análise de velocidade sísmica estendendo o caso 2D de otimização do semblance ao caso 3D, estender o presente estudo para o caso baseado na teoria do raio imagem com a finalidade de produzir um mapa continuo de velocidades para toda a seção sísmica de forma automática.

Avaliação de algoritmos para conversão de modelos de velocidade de tempo para profundidade

Ainda hoje, a migração em tempo é o processo de imageamento substancialmente empregado na indústria do petróleo. Tal popularidade é devida ao seu alto grau de eficiência e robustez, além de sua habilidade em focalizar refletores nos mais variados ambientes geológicos. Entretanto, em áreas de alta complexidade geológica a migração em tempo falha de tal forma que a migração em profundidade e um campo de velocidade em profundidade são indispensáveis. Esse campo é geralmente obtido através de processos tomográficos partindo de um campo de velocidade inicial. A conversão de campos de velocidade de tempo para profundidade é uma forma rápida de se obter um campo inicial mais consistente geologicamente para tais processos. Alguns algoritmos de conversão tempo-profundidade recentemente desenvolvidos baseados no traçamento de raios-imagem são revistos e um algoritmo alternativo baseado na propagação da frente de onda-imagem é proposto. Os algoritmos são aplicados a dados sintéticos bidimensionais e avaliados de acordo com suas eficiência e acurácia, destacando suas vantagens, desvantagens e limitações na obtenção de campos de velocidade em profundidade.

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.

Um estudo da estabilidade mecânica de poços a partir de perfis geofísicos em formações pouco consolidadas

A análise da estabilidade mecânica de um poço pode ser feita a partir do cálculo de parâmetros elásticos da formação utilizando a densidade do meio e as velocidades de propagação das ondas compressional e cisalhante na formação rochosa, os quais podem ser obtidos de perfis geofísicos do poço. Em formações sedimentares pouco consolidadas as ferramentas de perfilagem sônica convencionais (monopolares) não conseguem registrar com acuidade a velocidade da onda cisalhante pois a primeira chegada dessa onda é camuflada pela chegada de outras ondas que podem ser mais rápidas que a onda cisalhante num poço perfurado neste tipo de formação. Medidas das velocidades sônicas são feitas em laboratório em amostras da formação, sob condições semelhantes às condições in situ, servindo como ajuste das velocidades registradas no poço pela ferramenta de perfilagem sônica. Para a análise de estabilidade da formação, perfis auxiliares são necessários como o perfil de porosidade, saturação de fluidos e perfis de composição mineralógica da formação rochosa. Exige-se ainda dados de testes de avaliação da formação e de condições do reservatório, mas que são comuns em poços de petróleo, como o teste de formação e os testes de pressurização do poço, tais como o teste de micro-fraturamento hidráulico ou o teste de absorção. A avaliação das tensões principais efetivas que atuam distante do poço, e que não são afetadas pela sua presença, é feita através da associação de um modelo de deformação elástica apropriado e o resultado do teste de pressurização disponível para o poço em estudo. Utilizando resultados clássicos da teoria da elasticidade geral pode-se calcular o campo de tensões modificado na vizinhança da parede do poço devido ao efeito da própria presença do poço ali perfurado e da diferença de pressão existente entre o interior do poço e a formação rochosa. A determinação das propriedades mecânicas da formação a partir das velocidades sônicas e a avaliação do campo de tensões assumindo um modelo elástico de deformação, supõem o meio rochoso no qual as ondas se propagam como elástico, homogêneo e isotrópico. Esta suposição representa a principal aproximação assumida pela metodologia descrita neste trabalho. De posse das propriedades mecânicas da formação e do campo de tensões que age na vizinhança do poço resta definir o critério segundo o qual a rocha sofre instabilidade mecânica quando submetida aquele campo de tensões. Isto permite determinar se, nas condições avaliadas do poço e da formação, haverá quebra da parede do poço por excesso de tensão e, se houver, qual a sua extensão. Assim o problema é como analisar o comportamento mecânico de um poço em uma formação pouco consolidada a partir de perfis geofísicos os quais podem ter problemas no registro das propriedades físicas do meio em formações deste tipo. A metodologia proposta é aplicada a dois intervalos de profundidade pertencentes a dois poços onde arenitos e folhelhos se intercalam e nos quais todos os dados necessários estão disponíveis. Os resultados obtidos mostram que, exceto quando outros mecanismos de quebra da parede do poço agem na formação, a metodologia proposta consegue com sucesso detectar zonas de ocorrência de instabilidade mecânica do poço provocadas por um campo de tensões que excede a resistência mecânica da formação.