可左右滑动选省市

煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法

更新时间:2024-11-03
煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法 专利申请类型:发明专利;
源自:北京高价值专利检索信息库;

专利名称:煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法

专利类型:发明专利

专利申请号:CN202410259337.8

专利申请(专利权)人:中国矿业大学(北京)
权利人地址:北京市海淀区学院路丁11号

专利发明(设计)人:徐超,王文静,王凯,郭琳,杨通,蔡永博,苑永旺,胡园园

专利摘要:本发明提供了煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法。该方法包括:基于晶体充填法,构建煤体致密固体基质下形态各异的纳米级单个孔隙模型;获取孔隙模型的结构参数和形态分布特征;采用蒙特卡洛‑结构优化‑分子动力学耦合方法,模拟煤体纳米孔隙模型中气体吸附过程,得到目标压力和温度下能量最低的平衡态吸附构型;对构型内部气体分子的赋存形式进行定量化表征;在平衡态吸附构型上增加外部气相空间,形成解吸模拟的初始构型,进而开展气体分子宏观解吸过程模拟。本发明实现了煤体真实纳米孔隙结构的构建,在分子尺度上还原了煤层气在煤体纳米孔隙结构中的吸附解吸过程,对煤层气产能预测以及煤与瓦斯突出防治具有重要指导意义。

