A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

Diep onder de grond, in het complexe netwerk van gesteenten die olie, gas en water vasthouden, bewegen vloeistoffen niet altijd in eenvoudige, voorspelbare lijnen. Het gesteente zelf is zelden uniform; het is een lappendeken van materialen met zeer verschillende vermogens om vloeistof door te laten. Sommige lagen zijn als open snelwegen, terwijl andere als dichte muren fungeren. Wanneer ingenieurs proberen te modelleren hoe deze vloeistoffen stromen, moeten ze rekening houden met deze plotselinge veranderingen in de eigenschappen van het gesteente. Bovendien kunnen de vloeistoffen zelf lastig zijn. Terwijl water gemakkelijk stroomt, gedragen veel industriële vloeistoffen—zoals de slurries die worden gebruikt in de papierproductie of de boorvloeistoffen die worden gebruikt in olieputten—zich anders. Ze zijn "niet-Newtoniaans", wat betekent dat hun dikte en weerstand tegen stroming veranderen afhankelijk van hoe hard ze worden geduwd. Wanneer deze complexe vloeistoffen een gesteente tegenkomen dat abrupt van karakter verandert, wordt de wiskunde die hun beweging beschrijft ongelooflijk moeilijk en loopt deze vaak vast op precies de punten waar de eigenschappen van het gesteente verschuiven.

Decennialang hebben wetenschappers vertrouwd op een standaardmethode om deze stromingen te berekenen, een methode die in essentie de eigenschappen van het gesteente tussen twee punten middelt. Deze aanpak werkt goed genoeg wanneer het gesteente relatief consistent is, maar faalt dramatisch wanneer het gesteente zeer gevarieerd is. In die gevallen onderschat de standaardmethode consequent hoeveel vloeistof er daadwerkelijk kan bewegen, wat leidt tot fouten die enorm kunnen zijn. Om dit op te lossen, hebben onderzoekers aan de University of Science and Technology of China een nieuwe computationele techniek ontwikkend. In plaats van te vertrouwen op eenvoudige gemiddelden, gebruikt deze nieuwe methode een nauwkeurige wiskundige beschrijving van hoe een vloeistof zich gedraagt precies op de rand waar twee verschillende soorten gesteente elkaar ontmoeten. Door het ware, scherpe gedrag van de vloeistof bij deze kritieke grenzen vast te leggen, biedt de nieuwe methode een veel duidelijker en nauwkeuriger beeld van de ondergrondse stroming, zelfs wanneer het gesteente extreem ongelijkmatig is en de vloeistof complex is.

De onderzoekers concentreerden zich op een specifiek type vloeistof dat bekend staat als een power-law vloeistof. Dit zijn vloeistoffen waarbij de relatie tussen de kracht die hen duwt en de snelheid waarmee zij bewegen een specifieke, niet-lineaire regel volgt. In de echte wereld beslaat dit een breed scala aan substanties, van dikke polymeeroplossingen tot bepaalde soorten ruwe olie. De uitdaging ontstaat wanneer deze vloeistoffen door poreuze media stromen—gesteenten vol met minuscule gaatjes—waar de permeabiliteit, of de gemak van doorstroming, plotseling van de ene plek naar de andere verandert. In dergelijke omgevingen verandert de druk van de vloeistof niet geleidelijk; in plaats daarvan creëert het een scherpe piek of een "singulariteit" precies op de grens tussen het gemakkelijk stromende en het moeilijk stromende gesteente. Traditionele computermodellen, die het gesteente verdelen in een raster van kleine driehoeken, hebben moeite met deze pieken. Ze gaan ervan uit dat de druk lineair verandert tussen punten, wat fysiek onjuist is bij deze grenzen. Als gevolg hiervan produceren de modellen stroomsnelheden die veel te laag zijn, en ze vereisen een enorme hoeveelheid minuscule rasterpunten om zelfs een ruwe benadering van de werkelijkheid te krijgen.

