Simulação Numérica de Escoamento Bifásico em reservatório de Petróleo Heterogêneos e Anisotrópicos utilizando um Método de Volumes Finitos “Verdadeiramente” Multidimensional com Aproximação de Alta Ordem

• Autor(a) principal: SOUZA, Márcio Rodrigo de Araújo
• Data de Publicação: 2015
• Tipo de documento: Tese
• Idioma: por
• Título da fonte: Repositório Institucional da UFPE
• dARK ID: ark:/64986/001300000nnv1
• Texto Completo: https://repositorio.ufpe.br/handle/123456789/17248
• topic: Anisotropia e heterogeneidade; Escoamento bifásico (óleo e água); Reservatórios de petróleo; MPFA; Multidimensional Limiting Process (MLP); Fluxo não convexo; Métodos verdadeiramente multidimensionais; Correção de entropia; Oil and water displacementsi; High order finite volume methods; Multipoint flux approximation methods; Multidimensional Limiting Process (MLP); Truly multidimensional methods; Entropy fix

Sob certas hipóteses simplificadoras, o modelo matemático que descreve o escoamento de água e óleo em reservatórios de petróleo pode ser representado por um sistema não linear de Equações Diferenciais Parciais composto por uma equação elíptica de pressão (fluxo) e uma equação hiperbólica de saturação (transporte). Devido a complexidades na modelagem de ambientes deposicionais, nos quais são incluídos camadas inclinadas, canais, falhas e poços inclinados, há uma dificuldade de se construir um modelo que represente adequadamente certas características dos reservatórios, especialmente quando malhas estruturadas são usadas (cartesianas ou corner point). Além disso, a modelagem do escoamento multifásico nessas estruturas geológicas incluem descontinuidades na variável e instabilidades no escoamento, associadas à elevadas razões de mobilidade e efeitos de orientação de malha. Isso representa um grande desafio do ponto de vista numérico. No presente trabalho, uma formulação fundamentada no Método de Volumes Finitos é estudada e proposta para discretizar as equações elíptica de pressão e hiperbólica de saturação. Para resolver a equação de pressão três formulações robustas, com aproximação dos fluxos por múltiplos pontos são estudadas. Essas formulações são abeis para lidar com tensores de permeabilidade completos e malhas poligonais arbitrárias, sendo portanto uma generalização de métodos mais tradicionais com aproximação do fluxo por apenas dois pontos. A discretização da equação de saturação é feita com duas abordagens com característica multidimensional. Em uma abordagem mais convencional, os fluxos numéricos são extrapolados diretamente nas superfícies de controle por uma aproximação de alta resolução no espaço (2ª a 4ª ordem) usando uma estratégia do tipo MUSCL. Uma estratégia baseada na Técnica de Mínimos Quadrados é usada para a reconstrução polinomial. Em uma segunda abordagem, uma variação de uma esquema numérico Verdadeiramente Multidimensional é proposto. Esse esquema diminui o efeito de orientação de malha, especialmente para malhas ortogonais, mesmo embora alguma falta de robustez possa ser observada pra malhas excessivamente distorcidas. Nesse tipo de formulação, os fluxos numéricos são calculados de uma forma multidimensional. Consiste em uma combinação convexa de valores de saturação ou fluxo fracionário, seguindo a orientação do escoamento através do domínio computacional. No entanto, a maioria dos esquemas numéricos achados na literatura tem aproximação apenas de primeira ordem no espaço e requer uma solução implícita de sistemas algébricos locais. Adicionalmente, no presente texto, uma forma modificada desses esquemas “Verdadeiramente” Multidimensionais é proposta em um contexto centrado na célula. Nesse caso, os fluxos numéricos multidimensionais são calculados explicitamente usando aproximações de alta ordem no espaço. Para o esquema proposto, a robustez e o caráter multidimensional também leva em conta a distorção da malha por meio de uma ponderação adaptativa. Essa ponderação regula a característica multidimensional da formulação de acordo com a distorção da malha. Claramente, os efeitos de orientação de malha são reduzidos. A supressão de oscilações espúrias, típicas de aproximações de alta ordem, são obtidas usando, pela primeira vez no contexto de simulação de reservatórios, uma estratégia de limitação multidimensional ou Multidimensional Limiting Process (MLP). Essa estratégia garante soluções monótonas e podem ser usadas em qualquer malha poligonal, sendo naturalmente aplicada em aproximações de ordem arbitrária. Por fim, de modo a garantir soluções convergentes, mesmo para problemas tipicamente não convexos, associados ao modelo de Buckley-Leverett, uma estratégia robusta de correção de entropia é empregada. O desempenho dessas formulações é verificado com a solução de problemas relevantes achados na literatura.