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

Sintesi Tecnica: Un Modello Transitorio Completamente Accoppiato Pozzo–Frattura per l'Analisi della Stabilità delle Fratture durante le Operazioni di Estrazione in Serbatoi di Gas Profondi

Definizione del Problema
La perdita di circolazione indotta dall'effetto "swab" (aspirazione) durante le operazioni di estrazione (tripping-out) rappresenta un rischio critico nei giacimenti di gas fratturati, particolarmente nelle formazioni profonde dove le finestre di pressione sono strette. Mentre la colonna di perforazione si muove verso l'alto, induce una riduzione transitoria della pressione di fondo pozzo (pressione di swab). Nelle formazioni fratturate, questo calo di pressione aumenta lo stress efficace sulle pareti delle fratture, promuovendo la chiusura delle fratture e potenzialmente innescando la perdita di fluido. Sebbene la ricerca esistente abbia fatto progressi nella modellazione della pressione transitoria del pozzo e nella meccanica delle fratture separatamente, è mancato un quadro unificato capace di catturare il feedback bidirezionale tra i transitori di pressione del pozzo e la deformazione dinamica della frattura. Gli approcci attuali spesso trattano le fratture come confini statici o utilizzano un accoppiamento sequenziale, fallendo nel considerare l'interazione in tempo reale in cui la deformazione della frattura altera la capacità di flusso, che a sua volta modifica la propagazione della pressione.

Metodologia
Gli autori hanno sviluppato un modello numerico transitorio completamente accoppiato che integra l'idraulica del pozzo con il flusso di fluido nella frattura e la deformazione elastica. Il sistema fisico è idealizzato come una frattura verticale che interseca un pozzo, dove il movimento ascendente della colonna di perforazione riduce la pressione del pozzo, inducendo la chiusura elastica della frattura.

• Equazioni di Governo: Il modello accoppia l'equazione di bilancio della massa del pozzo (che tiene conto della compressibilità del fluido e dello spostamento del pistone) con l'equazione del flusso di fluido nella frattura (governata dalla legge cubica per il flusso laminare tra piastre parallele) e il modello costitutivo della deformazione della frattura (che mette in relazione la larghezza con lo stress efficace tramite la rigidezza normale della frattura).
• Meccanismo di Accoppiamento: Il sistema viene risolto tramite uno schema numerico completamente implicito. All'interfaccia pozzo–frattura, il modello impone la continuità della pressione e la continuità del flusso di massa. La portata proveniente dalla frattura agisce come termine sorgente/pozzo nel bilancio di massa del pozzo, mentre la pressione del pozzo funge da condizione al contorno per la diffusione della pressione nella frattura.
• Soluzione Numerica: Viene impiegato uno schema di discretizzazione temporale di Eulero completamente implicito (backward) per garantire la stabilità. Il risultante sistema non lineare di equazioni viene risolto iterativamente utilizzando il metodo di Newton–Raphson. La matrice Jacobiana è costruita analiticamente per gestire il forte accoppiamento tra pressione e larghezza della frattura.
• Validazione: Il modello è stato validato rispetto ai dati di campo di un pozzo di gas shale profondo (Pozzo Gao-X) nel Bacino del Sichuan, mirato alla Formazione Dengying. I risultati della simulazione hanno mostrato una buona concordanza con le pressioni di fondo misurate, con errori relativi compresi tra il 5,7% e il 6,5%.

Contributi Chiave e Risultati
Lo studio fornisce una base quantitativa per comprendere come i parametri operativi e le proprietà intrinseche della formazione interagiscano per influenzare la stabilità delle fratture.

