Morfodinâmica costeira e o uso da orla oceânica de Salinópolis (Nordeste do Pará)

Os aspectos morfodinâmicos relacionados à erosão ou acresção da linha de costa são alguns dos assuntos analisados na gestão das zonas costeiras que vêem sendo tratada em todo mundo no sentido de monitorar e proteger essas zonas. Esta tese objetiva analisar o comportamento da morfodinâmica costeira de Salinópolis, relacionando-o ao uso da orla oceânica. A área de estudo foi compartimentada em três setores: Oeste (praias da Corvina e do Maçarico), Central (praia do Farol Velho) e Leste (praia do Atalaia). A metodologia consistiu na: (a) aquisição e tratamento de imagens multitemporais (1988-2001-2013) do satélite Landsat 5 TM, 7 ETM e 8 OLI; (b) aplicação de entrevistas/questionários com banhistas, (c) aquisição de dados de campo durante as estações chuvosa (26, 27, 28/04/2013) e menos chuvosa (04, 05, 06/10/2013); e (d) análise laboratorial para o tratamento dos dados adquiridos em campo (topografia das praias estudadas, amostragem de sedimentos superficiais das mesmas e com o uso de armadilhas, e medições oceanográficas de ondas, marés, correntes e turbidez). Foram feitas as representações gráficas dos perfis topográficos das praias, calculados os parâmetros estatísticos granulométricos de Folk & Ward (1957), as taxas do transporte sedimentar nas praias e os parâmetros morfométricos de Short & Hesp (1982), estes últimos foram calculados com o intuito de relacioná-los aos estados morfodinâmicos de praias propostos por Wright & Short (1984) e Masselink & Short (1993). Para a classificação da costa oceânica de Salinópolis em termos de uso e ocupação foi utilizado o decreto nº 5.300 de 7 de dezembro de 2004. A partir das pesquisas sobre a urbanização na costa e das obras situadas nos ambientes costeiros foi utilizada uma matriz proposta por Farinaccio & Tessler (2010) que lista uma série de impactos ambientais, e o quadro de geoindicadores do comportamento da linha de costa proposto por Bush et al. (1999), para a identificação de locais com vulnerabilidade à erosão ou acresção. Para as condições oceanográficas em cada praia e periculosidade ao banho nas mesmas, foram integralizados os dados de ondas, de correntes, de morfodinâmica praial e questionários aplicados com banhistas. Atualmente, a orla oceânica de Salinópolis possui diferentes características quanto à utilização e conservação, abrangendo desde a tipologia de orlas naturais (Classe A) até orlas com urbanização consolidada (Classe C). A primeira ocorre nos extremos da área de estudo e, a segunda, na região da sede municipal. Quatro tipos de praias foram identificados segundo a exposição marítima e o grau das condições oceanográficas: tipo 1 (Maçarico), tipo 2 (Corvina), tipo 3 (Farol Velho) e tipo 4 (Atalaia). O trecho de costa com maiores impactos ambientais e com elevada erosão costeira localiza-se na praia do Farol Velho. O grau de periculosidade ao banho foi de 4 (praia do Maçarico) a 7 (praia do Atalaia) – médio a alto grau de risco. As praias de Salinópolis apresentam declives suaves (< 1,5°), grandes variações na linha de costa entre as estações do ano (9,6 a 88, 4 m) e volume sedimentar variável dependendo do grau de exposição das praias ao oceano aberto. Predominou o estado morfodinâmico dissipativo (Ω>5,5) para estas praias, mas com ocorrência do estado de banco e calha longitudinais (4,7<Ω<5,5) no setor oeste. As macromarés na área de estudo apresentaram altura máxima de 5,3 m (Setor Central, durante a estação menos chuvosa) e mínima de 4 m no mesmo setor, durante a estação chuvosa. As correntes longitudinais foram mais intensas no setor leste (>0,45 m/s) durante as duas estaçoes do ano. As alturas de ondas foram também maiores no setor leste (máximo de 1,05 m durante a maré enchente na estação menos chuvosa) e os períodos de ondas foram mais curtos (<4,5 s) no setor oeste. A média granulométrica obtida dos sedimentos coletados na face praial apresentou escala mais freqüente entre 2,6 a 2,8 phi, indicando a predominância de areia fina. O grau de seleção predominante dos sedimentos foi de 0,2 a 0,5 phi (muito bem selecionados e bem selecionados), e da assimetria foi de positiva (0,10 a 0,30) e de aproximadamente simétrica (-0,10 a 0,10). O grau de curtose variou desde muito platicúrtica (<0,67) a muito leptocúrtica (1,50 a 3,00). Foram observados eventos de acresção sedimentar da estação chuvosa a menos chuvosa. De 22/07/1988 a 28/08/2013 (25 anos) também houve predomínio de acresção, onde o avanço médio linear da linha de costa foi de 190,26 m. O recuo médio linear obtido para toda área de estudo foi de -42,25 m. Áreas com maior erosão são pontuais: divisas das praias da Corvina e Maçarico, e Farol Velho e Atalaia. Os traps portáteis indicaram uma maior quantidade de sedimentos transportados longitudinalmente na estação menos chuvosa (Mín. 280 g/m³: enchente, setor oeste; Máx. 1098 g/m³: vazante, setor leste). Nos traps de espraiamento, o balanço entre a quantidade de sedimentos entrando e saindo nas praias foi menor no setor central (Mín. 80 g/m³: vazante, estação menos chuvosa; Máx. 690 g/m³: enchente, estação menos chuvosa). A circulação costeira sedimentar é proveniente, principalmente, do efeito das marés, com direção governada pela enchente e vazante dos rios que atravessam a costa. Os dados indicam o transporte longitudinal de sedimentos da ilha de Atalaia e rio Sampaio para o setor oeste e as margens das faixas praiais.

