Oppervlakte-eigenschappen sturen olie transport in nanokanalen: inzichten uit moleculaire dynamica simulaties

Samenvatting

Deze studie gebruikt moleculaire dynamica simulaties om te verduidelijken hoe oppervlakte-interacties olie transport in nanokanalen beïnvloeden. We vinden dat kanaalbreedte, olie polariteit, materiaalsamenstelling en oppervlakteruwing de center-of-mass (COM) verplaatsing en diffusie van olie moleculen aanzienlijk beïnvloeden. In een 6-nm kanaal overschrijdt de COM verplaatsing die in een 2-nm kanaal met meer dan 30 keer, en de diffusiecoëfficiënt in het kanaalmidden is bijna het dubbele van die bij het oppervlak. Polaire olie componenten ervaren sterke elektrostatische interacties met het silice substraat, wat hun mobiliteit dramatisch reduceert ten opzichte van apolaire octadecaan. Goud substraten leggen een veel hogere interactie-energie op, bijna immobiliserend de olie, terwijl calciet matige verplaatsing toestaat. Tot slot wijzigen nanoschaal caviteiten lokaal transport: ondiepe caviteiten versterken verplaatsing licht door het kanaal effectief te verruimen, terwijl diepere caviteiten stroom belemmeren door verhoogde adsorptie. Deze bevindingen verhelderen de mechanismen die olie migratie in begrensde geometrieën sturen en informeren strategieën voor efficiënte shale-olie winning.

Achtergrond

Met de stijgende wereldwijde energievraag en het uitputten van conventionele bronnen, heeft onconventionele shale-olie intense onderzoeksinteresse getrokken vanwege de enorme reserves. Ondanks dat shale-olie reserves conventionele olie met meer dan drie keer overstijgen, blijft de exploiteerbare fractie beperkt door insluiting in nanoscoop poriën van 2 tot 100 nm. Deze insluiting induceert uitgesproken oppervlakte-effecten die vloeistofeigenschappen – zoals dichtheid, benettingsbaarheid en diffusie – veranderen, leidend tot transportgedrag dat merkbaar afwijkt van macroscopische systemen. Bijvoorbeeld stroomt water sneller in koolstofnanobuizen dan in equivalente macroscopische buizen, en het hydrofiele/hydrofobe karakter van koolstofnanobuizen verschuift met diameter. Eerder moleculaire dynamica werk toonde aan dat afschuifspanning en poriegrootte water- en olie-dynamica in koolstof- en silice nanoslits beïnvloeden. Systematische studies over hoe oppervlaktechemie en topografie shale-olie transport sturen, zijn echter schaars. Het begrijpen van deze effecten is cruciaal voor het optimaliseren van winningstechnologieën en heeft bredere implicaties voor membraan scheiding, biomoleculaire transport en energieopslag apparaten.

Methoden

Alle simulaties werden uitgevoerd met de Discover module van Materials Studio (Accelrys Inc.) met het COMPASS all-atom krachtveld, dat gebonden, cross-term en non-gebonden interacties in gecondenseerde fasen nauwkeurig vastlegt:

$$ {E}_{\mathrm{total}}={E}_{\mathrm{valence}}+{E}_{\mathrm{cross}\text{-}term}+{E}_{\mathrm{nonbond}} \tag{1} $$

Volledige uitdrukkingen voor elke term zijn gegeven in de oorspronkelijke publicatie. Van der Waals krachten werden gemodelleerd met een Lennard-Jones potentieel (afsnij 15,5 Å), en elektrostatica werden behandeld via de Ewald methode. Simulaties liepen in het NVT ensemble bij 298 K met de Andersen thermostaat. Periodieke randvoorwaarden werden in alle richtingen toegepast. Data werden elke 5 ps gesampled.

Het silice substraat, vertegenwoordigend het dominante mineraal in shale, was geconstrueerd vanuit het (001) oppervlak en gehydroxyleerd om geologische condities na te bootsen. Twee tegenoverliggende oppervlakken, elk 1,5 × 7 × 0,85 nm³, creëerden een spleetkanaal met breedtes van 2, 4 en 6 nm. Octadecaan, pyridine en fenol – gangbare shale-olie constituenten – werden geïntroduceerd om een bulk dichtheid van 0,8 g cm⁻³ te bereiken (40 moleculen in het 2-nm kanaal; 407, 344, 80 en 120 moleculen voor pyridine, fenol en octadecaan in de 4- en 6-nm kanalen, respectievelijk). Na energieminimalisatie en een 500-ps equilibratie, werd een constante volumekracht equivalent aan 3,1 × 10⁻¹⁴ N per atoom toegepast langs de kanaalas om de moleculen te drijven, een standaard aanpak in MD transport studies.

Resultaten en discussie

Effect van kanaalbreedte

Onder de opgelegde kracht nam het aantal moleculen dat het kanaal doorkruiste stabiel toe (Supporting Information, Fig. S1). Figuur 2 toont dat wijdere kanalen grotere COM verplaatsingen opleveren na 2 ns. De COM positie is gedefinieerd als

