A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

Profundamente no subsolo, na complexa rede de rochas que contêm petróleo, gás e água, os fluidos nem sempre se movem em linhas simples e previsíveis. A própria rocha raramente é uniforme; é um mosaico de materiais com capacidades vastamente diferentes de permitir a passagem de fluidos. Algumas camadas são como rodovias abertas, enquanto outras são como paredes densas. Quando os engenheiros tentam modelar como esses fluidos fluem, eles devem levar em conta essas mudanças súbitas nas propriedades da rocha. Além disso, os próprios fluidos podem ser complicados. Embora a água flua facilmente, muitos fluidos industriais — como as polpas usadas na fabricação de papel ou as lamas de perfuração usadas em poços de petróleo — comportam-se de maneira diferente. Eles são "não-newtonianos", o que significa que sua viscosidade e resistência ao fluxo mudam dependendo de quão fortemente são pressionados. Quando esses fluidos complexos encontram uma rocha que muda seu caráter abruptamente, a matemática que descreve seu movimento torna-se incrivelmente difícil, muitas vezes falhando justamente nos pontos onde as propriedades da rocha mudam.

Por décadas, cientistas confiaram em uma forma padrão de calcular esses fluxos, um método que essencialmente faz a média das propriedades da rocha entre dois pontos. Essa abordagem funciona bem o suficiente quando a rocha é relativamente consistente, mas falha dramaticamente quando a rocha é altamente variada. Nesses casos, o método padrão subestima consistentemente o quanto de fluido pode realmente se mover, levando a erros que podem ser massivos. Para resolver isso, pesquisadores da Universidade de Ciência e Tecnologia da China desenvolveram uma nova técnica computacional. Em vez de depender de médias simples, este novo método utiliza uma descrição matemática precisa de como o fluido se comporta exatamente na borda onde dois tipos diferentes de rocha se encontram. Ao capturar o comportamento real e nítido do fluido nessas fronteiras críticas, o novo método fornece uma imagem muito mais clara e precisa do fluxo subterrâneo, mesmo quando a rocha é extremamente irregular e o fluido é complexo.

Os pesquisadores focaram em um tipo específico de fluido conhecido como fluido de lei de potência. Estes são fluidos onde a relação entre a força que os empurra e a velocidade com que se movem segue uma regra não linear específica. No mundo real, isso abrange uma ampla gama de substâncias, desde soluções de polímeros espessas até certos tipos de petróleo bruto. O desafio surge quando esses fluidos fluem através de meios porosos — rochas cheias de pequenos buracos — onde a permeabilidade, ou a facilidade de fluxo, muda subitamente de um ponto para o outro. Em tais ambientes, a pressão do fluido não muda suavemente; em vez disso, cria um pico agudo ou uma "singularidade" exatamente na fronteira entre a rocha de fluxo fácil e a de fluxo difícil. Os modelos computacionais tradicionais, que dividem a rocha em uma grade de pequenos triângulos, têm dificuldade com esses picos. Eles assumem que a pressão muda em uma linha reta entre os pontos, o que é fisicamente incorreto nessas fronteiras. Como resultado, os modelos produzem taxas de fluxo muito baixas e exigem um número enorme de pontos de grade minúsculos para obter sequer uma aproximação grosseira da verdade.

Para corrigir isso, a equipe desenvolveu um Método Analítico Finito. A ideia central é parar de adivinhar como o fluido se comporta na fronteira e, em vez disso, usar uma descrição matemática exata desse comportamento. Os pesquisadores perceberam que, perto do ponto onde as propriedades da rocha mudam, a pressão do fluido segue um padrão curvo específico, subindo ou descendo de uma forma que depende da natureza de lei de potência do fluido. Eles calcularam a forma exata dessa curva, conhecida como autofunção, para a parte mais crítica do fluxo. Esse cálculo envolveu a resolução de um conjunto complexo de equações que descrevem como o fluido se move em uma área em forma de cunha onde as propriedades da rocha diferem. Ao encontrar a solução de menor ordem e mais significativa dessas equações, eles capturaram a forma dominante como o fluido se comporta justamente no ponto problemático.

Uma vez compreendido esse comportamento local, os pesquisadores construíram um novo esquema numérico em torno dele. Em vez de usar o antigo e falho método de média para conectar os pontos da grade, eles usaram sua nova solução analítica para calcular o fluxo através das fronteiras de cada célula triangular em sua grade. Isso permitiu que o computador "soubesse" exatamente como o fluido deveria se comportar ao ser espremido pelo estreito gargalo entre zonas de alta e baixa permeabilidade. O resultado é uma simulação que respeita a física do problema nos pontos mais difíceis, em vez de suavizá-los com uma aproximação.

A equipe testou este novo método contra vários cenários desafiadores para ver como ele se sairia. Começaram com um padrão clássico de "tabuleiro de xadrez", onde rochas de alta permeabilidade e baixa permeabilidade eram arranjadas em quadrados alternados. Nesta configuração, as propriedades da rocha mudavam por um fator de cem. Quando compararam o novo método ao tradicional, a diferença foi gritante. O método tradicional subestimava o fluxo por margens enormes, às vezes errando o valor real por quase noventa por cento em uma grade grossa. Em contraste, o novo método permaneceu notavelmente próximo da resposta correta, mesmo quando a grade não era muito fina. À medida que tornavam a grade mais fina, o novo método convergia para a resposta correta muito rapidamente, enquanto o método tradicional lutava para acompanhar, mesmo com grades extremamente finas.

Eles também testaram o método em um domínio hexagonal regular e em um domínio com uma distribuição log-normal aleatória de permeabilidade, que mimetiza a natureza caótica dos reservatórios subterrâneos reais. Em todos os casos, o novo método provou ser superior. Ele lidou com as variações extremas nas propriedades da rocha com facilidade e permaneceu preciso, independentemente do tipo específico de fluido de lei de potência sendo simulado. Uma das descobertas mais impressionantes foi que a precisão do novo método não dependia fortemente das propriedades específicas do fluido ou da rocha. Quer o fluido fosse muito espesso ou muito fino, ou se o contraste da rocha fosse moderado ou extremo, o método entregava resultados consistentes. Os pesquisadores descobriram que um nível moderado de refinamento de grade era suficiente para alcançar alta precisão, com erros caindo para menos de oito por cento mesmo nos piores cenários.

Este trabalho representa um passo significativo no avanço de como simulamos o fluxo de fluidos em ambientes subterrâneos complexos. Ao substituir uma técnica de média falha por uma solução que respeita a verdadeira física da fronteira, os pesquisadores criaram uma ferramenta que é tanto mais precisa quanto mais eficiente. Para indústrias que dependem da compreensão de como os fluidos se movem através da terra, da produção de energia à gestão ambiental, isso significa que as simulações podem ser mais confiáveis e executadas mais rapidamente. O método não exige o poder computacional massivo que era necessário anteriormente para obter uma resposta decente, tornando a modelagem de alta fidelidade de sistemas heterogêneos complexos muito mais acessível. Os achados sugerem que, ao olhar de perto para os pontos singulares onde o mundo muda, podemos construir uma melhor compreensão do todo.