Nesta tese apresentamos novas formulações de elementos finitos mistas estabilizadas para o escoamento de Darcy, e para o escoamento acoplado de Stokes-Darcy. A aplicação do Método de Galerkin às formulações clássicas do sistema de Darcy é revista, e são comentadas algumas estabilizações existentes na literatura. A partir da combinação consistente de resíduos de mínimos quadrados da lei de Darcy, da equação de conservação de massa e do rotacional de lei de Darcy, são desenvolvidas novas formulações para a aproximação do escoamento no meio poroso, que possuem características de estabilidade superiores às existentes. É apresentada uma nova abordagem que permite o emprego destas formulações em meios porosos com interface de descontinuidade da condutividade hidráulica. Combinando estas novas formulações estabilizadas para o escoamento de Darcy, com formulações de Petrov-Galerkin para o escoamento de Stokes, são desenvolvidos métodos estabilizados para o sistema acoplado, que permitem o emprego de interpolações lagrangianas, de mesma ordem inclusive, para a velocidade e para a pressão, tanto no domínio de Stokes quanto no de Darcy, com interpolações contínuas ou descontínuas para a pressão. / In this thesis we propose new stabilized finite element methods for Darcy flow and for the coupled Stokes-Darcy flow. We review the Galerkin method applied to the classical variational formulations of the Darcy's system, and comment on some stabilization techniques proposed in the literature. Through the consistent combination of least square residuals of the Darcy's law, the mass balance equation and the curl of Darcy's law, we develop new finite element formulations to porous media flow, that have higher stability than the usual stabilizations. We present a methodology to incorporate restrictions that allows the application of these formulations to flow in heterogeneous porous media with an interface of discontinuity of the hydraulic conductivity tensor. By combining the new stabilized formulations for Darcy flow with Petrov-Galerkin formulations to Stokes flow, we develop stabilized methods for the coupled Stokes-Darcy flow that employ Lagrangian finite element spaces in both domains, continuous for the velocity and continuous or discontinuous for the pressure. Equal-order interpolations can be adopted for both velocity and pressure spaces.

Neste trabalho desenvolvemos a modelagem computacional dos fenômenos inerentes ao acoplamento hidromecânico que governam o escoamento de dois fluidos imiscíveis em uma matriz porosa heterogênea e deformável. As equações do modelo na escala de Darcy são decompostas em dois subsistemas associados a poromecânica e ao transporte dos fluidos. Neste contexto novos métodos numéricos são propostos para a computação da velocidade de Darcy baseados em técnicas de pós-processamento de Petrov-Galerkin em conjunção com métodos localmente conservativos para a equação hiperbólica de transporte da saturação aliados à técnica de decomposição de operadores para a computação dos efeitos da evolução temporal da porosidade induzidos pela deformação da matriz porosa sobre o transporte. As simulações numéricas do modelo resultante nos permite identificar diferentes regimes do acoplamento hidromecânico durante o processo de extração secundária de petróleo. Em particular analisamos a influência da razão de viscosidade entre os fluidos e da heterogeneidade da matriz porosa sobre os diferentes regimes de acoplamento. Dentre os vários fenômenos capturados no modelo damos particular ênfase ao surgimento de um processo de compactação retardada do reservatório devido à inundação de água onde, em particular, ilustramos seu efeito sobre as curvas de produção de petróleo. No contexto da modelagem estocástica analisamos também os efeitos das heterogeneidades e incertezas presentes nos coeficientes de permeabilidade e das constantes elásticas do meio poroso sobre os diferentes regimes do acoplamento geomecânico. A caracterização destes regimes governados pela razão de viscosidade e pelo coeficiente de variação das heterogeneidades é ilustrada por meio de simulações numéricas do processo de extração secundária de um reservatório sujeito ao peso das camadas superiores. / In this work we develop the computational modeling of the hydromechanical couplings which govern two-phase flow in a heterogeneous poroelastic media. At the Darcy scale the governing equations are decomposed in two subsystems associated with the poromechanics and hydrodynamics. In this context new numerical methods are proposed for the computation of the Darcy velocity based on Petrov-Galerkin post processing techiniques in conjunction with locally conservative methods for the hyperbolic transport equation and for water saturation together with an operator splitting technique for the computation of the transient porosity effect upon the transport equation. Numerical simulations allow us to identify different regimes of hydromechanics coupling during the secondary of petroleum withdrawal.In particular, we analyze the influence of the viscosity ratio and strength of heterogeneity upon the various coupling regimes. Among the many phenomena captured by the model we give particular emphasis on the delayed compaction of the reservoir due to the water flooding illustrating its effects upon the oil production curves. Within the framework of the stochastic modelinf, we analyze the effects of the heterogeneity and the uncertainty in permeability and elastic coeddicientsupon the hydromechanical coupling regimes. The characterization of these regimes governed by yhe viscosity ratio and by the strength of heterogeneity is illustrated in numerical simulations of secondary oil withdrawal from a reservoir subject to the overburden due the weight of the overlaying formations.