$$ {z}_{\mathrm{COM}}={\displaystyle \sum_i\frac{m_i}{M}{r}_i} \tag{5} $$

In het 2-nm kanaal bewoog de COM slechts 0,85 nm, terwijl in het 6-nm kanaal de verplaatsing 30 nm overschreed. Dit uitgesproken verschil ontstaat door sterkere adsorptie in smalle poriën, gekwantificeerd door de gemiddelde interactie-energie

$$ {E}_{\mathrm{average}\ \mathrm{interaction}}=\frac{E_{\mathrm{total}}-(E_{\mathrm{oil}}+E_{\mathrm{substrate}})}{N} \tag{6} $$

Figuur 3 toont een duidelijke inverse correlatie: hogere interactie-energie leidt tot verminderde verplaatsing, met een drievoudige energietoename die een >30-voudige daling in COM verschuiving veroorzaakt. Interactie-effecten nemen af naarmate het kanaal wijder wordt, wat de grootte-afhankelijke aard van nanoschaal transport benadrukt.

Gelaagde dichtheidsprofielen (Fig. 2) onthullen een 5-Å dikke interfaciale laag angrenzend aan de wanden. Octadecaan oriëntatie in deze lagen is voornamelijk parallel aan het oppervlak (θ ≈ 80°–90°), terwijl moleculen in het kanaalmidden willekeurige oriëntaties vertonen. Diffusiecoëfficiënten, afgeleid van gemiddelde kwadratische verplaatsingen (Eq. 7), variëren met radiale positie: ze zijn hoogst in het midden en dalen naar de wand (Fig. 5a). Het 6-nm kanaal vertoont een steile parabolische profiel, terwijl het 2-nm kanaal bijna uniforme, lage diffusie toont, gelijkend op een zuiger-achtige front. Deze trends kunnen toegeschreven worden aan de concurrentie tussen wandadsorptie (die beweging in de nabije-oppervlakte "schaduw" regio onderdrukt) en intermoleculaire impulsoverdracht.

Effect van polariteit

Simulaties in een 4-nm silice kanaal met vergelijkbare aantallen fenol, pyridine en octadecaan illustreren de rol van moleculaire polariteit. Ondanks identieke krachten is octadecaan's COM verplaatsing ~16× groter dan die van fenol of pyridine (Fig. 8). Totale interactie-energieën (vdW + elektrostatisch) zijn aanzienlijk hoger voor de polaire soorten, met elektrostaticiteit die ongeveer de helft van het totaal bijdraagt (Fig. 8). Dipolmomenten berekend via first-principle methoden bevestigen dat octadecaan (0,032 Debye) virtueel apolaire is, terwijl fenol (1,306 Debye) en pyridine (2,245 Debye) significante dipolen bezitten. Gevolgd ervaren polaire moleculen sterkere binding aan het silice oppervlak, wat mobiliteit reduceert.

Effect van materiaalsamenstelling

Vergelijkingen tussen silice, goud en calciet kanalen (Fig. 9) tonen aan dat goud octadecaan sterk immobiliseert, terwijl calciet matige verplaatsing toestaat. Figuur 10 geeft aan dat de COM verschuiving correleert met de gemiddelde interactie-energie: silice vertoont de laagste energie en hoogste verplaatsing, goud het tegenovergestelde, en calciet intermediër. De gelijke energiëen tussen silice en calciet maar verschillende verplaatsingen suggereren dat oppervlakte atoomtype en textuur ook transport beïnvloeden.

Effect van oppervlakteruwing

Het introduceren van nanoschaal caviteiten (3 Å of 6 Å diep) in het silice oppervlak wijzigt lokaal stroom. Caviteiten herbergen vastzittende olie moleculen, die de snelheid van aangrenzende vloeistofflagen reduceren (Fig. 11). COM verplaatsingen voor 3-Å en 6-Å caviteiten zijn 3,95 Å en 3,07 Å, respectievelijk – licht hoger dan de vlakke-oppervlakte waarde (3,17 Å) voor de ondiepe caviteit, maar lager voor de diepere caviteit. De verbetering met ondiepe caviteiten resulteert waarschijnlijk uit een effectieve verbreiding van het kanaal, terwijl diepere caviteiten sterkere adsorptie opleggen en stroom hinderen.

Conclusies

Deze moleculaire-dynamica studie toont aan dat nano-kanaalbreedte, olie polariteit, substraatmateriaal en oppervlakteruwing collectief olie transport dicteren. Het 6-nm kanaal biedt een >30-voudig grotere COM verplaatsing dan het 2-nm kanaal, en de diffusiecoëfficiënt in het kanaalmidden bijna verdubbelt die bij de wand. Polaire moleculen en hoge-interactie substraten onderdrukken mobiliteit aanzienlijk, terwijl ondiepe oppervlakte caviteiten transport matig kunnen verbeteren. Deze inzichten zijn essentieel voor het ontwerpen van nano-gestructureerde winningsystemen en voor bredere toepassingen waar vloeistoftransport in begrensde geometrieën cruciaal is.