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

Technische Samenvatting: Een volledig gekoppeld transiënt putboring–fractuurmodel voor het analyseren van fractuurstabiliteit tijdens het uittrekken van boorverbindingen in diepe gasreservoirs

Probleemstelling
Verlies van circulatie door het swab-effect tijdens het uittrekken van de boorverbinding vormt een kritiek risico bij het boren in gebroken gasreservoirs, met name in diepe formaties waar de drukvensters nauw zijn. Terwijl de boorstreng omhoog beweegt, induceert dit een transiënte reductie van de bodemdruk (swab-druk). In gebroken formaties verhoogt deze drukval de effectieve spanning op de fractuurwanden, wat fractuursluiting bevordert en potentieel fluidverlies triggert. Hoewel bestaand onderzoek de transiënte wellbore-drukmodellering en fractuurmechanica afzonderlijk heeft gevorderd, ontbrak een verenigd kader dat de bidirectionele feedback tussen transiënte wellbore-druk en dynamische fractuurdeformatie vastlegt. Huidige benaderingen behandelen fracturen vaak als statische grenzen of maken gebruik van sequentiële koppeling, waardoor ze er niet in slagen de real-time interactie te vatten waarbij fractuurdeformatie de doorvoercapaciteit verandert, wat op zijn beurt de drukpropagatie wijzigt.

Methodologie
De auteurs hebben een volledig gekoppeld, transiënt numeriek model ontwikkeld dat wellbore-hydraulica integreert met fractuurvloeistofstroom en elastische deformatie. Het fysieke systeem wordt geïdealiseerd als een verticale fractuur die een wellbore kruist, waarbij de opwaartse beweging van de boorstreng de wellbore-druk vermindert, wat een elastische sluiting van de fractuur induceert.

• Sturende Vergelijkingen: Het model koppelt de wellbore massa-balansvergelijking (rekening houdend met vloeistofcompressibiliteit en zuigerverplaatsing) met de fractuur vloeistofstroomvergelijking (geregeerd door de kubische wet voor laminaire stroming tussen parallelle platen) en het fractuurdeformatie constitutieve model (dat de breedte relateert aan de effectieve spanning via de normale fractuurstijfheid).
• Koppelingsmechanisme: Het systeem wordt opgelost via een volledig impliciet numeriek schema. Op het interface tussen de wellbore en de fractuur handhaaft het model continuïteit van druk en continuïteit van massavloei. De debiet uit de fractuur fungeert als bron/putterm in de wellbore massa-balans, terwijl de wellbore-druk dient als de randvoorwaarde voor de fractuur-drukdiffusie.
• Numerieke Oplossing: Een volledig impliciet (backward) Euler-schema wordt gebruikt voor temporele discretisatie om stabiliteit te garanderen. Het resulterende niet-lineaire stelsel van vergelijkingen wordt iteratief opgelost met behulp van de Newton–Raphson methode. De Jacobiaan-matrix wordt analytisch geconstrueerd om de sterke koppeling tussen druk en fractuurbreedte te beheersen.
• Validatie: Het model is gevalideerd tegen veldgegevens van een diepe schaliegasput (Well Gao-X) in het Sichuan-bekken, gericht op de Dengying-formatie. De simulatieresultaten vertoonden een goede overeenstemming met de gemeten bodemdrukken, met relatieve fouten variërend van 5,7% tot 6,5%.

Belangrijkste Bijdragen en Resultaten
De studie biedt een kwantitatieve basis voor het begrijpen van hoe operationele parameters en intrinsieke formatie-eigenschappen interageren om de fractuurstabiliteit te beïnvloeden.

