Les propriétés de surface régissent le transport d'huile dans les nano-canaux : insights de la dynamique moléculaire

Résumé

Cette étude utilise des simulations de dynamique moléculaire pour élucider comment les interactions de surface influencent le transport d'huile dans les nano-canaux. Nous constatons que la largeur du canal, la polarité de l'huile, la composition du matériau et la rugosité de surface affectent tous de manière marquée le déplacement du centre de masse (CDM) et la diffusion des molécules d'huile. Dans un canal de 6 nm, le déplacement du CDM dépasse celui d'un canal de 2 nm de plus de 30 fois, et le coefficient de diffusion au centre du canal est près du double de celui près de la surface. Les composants polaires de l'huile subissent des interactions électrostatiques plus fortes avec le substrat de silice, réduisant drastiquement leur mobilité par rapport à l'octadécane non polaire. De plus, les substrats en or imposent une énergie d'interaction bien plus élevée, immobilisant presque l'huile, tandis que la calcite permet un déplacement modéré. Enfin, les cavités à l'échelle nanométrique modifient le transport local : les cavités peu profondes augmentent légèrement le déplacement en élargissant le canal, tandis que les cavités plus profondes entravent l'écoulement en raison d'une adsorption accrue. Ces résultats clarifient les mécanismes régissant la migration d'huile dans les géométries confinées et orientent les stratégies pour une extraction efficace du pétrole de schiste.

Contexte

Avec la hausse de la demande énergétique mondiale et l'épuisement des ressources conventionnelles, le pétrole de schiste non conventionnel a suscité un intense intérêt de recherche en raison de ses vastes réserves [1,2]. Bien que les ressources de pétrole de schiste dépassent le pétrole conventionnel de plus de trois fois [2], la fraction exploitable reste limitée en raison de son confinement dans des pores nanoscopiques allant de 2 à 100 nm [4,5]. Ce confinement induit des effets de surface prononcés qui altèrent les propriétés des fluides — telles que la densité, la mouillabilité et la diffusion — conduisant à des comportements de transport très différents des systèmes macroscopiques [10–12]. Par exemple, l'eau s'écoule plus rapidement dans les nanotubes de carbone que dans des tubes macroscopiques équivalents [6,7], et le caractère hydrophile/hydrophobe des nanotubes de carbone varie avec le diamètre [9]. Des travaux antérieurs en dynamique moléculaire ont montré que la contrainte de cisaillement et la taille des pores influencent la dynamique de l'eau et de l'huile dans des nano-fentes de carbone et de silice [13–15]. Cependant, les études systématiques sur la manière dont la chimie et la topographie de surface régissent le transport du pétrole de schiste sont rares. Comprendre ces effets est crucial pour optimiser les technologies d'extraction et a des implications plus larges pour les séparations membranaires, le transport biomoléculaire et les dispositifs de stockage d'énergie [16–20].

Méthodes

Toutes les simulations ont été effectuées à l'aide du module Discover de Materials Studio (Accelrys Inc.) avec le champ de force COMPASS tout-atome, qui capture avec précision les interactions liées, les termes croisés et les interactions non liées dans les phases condensées :

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

Les expressions complètes pour chaque terme sont données dans la publication originale. Les forces de van der Waals ont été modélisées avec un potentiel de Lennard-Jones (coupure 15,5 Å), et l'électrostatique a été traitée via la méthode d'Ewald. Les simulations ont été menées dans l'ensemble NVT à 298 K en utilisant le thermostat d'Andersen. Des conditions aux limites périodiques ont été appliquées dans toutes les directions. Les données étaient échantillonnées toutes les 5 ps.

Le substrat de silice, représentant le minéral dominant dans le schiste, a été construit à partir de la surface (001) et hydroxylé pour imiter les conditions géologiques [26]. Deux surfaces opposées, chacune de 1,5 × 7 × 0,85 nm³, créaient un canal fendu de largeurs 2, 4 et 6 nm. L'octadécane, la pyridine et le phénol — constituants courants du pétrole de schiste — ont été introduits pour atteindre une densité de masse de 0,8 g cm⁻³ (40 molécules dans le canal de 2 nm ; 407, 344, 80 et 120 molécules pour la pyridine, le phénol et l'octadécane dans les canaux de 4 et 6 nm, respectivement). Après minimisation d'énergie et une équilibration de 500 ps, une force de volume constante équivalente à 3,1 × 10⁻¹⁴ N par atome a été appliquée le long de l'axe du canal pour entraîner les molécules, une approche standard dans les études de transport MD [27–29].

Résultats et Discussion

Effet de la largeur du canal

Sous la force imposée, le nombre de molécules traversant le canal augmentait régulièrement (Information supplémentaire, Fig. S1). La Figure 2 montre que les canaux plus larges produisent des déplacements du CDM plus importants après 2 ns. La position du CDM est définie comme :

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

