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

Resumen Técnico: Un modelo transitorio totalmente acoplado de pozo y fractura para analizar la estabilidad de las fracturas durante el desembarque en yacimientos de gas profundo

Planteamiento del problema
La pérdida de circulación inducida por el efecto de succión (swab effect) durante las operaciones de desembarque (tripping-out) representa un riesgo crítico en la perforación de yacimientos de gas fracturados, particularmente en formaciones profundas donde las ventanas de presión son estrechas. A medida que la sarta de perforación se desplaza hacia arriba, se induce una reducción transitoria de la presión en el fondo del pozo (presión de succión). En formaciones fracturadas, esta caída de presión aumenta el esfuerzo efectivo sobre las paredes de la fractura, promoviendo el cierre de la fractura y potencialmente desencadenando la pérdida de fluidos. Si bien la investigación existente ha avanzado en el modelado de presiones transitorias en el pozo y en la mecánica de fractura de forma separada, ha faltado un marco unificado que capture la retroalimentación bidireccional entre las transitorias de presión del pozo y la deformación dinámica de la fractura. Los enfoques actuales suelen tratar las fracturas como límites estáticos o utilizan un acoplamiento secuencial, fallando al no contabilizar la interacción en tiempo real donde la deformación de la fractura altera la capacidad de flujo, lo que a su vez modifica la propagación de la presión.

Metodología
Los autores desarrollaron un modelo numérico transitorio totalmente acoplado que integra la hidráulica del pozo con el flujo de fluidos en la fractura y la deformación elástica. El sistema físico se idealiza como una fractura vertical que intersecta un pozo, donde el movimiento ascendente de la sarta de perforación reduce la presión del pozo.

• Ecuaciones gobernantes: El modelo acopla la ecuación de balance de masa del pozo (considerando la compresibilidad del fluido y el desplazamiento del pistón) con la ecuación de flujo de fluidos en la fractura (regida por la ley cúbica para el flujo laminar entre placas paralelas) y el modelo de constitución de la deformación de la fractura (que relaciona el ancho con el esfuerzo efectivo mediante la rigidez normal de la fractura).
• Mecanismo de acoplamiento: El sistema se resuelve mediante un esquema numérico totalmente implícito. En la interfaz pozo-fractura, el modelo impone la continuidad de la presión y la continuidad del flujo de masa. La tasa de flujo desde la fractura actúa como un término de fuente/sumidero en el balance de masa del pozo, mientras que la presión del pozo sirve como la condición de contorno para la difusión de la presión de la fractura.
• Solución numérica: Se emplea un esquema de Euler totalmente implícito (hacia atrás) para la discretización temporal para asegurar la estabilidad. El sistema de ecuaciones no lineales resultante se resuelve iterativamente utilizando el método de Newton-Raphson. La matriz Jacobiana se construye analíticamente para manejar el fuerte acoplamiento entre la presión y el ancho de la fractura.
• Validación: El modelo fue validado con datos de campo de un pozo de gas de esquisto profundo (Pozo Gao-X) en la Cuenca de Sichuan, dirigido a la Formación Dengying. Los resultados de la simulación mostraron una buena concordancia con las presiones medidas en el fondo del pozo, con errores relativos que oscilaron entre el 5.7% y el 6.5%.

Contribuciones clave y resultados
El estudio proporciona una base cuantitativa para comprender cómo los parámetros operativos y las propiedades intrínsecas de la formación interactúan para influir en la estabilidad de la fractura.

1. Dominio de la velocidad de desembarque: El análisis paramétrico identifica la velocidad de desembarque como la principal variable de control operativo. Aumentar la velocidad de 0.3 m/s a 1.5 m/s provoca una caída de presión en el pozo de 7.6 MPa, lo que resulta en una reducción del 78% en el ancho de la fractura (de 0.36 mm a 0.08 mm). Se identificó un umbral de velocidad crítica de aproximadamente 1.4 m/s, más allá del cual el cierre de la fractura se acelera notablemente.
2. Sensibilidad dependiente de la profundidad: El riesgo de cierre de la fractura depende fuertmente de la profundidad. A una velocidad constante de 0.3 m/s, la amplitud de cierre a 5200 m de profundidad es 3.2 veces mayor que a 1200 m. Esto se atribuye a una mayor presión hidrostática inicial y a una mayor atenuación de las ondas de presión a lo largo de trayectorias de propagación más largas en secciones más profundas.
3. Impacto del diseño de la sarta de perforación: La relación anular de la tubería de perforación (relación entre el diámetro exterior de la sarta y el diámetro interno del pozo) intensifica significamente el "efecto pistón". Aumentar esta relación de 0.65 a 0.68 reduce el ancho de la fractura en un 18%. En contraste, la relación anular del ensamblaje de fondo (BHA) tiene un efecto marginal.
4. Papel de la viscosidad del fluido: La viscosidad del fluido de perforación desempeña un papel secundario y menor. Aumentar la viscosidad de 0.025 a 0.035 Pa·s resulta en solo una reducción del 5.6% en el ancho de la fractura, lo que indica que la optimización de la viscosidad por sí sola es insuficiente para mitigar los riesgos de pérdida de circulación en comparación con el control de la velocidad y el diseño anular.
5. Propiedad intrínseca de la formación: La rigidez normal de la fractura se identifica como la propiedad intrínseca gobernante. Duplicar la rigidez de 50 a 100 GPa/m aumenta el ancho de la fractura en un 19% bajo la misma caída de presión. Esto explica por qué diferentes formaciones exhiben respuestas de estabilidad divergentes bajo las mismas operaciones de perforación.

Significancia e implicaciones de ingeniería
El artículo sostiene que su principal contribución no radica en la invención de nuevas ecuaciones físicas, sino en la integración exitosa de la mecánica de pozos y de fractura establecida en un marco transitorio totalmente acoplado. Este enfoque captura los mecanismos de retroalimentación bidireccional que los modelos desacoplados o acoplados secuencialmente omiten.

El estudio traduce estos conocimientos mecánicos en directrices de ingeniería específicas para yacimientos de gas profundo:

• Protocolos específicos por profundidad: Implementar estrategias de desembarque segmentadas con límites de velocidad reducidos (por ejemplo, ≤0.3 m/s para profundidades >4000 m, y ≤0.2 m/s para profundidades >5000 m) al aproximarse a intervalos fracturados.
• Optimización del diseño: Priorizar la optimización de la relación anular de la tubería de perforación (apuntando a ≤0.65) para minimizar el efecto pistón, en lugar de depender únicamente de ajustes en la viscosidad del fluido.
• Caracterización previa a la perforación: Incorporar la estimación de la rigidez normal de la fractura (vía DFIT, registros acústicos o datos de núcleos) en las evaluaciones geomecánicas para predecir los riesgos de estabilidad específicos de la formación.
• Gestión de riesgos: El marco propuesto respalda una transición de prácticas empíricas de "talla única" hacia una gestión de riesgos predictiva y específica para cada formación, permitiendo una evaluación más confiable del potencial de pérdida de circulación en operaciones complejas de pozos de gas.

Los autores reconocen limitaciones, incluyendo la suposición de fluidos Newtonianos y fracturas planas aisladas, señalando que trabajos futuros deberán abordar la reología no newtoniana, redes de fractura complejas, efectos térmicos y flujo multifásico para refinar aún más el modelo para aplicaciones en yacimientos de gas profundo.