1. Dominantie van de Boorsnelheid: Parametrische analyse identificeert boorsnelheid als de primaire operationele controlevariabele. Het verhogen van de snelheid van 0,3 m/s naar 1,5 m/s veroorzaakt een wellbore-drukval van 7,6 MPa, wat resulteert in een 78% reductie van de fractuurbreedte (van 0,36 mm naar 0,08 mm). Een kritische snelheidsdrempel van ongeveer 1,4 m/s werd geïdentificeerd, waarbij de fractuursluiting merkbaar versnelt.
2. Diepte-afhankelijke Sensitiviteit: Het risico op fractuursluiting is sterk diepte-afhankelijk. Bij een constante snelheid van 0,3 m/s is de sluitingsamplitude op 5200 m diepte 3,2 keer groter dan op 1200 m. Dit wordt toegeschreven aan de hogere initiële hydrostatische druk en de grotere attenuatie van drukgolven over langere propagatiepaden in diepere secties.
3. Impact van Boorstrengontwerp: De annulusverhouding van de boorbuis (de ratio van de buitenste diameter van de boorstreng tot de binnenste diameter van de wellbore) intensiveert het "zuigereffect" aanzienlijk. Het verhogen van deze ratio van 0,65 naar 0,68 vermindert de fractuurbreedte met 18%. In tegen tegenstelling hiertoe heeft de annulusverhouding van de BHA een marginaal effect.
4. Rol van Vloeistofviscositeit: De viscositeit van de boorvloeistof speelt een kleine, secundaire rol. Het verhogen van de viscositeit van 0,025 naar 0,035 Pa·s resulteert in slechts een 5,6% reductie van de fractuurbreedte, wat aangeeft dat optimalisatie van de viscositeit alleen onvoldoende is om het risico op circulatieverlies te beperken vergeleken met de controle op snelheid en annulusontwerp.
5. Intrinsieke Formatie-eigenschap: De normale fractuurstijfheid wordt geïdentificeerd als de leidende intrinsieke eigenschap. Het verdubbelen van de stijfheid van 50 naar 100 GPa/m verhoogt de fractuurbreedte met 19% onder dezelfde drukval. Dit verklaart waarom verschillende formaties uiteenlopende stabiliteitsreacties vertonen onder identieke booroperaties.

Betekenis en Engineering-implicaties
Het artikel stelt dat de primaire bijdrage niet ligt in het uitvinden van nieuwe fysieke vergelijkingen, maar in de succesvolle integratie van gevestigde wellbore- en fractuurmechanica in een volledig gekoppeld, transiënt kader. Deze benadering legt de bidirectionele feedbackmechanismen vast die eerdere ontkoppelde of sequentieel gekoppelde modellen negeren.

De studie vertaalt deze mechanistische inzichten naar specifieke engineering-richtlijnen voor diepe gasreservoirs:

• Diepte-specifieke Protocollen: Implementeer gesegmenteerde uittrekkingsstrategieën met gereduceerde snelheidslimieten (bijv. ≤0,3 m/s voor dieptes >4000 m, en ≤0,2 m/s voor dieptes >5000 m) wanneer men fractie-intervallen nadert.
• Ontwerpoptimalisatie: Geef prioriteit aan het optimaliseren van de annulusverhouding van de boorbuis (gericht op ≤0,65) om het zuigereffect te minimaliseren, in plaats van uitsluitend te vertrouwen op aanpassingen van de viscositeit van de boorvloeistof.
• Pre-Drill Karakterisering: Integreer de schatting van de normale fractuurstijfheid (via DFIT, akoestische logs of kerngegevens) in geomechanische beoordelingen om formatie-specifieke stabiliteitsrisico's te voorspellen.
• Risicobeheer: Het voorgestelde kader ondersteunt een verschuiving van empirische, "one-size-fits-all" praktijken naar voorspellend, formatie-specifiek risicobeheer, wat een betrouwbaardere beoordeling van het potentieel voor circulatieverlies mogelijk maakt in complexe gasput-operaties.

De auteurs erkennen beperkingen, waaronder de aanname van Newtoniaanse vloeistoffen en geïsoleerde planaire fracturen, en merken op dat toekomstig werk de non-Newtoniaanse rheologie, complexe fractuurnetwerken, thermische effecten en meerfasige stroming moet adresseren om het model verder te verfijnen voor toepassingen in diepe gasreservoirs.