Aplicação de deconvolução homomórfica a dados sísmicos

Um registro sísmico é frequentemente representado como a convolução de um pulso-fonte com a resposta do meio ao impulso, relacionada ao caminho da propagação. O processo de separação destes dois componentes da convolução é denominado deconvolução. Existe uma variedade de aproximações para o desenvolvimento de uma deconvolução. Uma das mais comuns é o uso da filtragem linear inversa, ou seja, o processamento do sinal composto, através de um filtro linear, cuja resposta de frequência é a recíproca da transformada de Fourier de um dos componentes do sinal. Obviamente, a fim de usarmos a filtragem inversa, tais componentes devem ser conhecidas ou estimadas. Neste trabalho, tratamos da aplicação a sinais sísmicos, de uma técnica de deconvolução não linear, proposta por Oppenheim (1965), a qual utiliza a teoria de uma classe de sistemas não lineares, que satisfazem um princípio generalizado de superposição, denominados de sistemas homomórficos. Tais sistemas são particularmente úteis na separação de sinais que estão combinados através da operação de convolução. O algoritmo da deconvolução homomórfica transforma o processo de convolução em uma superposição aditiva de seus componentes, com o resultado de que partes simples podem ser separadas mais facilmente. Esta classe de técnicas de filtragem representa uma generalização dos problemas de filtragem linear. O presente método oferece a considerável vantagem de que não é necessário fazer qualquer suposição prévia sobre a natureza do pulso sísmico fonte, ou da resposta do meio ao impulso, não requerendo assim, as considerações usuais de que o pulso seja de fase-mínima e que a distribuição dos impulsos seja aleatória, embora a qualidade dos resultados obtidos pela análise homomórfica seja muito sensível à razão sinal/ruído, como demonstrado.

Avaliação do efeito de janela e descoloração nos filtros Wiener-Hopf

A presente dissertação consta de estudos sobre deconvolução sísmica, onde buscamos otimizar desempenhos na operação de suavização, na resolução da estimativa da distribuição dos coeficientes de reflexão e na recuperação do pulso-fonte. Os filtros estudados são monocanais, e as formulações consideram o sismograma como o resultado de um processo estocástico estacionário, e onde demonstramos os efeitos de janelas e de descoloração. O principio aplicado é o da minimização da variância dos desvios entre o valor obtido e o desejado, resultando no sistema de equações normais Wiener-Hopf cuja solução é o vetor dos coeficientes do filtro para ser aplicado numa convolução. O filtro de deconvolução ao impulso é desenhado considerando a distribuição dos coeficientes de reflexão como uma série branca. O operador comprime bem os eventos sísmicos a impulsos, e o seu inverso é uma boa aproximação do pulso-fonte. O janelamento e a descoloração melhoram o resultado deste filtro. O filtro de deconvolução aos impulsos é desenhado utilizando a distribuição dos coeficientes de reflexão. As propriedades estatísticas da distribuição dos coeficientes de reflexão tem efeito no operador e em seu desempenho. Janela na autocorrelação degrada a saída, e a melhora é obtida quando ela é aplicada no operador deconvolucional. A transformada de Hilbert não segue o princípio dos mínimos-quadrados, e produz bons resultados na recuperação do pulso-fonte sob a premissa de fase-mínima. O inverso do pulso-fonte recuperado comprime bem os eventos sísmicos a impulsos. Quando o traço contém ruído aditivo, os resultados obtidos com auxilio da transformada de Hilbert são melhores do que os obtidos com o filtro de deconvolução ao impulso. O filtro de suavização suprime ruído presente no traço sísmico em função da magnitude do parâmetro de descoloração utilizado. A utilização dos traços suavizados melhora o desempenho da deconvolução ao impulso. A descoloração dupla gera melhores resultados do que a descoloração simples. O filtro casado é obtido através da maximização de uma função sinal/ruído. Os resultados obtidos na estimativa da distribuição dos coeficientes de reflexão com o filtro casado possuem melhor resolução do que o filtro de suavização.

Comparação de desempenho da deconvolução preditiva multicanal e da filtragem f-k na atenuação de múltiplas do fundo do mar