Made available in DSpace on 2015-03-04T18:50:41Z (GMT). No. of bitstreams: 1 Apresentacao.pdf: 81314 bytes, checksum: 1c3657984b3fd3891b1616bbf96f1250 (MD5) Previous issue date: 2006-11-24 / Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / In this thesis we propose new stabilized finite element methods for Darcy flow and for the coupled Stokes-Darcy flow. We review the Galerkin method applied to the classical variational formulations of the Darcy's system, and comment on some stabilization techniques proposed in the literature. Through the consistent combination of least square residuals of the Darcy's law, the mass balance equation and the curl of Darcy's law, we develop new finite element formulations to porous media flow, that have higher stability than the usual stabilizations. We present a methodology to incorporate restrictions that allows the application of these formulations to flow in heterogeneous porous media with an interface of discontinuity of the hydraulic conductivity tensor. By combining the new stabilized formulations for Darcy flow with Petrov-Galerkin formulations to Stokes flow, we develop stabilized methods for the coupled Stokes-Darcy flow that employ Lagrangian finite element spaces in both domains, continuous for the velocity and continuous or discontinuous for the pressure. Equal-order interpolations can be adopted for both velocity and pressure spaces. / Nesta tese apresentamos novas formulações de elementos finitos mistas estabilizadas para o escoamento de Darcy, e para o escoamento acoplado de Stokes-Darcy. A aplicação do Método de Galerkin às formulações clássicas do sistema de Darcy é revista, e são comentadas algumas estabilizações existentes na literatura. A partir da combinação consistente de resíduos de mínimos quadrados da lei de Darcy, da equação de conservação de massa e do rotacional de lei de Darcy, são desenvolvidas novas formulações para a aproximação do escoamento no meio poroso, que possuem características de estabilidade superiores às existentes. É apresentada uma nova abordagem que permite o emprego destas formulações em meios porosos com interface de descontinuidade da condutividade hidráulica. Combinando estas novas formulações estabilizadas para o escoamento de Darcy, com formulações de Petrov-Galerkin para o escoamento de Stokes, são desenvolvidos métodos estabilizados para o sistema acoplado, que permitem o emprego de interpolações lagrangianas, de mesma ordem inclusive, para a velocidade e para a pressão, tanto no domínio de Stokes quanto no de Darcy, com interpolações contínuas ou descontínuas para a pressão.

