Las Propiedades de Superficie Gobiernan el Transporte de Petróleo en Nano-canales: Perspectivas de Dinámica Molecular

Resumen

Este estudio emplea simulaciones de dinámica molecular para elucidar cómo las interacciones de superficie influyen en el transporte de petróleo en nano-canales. Encontramos que el ancho del canal, polaridad del petróleo, composición del material y rugosidad superficial afectan marcadamente el desplazamiento del centro de masa (COM) y difusión de moléculas de petróleo. En un canal de 6 nm, el desplazamiento COM excede al de un canal de 2 nm en más de 30 veces, y el coeficiente de difusión en el centro del canal es casi el doble del valor cerca de la superficie. Componentes polares de petróleo experimentan interacciones electrostáticas más fuertes con sustrato de sílice, reduciendo dramáticamente su movilidad comparado con octadecano no polar. Además, sustratos de oro imponen energía de interacción mucho mayor, casi inmovilizando el petróleo, mientras calcita permite desplazamiento moderado. Finalmente, cavidades nanoescalares modifican transporte local: cavidades poco profundas mejoran ligeramente desplazamiento al ensanchar el canal, mientras cavidades más profundas impiden flujo debido a mayor adsorción. Estos hallazgos clarifican mecanismos que gobiernan migración de petróleo en geometrías confinadas e informan estrategias para extracción eficiente de petróleo de lutitas.

Antecedentes

Con demanda energética global aumentando y recursos convencionales agotándose, el petróleo de lutitas no convencional ha atraído intenso interés investigativo por sus vastas reservas. A pesar de que recursos de petróleo de lutitas superan al petróleo convencional en más de tres veces, la fracción explotable permanece limitada debido a su confinamiento dentro de poros nanoscópicos de 2 a 100 nm. Este confinamiento induce efectos de superficie pronunciados que alteran propiedades de fluidos —densidad, mojabilidad y difusión— llevando a comportamientos de transporte que difieren marcadamente de sistemas macroscópicos. Por ejemplo, agua fluye más rápido en nanotubos de carbono que en tubos macroscópicos equivalentes, y el carácter hidrofílico/hidrofóbico de nanotubos de carbono cambia con diámetro. Trabajos previos de dinámica molecular han mostrado que esfuerzo de cizallamiento y tamaño de poro influyen dinámica de agua y petróleo en nanorranuras de carbono y sílice. Sin embargo, estudios sistemáticos sobre cómo la química y topografía de superficie gobiernan transporte de petróleo de lutitas son escasos. Entender estos efectos es crucial para optimizar tecnologías de extracción y tiene implicaciones más amplias para separaciones por membrana, transporte biomolecular y dispositivos de almacenamiento de energía.

Métodos

Todas las simulaciones se realizaron usando el módulo Discover de Materials Studio (Accelrys Inc.) con el campo de fuerza COMPASS all-atom, que captura con precisión interacciones enlazadas, de término cruzado y no enlazadas en fases condensadas:

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

Expresiones completas para cada término se dan en la publicación original. Fuerzas de van der Waals se modelaron con potencial Lennard-Jones (corte 15,5 Å), y electrostática se trató vía método Ewald. Simulaciones corrieron en ensamble NVT a 298 K usando termostato Andersen. Límites periódicos se aplicaron en todas direcciones. Datos se muestrearon cada 5 ps.

El sustrato de sílice, representando mineral dominante en lutitas, se construyó desde superficie (001) e hidroxilado para imitar condiciones geológicas. Dos superficies opuestas, cada una 1,5 × 7 × 0,85 nm³, crearon canal de hendidura con anchos de 2, 4 y 6 nm. Octadecano, piridina y fenol —constituyentes comunes de petróleo de lutitas— se introdujeron para lograr densidad masa 0,8 g cm⁻³ (40 moléculas en canal 2 nm; 407, 344, 80 y 120 moléculas para piridina, fenol y octadecano en canales 4 y 6 nm, respectivamente). Tras minimización de energía y equilibración de 500 ps, se aplicó fuerza de cuerpo constante equivalente a 3,1 × 10⁻¹⁴ N por átomo a lo largo del eje del canal para impulsar moléculas, enfoque estándar en estudios de transporte MD.

Resultados y Discusión

Efecto del Ancho del Canal