1. Dominanza della Velocità di Estrazione: L'analisi parametrica identifica la velocità di estrazione come la principale variabile di controllo operativo. L'aumento della velocità da 0,3 m/s a 1,5 m/s causa una caduta di pressione nel pozzo di 7,6 MPa, determinando una riduzione del 78% della larghezza della frattura (da 0,36 mm a 0,08 mm). È stata identificata una soglia di velocità critica di circa 1,4 m/s, oltre la quale la chiusura della frattura accelera marcatamente.
2. Sensibilità Dipendente dalla Profondità: Il rischio di chiusura delle fratture dipende fortemente dalla profondità. A una velocità costante di 0,3 m/s, l'ampiezza della chiusura a 5200 m di profondità è 3,2 volte maggiore rispetto a quella a 1200 m. Ciò è attribuito alla maggiore pressione idrostatica iniziale e al maggiore attenuamento delle onde di pressione su percorsi di propagazione più lunghi nelle sezioni più profonde.
3. Impatto della Progettazione della Colonna di Perforazione: Il rapporto anulare della colonna di perforazione (rapporto tra il diametro esterno della colonna e il diametro interno del pozzo) intensifica significativamente l'effetto "pistone". L'aumento di questo rapporto da 0,6 5 a 0,68 riduce la larghezza della frattura del 18%. Al contrario, il rapporto anulare del BHA (Bottom Hole Assembly) ha un effetto marginale.
4. Ruolo della Viscosità del Fluido: La viscosità del fluido di perforazione gioca un ruolo secondario e minore. L'aumento della viscosità da 0,025 a 0,035 Pa·s risulta in una riduzione della larghezza della frattura di soli il 5,6%, indicando che l'ottimizzazione della sola viscosità è insufficiente per mitigare i rischi di perdita di circolazione rispetto al controllo della velocità e del design anulare.
5. Proprietà Intrinseca della Formazione: La rigidezza normale della frattura è identificata come la proprietà intrinseca determinante. Raddoppiare la rigidezza da 50 a 100 GPa/m aumenta la larghezza della frattura del 19% sotto lo stesso calo di pressione. Ciò spiega perché diverse formazioni esibiscono risposte di stabilità divergenti sotto identiche operazioni di perforazione.

Significato e Implicazioni Ingegneristiche
L'articolo sostiene che il suo contributo primario non risieda nell'invenzione di nuove equazioni fisiche, ma nel successo dell'integrazione di meccaniche consolidate del pozzo e della frattura in un quadro transitorio completamente accoppiato. Questo approccio cattura i meccanismi di feedback bidirezionale che i modelli precedentemente disaccoppiati o accoppiati sequenzialmente trascurano.

Lo studio traduce queste intuizioni meccanicistiche in specifiche linee guida ingegneristiche per i giacimenti di gas profondi:

• Protocolli Specifici per Profondità: Implementare strategie di estrazione segmentate con limiti di velocità ridotti (ad esempio, ≤0,3 m/s per profondità >4000 m, e ≤0,2 m/s per profondità >5000 m) quando ci si avvicina a intervalli fratturati.
• Ottimizzazione del Design: Dare priorità all'ottimizzazione del rapporto anulare della colonna di perforazione (puntando a ≤0,65) per minimizzare l'effetto pistone, piuttosto che fare affidamento esclusivamente sugli aggiustamenti della viscosità del fluido.
• Caratterizzazione Pre-Drilling: Incorporare la stima della rigidezza normale della frattura (tramite DFIT, log acustici o dati di carotaggio) nelle valutazioni geomeccaniche per prevedere i rischi di stabilità specifici della formazione.
• Gestione del Rischio: Il quadro proposto supporta il passaggio da pratiche empiriche "one-size-fits-all" a una gestione del rischio predittiva e specifica per la formazione, consentendo una valutazione più affidabile del potenziale di perdita di circolazione in operazioni complesse di pozzi di gas.

Gli autori riconoscono i limiti, tra cui l'assunzione di fluidi Newtoniani e di fratture planari isolate, notando che i lavori futuri dovrebbero affrontare la reologia non newtoniana, le reti di fratture complesse, gli effetti termici e il flusso multifase per perfezionare ulteriormente il modello per le applicazioni nei giacimenti di gas profondi.