Estágio realizado na Procalçado - Produtora de Componentes para Calçado, S. A. - e orientado pelo Diretor Industrial da Empresa Rui Jorge Santos Russo / Tese de mestrado integrado. Engenharia do Ambiente (Área de Especialização de Projeto). Faculdade de Engenharia. Universidade do Porto. 2012

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / No presente trabalho, um algoritmo algébrico sequencial é utilizado para descrever a propagação de ondas acústicas ao longo de uma barra e utilizado na identificação de danos. Esse algoritmo é validado com base nos conficientes de sensibilidade dos ecos correspondentes aos diferentes cenários de danos apresentados. Na formulação do problema de identificação de dano, o campo de impedância generalizada, que minimiza o funcional definido como a distância entre o eco calculado e o eco experimental sintético é procurado. Os tempos de percurso da resposta, obtidos a partir de experimentos numéricos, são utilizados para identificar a posição, intensidade e forma do dano. Para simular dados corrompidos, diferentes níveis de ruído - variando de 30 a 0 dB - são introduzidos. O processo de identificação foi avaliado com os seguintes métodos de otimização: Otimização por Enxame de Partículas (PSO); Luus-Jaakola (LJ); Algoritmo de Colisão de Partículas (PCA); Algoritmos Genéticos (GA) e Recozimento Simulado (SA); e a hibridização desses métodos com o método determinístico de Levenberg-Marquardt. É mostrado que o processo de identificação de dano construído sobre a abordagem de propagação de ondas acústicas foi bem sucedido, mesmo para dados ruidosos altamente corrompidos. Os resultados dos casos testes são apresentados e algumas observações sobre as vantagens dos métodos determinísticos e estocásticos e sua combinação também são relatados. / In the present work, a sequential algorithm is used for describing the acoustic wave propagation along a bar and applied for damage identification purposes. The algorithm is validated based on the sensitivity coefficients of the corresponding echoes to the adressed damage scenarios. In the formulation of the damage identification problem, the generalized impedance field, that minimizes the functional defined as the distance between the calculated echo and the synthetic experimental one is sought. Time history responses, obtained from pulse-echo experiments, are used to identify damage position, severity and shape. In oder to account for noise corrupted data, different levels of signal to noise ratio - varying from 30 to 0 dB - are introduced. In the identification procedure the following optimization methods were applied: Particle Swarm Optimization (PSO); Luus-Jaakola (LJ); Particle Collision Algorithm (PCA); Genetic Algorithms (GA); and Simmulated Annealing (SA): and the hybridization of these methods with the deterministic Levenberg-Marquardt method. It is shown that the damage identification procedure built on the acoustic wave propagation approach was successful, even for highly corrupted noisy data. Test case results are presented and a few comments on the advantages of deterministic and stochastic methods and their combination are also reported.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / No presente trabalho, um algoritmo algébrico sequencial é utilizado para descrever a propagação de ondas acústicas ao longo de uma barra e utilizado na identificação de danos. Esse algoritmo é validado com base nos conficientes de sensibilidade dos ecos correspondentes aos diferentes cenários de danos apresentados. Na formulação do problema de identificação de dano, o campo de impedância generalizada, que minimiza o funcional definido como a distância entre o eco calculado e o eco experimental sintético é procurado. Os tempos de percurso da resposta, obtidos a partir de experimentos numéricos, são utilizados para identificar a posição, intensidade e forma do dano. Para simular dados corrompidos, diferentes níveis de ruído - variando de 30 a 0 dB - são introduzidos. O processo de identificação foi avaliado com os seguintes métodos de otimização: Otimização por Enxame de Partículas (PSO); Luus-Jaakola (LJ); Algoritmo de Colisão de Partículas (PCA); Algoritmos Genéticos (GA) e Recozimento Simulado (SA); e a hibridização desses métodos com o método determinístico de Levenberg-Marquardt. É mostrado que o processo de identificação de dano construído sobre a abordagem de propagação de ondas acústicas foi bem sucedido, mesmo para dados ruidosos altamente corrompidos. Os resultados dos casos testes são apresentados e algumas observações sobre as vantagens dos métodos determinísticos e estocásticos e sua combinação também são relatados. / In the present work, a sequential algorithm is used for describing the acoustic wave propagation along a bar and applied for damage identification purposes. The algorithm is validated based on the sensitivity coefficients of the corresponding echoes to the adressed damage scenarios. In the formulation of the damage identification problem, the generalized impedance field, that minimizes the functional defined as the distance between the calculated echo and the synthetic experimental one is sought. Time history responses, obtained from pulse-echo experiments, are used to identify damage position, severity and shape. In oder to account for noise corrupted data, different levels of signal to noise ratio - varying from 30 to 0 dB - are introduced. In the identification procedure the following optimization methods were applied: Particle Swarm Optimization (PSO); Luus-Jaakola (LJ); Particle Collision Algorithm (PCA); Genetic Algorithms (GA); and Simmulated Annealing (SA): and the hybridization of these methods with the deterministic Levenberg-Marquardt method. It is shown that the damage identification procedure built on the acoustic wave propagation approach was successful, even for highly corrupted noisy data. Test case results are presented and a few comments on the advantages of deterministic and stochastic methods and their combination are also reported.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.

Formulamos um modelo simplificado para o estudo do processo de injeção de solvente em reservatórios de petróleo, onde o fluido injetado (um ácido) tem a capacidade de dissolver parcialmente a matriz sólida. Como hipóteses principais, consideramos que o solvente e o soluto (componente químico que constitui o meio poroso) são espécies totalmente miscíveis, a viscosidade da mistura solvente + soluto não varia com a concentração de soluto, há significativa transferência de massa entre as fases e a permeabilidade do meio poroso varia linearmente com a porosidade. O modelo é formado por duas Equações Diferenciais Parciais, uma do tipo Convecção-Difusão a outra é do tipo Convecção-Reação. Para resolução numérica, desenvolvemos uma metodologia que denominamos de EPEC (Explícita Porosidade e Explícita Concentração). Tal metodologia se baseia em um limitador de fluxo do tipo TVD e em diferenças finitas centradas de segunda ordem. Em adição, o EPEC emprega uma técnica de separação de operadores. Deste modo, em cada passo de tempo, realizamos inicialmente o cálculo explícito da porosidade seguido do cálculo explícito da concentração do solvente. Assim, obtemos um desacoplamento natural das equações que descrevem o problema. Resultados de simulações são apresentados para um meio poroso bidimensional, após sessenta dias de injeção de solvente. / We formulate a simplified Model to study the process of solvent injection in petroleum Reservoir, where the injected fluid (an acid) can partially dissolve a solid matrix. As prime hypotheses, we considered that solvent an soluble component are completely mixed, the viscosity of the fluid does not vary with the concentration of the soluble component, theres significant transfer of mass between the parts and, the permeability of media porous changes linearly with porosity. The model is formed by two Partial Differential Equation, one is convection-diffusion type and another is a convection-reaction type. The Numerical Resolution weve developed a method called EPEC (Explicit Porosity Explicit Concentration). Such methodology is based upon a Limiting of Flow of TVD type and, used Centered Finite Differences of second order. In addition, the EPEC use a operators separation technique. This way, every time, first we clearly calculate the porosity and then the concentration of solvent is calculated. Thus we obtain a natural decoupling of the equations that describe the problem. Simulation results are presented to a two dimensional media porous after sixty days of solvent injection.