A fully coupled transient wellbore–fracture model for analyzing fracture stability during tripping out deep gas reservoirs

Resumo Técnico: Um Modelo Transiente Totalmente Acoplado de Poço–Fratura para Analisar a Estabilidade de Fraturas Durante o Desmonte de Reservatórios de Gás Profundos

Definição do Problema
A perda de circulação induzida pelo efeito de sucção (swab effect) durante as operações de desmonte (tripping-out) representa um risco crítico em perfurações de reservatórios de gás fraturados, particularmente em formações profundas onde as janelas de pressão são estreitas. À medida que a coluna de perfuração se move para cima, induz-se uma redução transiente na pressão do fundo do poço (pressão de swab). Em formações fraturadas, essa queda de pressão aumenta a tensão efetiva nas paredes das fraturas, promovendo o fechamento das fraturas e potencialmente desencadeando a perda de fluido. Embora pesquisas existentes tenham avançado nos modelos de pressão transiente do poço e na mecânica de fraturas separadamente, faltava um framework unificado que capturasse o feedback bidirecional entre os transientes de pressão do poço e a deformação dinâmica da fratura. As abordagens atuais frequentemente tratam as fraturas como limites estáticos ou utilizam acoplamento sequencial, falhando em considerar a interação em tempo real onde a deformação da fratura altera a capacidade de fluxo, o que, por sua vez, modifica a propagação da pressão.

Metodologia
Os autores desenvolveram um modelo numérico transiente totalmente acoplado, integrando a hidráulica do poço com o fluxo de fluido e a deformação elástica da fratura. O sistema físico é idealizado como uma fratura vertical que intercepta um poço, onde o movimento ascendente da coluna de perfuração reduz a pressão do poço, induzindo o fechamento elástico da fratura.

• Equações Governantes: O modelo acopla a equação de balanço de massa do poço (considerando a compressibilidade do fluido e o deslocamento do pistão) com a equação de fluxo de fluido da fratura (regida pela lei cúbica para fluxo laminar entre placas paralelas) e o modelo constitutivo de deformação da fratura (relacionando a largura à tensão efetiva via rigidez normal da fratura).
• Mecanismo de Acoplamento: O sistema é resolvido através de um esquema numérico totalmente implícito. Na interface poço–fratura, o modelo impõe continuidade de pressão e continuidade de fluxo de massa. A taxa de fluxo proveniente da fratura atua como um termo de fonte/sumidouro no balanço de massa do poço, enquanto a pressão do poço serve como a condição de contorno para a difusão de pressão da fratura.
• Solução Numérica: Um esquema de Euler totalmente implícito (retroativo) é empregado para a discretização temporal para garantir a estabilidade. O sistema de equações não lineares resultante é resolvido iterativamente utilizando o método de Newton–Raphson. A matriz Jacobiana é construída analiticamente para lidar com o forte acoplamento entre pressão e largura da fratura.
• Validação: O modelo foi validado contra dados de campo de um poço de gás de xisto profundo (Poço Gao-X) na Bacia de Sichuan, visando a Formação Dengying. Os resultados da simulação mostraram boa concordância com as pressões medidas no fundo do poço, com erros relativos variando de 5,7% a 6,5%.

Principais Contribuições e Resultados
O estudo fornece uma base quantitativa para entender como os parâmetros operacionais e as propriedades intrínsecas da formação interagem para influenciar a estabilidade da fratura.

1. Dominância da Velocidade de Desmonte: A análise paramétrica identifica a velocidade de desmonte como a principal variável de controle operacional. O aumento da velocidade de 0,3 m/s para 1,5 m/s causa uma queda de pressão no poço de 7,6 MPa, resultando em uma redução de 78% na largura da fratura (de 0,36 mm para 0,08 mm). Um limiar de velocidade crítica de aproximadamente 1,4 m/s foi identificado, além do qual o fechamento da fratura acelera marcadamente.
2. Sensibilidade Dependente da Profundidade: O risco de fechamento de fratura é fortemente dependente da profundidade. A uma velocidade constante de 0,3 m/s, a amplitude de fechamento a 5200 m de profundidade é 3,2 vezes maior do que a 1200 m. Isso é atribuído à maior pressão hidrostática inicial e à maior atenuação das ondas de pressão ao longo de caminhos de propagação mais longos em seções mais profundas.
3. Impacto do Design da Coluna de Perfuração: A razão anular da coluna de perfuração (razão entre o diâmetro externo da coluna e o diâmetro interno do poço) intensifica significativamente o "efeito pistão". Aumentar essa razão de 0,65 para 0,68 reduz a largura da fratura em 18%. Em contraste, a razão anular do BHA (Bottom Hole Assembly) tem um efeito marginal.
4. Papel da Viscosidade do Fluido: A viscosidade do fluido de perfuração desempenha um papel secundário e menor. O aumento da viscosidade de 0,025 para 0,035 Pa·s resulta em apenas uma redução de 5,6% na largura da fratura, indicando que a otimização da viscosidade isoladamente é insuficiente para mitigar os riscos de perda de circulação em comparação ao controle da velocidade e do design anular.
5. Propriedade Intrínseca da Formação: A rigidez normal da fratura é identificada como a propriedade intrínseca governante. Dobrar a rigidez de 50 para 100 GPa/m aumenta a largura da fratura em 19% sob a mesma queda de pressão. Isso explica por que diferentes formações exibem respostas de estabilidade divergentes sob as mesmas operações de perfuração.

Significância e Implicações de Engenharia
O artigo afirma que sua principal contribuição não reside na invenção de novas equações físicas, mas na integração bem-sucedida da mecânica de poço e de fratura em um framework transiente totalmente acoplado. Esta abordagem captura mecanismos de feedback bidirecional que modelos desacoplados ou sequencialmente acoplados negligenciam.

O estudo traduz esses insights mecanísticos em diretrizes de engenharia específicas para reservatórios de gás profundos:

• Protocolos Específicos por Profundidade: Implementar estratégias de desmonte segmentadas com limites de velocidade reduzidos (ex: ≤0,3 m/s para profundidades >4000 m, e ≤0,2 m/s para profundidades >5000 m) ao se aproximar de intervalos fraturados.
• Otimização de Design: Priorizar a otimização da razão anular da coluna de perfuração (visando ≤0,65) para minimizar o efeito pistão, em vez de depender apenas de ajustes de viscosidade do fluido.
• Caracterização Pré-Perfuração: Incorporar a estimativa da rigidez normal da fratura (via DFIT, logs acústicos ou dados de testemunho) nas avaliações geomecânicas para prever riscos de estabilidade específicos da formação.
• Gestão de Risco: O framework proposto apoia uma mudança de práticas empíricas de "tamanho único" para uma gestão de risco preditiva e específica para cada formação, permitindo uma avaliação mais confiável do potencial de perda de circulação em operações complexas de poços de gás.

Os autores reconhecem limitações, incluindo a suposição de fluidos Newtonianos e fraturas planas isoladas, observando que trabalhos futuros devem abordar reologia não-newtoniana, redes de fraturas complexas, efeitos térmicos e fluxo multifásico para refinar ainda mais o modelo para aplicações em reservatórios de gás profundos.