藏在分子里的油藏密码：分子模拟如何揭示多组分烃类的气液相行为？
藏在分子里的油藏密码：分子模拟如何揭示多组分烃类的气液相行为？| MDPI Methane
论文标题：Phase Behaviour of Multicomponent Mixtures of Hydrocarbons: MD Simulation
论文链接： https://doi.org/10.3390/methane4040024
期刊名： Methane
期刊主页： https://www.mdpi.com/journal/methane
看不见、摸不着的油气在地下多孔岩石中究竟如何分布？当压力降低时，复杂烃类混合物中的轻重组分会怎样“各奔东西”？这一切的答案，都藏在分子的微观世界里。一项发表于Methane期刊的研究，利用分子动力学模拟技术，在计算机中复现了地下油气藏的真实相态变化，为我们揭示了一个充满细节的微观世界。
1. 研究背景
在油气开采的实践中，准确预测地层流体的相态变化至关重要。当钻井深入地下，压力随着开采而逐渐降低，原本处于液态的烃类混合物会分离出气体——这个过程被称为“液气两相分离”。
对于甲烷、乙烷等少数几种成分组成的简单体系，工程师们可以利用经典热力学模型进行预测。然而，真实油气藏的烃类流体是 包含数十甚至上百种组分 的多组分复杂混合物， 轻组分和重组分在气液两相中的分布极不均匀 。现有的模型在面对如此复杂的体系时，准确度大打折扣。
此外，地下储层广泛存在纳米级孔隙，在纳米受限空间内，流体分子与孔壁之间存在着复杂的相互作用，使得流体的相行为与体相（Bulk Phase）状态截然不同。 如何从分子层面理解多组分烃类混合物在宏观体相中的相变规律，成为提升油气资源评估与开发效率的关键。
分子动力学（MD）模拟是一种在原子/分子尺度上通过求解牛顿运动方程来研究物质微观演化过程的计算方法。相比于昂贵且难以在纳米尺度开展的物理实验，它提供了一种高效、可重复且可视化的研究手段。
2. 研究方法
为了研究多组分烃类混合物的气液相变行为，该研究利用分子动力学模拟技术，构建了一个包含 8种典型烃类组分的模型流体 ，其组成与 挥发油藏 的流体高度相似。
8组分模型流体的摩尔组成如下：
甲烷 (CH 4 )
30%
乙烷 (C 2 H 6 )
10%
正丁烷 (n-C 4 H 10 )
10%
正己烷 (n-C 6 H 14 )
10%
正辛烷 (n-C 8 H 18 )
10%
十二烷 (C 12 H 26 )
10%
十六烷 (C 16 H 34 )
10%
二十四烷 (C 24 H 50 )
10%
研究人员为C 12 、C 16 和C 24 这三种高分子量组分构建了支链分子模型，以更真实地反映原油组分的化学结构。
在恒定温度375.15K（约102 °C）下，通过改变模拟盒子的大小来控制体系密度，进而调控压力。他们系统性地研究了在 不同形状的气液界面 （平面、球形和圆柱形）条件下，气液两相的 相平衡与组分再分布 特征
3．研究结果
模拟结果生动地展示了多组分混合物在压力降低时发生相分离的微观过程。
核心发现 1：甲烷是气泡形成的主要驱动力
当压力降低至约 7 MPa 的泡点时，体系中开始形成球形气泡。根据图1分析表明， 甲烷分子是气泡的主要“建造者” 。气泡内部的气相几乎全部由甲烷组成，伴随少量的乙烷和丁烷，而重质组分则几乎全部留在液相中。
图1. 八组分混合物在恒温下的相图。低压区间表示两种共存相。数据为球形气泡。为了便于比较，图中还列出了甲烷、乙烷、正丁烷、正己烷和正辛烷的气液分离临界点（请注意，这些是实验数据，适用于热力学极限，而我们的有限尺寸系统并非如此）。空心圆圈表示球形气泡内部低密度相消失的情况。.
核心发现 2：气液界面处的组分富集现象
在气液界面层中，甲烷和乙烷的摩尔密度达到峰值，甚至高于它们在气泡内部和液体本体中的密度。这意味着， 界面区域对这些轻质组分有着独特的“吸附”或富集作用 ，揭示了气液传质过程中的关键微观机制。
图2. (a) 流体组分的摩尔密度（左侧纵轴，彩色曲线）与混合物的平均质量密度（右侧纵轴）。横轴对应于体素（即局部微小子体积）的累积体积分数；这些体素具有特定的混合物平均质量密度（数值可读自右侧纵轴）及混合物组分（由各组分的摩尔密度表示，数值可读自左侧纵轴）。图中的水平线段表示具有特定密度和组分的体素所占的体积分数区间。在本图的另外两个子图 (b) 和 (c) 中，混合物组分分别通过摩尔组分和质量密度进行表示。(b) 流体的摩尔组分（左侧纵轴，彩色曲线）与混合物的平均质量密度（右侧纵轴）。(c) 流体各组分的平均质量密度（左侧纵轴，彩色曲线）与混合物的平均质量密度（右侧纵轴；在此图中，该轴与左侧纵轴数值一致）。本图 (a–c) 所示数据均取自以下工况：体系尺寸为 26×25×25 nm³，内含一个直径约为 7 nm 的球形气泡；体系压力为 5.15 MPa，平均密度为 524.68 kg/m³。
核心发现 3：界面形状对相平衡的影响
通过对比球形、圆柱形和平面气液界面在不同压力下的密度变化，研究团队发现，界面曲率（即气泡的大小和形状）对相平衡状态有显著影响，这与经典的开尔文方程（描述弯曲液面饱和蒸气压变化的方程）所描述的现象定性一致。
表1. 主要模拟结果汇总表