Made available in DSpace on 2015-03-04T18:50:55Z (GMT). No. of bitstreams: 1 TESE.pdf: 2353628 bytes, checksum: 5452725da7b8af8e0ba71da5e8c21c09 (MD5) Previous issue date: 2007-12-17 / Conselho Nacional de Desenvolvimento Cientifico e Tecnologico / In this work we develop the computational modeling of the hydromechanical couplings which govern two-phase flow in a heterogeneous poroelastic media. At the Darcy scale the governing equations are decomposed in two subsystems associated with the poromechanics and hydrodynamics. In this context new numerical methods are proposed for the computation of the Darcy velocity based on Petrov-Galerkin post processing techiniques in conjunction with locally conservative methods for the hyperbolic transport equation and for water saturation together with an operator splitting technique for the computation of the transient porosity effect upon the transport equation. Numerical simulations allow us to identify different regimes of hydromechanics coupling during the secondary of petroleum withdrawal.In particular, we analyze the influence of the viscosity ratio and strength of heterogeneity upon the various coupling regimes. Among the many phenomena captured by the model we give particular emphasis on the delayed compaction of the reservoir due to the water flooding illustrating its effects upon the oil production curves. Within the framework of the stochastic modelinf, we analyze the effects of the heterogeneity and the uncertainty in permeability and elastic coeddicientsupon the hydromechanical coupling regimes. The characterization of these regimes governed by yhe viscosity ratio and by the strength of heterogeneity is illustrated in numerical simulations of secondary oil withdrawal from a reservoir subject to the overburden due the weight of the overlaying formations. / Neste trabalho desenvolvemos a modelagem computacional dos fenômenos inerentes ao acoplamento hidromecânico que governam o escoamento de dois fluidos imiscíveis em uma matriz porosa heterogênea e deformável. As equações do modelo na escala de Darcy são decompostas em dois subsistemas associados a poromecânica e ao transporte dos fluidos. Neste contexto novos métodos numéricos são propostos para a computação da velocidade de Darcy baseados em técnicas de pós-processamento de Petrov-Galerkin em conjunção com métodos localmente conservativos para a equação hiperbólica de transporte da saturação aliados à técnica de decomposição de operadores para a computação dos efeitos da evolução temporal da porosidade induzidos pela deformação da matriz porosa sobre o transporte. As simulações numéricas do modelo resultante nos permite identificar diferentes regimes do acoplamento hidromecânico durante o processo de extração secundária de petróleo. Em particular analisamos a influência da razão de viscosidade entre os fluidos e da heterogeneidade da matriz porosa sobre os diferentes regimes de acoplamento. Dentre os vários fenômenos capturados no modelo damos particular ênfase ao surgimento de um processo de compactação retardada do reservatório devido à inundação de água onde, em particular, ilustramos seu efeito sobre as curvas de produção de petróleo. No contexto da modelagem estocástica analisamos também os efeitos das heterogeneidades e incertezas presentes nos coeficientes de permeabilidade e das constantes elásticas do meio poroso sobre os diferentes regimes do acoplamento geomecânico. A caracterização destes regimes governados pela razão de viscosidade e pelo coeficiente de variação das heterogeneidades é ilustrada por meio de simulações numéricas do processo de extração secundária de um reservatório sujeito ao peso das camadas superiores.

Estudamos problemas de Controle Ótimo de Sistemas de Parâmetros Distribuídos. O modelo consiste de uma equação diferencial parcial do tipo parabólica que modela o transporte de um poluente em um meio viscoso e incompressível, mais condições de contorno e condição inicial; ou seja, em nosso modelo estamos considerando além do termo difusivo, o termo convectivo. A modelagem matemática desenvolvida permite calcular a concentração do poluente que é despejada numa região do espaço, de forma que no tempo t=T, a concentração do poluente esteja o mais próximo possível da concentração máxima aceitável no meio. Mostramos a existência e unicidade de solução do sistema de estado e de um único controle ótimo sobre um conjunto de funções admissíveis. Caracterizamos o controle ótimo de forma a obter um sistema de otimalidade que permita o cálculo numérico do problema e mostramos a convergência do método. Provamos que o funcional de custo ótimo j(b)=J(u(b,t)) é diferenciável com derivadas parciais contínuas sobre um subconjunto convexo fechado U_ad do espaço de controles U. Como aplicação, estudamos o caso da contaminação de um rio/lagoa pelo mercúrio em água parada e em movimento.O objetivo maior desta parte do trabalho é minimizar, através dos controles, os efeitos causados pelos detritos do agente poluente. Na segunda parte deste trabalho, estudamos o problema de contato termoelástico que descreve as vibrações dinâmicas de uma viga que está engastada num de seus extremos enquanto o outro está livre e pode se mover entre dois obstáculos rígidos. O principal resultado desta parte do trabalho é mostrar a existência de solução do problema e que a energia associada ao sistema decai exponencialmente a zero quando o tempo vai para o infinito. Finalmente, em todos os casos, graficamos a solução do problema e analisamos os resultados.

