一、行业背景与模拟价值页岩油作为非常规油气资源的核心，其高效开发是能源领域的研究热点。CO₂驱替技术因能同时实现页岩油采收率提升与碳封存，成为行业重点方向。蒙脱土是页岩储层的核心黏土矿物，其纳米级孔道结构、表面性质直接影响油相运移和CO₂驱替效率。分子动力学（MD）模拟能从原子尺度揭示蒙脱石壁面-页岩油-CO₂的相互作用机理、流体运移规律（如壁面滑移、流速剖面），为现场驱替方案优化提供微观依据。本文基于LAMMPS平台，完整拆解蒙脱石壁面页岩油CO₂驱替的分子模拟全流程（建模→组装→平衡→驱替→后处理），涵盖核心代码、参数设计与商业化落地要点，可直接复用于科研项目、技术服务与成果转化。二、模拟基础条件与核心目标2.1 核心模拟参数（Real单位体系）参数项取值/说明温度/压强373K、38MPa（地层实际工况）油相组成正辛烷（147个）+3-甲基屈（70个），饱和烃:芳香烃=2.1:1空间尺寸蒙脱石尺寸157×36.23×11.929 Å，孔径宽度5nm，模拟盒70×37×75 Å力场体系蒙脱石（CLAYFF）+页岩油（OPLS-AA/L）+CO₂（TraPPE）核心观测指标驱替过程流速剖面、壁面滑移长度、油-CO₂/蒙脱石相互作用能、MSD（均方位移）、质心位移2.2 模拟核心目标揭示蒙脱石壁面对页岩油/CO₂的吸附特性与运移阻力；量化CO₂驱替过程中油相的运移规律（流速、扩散）；计算壁面滑移长度，优化储层改造方案；评估CO₂-油相的相互作用对驱替效率的影响。三、全流程建模详解（核心代码+逻辑）3.1 蒙脱石壁面建模：从基础文件到拓扑修正蒙脱石的原子级建模是整个体系的核心，需保证壁面结构与力场参数的准确性：基础文件复用：基于预生成的石英壁面文件150mmt_base.data（含4890个原子、11种原子类型，CGCMM力场），包含Si/Al/Mg/Na/O/H等原子的质量、LJ参数、电荷与坐标；壁面翻转处理：通过LAMMPS旋转指令修正壁面朝向，保证孔道空间匹配：units real boundary p p p atom_style full pair_style lj/cut/coul/long 2.57675 read_data 150mmt_base.data # 计算重心并绕Y轴旋转180° variable x0 equal xcm(all,x) variable y0 equal xcm(all,y) variable z0 equal xcm(all,z) displace_atoms all rotate ${x0} ${y0} ${z0} 0 1 0 180 write_data 150mmt.data力场关键参数：蒙脱石各原子的LJ势参数（如Na：ε=0.130100，σ=2.350013）需严格匹配CGCMM力场，保证壁面与流体的相互作用精度。3.2 页岩油分子建模：Packmol+Moltemplate双工具协同页岩油为混合体系（正辛烷+3-甲基屈），需先构建原子坐标，再补全拓扑/力场信息：Packmol构建混合油相坐标：编写model.inp定义分子数量与空间范围，生成包含147个正辛烷+70个3-甲基屈的model.pdb：tolerance 1 add_box_sides 1 filetype pdb output model.pdb structure c8.pdb # 正辛烷 **** structure 3ji.pdb # 3-甲基屈 *** end structure seed -1执行指令：packmol model.inp；Moltemplate补全拓扑信息：编写model.lt导入分子拓扑，生成LAMMPS可识别的model.data（含键/角/二面角）：import "c8.lt" import "3ji.lt" *** }执行指令：moltemplate.sh -pdb model.pdb model.lt；力场适配：油相采用OPLS-AA/L力场，电荷由1.2xCM5方法计算，model.in.settings中定义键/角/二面角的谐波势/OPLS势参数（如C-C键：k=268.0 kcal/(mol·Å²)，r0=1.53 Å）。3.3 挡板（活塞）建模：驱替动力的物理实现挡板（活塞）用于模拟CO₂驱替的“推进力”，分为左右挡板，通过LAMMPS直接创建原子区域：# 左挡板示例（right.data仅修改z范围为0-49） dimension 3 units real atom_style full boundary p p p region box block 0 1 0 36.23 0 74 units box *** run 1 write_data left.data执行指令：mpiexec -np 180 lmp_mpi -in makepistl.lmp（多核心加速）。3.4 CO₂分子建模：TraPPE力场的拓扑定义CO₂采用TraPPE模型，编写CO2.txt定义分子拓扑（3原子、2键、1角）与电荷（C：+0.7，O：-0.35）：3 atoms 2 bonds 1 angles Coords 1 0.0 0.0 0.0 2 -1.16 0.0 0.0 3 1.16 0.0 0.0 Types 1 1 # C 2 2 # O 3 2 # O Charges 1 0.7 2 -0.35 3 -0.35 Bonds 1 1 1 2 2 1 1 3 Angles 1 1 2 1 3四、多组分体系组装：尺寸匹配与参数整合多组分（油+蒙脱石+CO₂+挡板）的组装是模拟成功的关键，核心是原子类型规划与空间尺寸匹配，编写combine.in实现：4.1 核心参数规划（避免类型冲突）组分原子类型数键类型数角类型数核心说明页岩油868OPLS-AA/L力场蒙脱石1100CGCMM力场CO₂211TraPPE力场挡板200纯C原子，无键相互作用总计2379需提前定义create_box参数4.2 体系组装核心代码log log_1mmol.log units real boundary p p p atom_style full special_bonds lj/coul 0.0 0.0 0.5 kspace_style pppm 1.0e-4 # 长程静电相互作用 bond_style harmonic angle_style harmonic dihedral_style opls improper_style cvff pair_style lj/cut/coul/long 10.0 pair_modify mix arithmetic tail yes # 定义模拟盒尺寸（适配蒙脱石+油+CO₂+挡板） region box block -70 100 0 37 0 74 units box # 创建盒子，预分配拓扑空间 create_box 23 box bond/types 7 angle/types 9 dihedral/types 4 improper/types 1 **** # 补全CO₂/挡板力场参数 pair_coeff 20 20 0.05365 2.8 # CO₂-C pair_coeff 21 21 0.15697 3.05044 # CO₂-O pair_coeff 22 22 0.0207