Bajo fuerza impuesta, el número de moléculas cruzando el canal aumentó constantemente (Información Suplementaria, Fig. S1). Figura 2 demuestra que canales más anchos producen mayores desplazamientos COM tras 2 ns. La posición COM se define como:

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

En canal 2 nm, COM se movió solo 0,85 nm, mientras en canal 6 nm el desplazamiento excedió 30 nm. Esta diferencia pronunciada surge de adsorción más fuerte en poros estrechos, cuantificada por energía de interacción promedio:

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

Figura 3 muestra correlación inversa clara: mayor energía de interacción conduce a desplazamiento reducido, con aumento triple de energía causando caída >30 veces en desplazamiento COM. Efectos de interacción disminuyen a medida que canal se ensancha, subrayando naturaleza tamaño-dependiente de transporte nanoescala.

Perfiles de densidad estratificados (Fig. 2) revelan capa interfacial de 5 Å adyacente a paredes. Orientación de octadecano dentro de estas capas es predominantemente paralela a superficie (θ ≈ 80°–90°), mientras moléculas en centro del canal exhiben orientaciones aleatorias. Coeficientes de difusión, derivados de desplazamientos cuadráticos medios (Ec. 7), varían con posición radial: son máximos en centro y disminuyen hacia pared (Fig. 5a). Canal 6 nm muestra perfil parabólico pronunciado, mientras canal 2 nm muestra difusión casi uniforme, baja, pareciendo frente tipo pistón. Estas tendencias se atribuyen a competencia entre adsorción en pared (suprime movimiento en región "sombra" cerca de superficie) y transferencia de momentum inter-molecular.

Efecto de Polaridad

Simulaciones en canal de sílice 4 nm con números comparables de fenol, piridina y octadecano ilustran papel de polaridad molecular. A pesar de fuerzas idénticas, desplazamiento COM de octadecano es ~16× mayor que fenol o piridina (Fig. 8). Energías de interacción totales (vdW + electrostáticas) son marcadamente mayores para especies polares, con electrostática contribuyendo aproximadamente la mitad del total (Fig. 8). Momentos dipolares calculados vía métodos de primeros principios confirman que octadecano (0,032 Debye) es virtualmente no polar, mientras fenol (1,306 Debye) y piridina (2,245 Debye) poseen dipolos significativos. Consecuentemente, moléculas polares experimentan unión más fuerte a superficie de sílice, reduciendo movilidad.

Efecto de Composición de Material

Comparaciones entre canales de sílice, oro y calcita (Fig. 9) muestran que oro inmoviliza fuertemente octadecano, mientras calcita permite desplazamiento moderado. Figura 10 indica que desplazamiento COM correlaciona con energía de interacción promedio: sílice exhibe menor energía y mayor desplazamiento, oro lo opuesto, y calcita intermedio. Energías similares entre sílice y calcita pero desplazamientos diferentes sugieren que tipo de átomo superficial y textura también influyen transporte.

Efecto de Rugosidad Superficial

Introducir cavidades nanoescalares (3 Å o 6 Å profundas) en superficie de sílice altera flujo local. Cavidades albergan moléculas de petróleo atrapadas, que reducen velocidad de capas de fluido adyacentes (Fig. 11). Desplazamientos COM para cavidades 3 Å y 6 Å son 3,95 Å y 3,07 Å, respectivamente —ligeramente mayores que valor superficie plana (3,17 Å) para cavidad poco profunda, pero menores para cavidad más profunda. Mejora con cavidades poco profundas probablemente resulta de ensanchamiento efectivo del canal, mientras cavidades más profundas imponen adsorción más fuerte y obstaculizan flujo.

Conclusiones

Este estudio de dinámica molecular demuestra que ancho de nano-canal, polaridad de petróleo, material de sustrato y rugosidad superficial dictan colectivamente transporte de petróleo. Canal 6 nm brinda >30 veces mayor desplazamiento COM que canal 2 nm, y coeficiente de difusión centro-canal casi duplica al cerca de pared. Moléculas polares y sustratos de alta interacción suprimen marcadamente movilidad, mientras cavidades superficiales poco profundas pueden mejorar modestamente transporte. Estos insights son esenciales para diseñar sistemas de extracción nano-estructurados y aplicaciones más amplias donde transporte de fluido en geometrías confinadas es crítico.

Abreviaturas

COM:

Centro de masa

COMPASS:

Potencial molecular optimizado para fases condensadas para estudios de simulación atomística

MD:

Dinámica molecular

vdW:

van der Waals