Made available in DSpace on 2015-03-04T18:50:42Z (GMT). No. of bitstreams: 1 Apresentacao.pdf: 169709 bytes, checksum: 574c46462e9fcf63aa87693349129a55 (MD5) Previous issue date: 2006-11-29 / Conselho Nacional de Desenvolvimento Cientifico e Tecnologico / Estudamos problemas de Controle Ótimo de Sistemas de Parâmetros Distribuídos. O modelo consiste de uma equação diferencial parcial do tipo parabólica que modela o transporte de um poluente em um meio viscoso e incompressível, mais condições de contorno e condição inicial; ou seja, em nosso modelo estamos considerando além do termo difusivo, o termo convectivo. A modelagem matemática desenvolvida permite calcular a concentração do poluente que é despejada numa região do espaço, de forma que no tempo t=T, a concentração do poluente esteja o mais próximo possível da concentração máxima aceitável no meio. Mostramos a existência e unicidade de solução do sistema de estado e de um único controle ótimo sobre um conjunto de funções admissíveis. Caracterizamos o controle ótimo de forma a obter um sistema de otimalidade que permita o cálculo numérico do problema e mostramos a convergência do método. Provamos que o funcional de custo ótimo j(b)=J(u(b,t)) é diferenciável com derivadas parciais contínuas sobre um subconjunto convexo fechado U_ad do espaço de controles U. Como aplicação, estudamos o caso da contaminação de um rio/lagoa pelo mercúrio em água parada e em movimento.O objetivo maior desta parte do trabalho é minimizar, através dos controles, os efeitos causados pelos detritos do agente poluente. Na segunda parte deste trabalho, estudamos o problema de contato termoelástico que descreve as vibrações dinâmicas de uma viga que está engastada num de seus extremos enquanto o outro está livre e pode se mover entre dois obstáculos rígidos. O principal resultado desta parte do trabalho é mostrar a existência de solução do problema e que a energia associada ao sistema decai exponencialmente a zero quando o tempo vai para o infinito. Finalmente, em todos os casos, graficamos a solução do problema e analisamos os resultados.

Técnicas de injeção de traçadores são bastante utilizadas nos estudos de escoamentos em meios porosos heterogêneos, principalmente em problemas relacionados à simulação numérica de escoamentos miscíveis em reservatórios de petróleo e à dispersão de contaminantes em aqüíferos. Neste trabalho apresentamos novos algoritmos para a aproximação numérica do problema de injeção de traçadores. Apresentaremos desenvolvimentos recentes do método Forward Integral-Tube Tracking (FIT) que foi originalmente apresentado em Aquino et al. (2007a). O FIT é um método lagrangeano localmente conservativo utilizado na resolução de problemas de transporte linear. Este método não faz o uso de soluções de problemas de Riemann e baseia-se na construção dos tubos integrais introduzidas em Douglas Jr. et al. (2000b). Além disso, ele possui excelente eficiência computacional e é virtualmente livre de difusão numérica. Resultados numéricos são apresentados com o objetivo de comparar a precisão das soluções fornecidas por novas implementações do método FIT na resolução do problema do traçador em reservatórios de petróleo. / The injection of tracers are used in the investigation of flows in heterogeneous porous media, in studies related to the simulation of miscible dispacements in petroleum reservoirs and the dispersion of contaminants in aquifers. In this work we present new algorithms for the numerical approximation of tracer injection problems. We discuss recent developments of the Forward Integral-Tube Tracking (FIT) scheme which was introduced in Aquino et al. (2007a). The FIT is a locally conservative lagrangian scheme for the approximation of the linear transport problems. This scheme does not use analytic solutions of Riemann problems and is based on the construction of the integral tubes introduced in Douglas Jr. et al. (2000b). The FIT scheme is computationally very eficient and is virtually free of numerical diffusion. Numerical results are presented to compare the accuracy of the solutions provided by new implementation of the FIT scheme for the injection of tracers in petroleum reservoirs.

