表面性质主导纳米通道油运移：分子动力学洞见

摘要

本研究采用分子动力学模拟阐明表面相互作用如何影响纳米通道内油运移。发现通道宽度、油极性、材料组成、表面粗糙度均显著影响油分子质心 (COM) 位移和扩散。6 nm 通道中 COM 位移比 2 nm 通道大 30 余倍，通道中心扩散系数近乎近壁面值的两倍。极性油组分与二氧化硅基底静电相互作用更强，迁移率远低于非极性十八烷。金基底施加极高相互作用能，几乎固定油分子，而方解石允许中等位移。纳米级空腔改变局部运移：浅空腔略微增强位移(有效拓宽通道)，深空腔因吸附增加而阻碍流动。这些发现阐明限域几何中油迁移机制，为高效页岩油采收提供策略依据。

背景

全球能源需求上升、常规资源枯竭，非常规页岩油因巨大储量受关注 [1,2]。尽管页岩油资源超常规油三倍以上 [2]，可采比例仍受限于其限域于 2-100 nm 纳米孔隙 [4,5]。这种限域诱导显著表面效应，改变流体性质——密度、润湿性、扩散——导致运移行为与宏观体系截然不同 [10-12]。例如，水在碳纳米管流动快于等效宏观管 [6,7]，碳纳米管亲疏水性随直径变化 [9]。先前分子动力学工作显示剪切应力和孔径影响碳和二氧化硅纳米狭缝中水油动力学 [13-15]。但系统研究表面化学和形貌如何主导页岩油运移者鲜见。理解这些效应对优化采收技术至关重要，对膜分离、生物分子运移、储能装置也有广泛意义 [16-20]。

方法

所有模拟使用 Materials Studio (Accelrys Inc.) Discover 模块，采用 COMPASS 全原子力场，精准捕捉凝聚相中键合、交叉项、非键合相互作用：

$$ E_{total} = E_{valence} + E_{cross-term} + E_{nonbond} \tag{1} $$

各项完整表达式见原文。范德华力用 Lennard-Jones 势 (截断 15.5 Å)，静电用 Ewald 方法。NVT 系综 298 K 下用 Andersen 温控器模拟。三向周期性边界。每 5 ps 采样一次数据。

二氧化硅基底模拟页岩主要矿物，由 (001) 面构建并羟基化模拟地质条件 [26]。两相对表面各 1.5 × 7 × 0.85 nm³ 形成狭缝通道，宽度 2、4、6 nm。十八烷、吡啶、苯酚——页岩油常见组分——按体相密度 0.8 g cm⁻³ 引入 (2 nm 通道 40 分子；4、6 nm 通道分别 407、344、80、120 分子)。能量最小化和 500 ps 平衡后，沿通道轴施加恒定体积力 (相当于 3.1×10⁻¹⁴ N/原子) 驱动分子，这是 MD 运移研究标准方法 [27-29]。

结果与讨论

通道宽度效应

施加力下穿越通道分子数稳步增加 (支持信息 Fig. S1)。图 2 显示 2 ns 后宽通道 COM 位移更大。COM 位置定义为 $$ z_{COM} = \sum_i \frac{m_i}{M} r_i \tag{5} $$。2 nm 通道 COM 仅移动 0.85 nm，6 nm 通道超 30 nm。显著差异源于窄孔更强吸附，量化为平均相互作用能 $$ E_{average interaction} = \frac{E_{total} - (E_{oil} + E_{substrate})}{N} \tag{6} $$。图 3 显示反比相关：相互作用能增高导致位移减少，能量三倍增导致 COM 位移下降 30 余倍。相互作用效应随通道变宽减弱，凸显纳米运移尺寸依赖性。

分层密度分布 (图 2) 显示壁面邻近 5 Å 界面层。十八烷在界面层取向多平行表面 (θ ≈ 80°-90°)，通道中心分子取向随机。扩散系数由均方位移导出 (Eq. 7)，随径向位置变化：中心最大、近壁减小 (图 5a)。6 nm 通道呈陡峭抛物线分布，2 nm 通道近乎均匀低扩散，类似活塞前沿。趋势归因于壁面吸附(抑制近表面"阴影区"运动)与分子间动量传递的竞争。

极性效应

4 nm 二氧化硅通道中以可比数量模拟苯酚、吡啶、十八烷说明分子极性作用。尽管力相同，十八烷 COM 位移比苯酚或吡啶大 ~16 倍 (图 8)。总相互作用能 (vdW + 静电) 极性物种显著更高，静电贡献约一半 (图 8)。第一性原理计算偶极矩证实十八烷 (0.032 Debye) 几乎非极性，苯酚 (1.306 Debye) 和吡啶 (2.245 Debye) 具显著偶极。因此极性分子与二氧化硅表面结合更强，迁移率降低。

材料组成效应

二氧化硅、金、方解石通道对比 (图 9) 显示金强固定十八烷，方解石允许中等位移。图 10 指出 COM 位移与平均相互作用能相关：二氧化硅能量最低、位移最大；金相反；方解石居中。二氧化硅和方解石能量相近但位移不同，提示表面原子类型和纹理也影响运移。

表面粗糙度效应

在二氧化硅表面引入纳米级空腔 (3 Å 或 6 Å 深) 改变局部流动。空腔捕获油分子，降低相邻流体层速度 (图 11)。3 Å 和 6 Å 空腔 COM 位移分别为 3.95 Å 和 3.07 Å——浅空腔略高于平面值 (3.17 Å)，深空腔则更低。浅空腔增强可能源于有效拓宽通道，深空腔因吸附增强而阻碍流动。

结论

本分子动力学研究表明纳米通道宽度、油极性、基底材料、表面粗糙度共同决定油运移。6 nm 通道 COM 位移比 2 nm 大 30 余倍，中心扩散系数近乎近壁两倍。极性分子和高相互作用基底显著抑制迁移，浅表面空腸可适度增强运移。这些见解对设计纳米结构采收系统及更广泛限域几何流体运移应用至关重要。