Analysis of Casing Loads in Coalbed Methane Horizontal Wells during Multi-Stage Hydraulic Fracturing

Resumo Técnico: Análise de Cargas de Revestimento em Poços Horizontais de Gás de Xisto de Carvão durante Fraturamento Hidráulico de Múltiplas Etapas

Definição do Problema
Os recursos de gás de xisto de carvão (CBM) profundos, particularmente nos da Bacia de Ordos, na China, apresentam um valor estratégico significativo, mas impõem desafios de engenharia únicos. Durante o fraturamento hidráulico de múltiplas etapas de poços horizontais profundos, a propagação da fratura induz uma redistribuição complexa de tensões que se sobrepõe ao campo de tensão in situ original. Isso resulta em aumentos anormais nas tensões do revestimento, elevando os riscos de escoamento, deformação e falha. Os modelos teóricos e gráficos de engenharia existentes, amplamente derivados de reservatórios de gás de xisto ou camadas de carvão rasas, sofrem de três limitações primárias quando aplicados ao CBM profundo: (1) análise de sensibilidade quantitativa insuficiente dos parâmetros operacionais fundamentais (ex: taxa de injeção, espaçamento de fratura); (2) aplicabilidade limitada devido à falha em considerar a alta tensão in situ, forte heterogeneidade e redes de fraturas complexas específicas do carvão profundo; e (3) refinamento inadequado do modelo, visto que os modelos analíticos tradicionais dependem de suposições simplificadoras em relação à atenuação da pressão líquida e aos intervalos de interferência de tensão que não se sustentam em ambientes de carvão profundo.

Metodologia
O estudo foca na camada de carvão da Formação Taiyuan no Campo de Gás Daniudi, Bacia de Ordos, empregando uma abordagem integrada de múltiplas etapas:

1. Modelagem Geomecânica 3D: Um modelo geomecânico tridimensional de alta precisão foi construído utilizando métodos geoestatísticos, integrando dados estratigráficos, curvas de perfuração, experimentos de mecânica de rochas e testes de tensão in situ. A confiabilidade do modelo foi verificada contra dados de monitoramento eletromagnético de ampla área, atingindo um grau de correspondência da área afetada pela fratura superior a 92%.
2. Modelo de Tensão Induzida Modificado: Os autores desenvolveram um modelo de previsão de tensão induzida composto modificando a solução analítica clássica para a tensão induzida de fratura única. Dois coeficientes de correção foram introduzidos:
   • Coeficiente de Correção da Pressão Líquida da Fratura ($A$): Para contabilizar o aumento do vazamento de fluido devido ao sistema de clivagem (cleat) nas camadas de carvão.
   • Coeficiente de Correção da Distância de Atenuação de Tensão ($B$): Para abordar a atenuação da propagação de tensão dependente da distância em meios de carvão heterogêneos.
3. Inversão de Parâmetros: Cinco algoritmos de otimização inteligente (Otimização por Enxame de Partículas, Algoritmo de Busca de Pardal, Simulated Annealing, Algoritmo Genético e Mínimos Quadrados Não Lineares) foram comparados para inverter os coeficientes de correção. O Algoritmo de Busca de Pardal (SSA) foi selecionado por sua superioridade em convergência e precisão.
4. Superposição de Múltiplas Etapas: Um modelo de tensão composta foi estabelecido acoplando escala de tratamento, taxa de injeção, espaçamento de fratura e número de etapas usando o princípio da superposição linear.
5. Simulação de Elementos Finitos e Metodologia de Superfície de Resposta (RSM): Um modelo 3D de elementos finitos do sistema revestimento–camada de cimento–formação foi construído no ABAQUS. Um delineamento Box–Behnken (BBD) com 25 esquemas de simulação numérica foi utilizado para desenvolver um modelo de regressão quadrática para a previsão da tensão do revestimento.
6. Validação de Campo: O modelo preditivo foi validado contra dados de campo do Caso de Poço 3, comparando as etapas de escoamento previstas com as variações reais do diâmetro interno do revestimento medidas por registros de calibração (caliper logs).

Principais Contribuições e Resultados

• Precisão e Otimização do Modelo: O algoritmo SSA alcançou o melhor desempenho de inversão de parâmetros ($R^2 = 0,9371$, RMSE = 0,3761). Os coeficientes de correção ideais foram determinados como $A = 0,7330$ e $B = 0,9238$. Estes valores indicam que a pressão líquida efetiva em camadas de carvão profundas é aproximadamente 26,7% menor do que as previsões do modelo PKN teórico devido ao vazamento de fluido, e a distância de propagação de tensão é ligeiramente atenuada em comparação com modelos homogêneos.
• Tendências de Evolução de Tensão: A tensão induzida composta aumenta rapidamente durante as etapas iniciais de fraturamento e estabiliza após aproximadamente a sétima etapa. Taxas de injeção mais altas e menores espaçamentos de fratura foram encontrados como fatores que aumentam significativamente a tensão composta de pico. Especificamente, uma taxa de injeção de 24 m³/min com espaçamento de 50 m criou as condições de maior risco.
• Modelo de Regressão e Sensibilidade: Um modelo de regressão quadrática foi estabelecido vinculando parâmetros operacionais à tensão equivalente máxima do revestimento. O modelo demonstrou alta significância estatística ($P < 0,0001$) e precisão de ajuste ($R^2 = 0,9826$). A análise de sensibilidade revelou a seguinte classificação de influência na tensão do revestimento:
  1. Número de etapas de fraturamento
  2. Taxa de injeção
  3. Espaçamento de fratura
  4. Escala de tratamento
• Validação de Campo: No caso de validação (Caso de Poço 3), o modelo previu o escoamento do revestimento na 6ª etapa de fraturamento, enquanto os dados de calibração de campo indicaram que uma mudança abrupta no diâmetro interno (escoamento) ocorreu após a 4ª etapa. Os autores atribuem esse desvio de duas etapas a fatores como o uso de pressão interna média pelo modelo (ignorando flutuações dinâmicas), a simplificação das fraturas como entidades únicas simétricas (ignorando redes complexas), a heterogeneidade espacial dos parâmetros da rocha e a exclusão dos efeitos de tensão térmica provenientes de fluidos de fraturamento de baixa temperatura.

Significância e Alegações
O artigo afirma fornecer suporte teórico e técnico para a otimização de parâmetros de fraturamento e proteção da integridade do revestimento em reservatórios de gás de xisto de carvão profundos. Ao abordar a heterogeneidade específica e as condições de tensão das camadas de carvão profundo, o estudo oferece uma ferramenta de previsão refinada que supera as limitações dos modelos baseados em xisto existentes. Os gráficos de previsão de tensão de revestimento e as equações de regressão geradas são apresentados como ferramentas práticas para engenheiros ajustarem as taxas de injeção e o espaçamento de fratura para mitigar os riscos de falha do revestimento. O estudo enfatiza que, embora o modelo apresente uma precisão razoável, ele serve como um guia para a otimização de parâmetros e não como um preditor perfeito de eventos dinâmicos de campo, reconhecendo as complexidades inerentes aos reservatórios de carvão profundo.