Dans le canal de 2 nm, le CDM n'a bougé que de 0,85 nm, alors que dans le canal de 6 nm, le déplacement dépassait 30 nm. Cette différence prononcée provient d'une adsorption plus forte dans les pores plus étroits, comme quantifié par l'énergie d'interaction moyenne :

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

La Figure 3 montre une corrélation inverse claire : une énergie d'interaction plus élevée entraîne un déplacement réduit, avec une augmentation triplement de l'énergie causant une chute de >30 fois du décalage du CDM. Les effets d'interaction diminuent à mesure que le canal s'élargit, soulignant la nature taille-dépendante du transport nanoscopique.

Les profils de densité en couches (Fig. 2) révèlent une couche interfaciale de 5 Å d'épaisseur adjacente aux parois. L'orientation de l'octadécane dans ces couches est majoritairement parallèle à la surface (θ ≈ 80°–90°), tandis que les molécules au centre du canal présentent des orientations aléatoires. Les coefficients de diffusion, dérivés des déplacements quadratiques moyens (Éq. 7), varient avec la position radiale : ils sont maximaux au centre et diminuent vers la paroi (Fig. 5a). Le canal de 6 nm affiche un profil parabolique marqué, tandis que le canal de 2 nm montre une diffusion quasi uniforme et faible, ressemblant à un front de type piston. Ces tendances s'expliquent par la compétition entre l'adsorption sur les parois (qui supprime le mouvement dans la région d'« ombre » proche de la surface) et le transfert de momentum intermoléculaire.

Effet de la polarité

Des simulations dans un canal de silice de 4 nm avec un nombre comparable de phénol, pyridine et octadécane illustrent le rôle de la polarité moléculaire. Malgré des forces identiques, le déplacement du CDM de l'octadécane est ~16 fois supérieur à celui du phénol ou de la pyridine (Fig. 8). Les énergies d'interaction totales (vdW + électrostatiques) sont nettement plus élevées pour les espèces polaires, l'électrostatique contribuant à environ la moitié du total (Fig. 8). Les moments dipolaires calculés par méthodes ab initio confirment que l'octadécane (0,032 Debye) est virtuellement non polaire, tandis que le phénol (1,306 Debye) et la pyridine (2,245 Debye) possèdent des dipoles significatifs. Par conséquent, les molécules polaires subissent une liaison plus forte à la surface de silice, réduisant leur mobilité.

Effet de la composition du matériau

Les comparaisons entre canaux de silice, d'or et de calcite (Fig. 9) montrent que l'or immobilise fortement l'octadécane, tandis que la calcite permet un déplacement modéré. La Figure 10 indique que le décalage du CDM corrèle avec l'énergie d'interaction moyenne : la silice présente l'énergie la plus faible et le déplacement le plus élevé, l'or l'inverse, et la calcite est intermédiaire. Les énergies similaires entre silice et calcite mais déplacements différents suggèrent que le type d'atome de surface et la texture influencent également le transport.

Effet de la rugosité de surface

L'introduction de cavités à l'échelle nanométrique (3 Å ou 6 Å de profondeur) dans la surface de silice modifie l'écoulement local. Les cavités piègent des molécules d'huile, ce qui réduit la vitesse des couches de fluide adjacentes (Fig. 11). Les déplacements du CDM pour les cavités de 3 Å et 6 Å sont de 3,95 Å et 3,07 Å, respectivement — légèrement supérieurs à la valeur de surface plane (3,17 Å) pour la cavité peu profonde, mais inférieurs pour la cavité plus profonde. L'amélioration avec les cavités peu profondes résulte probablement d'un élargissement effectif du canal, tandis que les cavités plus profondes imposent une adsorption plus forte et entravent l'écoulement.

Conclusions

Cette étude de dynamique moléculaire démontre que la largeur du nano-canal, la polarité de l'huile, le matériau du substrat et la rugosité de surface dictent collectivement le transport d'huile. Le canal de 6 nm offre un déplacement du CDM >30 fois supérieur à celui du canal de 2 nm, et le coefficient de diffusion au centre du canal double presque celui près de la paroi. Les molécules polaires et les substrats à forte interaction suppriment nettement la mobilité, tandis que les cavités de surface peu profondes peuvent modestement améliorer le transport. Ces insights sont essentiels pour concevoir des systèmes d'extraction nano-structurés et pour des applications plus larges où le transport de fluides dans des géométries confinées est critique.

Abréviations

CDM :

Centre de masse

COMPASS :

Potentiel moléculaire optimisé pour la phase condensée pour études de simulation atomistique

DM :

Dynamique moléculaire

vdW :

van der Waals