The seismic reflection is used on a large scale in oil exploration. In case of marine acquisition the presence of high impedance contrast at the interfaces water/air generates multiple reflections events. Such multiple events can mask primary events; thus from the interpretational viewpoint it is necessary to mitigate the multiples. In this manuscript we compare two methods of multiple attenuation: the predictive multichannel deconvolution (DPM) and the F-K filtering (FKF). DPM is based in the periodicity of the multiples while FKF is based in multiples and primaries splitting in F-K domain. DPM and FKF were applied in common-offset and CDP gathers, respectively. DPM is quite sensible to the correct identification of the period and size of the filter while FKF is quite sensible to an adequate choice of the velocity in order to split multiples and primaries events in the F-K domain. DPM is a method that is designed to act over a specific event. So, when the parameters are well selected, DPM is very efficient in removing the specified multiple. Then DPM can be optimized by applying it several times, each time with a different parameterization. A deficiency of DPM occurs when a multiple is superposed to a primary event: in this situation, DPM can attenuate also the primary event. On the other hand, FKF presents almost the same performance to all multiples that are localized in the same sector of the F-K domain. The two methods can be combined in order to take advantage of their associated potentials. In this situation, DPM is firstly applied, with a focus in the sea bed multiples. Then FKF is applied in order to attenuate the remaining multiples

Atenuação de múltiplas e compressão do pulso fonte em dados de sísmica de reflexão utilizando o filtro Kalman-Bucy

