一种分子模拟不同温压条件下吸附初始模型快速建模方法技术

本发明专利技术涉及一种分子模拟不同温压条件下吸附初始模型快速建模方法，包括：步骤一、采用Materials Studio软件构建狭缝结构和吸附质粒子源，生成狭缝和吸附质粒子源坐标文件、拓扑文件，在力场参数文件Par.inp中定义原子的非键相互作用参数；步骤二、在计算控制文件in.conf中写入算法类型“GEMC NPT”；明确模拟目标温度、目标压力所在文本行数，将所有文件放在文件夹model中作为模板模型，使用开源软件GOMC进行模拟；步骤三、设置温度、压力数组，复制模板模型并重命名为目标温度压力$T.$P；步骤四、通过Shell语言修改目标模型的温度、压力，遍历温度、压力数组，批量生成各温度、压力下的吸附模型。本发明专利技术精度高，简便适用，不仅大幅提升批量建模效率，还能避免人为操作出错的可能性。可能性。可能性。

【技术实现步骤摘要】
一种分子模拟不同温压条件下吸附初始模型快速建模方法


[0001]本专利技术涉及油气田开发领域，具体涉及微观吸附模拟技术中分子模拟吸附模型批量建模方法。

技术介绍

[0002]随着能源结构的调整，页岩气发挥着越来越重要的作用。页岩气主要是指以游离、吸附、溶解状态存在于暗色泥页岩中的非常规天然气。其中，吸附气所占比例为20％～85％，主要吸附于孔径小于10nm的纳米孔隙中，其吸附能力与有机质含量、矿物成分、储层温度、压力、孔径大小等因素有关。目前，研究页岩气吸附的主要方法为物理模拟实验和分子模拟方法。物理模拟实验可获取页岩气在不同条件下的吸附量以及页岩储层条件、有机质性质和水分等因素对吸附的影响，但其存在明显的不足：成本高、高温高压限制、实验风险和安全隐患，且无法系统深入地揭示页岩气微观吸附机理。物理模拟实验的这些不足之处恰好是分子模拟技术的优势。分子模拟采用计算机辅助实验技术，以统计力学和量子力学为基础，基于原子水平的分子模型来模拟计算分子的结构与行为，以及体系的各种物理化学性质，其优势在于成本低、安全性高、具有预测性、能从原子和分子水平上直观再现物质的微观性质和各种物理化学过程等。
[0003]分子模拟主要分为分子力学(molecular mechanics,MM)、蒙特卡洛(Monte Carlo,MC)以及分子动力学(molecular dynamics,MD)模拟等。其中，蒙特卡洛方法常用来研究吸附。目前，页岩气吸附模拟采用巨正则蒙特卡洛方法研究石墨片层、干酪根、石英、伊利石、蒙脱石和高岭石等单结构模型以及混合结构模型中吸附质的吸附特性[Wang T,Tian S,Li G,et al.Molecular simulation of gas adsorption in shale nanopores:A critical review[J].Renewable and Sustainable Energy Reviews,2021,149:111391]。除了涉及不同矿物组成外，研究还涉及不同形状、不同孔径的孔隙模型，通过计算得到不同温度、压力条件下的吸附量，例如文献[李晶辉,韩鑫,黄思婧,等.页岩干酪根吸附规律的分子模拟研究[J].油气藏评价与开发,2022,12(03):455
‑
461]涵盖了8个孔径，各孔径下涵盖5个温度点、6个压力点，共计240组，具体地，模拟的孔径设置分别为1、2、3、4、5、10、15、20nm；温度介于298～378K；压力介于5～30MPa。类似的，文献[王晓琦,翟增强,金旭,等.地层条件下页岩有机质孔隙内CO2与CH4竞争吸附的分子模拟[J].石油勘探与开发,2016,43(05):772
‑
779]模拟了页岩气在343K、368K、398K、428K共4个温度下的吸附等温线，各温度下设置27个压力点，范围0～80MPa，共计108组。在使用蒙特卡洛方法研究吸附时，需要建立不同矿物结构模型以及各矿物不同孔径模型，提交计算时需要逐个修改温度和压力(逸度)，少则设置几十组，多则需要修改成百上千组。如果每个模型的温度、压力(逸度)都是人为操作修改，则需要不断地复制基本模型、打开模型、修改模型的温度、压力，工作量巨大，耗时费力，而且容易出现错误。因此，急需一种快速批量建立吸附初始模型的方法，实现各吸附模型中不同温度压力条件的修改，为后续对不同矿物以及不同孔径模型下模拟吸附提供模型基础。

技术实现思路