Om dit op te lossen, ontwikkelde het team een Finite Analytic Method. Het kernidee is om te stoppen met het raden van hoe de vloeistof zich bij de grens gedraagt en in plaats daarvan een bekende, exacte wiskundige beschrijving van dat gedrag te gebruiken. De onderzoekers realiseerden zich dat de vloeistofdruk nabij het punt waar de gesteentseigenschappen veranderen, een specifiek gebogen patroon volgt, waarbij de druk stijgt of daalt op een manier die afhangt van de power-law natuur van de vloeistof. Ze berekenden de exacte vorm van deze curve, bekend als een eigenfunctie, voor het meest kritieke deel van de stroming. Deze berekening omvatte het oplossen van een complexe reeks vergelijkingen die beschrijven hoe de vloeistof beweegt in een wigvormig gebied waar de gesteentseigenschappen verschillen. Door de kleinste en meest significante oplossing van deze vergelijkingen te vinden, legden ze het dominante gedrag van de vloeistof vast, precies op de plek waar de problemen zich voordoen.

Zodra dit lokale gedrag begrepen was, bouwden de onderzoekers een nieuw numeriek schema rondom dit principe. In plaats van de oude, gebrekkige middelingsmethode te gebruiken om de rasterpunten te verbinden, gebruikten ze hun nieuwe analytische oplossing om de stroming over de grenzen van elke driehoekige cel in hun raster te berekenen. Dit stelde de computer in staat om precies te "weten" hoe de vloeistof zich zou gedragen terwijl deze door de nauwe hals tussen zones met een hoge en lage permeabiliteit werd geperst. Het resultaat is een simulatie die de fysica van het probleem respecteert op de meest moeilijke punten, in plaats van deze te verzachten met een benadering.

Het team testte deze nieuwe methode tegen verschillende uitdagende scenario's om te zien hoe goed deze presteerde. Ze begonnen met een klassiek "schaakbordpatroon", waarbij gesteenten met een hoge en lage permeabiliteit werden gerangschikt in afwisselende vierkanten. In deze opstelling veranderden de gesteentseigenschappen met een factor honderd. Wanneer ze hun nieuwe methode vergeleken met de traditionele methode, was het verschil treffend. De traditionele methode onderschatte de stroming met enorme marges, waarbij de werkelijke waarde soms met bijna negentig procent werd gemist op een grof raster. In contrast hiermee bleef de nieuwe methode opmerkelijk dicht bij het juiste antwoord, zelfs wanneer het raster niet erg fijn was. Naarmate ze het raster verfijnden, convergeerde de nieuwe methode zeer snel naar het juiste antwoord, terwijl de traditionele methode er moeite mee had om in te halen, zelfs met extreem fijne rasters.

Ze testten de methode ook op een regelmatig hexagonaal domein en op een domein met een willekeurige, log-normale verdeling van permeabiliteit, wat de chaotische aard van echte ondergrondse reservoirs nabootst. In elk geval bleek de nieuwe methode superieur. Het ging moeiteloos om met de extreme variaties in gesteentseigenschappen en bleef nauwkeurig, ongeacht het specifieke type power-law vloeistof dat werd gesimuleerd. Een van de meest opvallende bevindingen was dat de nauwkeurigheid van de nieuwe methode niet zwaar afhankelijk was van de specifieke eigenschappen van de vloeistof of het gesteente. Of de vloeistof nu zeer dik of zeer dun was, of of het contrast in het gesteente mild of extreem was, de methode leverde consistente resultaten. De onderzoekers ontdekten dat een matig niveau van rasterverfijning voldoende was om een hoge nauwkeurigheid te bereiken, waarbij de fouten zelfs in de slechtste scenario's onder de acht procent bleven.

Dit werk vertegenwoordigt een belangrijke stap voorwaarts in hoe we de stroming van vloeistoffen in complexe ondergrondse omgevingen simuleren. Door een gebrekkige middelingsmethode te vervangen door een oplossing die de ware fysica van de grens respecteert, hebben de onderzoekers een instrument gecreëerd dat zowel nauwkeuriger als efficiënter is. Voor industrieën die afhankelijk zijn van het begrijpen van hoe vloeistoffen door de aarde bewegen, van energieproductie tot milieubeheer, betekent dit dat simulaties beter te vertrouwen zijn en sneller kunnen worden uitgevoerd. De methode vereist niet de enorme rekenkracht die voorheen nodig was om een redelijk antwoord te krijgen, waardoor high-fidelity modellering van complexe, heterogene systemen veel toegankelijker wordt. De bevindingen suggereren dat door nauwkeurig te kijken naar de singulariteiten waar de wereld verandert, we een beter begrip van het geheel kunnen opbouwen.