O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função b_i(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas.

Avaliação da vulnerabilidade sísmica de edifícios escolares

A grande motivação para este trabalho resultou na avaliação da vulnerabilidade sísmica das escolas básicas e secundárias que fazem parte integrante do parque escolar de Portugal Continental. Apesar de até ao momento apenas se terem estudado a vulnerabilidade de algumas escolas em algumas zonas particulares do nosso país, para este trabalho de investigação tivemos uma ambição muito maior e obviamente fizemos um esforço muito maior. Estabelecemos assim como meta, a avaliação de todo o parque escolar. Embora todo o parque escolar possua na ordem das três centenas de escolas em todo o território nacional e sendo este projeto de reabilitação, um projeto com a duração de 2007 a 2015. Em 2011, por motivos da crise económica, todo o projeto congelou sendo reabilitadas até à data cerca de apenas um terço das escolas. Esta quantidade foi o número de escolas que avaliamos em todo o país. As escolas, sendo edifícios públicos com uma importância fundamental, tanto pela elevada concentração de jovens, como pela função essencial como centros de aprendizagem para as gerações vindouras, como também pela ameaça que representam na eventualidade dum cenário sísmico pela enorme densidade de utilizadores, e pela vantagem de nesse cenário de catástrofe a importância estrutural ser superior em relação à maior parte dos edifícios correntes, devidamente demonstrado pelos argumentos enumerados, consequentemente as escolas podem servir como instalações de proteção civil perante uma catástrofe sísmica para apoio das populações circundantes afetadas. Portanto para cada uma das escolas cordialmente fornecidas pelo Parque Escolar, E.P.E., foi feito um estudo exaustivo e muito individual de cada uma das escolas, onde cada análise foi desenvolvida por uma metodologia simplificada, sendo cada análise sempre individual, e nunca tendo sida aplicada em série, este fator melhora substancialmente a eficácia da avaliação para a quantificação das vulnerabilidades e da determinação do grau de dano e das frações de perda para os requisitos fundamentais de limitação de danos, de não colapso, e de colapso iminente que correspondem a ações com períodos de retorno de 95, 475 e 975 anos. Este trabalho é fundamental para as entidades competentes terem a consciência da vulnerabilidade das escolas secundárias, para puderem atuar a nível estrutural e diminuir assim a vulnerabilidade sísmica, e mesmo que por impossibilidade económica o poder governamental não intervenha, então pode e principalmente deve, elaborar planos de emergência tanto com engenheiros civis qualificados como com a total colaboração das corporações de bombeiros que fazem parte das forças de operação e socorro da Autoridade Nacional de Proteção Civil (ANPC).

Avaliação da resposta sísmica de edifícios sitos em solos brandos estratificados

O Eurocódigo 8 exige um estudo mais aprofundado da resposta dinâmica de solos brandos classificados como tipo S1, tendo em vista o dimensionamento de edifícios. Neste trabalho apresenta-se uma metodologia visando a avaliação da segurança, ou o dimensionamento, de edifícios construídos sobre solos brandos estratificados. A metodologia proposta baseia-se nas velocidades das ondas de corte VS, obtidas de estudos geofísicos/geotécnicos, realizados in situ, tal como é proposto na parte 5 do EC8. A amplificação dinâmica das vibrações, em resultado dos contrastes das impedâncias dos diversos estratos, é obtida admitindo a propagação unidimensional das ondas sísmicas VS, na vertical, a partir do substrato rochoso. A resposta superficial é quantificada com base na análise não linear da coluna de solo sobre o substrato rochoso. A acção sísmica é caracterizada com base em vários acelerogramas “strong motion”, registados noutros países, em afloramentos rochosos, multiplicados por um factor de escala. Os factores de escala são obtidos de modo a que o espectro de resposta médio se ajuste ao espectro objectivo do EC8 (solo tipo A). A análise sísmica da estrutura do edifício é realizada com o espectro de resposta médio, que é obtido tendo em conta as funções de transferência determinadas para o terreno em estudo. Aplica-se a metodologia proposta a um caso de estudo no Algarve.

Parâmetros que condicionam os valores da casualidade sísmica

As dispersões nos resultados das análises da casualidade sísmica têm originado alguma prudência na sua utilização em aplicações práticas de minimização de risco sísmico. Este trabalho visa a verificação, de forma sistemática e crítica, da influência de alguns parâmetros da análise, na credibilidade dos resultados, nomeadamente, o tipo de modelo probabilista adoptado, a delimitação das zonas geradoras de eventos sísmicos, a definição das leis de atenuação e a respectiva adequação à sismicidade verificada e às evidências sismotectónicas, no contexto da cidade de Lisboa.

Análise comparativa da casualidade sísmica do Algarve

In this paper we compared the seismic hazard of Algarve obtained with two different models: a point-source model (proposed by McGuire) and a fault-rupture model (proposed by Der Kiureghian and Ang).

Influência da modelação de elementos não estruturais no dimensionamento à acção sísmica

Realizou-se um estudo paramétrico sobre a influência de elementos não estruturais, executados em betão, no dimensionamento de elementos estruturais de betão armado, face às acções sísmicas regulamentares.

Nova metodologia de execução do ensaio SCPTu com avaliação da influência da fonte sísmica quanto à sua dimensão e afastamento

O objectivo deste trabalho consiste na determinação das Ondas de Corte, onda S, através de uma nova metodologia de ensaio do piezocone sísmico, SCPTu. A nova metodologia adoptada de aquisição de dados, consiste na aquisição da onde de corte no movimento ascendente de remoção de varas, após execução do já conhecido ensaio CPTu, ao invés da tradicional aquisição, no movimento descendente, previamente ao acréscimo de nova vara. Além da introdução de uma nova metodologia de aquisição sísmica, tentou-se avaliar a influência da dimensão da fonte geradora da onda de corte e a influência do afastamento das mesmas ao ponto de cravação, do ensaio de penetração estática mencionado. Na avaliação das influências foram usadas três fontes de madeira de secção rectangular, com uma chapa de ferro a topo e diferenciadas entre si pelas dimensões de 1m, 2m e 3m e afastamentos até 4m.

Identificação modal e simulação da resposta sísmica de barragens abóbada

Trabalho Final de Mestrado elaborado no Laboratório Nacional de Engenharia Civil (LNEC) para a obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de cooperação ente o ISEL e o LNEC

Verificação da segurança à ação sísmica da estrutura de um colégio em Ponta Delgada utilizando o Eurocódigo 8 - parte 3

Dissertação para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização em Estruturas

Análise sísmica e proposta de reforço de uma estrutura do séc. XIX sujeita a uma reabilitação parcial

Trabalho final de Mestrado para a obtenção do grau de mestre em Engenharia Civil

Modelação do comportamento dinâmico de sistemas barragem-fundação-albufeira: análise sísmica de uma barragem de abóbadas múltiplas

Trabalho Final de Mestrado elaborado no Laboratório Nacional de Engenharia Civil (LNEC) para a obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de cooperação entre o ISEL e o LNEC

Verificação da segurança sísmica de um edifício existente em betão armado

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil na Área de Especialização de Estruturas