Estudo da técnica de enchimento de chanfro do aço naval pelo processo GMAW-CW

Nesse trabalho estudou-se a viabilidade operacional e as características econômicas e geométrica da técnica do processo de soldagem GMAW-CW (alimentação adicional de um arame frio) em comparação Gas Metal Arc Welding (GMAW). O sistema de alimentação de arame frio foi projetado e adaptado a pistola de soldagem MIG/MAG. Utilizou-se uma fonte eletrônica de múltiplos processos ajustada em tensão constante e CC+, a proteção gasosa foi uma mistura 75%Ar+25%CO² e CO² comercialmente puro. O arame utilizado foi o da classe ER70S-6 com diametro de 1,2 mm para o arame eletrodo e 1,0 mm para o arame frio, os dois arames foram alimentados em cabeçotes independentes. As variáveis operacionais de entrada foram: a velocidade de alimentação de arame energizado, em três níveis, 4, 6 e 8 m/min e a velocidade de alimentação do arame frio em 50%, 60% e 70% da velocidade de alimentação do arame energizado. As soldagens foram automatizadas em simples deposição no sentido empurrando e o posicionamento do arame frio, em um único nível, Tandem em chanfro em “U” de chapas de aço ASTM 1020. As variáveis de resposta utilizadas foram: inspeção superficial dos cord��es; análise da geometria (largura, penetração, reforço e diluição) da solda e econômicas (taxa de fusão, taxa de deposição, rendimento e custo operacional). Os resultados indicaram que para a análise superficial, com o uso do gás Ar25%CO² a superfície dos cord��es mostraram-se mais homogêneas em relação ao CO² e com menor índice de salpicagem, para a análise das características econômicas, o processo GMAW-CW sempre foi superior ao processo convencional, quanto aos custos operacionais o processo convencional mostrou-se menor, porém não houve o preenchimento do chanfro, o que ocorreu com a utilização do processo GMAW-CW.

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.

Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.