Técnicas de injeção de traçadores são bastante utilizadas nos estudos de escoamentos em meios porosos heterogêneos, principalmente em problemas relacionados à simulação numérica de escoamentos miscíveis em reservatórios de petróleo e à dispersão de contaminantes em aqüíferos. Neste trabalho apresentamos novos algoritmos para a aproximação numérica do problema de injeção de traçadores. Apresentaremos desenvolvimentos recentes do método Forward Integral-Tube Tracking (FIT) que foi originalmente apresentado em Aquino et al. (2007a). O FIT é um método lagrangeano localmente conservativo utilizado na resolução de problemas de transporte linear. Este método não faz o uso de soluções de problemas de Riemann e baseia-se na construção dos tubos integrais introduzidas em Douglas Jr. et al. (2000b). Além disso, ele possui excelente eficiência computacional e é virtualmente livre de difusão numérica. Resultados numéricos são apresentados com o objetivo de comparar a precisão das soluções fornecidas por novas implementações do método FIT na resolução do problema do traçador em reservatórios de petróleo. / The injection of tracers are used in the investigation of flows in heterogeneous porous media, in studies related to the simulation of miscible dispacements in petroleum reservoirs and the dispersion of contaminants in aquifers. In this work we present new algorithms for the numerical approximation of tracer injection problems. We discuss recent developments of the Forward Integral-Tube Tracking (FIT) scheme which was introduced in Aquino et al. (2007a). The FIT is a locally conservative lagrangian scheme for the approximation of the linear transport problems. This scheme does not use analytic solutions of Riemann problems and is based on the construction of the integral tubes introduced in Douglas Jr. et al. (2000b). The FIT scheme is computationally very eficient and is virtually free of numerical diffusion. Numerical results are presented to compare the accuracy of the solutions provided by new implementation of the FIT scheme for the injection of tracers in petroleum reservoirs.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho é apresentado um novo método acurado com passo de tempo fracionário, baseado em uma técnica de decomposição de operadores, para a solução numérica de um sistema governante de equações diferenciais parciais que modela escoamento trifásico água-gás-óleo imiscível em reservatórios de petróleo heterogêneos no qual os efeitos de compressibilidade do gás não foram levados em conta. A técnica de decomposição de operadores em dois níveis permite o uso de passos de tempo distintos para os três problemas definidos pelo procedimento de decomposição: convecção, difusão e pressão-velocidade. Um sistema hiperbólico de leis de conservação que modela o transporte convectivo das fases fluidas é aproximado por um esquema central de diferenças finitas explícito, conservativo, não oscilatório e de segunda ordem. Este esquema é combinado com elementos finitos mistos, localmente conservativos, para a aproximação numérica dos sistemas de equações parabólico e elíptico associados aos problemas de transporte difusivo e de pressão-velocidade, respectivamente. O operador temporal associado ao sistema parabólico é resolvido fazendo-se uso de uma estratégia implícita de solução (Backward Euler). O modelo matemático para escoamento trifásico considerado neste trabalho leva em conta as forças de capilaridade e expressões gerais para as funções de permeabilidade relativa, campos variáveis de porosidade e de permeabilidade e os efeitos da gravidade. A escolha de expressões gerais para as funções de permeabilidade relativa pode levar à perda de hiperbolicidade escrita e, desta maneira, à existência de uma região elíptica ou de pontos umbílicos para o sistema não linear de leis de conservação hiperbólicas que descreve o transporte convectivo das fases fluidas. Como consequência, a perda de hiperbolicidade pode levar à existência de choques não clássicos (também chamados de choques transicionais ou choques subcompressivos) nas soluções de escoamentos trifásicos. O novo procedimento numérico foi usado para investigar a existência e a estabilidade de choques não clássicos, com respeito ao fenômeno de fingering viscoso, em problemas de escoamentos trifásicos bidimensionais em reservatórios heterogêneos, estendendo deste modo resultados disponíveis na literatura para problemas de escoamentos trifásicos unidimensionais. Experimentos numéricos, incluindo o estudo de estratégias de injeção alternada de água e gás (Water-Alternating-Gas (WAG)), indicam que o novo procedimento numérico proposto conduz com eficiência computacional a resultados numéricos com precisão. Perspectivas para trabalhos de pesquisa futuros são também discutidas, tomando como base os desenvolvimentos reportados nesta tese. / We present a new, accurate fractional time-step method based on an operator splitting technique for the numerical solution of a system of partial differential equations modeling three-phase immiscible water-gas-oil flow problems in heterogeneous petroleum reservoirs in which the compressibility effects of the gas was not take into account. A two-level operator splitting technique allows for the use of distinct time steps for the three problems defined by the splitting procedure: convection, diffusion and pressure-velocity. A system of hyperbolic conservation laws modelling the convective transport of the fluid phases is approximated by a high resolution, nonoscillatory, second-order, conservative central difference scheme in the convection step. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the parabolic and elliptic problems associated with the diffusive transport of fluid phases and the pressure-velocity problem, respectively. The time discretization of the parabolic problem is performed by means of the implicit backward Euler method. The mathematical model for the three-phase flow considered in this work takes into account capillary forces and general expressions for the relative permeability functions, variable porosity and permeability fields, and the effect of gravity. The choice of general expressions for the relative permeability functions may lead to the loss of strict hyperbolicity and, therefore, to the existence of an elliptic region of umbilic points for the systems of nonlinear hyperbolic conservation laws describing the convective transport of the fluid phases. As a consequence, the loss of hyperbolicity may lead to the existence of nonclassical shocks (also called transitional shocks or undercompressive shocks) in three-phase flow solutions. The numerical procedure was used in an investigation of the existence and stability of nonclassical shocks with respect to viscous fingering in heterogeneous two-dimensional flows, thereby extending previous results for one-dimensional three-phase flow available in the literature. Numerical experiments, including the study of Water-Alternating-Gas (WAG) injection strategies, indicate that the proposed new numerical procedure leads to computational efficiency and accurate numerical results. Directions for further research are also discussed, based on the developments reported in this thesis.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho é apresentado um novo método acurado com passo de tempo fracionário, baseado em uma técnica de decomposição de operadores, para a solução numérica de um sistema governante de equações diferenciais parciais que modela escoamento trifásico água-gás-óleo imiscível em reservatórios de petróleo heterogêneos no qual os efeitos de compressibilidade do gás não foram levados em conta. A técnica de decomposição de operadores em dois níveis permite o uso de passos de tempo distintos para os três problemas definidos pelo procedimento de decomposição: convecção, difusão e pressão-velocidade. Um sistema hiperbólico de leis de conservação que modela o transporte convectivo das fases fluidas é aproximado por um esquema central de diferenças finitas explícito, conservativo, não oscilatório e de segunda ordem. Este esquema é combinado com elementos finitos mistos, localmente conservativos, para a aproximação numérica dos sistemas de equações parabólico e elíptico associados aos problemas de transporte difusivo e de pressão-velocidade, respectivamente. O operador temporal associado ao sistema parabólico é resolvido fazendo-se uso de uma estratégia implícita de solução (Backward Euler). O modelo matemático para escoamento trifásico considerado neste trabalho leva em conta as forças de capilaridade e expressões gerais para as funções de permeabilidade relativa, campos variáveis de porosidade e de permeabilidade e os efeitos da gravidade. A escolha de expressões gerais para as funções de permeabilidade relativa pode levar à perda de hiperbolicidade escrita e, desta maneira, à existência de uma região elíptica ou de pontos umbílicos para o sistema não linear de leis de conservação hiperbólicas que descreve o transporte convectivo das fases fluidas. Como consequência, a perda de hiperbolicidade pode levar à existência de choques não clássicos (também chamados de choques transicionais ou choques subcompressivos) nas soluções de escoamentos trifásicos. O novo procedimento numérico foi usado para investigar a existência e a estabilidade de choques não clássicos, com respeito ao fenômeno de fingering viscoso, em problemas de escoamentos trifásicos bidimensionais em reservatórios heterogêneos, estendendo deste modo resultados disponíveis na literatura para problemas de escoamentos trifásicos unidimensionais. Experimentos numéricos, incluindo o estudo de estratégias de injeção alternada de água e gás (Water-Alternating-Gas (WAG)), indicam que o novo procedimento numérico proposto conduz com eficiência computacional a resultados numéricos com precisão. Perspectivas para trabalhos de pesquisa futuros são também discutidas, tomando como base os desenvolvimentos reportados nesta tese. / We present a new, accurate fractional time-step method based on an operator splitting technique for the numerical solution of a system of partial differential equations modeling three-phase immiscible water-gas-oil flow problems in heterogeneous petroleum reservoirs in which the compressibility effects of the gas was not take into account. A two-level operator splitting technique allows for the use of distinct time steps for the three problems defined by the splitting procedure: convection, diffusion and pressure-velocity. A system of hyperbolic conservation laws modelling the convective transport of the fluid phases is approximated by a high resolution, nonoscillatory, second-order, conservative central difference scheme in the convection step. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the parabolic and elliptic problems associated with the diffusive transport of fluid phases and the pressure-velocity problem, respectively. The time discretization of the parabolic problem is performed by means of the implicit backward Euler method. The mathematical model for the three-phase flow considered in this work takes into account capillary forces and general expressions for the relative permeability functions, variable porosity and permeability fields, and the effect of gravity. The choice of general expressions for the relative permeability functions may lead to the loss of strict hyperbolicity and, therefore, to the existence of an elliptic region of umbilic points for the systems of nonlinear hyperbolic conservation laws describing the convective transport of the fluid phases. As a consequence, the loss of hyperbolicity may lead to the existence of nonclassical shocks (also called transitional shocks or undercompressive shocks) in three-phase flow solutions. The numerical procedure was used in an investigation of the existence and stability of nonclassical shocks with respect to viscous fingering in heterogeneous two-dimensional flows, thereby extending previous results for one-dimensional three-phase flow available in the literature. Numerical experiments, including the study of Water-Alternating-Gas (WAG) injection strategies, indicate that the proposed new numerical procedure leads to computational efficiency and accurate numerical results. Directions for further research are also discussed, based on the developments reported in this thesis.