主权利要求:
1.煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法,其特征在于,包括如下五个步骤:S1构建煤体内不同形态的纳米级单个孔隙模型:运用晶体充填法,构建煤体致密固体基质下形态各异的纳米级单个孔隙模型,在与构型平面垂直和平行的方向上均使用周期性边界,形成初始化的孔隙模拟空间;
S2获取孔隙模型的结构参数和形态分布特征:以氦气分子动力学直径大小的球形探针为标准,计算孔隙模型的体积和表面积大小;沿XZ平面方向制作孔隙形态切片,获取所建孔隙模型的形态分布特征;
S3孔隙模型内气体吸附过程模拟:将S1中构建的孔隙模型作为气体吸附的初始构型,获取气体逸度与压力的对应关系,以逸度为输入参数,采用蒙特卡洛‑结构优化‑分子动力学耦合方法,模拟单个孔隙模型在给定条件下对气体的物理吸附过程,得到目标压力值和温度值下能量最低的平衡态吸附构型和等温吸附曲线;
S4吸附平衡后气体分子赋存形式的定量化表征方法:利用原子组定义法则,在平衡态吸附构型中提取孔壁附近某一特定范围内的气体分子,定义为紧密束缚的吸附层,剩余气体分子定义为游离扩散的自由层,分别计算各分子层所占比例,实现不同赋存形式气体分子的量化区分;
S5孔隙模型内气体分子宏观解吸过程模拟:在S3所得平衡态吸附构型的Z方向上插入真空层,使得孔隙结构与外部气相空间相连;随后,应用分子动力学方法模拟煤体孔隙模型内气体分子瞬时解吸并扩散至气相空间的全过程,揭示实际工程应用中常见的宏观解吸现象;
在步骤S1中,所述的煤体纳米级单个孔隙模型的构建过程包括如下步骤:
S11以金刚石晶体作为内置充填材料的基础构型,将其聚合成不同的几何形状;
S12确定周期性晶胞盒子的尺寸(X、Y、Z),将聚合成的金刚石晶体构型固定于晶胞盒子的目标位置上;
S13基于受限空间充填‑分子动力学耦合方法,循环“煤大分子充填‑分子动力学方法平衡弛豫‑煤大分子充填”的过程,将煤大分子模型随机插入至晶胞盒子的剩余空间内,当晶胞盒子内煤大分子密度分布达到预期密度值时,停止循环过程;对完成充填的晶胞盒子依次进行几何结构优化、退火模拟和动力学弛豫,完成系统的能量和结构优化;
S14移除晶胞盒子中内置的金刚石晶体,得到目标孔径和形态的煤体纳米级单个孔隙模型;
在步骤S4中,所述的气体分子赋存形式的定量化表征方法包括如下步骤:
S41在平衡态吸附构型中,基于原子组定义法则依次选定孔壁附近0.1nm以内和0.1nm以外的气体分子,分别定义为紧密束缚的吸附层分子Qa和游离扩散的自由层分子Qf;
S42依据吸附平衡状态下孔隙模型中气体的绝对吸附量Q,得到不同赋存形式的吸附层分子Qa和自由层分子Qf占绝对吸附量Q的比例值na和nf,记为计算值;
S43基于绝对吸附量Q和S2中获取的孔隙结构参数,计算气体的过剩吸附量Qe,将过剩吸附量Qe与绝对吸附量Q的比值定义为吸附层分子比例的理论值ne,(1‑ne)定义为自由层分子比例的理论值,一一对比计算值和理论值,当误差范围控制在10%以内时,证明原子组定义法则中所取间距大小合理。
2.根据权利要求1所述的煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法,其特征在于,在步骤S5中,所述的气体分子宏观解吸过程模拟包括如下步骤:S51在平衡态吸附构型的Z方向上插入真空层,模拟吸附构型外部的气相空间,为平衡态吸附构型中气体分子的解吸扩散行为提供空间条件,进而形成气体分子宏观解吸模拟的初始构型;
S52对解吸模拟的初始构型进行总时长为100皮秒(ps)的动力学模拟,设定系综为正则系综(NVT),观察不同时段下孔隙模型中气体分子向低浓度气相空间内解吸扩散的全过程;
S53分析煤体纳米孔隙结构中气体分子的解吸特征。 说明书 : 煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法技术领域[0001] 本发明属于煤层气勘探开发技术领域,具体涉及煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法。背景技术[0002] 随着煤炭资源开发强度和深度的不断增加,煤层气的勘探潜力得以突破,推动着煤炭精准开采和能源供需跃进目标的早日实现。然而,煤体渗透率低、孔隙度小、地质构造复杂等因素使得煤层气的开发和对其储运特性的理解面临巨大的挑战。因此,煤层气在煤体内的赋存特性和运移规律依旧是目前煤层气藏开发研究中的热点和难点。[0003] 煤是一种天然的双重多孔介质,具有复杂多样的多尺度孔裂隙网络,限制着煤体储气含量、解吸速率及扩散流动等发展过程。其中,孔隙形态是重要的孔隙结构特征参数,对煤层气的富集程度、赋存状态和运移规律的影响不可忽视。虽然已有学者围绕孔隙形态在煤层气解吸和扩散过程中的重要作用开展相关研究,应用分子模拟手段从微观角度揭示煤层气的赋存机理,以刚性碳原子、石墨烯或碳纳米管为基础构型构建有机质孔隙模型,运用巨正则蒙特卡罗及分子动力学方法,提出了孔隙结构对气体解吸扩散规律影响的相似模拟研究,但发现这些研究中所建孔隙模型的壁面材料对气体吸附解吸特征影响较大,同时无法契合煤层气的实际储存环境。由此可知,对于不同形态孔隙结构中气体分子微观吸附解吸过程的探索不够深入,尚未构建真实的煤体孔隙结构模型,难以揭示煤体孔隙形态与煤层气储运特性之间的相应关系。[0004] 因此,亟待从微观分子层面提出一种煤体纳米级孔隙结构构建方法及气体吸附解吸过程的模拟方法,揭示煤体孔隙形态差异对煤层气吸附解吸特性的影响,进而解决煤体内不同赋存形式气体定量化表征难题,努力实现煤层气产能精准预测和煤与瓦斯突出防治的最终目标。发明内容[0005] 针对上述现有技术存在的问题,本发明提供了一种煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法,尤其涉及一种基于蒙特卡洛和分子动力学方法模拟煤体不同形态孔隙模型中气体吸附解吸过程的新方法。[0006] 为解决上述问题,本发明提供了如下的技术方案:[0007] 本方案提供的一种煤体纳米孔隙结构中气体吸附解吸过程的分子模拟方法,包括如下五个步骤:[0008] S1构建煤体内不同形态的纳米级单个孔隙模型:运用晶体充填法,构建煤体致密固体基质下形态各异的纳米级单个孔隙模型,在与构型平面垂直和平行的方向上均使用周期性边界,形成初始化的孔隙模拟空间;[0009] S2获取孔隙模型的结构参数和形态分布特征:以氦气分子动力学直径大小的球形探针为标准,计算孔隙模型的体积和表面积大小;沿XZ平面方向制作孔隙形态切片,获取所建孔隙模型的形态分布特征;[0010] S3孔隙模型内气体吸附过程模拟:将S1中构建的孔隙模型作为气体吸附的初始构型,获取气体逸度与压力的对应关系,以逸度为输入参数,采用蒙特卡洛‑结构优化‑分子动力学耦合方法,模拟单个孔隙模型在给定条件下对气体的物理吸附过程,得到目标压力值和温度值下能量最低的平衡态吸附构型和等温吸附曲线;[0011] S4吸附平衡后气体分子赋存形式的定量化表征方法:利用原子组定义法则,在平衡态吸附构型中提取孔壁附近某一特定范围内的气体分子,定义为紧密束缚的吸附层,剩余分子定义为游离扩散的自由层,分别计算各分子层所占比例,实现不同赋存形式气体分子的量化区分;[0012] S5孔隙模型内气体分子宏观解吸过程模拟:在S3所得平衡态吸附构型的Z方向上插入真空层,使得孔隙结构与外部气相空间相连;随后,应用分子动力学方法模拟煤体孔隙模型内气体分子瞬时解吸并扩散至气相空间的全过程,揭示实际工程应用中常见的宏观解吸现象。[0013] 在步骤S1中,所述的煤体纳米级单个孔隙模型的构建过程包括如下步骤:[0014] S11以金刚石晶体作为内置充填材料的基础构型,将其聚合成不同的几何形状;[0015] S12确定周期性晶胞盒子的尺寸(X、Y、Z),将聚合成的金刚石晶体构型固定于晶胞盒子的目标位置上;[0016] S13基于受限空间充填‑分子动力学耦合方法,循环“煤大分子充填‑分子动力学方法平衡弛豫‑煤大分子充填”的过程,将煤大分子模型随机插入至晶胞盒子的剩余空间内,当晶胞盒子内煤大分子密度分布达到预期密度值时,停止循环过程;对完成充填的晶胞盒子依次进行几何结构优化、退火模拟和动力学弛豫,完成系统的能量和结构优化;[0017] S14移除晶胞盒子中内置的金刚石晶体,得到目标孔径和形态的煤体纳米级单个孔隙模型。[0018] 在步骤S4中,所述的气体分子赋存形式的定量化表征方法包括如下步骤:[0019] S41在平衡态吸附构型中,基于原子组定义法则依次选定孔壁附近0.1nm以内和0.1nm以外的气体分子,分别定义为紧密束缚的吸附层分子Qa和游离扩散的自由层分子Qf;[0020] S42依据吸附平衡状态下孔隙模型中气体的绝对吸附量Q,得到不同赋存形式的吸附层分子Qa和自由层分子Qf占绝对吸附量Q的比例值na和nf,记为计算值;[0021] S43基于绝对吸附量Q和S2中获取的孔隙结构参数,计算气体的过剩吸附量Qe,具体公式如下:[0022] Qe=Q‑NAPVmp/RT(1)[0023] 式中:Qe——过剩吸附量,cm3/g;[0024] Q——绝对吸附量,cm3/g;[0025] NA——阿伏伽德罗常数,6.02×1023;[0026] Vmp——孔隙体积,cm3;[0027] R——通用气体常数,8.31J·mol‑1·K‑1;[0028] T——温度,K;[0029] 将过剩吸附量Qe与绝对吸附量Q的比值定义为吸附层分子比例的理论值ne,(1‑ne)定义为自由层分子比例的理论值,一一对比计算值和理论值,当误差范围控制在10%以内时,证明原子组定义法则中所取间距大小合理。[0030] 在步骤S5中,所述的气体分子宏观解吸过程模拟包括如下步骤:[0031] S51在平衡态吸附构型的Z方向上插入真空层,模拟吸附构型外部的气相空间,为平衡态吸附构型中气体分子的解吸扩散行为提供空间条件,进而形成气体分子宏观解吸模拟的初始构型;[0032] S52对初始构型进行总时长为100皮秒(ps)的动力学模拟,设定系综为正则系综(NVT),观察不同时段下孔隙模型中气体分子向低浓度气相空间内解吸扩散的全过程;[0033] S53分析煤体纳米孔隙结构中气体分子的解吸特征。[0034] 由于采用上述的技术方案,本发明专利的有益效果是:[0035] (1)本发明提出了基于煤大分子模型的煤体孔隙模型的构建方法,实现了煤体中不同形态纳米级孔隙结构的构建,还原了煤基质内孔隙结构的真实形貌,考虑了孔隙形态差异引起的煤层气储量差异,揭示了煤层气在真实煤体纳米级孔隙中的复杂赋存行为,使得从分子层面剖析煤体内部复杂孔隙网络与煤层气赋存微观机理之间的关系成为可能,为煤层气高效开采提供新思路。[0036] (2)本发明提出了一种煤体中不同赋存形式气体定量表征的计算方法,使得吸附态和游离态气体之间的划分更加便捷,实现了从分子角度定量化表征煤层气赋存形式的目标。[0037] (3)本发明提出了一种气体吸附解吸全过程的分子模拟方法,在分子尺度上还原了气体在煤体纳米孔隙结构中的吸附解吸过程,揭示了煤层气快速解吸现象的微观机理,对煤层气产能预测以及煤与瓦斯突出防治具有重要指导意义。附图说明[0038] 图1为煤体纳米孔隙结构中气体吸附解吸模拟过程的实施流程图;[0039] 图2为煤体纳米级孔隙模型的构建流程图;[0040] 图3为单个孔隙模型中沿XZ平面方向的孔隙形态切片;[0041] 图4为煤体纳米孔隙模型吸附甲烷时绝对吸附量随压力变化的等温吸附曲线;[0042] 图5为吸附平衡后孔隙模型内不同赋存形式甲烷分子的划分示意图;[0043] 图6为解吸模拟的初始构型和甲烷分子解吸过程示意图(以墨水瓶形孔隙为例);[0044] 图7为解吸过程中体系压力和体系能量随解吸时间变化的解吸曲线;具体实施方式[0045] 为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。[0046] 本实施例以煤体狭缝形、圆柱形、圆锥形和墨水瓶形孔隙模型的构建和甲烷在所建孔隙模型中的吸附解吸模拟过程为例对本发明的具体实施方式进行说明。[0047] 按照图1所示流程:[0048] S1构建煤体内不同形态的纳米级单个孔隙模型:[0049] 运用晶体充填法,构建煤体致密固体基质下狭缝形、圆柱形、圆锥形和墨水瓶形的纳米级单个孔隙模型,在与构型平面垂直和平行的方向上均使用周期性边界,形成初始化的孔隙模拟空间;[0050] S2获取孔隙模型的结构参数和形态分布特征:[0051] 以氦气分子动力学直径大小(0.26nm)的球形探针为标准,计算孔隙模型的体积和表面积大小;沿XZ平面方向制作孔隙形态切片,获取所建孔隙模型的形态分布特征;[0052] 表1孔隙模型的结构参数[0053] 3 2孔隙形状 体积/nm 表面积/nm狭缝形 13.26 41.65圆柱形 8.90 28.48圆锥形 3.59 17.05墨水瓶形 8.26 27.05[0054] S3孔隙模型内气体吸附过程模拟:[0055] 将S1中构建的孔隙模型作为气体吸附的初始构型,获取气体逸度与压力的对应关系,以逸度为输入参数,采用蒙特卡洛‑结构优化‑分子动力学耦合方法,模拟单个孔隙模型在给定条件下对气体的物理吸附过程,得到目标压力值和温度值下能量最低的平衡态吸附构型和等温吸附曲线;[0056] S4吸附平衡后气体分子赋存形式的定量化表征方法:[0057] 利用原子组定义法则,在平衡态吸附构型中提取孔壁附近某一特定范围内的气体分子,定义为紧密束缚的吸附层,剩余分子定义为游离扩散的自由层,分别计算各分子层所占比例,实现不同赋存形式气体分子的量化区分;[0058] S5孔隙模型内气体分子宏观解吸过程模拟:[0059] 在S3所得平衡态吸附构型的Z方向上插入真空层,使得孔隙结构与外部气相空间相连;随后,应用分子动力学方法模拟煤体孔隙模型内气体分子瞬时解吸并扩散至气相空间的全过程,揭示实际工程应用中常见的宏观解吸现象。[0060] 表2吸附解吸模拟过程的参数设置[0061] 模拟参数 设置 模拟参数 设置6方法 Metropolis 平衡步数 1×106精度控制 Customised 生产步数 1×10力场 COMPASSII 样本间隔 25电场 Forcefieldassigned 系综 NVT静电相互作用 Ewald&Group 时间步长 1fs非键相互作用 Atombased 逸度步数 10截断类型 0.001kcal/mol 温度 303K截断距离 1.25nm 压力 0.1~10MPa[0062] 参阅图2所示,在步骤S1中,所述的煤体纳米级单个孔隙模型的构建过程包括如下步骤:[0063] S11以金刚石晶体作为内置充填材料的基础构型,将其聚合成狭缝形(长、宽、高设置为2nm)、圆柱形(直径和高设为2nm)、圆锥形(内切圆半径和高设为2nm)和墨水瓶形(球形部分直径设为2nm、孔口直径设为1nm)的几何构型;[0064] S12确定周期性晶胞盒子的尺寸,其中狭缝形、圆柱形和圆锥形孔隙模型的晶胞盒子尺寸为:X=Y=5nm,Z=4nm;墨水瓶形孔隙模型:X=Y=3nm,Z=5nm,将聚合成的金刚石晶体构型固定于晶胞盒子的目标位置上;[0065] S13基于受限空间充填‑分子动力学耦合方法,循环“煤大分子充填‑分子动力学方法平衡弛豫‑煤大分子充填”的过程,将煤大分子模型(C196H132O10N2)随机插入至晶胞盒子的剩余空间内,当晶胞盒子内煤大分子密度分布达到预期密度值时,停止循环过程;对完成充填的晶胞盒子依次进行几何结构优化、退火模拟和动力学弛豫,完成系统的能量和结构优化;[0066] S14移除晶胞盒子中内置的金刚石晶体,得到煤体中孔径为2nm的单个狭缝形、圆柱形、圆锥形和墨水瓶形孔隙结构模型。[0067] 在步骤S4中,所述的气体分子赋存形式的定量化表征方法包括如下步骤:[0068] S41在平衡态吸附构型中,基于原子组定义法则依次选定孔壁附近0.1nm以内和0.1nm以外的气体分子,分别定义为紧密束缚的吸附层分子Qa和游离扩散的自由层分子Qf;[0069] S42依据吸附平衡状态下孔隙模型中气体的绝对吸附量Q,得到不同赋存形式的吸附层分子Qa和自由层分子Qf占绝对吸附量Q的比例值na和nf,记为计算值;[0070] S43基于绝对吸附量Q和S2中获取的孔隙结构参数,计算气体的过剩吸附量Qe,具体公式如下:[0071] Qe=Q‑NAPVmp/RT(1)[0072] 式中:Qe——过剩吸附量,cm3/g;[0073] Q——绝对吸附量,cm3/g;[0074] NA——阿伏伽德罗常数,6.02×1023;[0075] Vmp——孔隙体积,cm3;[0076] R——通用气体常数,8.31J·mol‑1·K‑1;[0077] T——温度,K;[0078] 将过剩吸附量Qe与绝对吸附量Q的比值定义为吸附层分子比例的理论值ne,(1‑ne)定义为自由层分子比例的理论值,一一对比计算值和理论值,当误差范围控制在10%以内时,证明原子组定义法则中所取间距大小合理。[0079] 在步骤S5中,所述的气体分子宏观解吸过程模拟包括如下步骤:[0080] S51固定平衡态吸附构型的X和Y不变,在Z方向上增加厚度为10nm的真空层,模拟吸附构型外部的气相空间,为平衡态吸附构型中气体分子的解吸扩散行为提供空间条件,进而形成气体分子宏观解吸模拟的初始构型;[0081] S52参阅图6所示,对解吸模拟的初始构型进行总时长为100皮秒(ps)的动力学模拟,设定系综为NVT,时间间隔为10ps,观察不同时段下孔隙模型中气体分子向低浓度气相空间内解吸扩散的全过程;[0082] S53参阅图7所示,分析煤体纳米孔隙结构中气体分子的解吸特征。[0083] 至此,本领域技术人员应认识到,虽本文已详尽示出和描述了本发明的示例性实施例,但是,在不脱离本发明精神和范围的情况下,仍然可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其他变型或修改。因此,本发明的范围应被理解和认定为覆盖了所有这些其他变型或修改。

专利地区:北京

专利申请日期:2024-03-07

专利公开日期:2024-09-03

专利公告号:CN118116480B


以上信息来自国家知识产权局,如信息有误请联系我方更正!
电话咨询
读内容
搜本页
回顶部