[0004]本专利技术的目的在于提供一种分子模拟不同温压条件下吸附初始模型快速建模方法，并对建模过程中的关键参数进行优化，该方法原理可靠、精度高，简便适用，不仅大幅度提升批量建模效率，减轻科研人员的工作强度，而且能够避免人为操作出错的可能性，克服了现有技术的缺陷和不足。
[0005]为达到以上技术目的，本专利技术采用以下技术方案。
[0006]一种分子模拟不同温压条件下吸附初始模型快速建模方法，依次包括以下步骤：
[0007]步骤一、建立吸附初始模型：
[0008]采用Materials Studio软件构建狭缝结构和吸附质粒子源，生成狭缝坐标文件BOX0.pdb、拓扑文件BOX0.psf和吸附质粒子源坐标文件BOX1.pdb、拓扑文件BOX1.psf，在力场参数文件Par.inp中定义原子的非键相互作用参数；
[0009]步骤二、建立模板模型：
[0010]在计算控制文件in.conf中写入算法类型“GEMC NPT”；明确模拟目标温度“Temperature”、目标压力“Pressure”所在文本行数，将所有文件放在文件夹model中作为模板模型，使用开源软件GOMC进行模拟；
[0011]步骤三、设置温度、压力数组，复制模板模型并重命名：
[0012]根据需要计算的目标温度、压力，确定温度计算区间、压力计算区间，并设置温度数组、压力数组；依次循环遍历温度数组、压力数组，通过Shell语言的cp命令复制模板模型文件夹，并重命名文件夹为目标温度压力$T.$P；
[0013]步骤四、通过Shell语言修改目标模型的温度、压力：
[0014]通过Shell语言的sed命令，将目标模型控制文件in.conf中的“Temperature”、“Pressure”行设置为目标温度$T、目标压力$P；循环重复步骤三和步骤四，遍历温度数组、压力数组，批量生成各温度、压力下的吸附模型。
[0015]进一步地，步骤一中，吸附质为甲烷，狭缝为石墨片层，甲烷采用TraPPE
‑
UA力场，石墨片层采用OPLS力场。
[0016]进一步的，更改吸附质并设置为循环变量，采用该方法，可以批量建立不同吸附质在给定孔径纳米孔中各温度、压力下的模型。
[0017]进一步的，更改纳米孔孔径并设置为循环变量，采用该方法，可以快速建立孔径、温度、压力正交的复合模型。
[0018]与现有技术相比，本专利技术具有以下有益效果：
[0019](1)本专利技术使用开源软件GPU Optimized Monte Carlo(GOMC)进行Gibbs系综蒙特卡洛模拟吸附，可直接进行温度压力设定，比巨正则系综方法减少了计算逸度的步骤。
[0020](2)本专利技术基于Shell语言修改各吸附模型中温度、压力，不仅大幅度提升了批量建模效率，减轻了科研人员的工作强度，而且避免了人为操作出错的可能性。
附图说明
[0021]构成本专利技术的一部分的说明书附图用来提供对本专利技术的进一步理解，与本专利技术的示意性实施例用于解释本专利技术，并不构成对本专利技术的不当限定。
[0022]图1是本专利技术的分子模拟吸附批量建模方法流程图。
[0023]图2是实施例模型吸附结果。
[0024]图3是各温度压力模型文件夹。
具体实施方式
[0025]下面根据附图和实施例进一步说明本专利技术。
[0026]应该指出，以下详细说明都是示例性的，旨在对本专利技术提供进一步的说明，以便于本本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种分子模拟不同温压条件下吸附初始模型快速建模方法，依次包括以下步骤：步骤一、建立吸附初始模型：采用Materials Studio软件构建狭缝结构和吸附质粒子源，生成狭缝坐标文件BOX0.pdb、拓扑文件BOX0.psf和吸附质粒子源坐标文件BOX1.pdb、拓扑文件BOX1.psf，在力场参数文件Par.inp中定义原子的非键相互作用参数；步骤二、建立模板模型：在计算控制文件in.conf中写入算法类型“GEMC NPT”；明确模拟目标温度“Temperature”、目标压力“Pressure”所在文本行数，将所有文件放在文件夹model中作为模板模型，使用开源软件GOMC进行模拟；步骤三、设置温度、压力数组，复制模板模型并重命名：根据需要计算的目标温度、压力，确定温度计算区间、压力计算区间，并设置温度数组、压力数组；依次循环遍历温度数组、压力数组，通过Shell语言的cp命令复制模板模型文件夹，并重命名文件夹为目标温度压力$T.$...

【专利技术属性】
技术研发人员：汪周华，赵建飞，余曹，郭平，刘煌，涂汉敏，胡义升，王烁石，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：