天然气地球科学 >

2025 , Vol. 36 >Issue 11: 2143 - 2153

DOI: https://doi.org/10.11764/j.issn.1672-1926.2025.06.008

非烃气：氢气

干酪根热演化过程中热解生氢分子机理

  • 安思润 , 1, 2 ,
  • 黄亮 , 1, 2 ,
  • 曲自硕 1, 2 ,
  • 杨哲 1, 2 ,
  • 徐侦耀 1, 2 ,
  • 陈秋桔 1, 2 ,
  • 冯鑫霓 1, 2 ,
  • 朱海燕 1, 2
展开
  • 1. 油气藏地质及开发工程全国重点实验室，成都理工大学，四川 成都 610059
  • 2. 成都理工大学能源学院（页岩气现代产业学院），四川 成都 610059
黄亮（1991-） ，男，四川遂宁人，博士，研究员，主要从事非常规油气和新能源开发研究. E-mail： .

安思润（2001-），男，四川江油人，硕士研究生，主要从事天然氢气勘探开发研究. E-mail： .

收稿日期: 2025-03-31

  修回日期: 2025-06-04

  网络出版日期: 2025-07-09

基金资助

国家自然科学基金(52204031)

收起

Molecular mechanisms of pyrolytic hydrogen production during thermal evolution of kerogen

  • Sirun AN , 1, 2 ,
  • Liang HUANG , 1, 2 ,
  • Zishuo QU 1, 2 ,
  • Zhe YANG 1, 2 ,
  • Zhenyao XU 1, 2 ,
  • Qiuju CHEN 1, 2 ,
  • Xinni FENG 1, 2 ,
  • Haiyan ZHU 1, 2
Expand
  • 1. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation，Chengdu University of Technology，Chengdu 610059，China
  • 2. College of Energy （College of Modern Shale Gas Industry），Chengdu University of Technology，Chengdu 610059，China

Received date: 2025-03-31

  Revised date: 2025-06-04

  Online published: 2025-07-09

Supported by

The National Natural Science Foundation of China(52204031)

Fold

摘要

烃源岩干酪根在热成熟演化过程中具有巨大的生氢潜力。目前干酪根热演化过程产氢特征不清，干酪根化学结构与孔隙结构对产氢能力的影响机制不明，水分对干酪根热解产氢的作用机理有待阐明。采用基于ReaxFF力场的反应分子动力学模拟方法，构建干燥和含水条件下干酪根单元结构与基质模型，开展未成熟干酪根升温热演化热解模拟，对比不同成熟度阶段干酪根热解生氢特征。结果表明：①热成熟演化过程中，低成熟干酪根热解产氢能力更强，氢气主要在其高温热解阶段产生；②脂肪族结构上氢原子脱落组合形成氢气是干酪根热演化产氢的主要模式；③水分通过氢源作用以及催化作用促进干酪根脂肪族结构热解产氢；④干酪根热演化产氢受化学结构和孔隙结构共同作用，化学结构影响大于孔隙结构影响。研究结果完善了干酪根热解产氢理论，可为天然氢藏勘探开发提供理论指导。

关键词： 干酪根; 热解; 天然氢气; ReaxFF; 分子模拟

本文引用格式

安思润 , 黄亮 , 曲自硕 , 杨哲 , 徐侦耀 , 陈秋桔 , 冯鑫霓 , 朱海燕 . 干酪根热演化过程中热解生氢分子机理[J]. 天然气地球科学, 2025 , 36(11) : 2143 -2153 . DOI: 10.11764/j.issn.1672-1926.2025.06.008

Abstract

Kerogen of source rock has great potential for hydrogen production during thermal maturation evolution. At present， the characteristics of hydrogen production during the thermal evolution of kerogen are unclear， the influence mechanisms of chemical structure and pore structure of kerogen on the hydrogen production capacity are unknown， and the mechanisms of the role of water on the pyrolytic hydrogen production of kerogen need to be elucidated. In this work， the unit structures and matrix models of kerogen under dry and water-bearing conditions were constructed， and the molecular dynamics simulation method based on the ReaxFF force field was used to conduct the pyrolysis simulation of immature kerogen at elevated temperatures and kerogen at different maturity stages. The results show that： （1） During the thermal maturation process， the lower matured kerogen is more capable of producing hydrogen， and hydrogen is mainly produced in the high-temperature stage； （2） The main mode of hydrogen production during kerogen thermal evolution is the combination of hydrogen atoms from the aliphatic structures； （3） Water promotes the pyrolysis of hydrogen in the aliphatic structure of kerogen through the role of hydrogen source and the catalytic effect； （4） Hydrogen production from thermal evolution of kerogen is affected by both chemical structure and pore structure， with chemical structure having a greater influence than pore structure. The results improve the theory of hydrogen production from pyrolysis of kerogen， which can provide theoretical guidance for the exploration and development of natural hydrogen reservoirs.

Key words： Kerogen; Pyrolysis; Natural hydrogen; ReaxFF; Molecular simulation

0 引言

当前，化石能源占全球能源结构比重约为80%［1-3］，化石能源的过度开发与消耗引发了严重的生态破坏及气候变化问题［4］。为此，中国政府提出了双碳目标战略部署，旨在逐步实现化石能源向可再生能源的转向，构建以新能源为核心的低碳能源系统［5-6］。氢气因其高能量密度、环境友好、安全可靠的特性被认为是具有广泛前景的低碳能源［7-9］。近年来，中国对氢能的需求逐年上涨且长期处于供不应求状态［10］。现有制氢技术面临制氢工艺成本高、制氢效率不稳定等挑战［11］。因此，急需探索行之有效的方法获取更多的氢气资源。
近年来，世界各地陆续发现地质条件下的天然氢藏。天然氢气的形成机理主要分为有机成因与无机成因［12］。目前大部分研究聚焦天然氢气的无机成因，认为蛇纹石化是天然氢气产生的主要方式［13-15］，但对于天然氢气的有机成因研究缺乏。作为烃源岩的主要有机质，干酪根热演化生烃过程中伴随着氢气的生成。选用干酪根代表有机质研究其热演化过程中的生氢特征有助于揭示天然氢气的有机成因机理。
干酪根室内热解实验是研究干酪根热解阶段和行为的直接手段［16-17］。但干酪根热解实验更多关注最终热解产物，在产物生成路径和微观机理解释方面仍存在挑战［18-22］。反应分子动力学（RMD）方法是一种追踪干酪根化学反应路径和阐明干酪根热解机理的有效方法［23-28］。近年来，大量学者基于Reactive Force Field （ReaxFF）力场结合RMD方法开展了干酪根热解分子模拟。LIU等［23］系统研究了干酪根热解过程中的断键机理。QIAN等［29］研究了油页岩干酪根热解过程的化学机理和产物分布。KUANG等［30］划分了绿河油页岩干酪根热解反应进程的主要阶段。ZHANG等［31］系统对比了干酪根直接热解和加氢热解过程的产物区别。尽管学者们已经针对干酪根热解过程和热解机理开展了相关研究，但H2在目前研究中主要作为烃类物质的副产物，其生成规律、影响因素以及分子成因尚未明确。另外，烃源岩长期被水饱和，目前水分在干酪根热解过程中对氢气产生的影响机理仍不清楚。
本文围绕干酪根热解生氢机理，采用基于Reactive Force Field力场的RMD方法，构建干酪根结构单元和基质模型，分别在干燥和含水条件下，开展未成熟干酪根升温热解模拟以及不同成熟度干酪根热解模拟，探究干酪根热演化过程中的热解生氢特征和微观机制。

1 模型和方法

1.1　干酪根热演化过程模拟方法

干酪根热演化过程的模拟通过2种方式（图1）体现。在第一种方式中，基于未成熟阶段的干酪根开展升温热解模拟，不同温度条件下的热解生氢结果差异可反映干酪根热演化进程的影响。在第二种方式中，通过对比4种不同热成熟阶段的干酪根热解生氢行为，可独立反映干酪根热成熟度的影响。本文通过干酪根结构单元体现化学结构差异，解析干酪根热解生氢机制；通过对比干酪根基质和结构单元的热解生氢差异，分析物理孔隙结构的影响；通过对比干燥和含水条件下的热解生氢差异，分析水分的影响。
图1 干酪根热演化过程模拟示意

Fig. 1 Simulation diagram of kerogen thermal evolution process

1.2　干酪根模型构建和验证

II型干酪根主要来源于藻类和浮游生物，在缺氧的深湖或海洋环境中沉积形成。本文研究选取了4种不同成熟度的II型干酪根（图1），用于研究热演化过程的影响。II型干酪根结构单元由UNGERER等［32］构建，其分子式分别为：C252H294N6O24S3（II-A）、C234H263N5O14S2（II-B）、C242H219N5O13S2（II-C）、C175H102N4O9S2 （II-D）。基于干酪根结构单元，结合笔者先前提出的模型构建流程［33-35］，通过几何优化、高温弛豫、退火模拟程序，构建上述4种成熟度干酪根的基质模型（图1）。由于实际烃源岩中干酪根是一种分子量高、组成成分复杂且没有固定化学结构的三维大分子混合物，无法获得其确定的分子结构，目前主要采用平均分子结构来近似模拟其化学结构，借此反映干酪根的典型性质。目前学者普遍采用平均化模型开展对干酪根热解行为的研究［36-39］，因此，笔者也选用常用的4种不同成熟度的干酪根平均化模型开展后续研究。
干酪根基质模型的物理性质如表1所示。干酪根基质模型密度在1.12~1.18 g/cm³范围内，与报道的实验结果接近［40］。干酪根基质模型孔隙度介于11.78%~18.04%之间，符合报道的有机质样品孔隙度范围［41］（1.95%~26.03%）。当干酪根由未成熟（II-A）逐渐演化至过成熟阶段（II-D），干酪根结构中的H/C原子比例逐渐减小。
表1 干酪根基质模型物理性质

Table 1 Physical properties of kerogen matrix model

干酪根类型 密度/（g/cm3） H/C 孔隙度/%
II-A 1.14 1.16 11.78
II-B 1.12 1.12 12.55
II-C 1.18 0.91 15.54
II-D 1.18 0.58 18.04
本文主要研究水分子对干酪根热演化生氢机制的影响，而水分子个数主要影响最后产氢数量，不会影响机理性结论，因此本文并未考虑实际烃源岩中水分含量波动的情况。对于水分子的处理，首先利用巨正则蒙特卡洛方法在干酪根结构单元模型和基质模型中分别吸附20和60个水分子，其最终水分含量为2%和5%，然后利用分子动力学方法开展水分平衡模拟，构建含水干酪根结构单元模型和基质模型。

1.3　模拟与计算细节

本研究中的模型构建和模拟基于Materials Studio软件完成。模型构建选用COMPASS III力场描述原子间相互作用，范德华作用基于atom-based 方法计算，截断半径设置为12.5 Å，静电作用基于Ewald方法计算。动力学模拟中温度由Nose方法控制，压力由Berendsen方法控制。
干酪根热解模拟基于ReaxFF 6.0力场进行。每个热解温度点的平衡时间为50 ps，生产时间为350 ps，模拟步长为0.5 fs。热解模拟的最后100帧用于定量分析热解结果。热解模拟之前，干酪根单元模型与基质模型在500 K（除非另有说明，本文统一采用 K为温度单位，其与 ℃的转换关系定义为：T（K）=T（℃）+273.15）下进行100 ps热平衡处理。通过计算原子间键长判断化学键的断裂或新产物的形成，从而定量分析H2的产生量。键长介于0.75~1.15 Å范围视为有键连接，超出此范围则认为化学键断裂，形成热解产物［34］。由于热解过程中原子电荷和编号保持不变，本研究结合H2的电荷信息和原子编号，追踪氢原子的来源，以阐明H2的形成机理。针对干酪根基质模型，采用Connolly提出的探针法［42］表征其在热解过程中的孔隙结构变化，计算精度选择Medium，截断半径为1.3 Å。该方法的主要思路是将固定半径为1.3 Å的氦气作为探针分子随机插入干酪根基质，让探针分子在干酪根原子表面上滚动并追踪其与表面接触点的运动轨迹以形成Connolly surface，从而计算出孔隙度以及比表面积。孔隙度的计算公式如式（1）所示：
ϕ = V 1 V 2 × 100 %
式中：ϕ为孔隙度，%；V1 为晶胞内部的自由孔隙体积，Å3；V2 为晶胞总体积，Å3。模型比表面积的计算公式如式（2）所示：
S 1 = S 2 ρ × V 2 × 10 4
式中：S 1为比表面积，m2/g；S 2为每个晶胞的孔隙表面积，Å2；ρ为晶胞内的干酪根密度，g/cm3。
在干酪根的地质热演化过程中时间尺度为百万年，在分子模拟中无法复现如此长时间的热演化过程。因此，本研究使用目前学者普遍采用的温度加速技术［23］，即在超高温下模拟从而减少模拟时间，虽无法克服与百万年时间尺度的差异，但结合Arrhenius的方法能尽量接近实际热解工程的时间尺度［43］。实验温度和模拟温度之间的对应关系是通过等效时率乘积［式（3）］来确定的，以确保模拟和实验中发生的事件数量相同。通过以下公式建立实验与模拟温度之间的关系，可以实现模拟时间尺度变化在实验中的准确反映。
f e x p t e x p = f s i n t s i n
式中： f e x p是实验频率，Hz； t e x p是实验时间，s； f s i n是模拟频率，Hz； t s i n是模拟时间，s。再结合Arrhenius方程得出式（4），它将实验温度和时间与模拟温度和时间联系起来：
T s i n = - E a k B l o g ( t e x p t s i n ) - E a k B T e x p - 1
式中： E a是活化势垒，eV； f是事件频率，Hz； A表示尝试频率，Hz； T表示加热温度，K； k B表示玻尔兹曼常数。

2 干酪根单元热演化过程中生氢机制

2.1　干酪根单元升温热解过程中生氢机制

图2为II-A型干酪根单元在升温热解过程的H2与H自由基生成量。II-A型干酪根的热解过程可分为2个阶段：低温热解阶段（800~1 200 K）和高温热解阶段（1 600~2 400 K）。低温阶段几乎不产生H2和H自由基，高温阶段H2和H自由基随热解温度增加呈持续上升趋势。低温阶段干酪根结构中C—H键的高键级和H自由基的强活性制约了H2和H自由基的产生。高温阶段，干酪根中大分子结构片段逐步分解，持续释放出羟基和H自由基，并通过自由基重组生成H2。当II-A型干酪根单元在含水条件下升温热解时，水分会对干酪根的热解行为产生较大影响。800 K时，含水和干燥条件下干酪根热解均未产生H2和H自由基。当温度升至1 200 K，含水条件下干酪根产生少量H2，但未产生H自由基。在1 200~2 400 K范围内，含水条件下干酪根产生的H2和H自由基数量显著大于干燥条件，水分有效促进了干酪根热解产氢。
图2 II-A型干酪根单元升温热解过程H2与H自由基生成量

Fig.2 The production of H2 and H free radicals in the pyrolysis process of type II-A kerogen unit

为探究有机成因的天然氢气是否具有成为天然氢藏的潜力，运用原子编号追踪的方法，对Ⅱ-A型干酪根在升温热演化过程中产生的每一个H₂分子逐一编号，并持续追踪其中H原子的去向，从而量化氢气的“留存比例”（即某温度区间起始时生成的H₂，在更高温度下仍保持为H₂的百分比）。当温度从1 200 K升至1 600 K时，含水条件下干酪根在1 200 K热解产生的4个H2中仅有1个H原子留存至1 600 K，此时的H2的留存比例仅为12.5%，说明低温下产生的H2更多又作为反应物重新参与反应。当温度从1 600 K升至2 000 K时，含水情况下1 600 K产生的38个H2中有30个H原子留存至2 000 K，H2留存比例达到了39.5%，而干燥情况下1 600 K产生的13个H2中也有10个H原子留存至2 000 K，留存比例同样达到了38.5%。同样的，温度从2 000 K升至2 400 K时，含水与干燥条件下的H2留存比例也分别达到了38.8%和41.6%，说明在高温热解阶段的氢气留存比例较高，有成为天然氢藏的潜力。
为阐释干酪根热解过程生氢机理，量化分析干酪根单元升温热解产生H2的H原子来源（图3）。1 200 K时，干燥条件下未出现H2分子，含水时产生的H数量较少，且全部来自于脂肪族结构，此时水分没有起到供H的氢源作用。随热解温度升高，含水和干燥条件下源自脂肪族和芳香环的H数量均显著增加。此外，含水条件下源自水分的H数量随温度升高而持续增长，表明水分在高温热解阶段具有氢源作用。含水条件下脂肪族H数量相比干燥条件的增加幅度始终大于水分的供H数量，表明水分通过催化作用有效促进了脂肪族结构的热解生氢。对比不同基团，水分和热解温度对脂肪族和芳香环的生氢能力具有显著影响，对吡咯基团生氢能力的影响较小。随后，通过追踪水分子与干酪根模型的变化，从微观的角度阐述水分子的氢源作用以及催化作用。由于水分子在1 200 K之前热解程度较低，氢源作用以及催化作用均不明显，因此选取1 600 K观察水分子的微观作用机制。水分子的氢源作用表现在两个方面，水分子热解产生的H自由基既可以与其他水分子产生的自由基结合，如H-593与H-573、H-583与H-540等形成H2，又可以和干酪根热解出的H自由基组合形成H2。由于1 600 K时干酪根主要表现为脂肪链断裂，因此水分子提供的H自由基大多与脂肪链产生的H自由基结合，如H-614（水）与H-397、H-609（水）与H-233、H-638（水）与H-400等结合形成H2，但也存在部分水分子的H自由基与芳香环上断裂的H自由基组合形成H2，如H-605（水）与H-153。上文通过对比干燥与含水情况下H2分子数量发现水分子对干酪根热解产氢行为具有促进作用，随后观察干燥与含水情况下的产物发现，含水情况下产生的额外H2中有许多是完全由干酪根脂肪链提供的H自由基组合形成，如H-379和H-160、H-112和H-69等。说明水分子的存在会催化干酪根脂肪链的断裂，提高干酪根热解产氢能力。
图3 II-A型干酪根单元升温热解产生H2的H原子来源

Fig.3 H atom source of H2 produced by pyrolysis of type II-A kerogen unit

图4展示了II-A型干酪根单元升温热解产生H2的H原子来源占比。在干燥条件下，生氢潜力随温度变化较小，绝大部分H2来源于脂肪族，仅小部分来源于芳香环和吡咯基团，这表明脂肪族结构是干燥条件下干酪根热解的主要氢源。随着温度升高，水分使得脂肪族的产氢占比呈现先下降后稳定的趋势，芳香环的产氢占比略有增加。此外，含水条件下水分自身的产氢能力介于12%~17%之间，表明水分在干酪根热解过程中不仅可以作为催化剂促进干酪根结构生氢，也能够以氢源作用供H参与生氢。
图4 II-A型干酪根单元升温热解产生H2的H原子来源占比

Fig.4 The ratio of H atom source of H2 produced by II-A type kerogen unit heating pyrolysis

干酪根中的氢原子以不同的电荷赋存于不同的结构基团中。本研究结合原子编号与电荷对热解生氢的H来源进行追踪，进而识别热解过程中H2的形成路径。图5为2 000 K时干酪根单元热解生氢的H2和H分布模式图。干酪根热解过程具有多种H2形成路径，其中最主要的方式是从脂肪族结构分离的H自由基组合形成氢气（图5①）。其次是电荷为0.053的H自由基与电荷为0.41的H自由基组合形成氢气（图5②）。对于II-A型干酪根，其电荷为0.41的H原子分别是来自脂肪链羧基的H-177和H-380、脂肪链羧基的H-179和H-416以及来自芳香环羟基的H-570和H-448，但仅有羟基的H-426、H-448以及羧基的H-177在II-A型干酪根升温过程中发生断裂，因此方式②形成的H2大部分是干酪根脂肪族H自由基与水分子产生的H自由基的结合。方式③则是不同水分子产生的H原子自由基在重新组分的方式形成H2。方式④形成的H2是干酪根脂肪族H自由基与芳香环上H自由基的组合。此外，还有极少部分H2以其他的方式组合，如2个芳香环上H自由基组合、芳香环上H自由基与水分H自由基组合、吡咯上H自由基与其他H自由基组合。
图5 II-A型干酪根单元热解生氢机理图（2 000 K）

Fig.5 Mechanism diagram of hydrogen generation from type II-A kerogen pyrolysis （2 000 K）

2.2　不同成熟度干酪根单元的生氢机制

干酪根的热演化成熟度变化会显著影响其化学结构和热解行为。图6展示了4种不同成熟度干酪根结构单元中氢原子的赋存形态分布。4种干酪根中氢原子均主要赋存于脂肪族碳原子上，随热演化成熟度增加，脂肪族碳原子上氢原子数量逐渐减小。此外，4种干酪根的芳香环上也赋存大量氢原子，随热演化成熟度增加，脂肪族碳原子上氢原子与芳香环碳原子上氢原子的数量差距减小。另外，干酪根脂肪族的羟基和羧基，芳香环的羟基以及吡咯基团上也赋存少量的氢原子，随热演化成熟度增加，这部分氢原子数量呈减少趋势。
图6 不同成熟度化学结构统计（结构单元）

Fig.6 Chemical structure statistics of different maturity （structural unit）

图7展示了4种干酪根单元热解生成的H2与H自由基数量。干燥与含水条件下，不同成熟度干酪根的热解产H2能力趋势相同，均为II-A>II-B>II-D>II-C，水分的催化作用及氢源作用使得不同成熟度干酪根的产H2能力差距增大。由于1 600 K处于高温热解的初级阶段，H自由基产生较少且H自由基的活性较高，易于产生二次反应。因此，不同成熟度的干酪根产生H自由基数量的差距较小。
图7 不同成熟度干酪根单元热解生成的H2与H自由基数量（1 600 K）

Fig.7 The number of H2 and H radicals produced by pyrolysis of kerogen units with different maturity（1 600 K）

图8为不同成熟度干酪根单元热解生成H2的H原子来源。在干燥条件下，热解生成H2的主要来源仍然是脂肪族上氢原子，随着成熟度的增加，干酪根中脂肪族结构数量减小，其提供的氢原子数量减小，但II-D型干酪根中脂肪族结构为H2提供的氢原子数量略有上升。由于1 600 K下热解程度较低，因此来自芳香环和吡咯上的氢原子数量变化不明显。水分对不同成熟度干酪根的热解生H2均有促进作用，其对脂肪族结构的影响大于芳香环结构。
图8 不同成熟度干酪根单元热解生成H2的H原子来源

Fig.8 The H atom source of H2 generated by pyrolysis of kerogen units with different maturity

图9对比了不同成熟度干酪根单元主要结构的产氢效率。产氢效率为干酪根结构基团热解生成的H占该结构基团上H原子总数的比例。相比低成熟II-B型和未成熟II-A型干酪根，过成熟II-D型和成熟II-C型干酪根的芳香环产氢效率更大。随着干酪根热成熟度增加，脂肪族的产氢效率逐渐降低，但过成熟II-D型干酪根脂肪族上68%的H原子参与热解生成H2，其热解效率远大于其他成熟度干酪根。
图9 不同成熟度干酪根单元主要结构的产氢效率（1 600 K）

Fig.9 Hydrogen production efficiency of the main structure of kerogen units with different maturity（1 600 K）

3 干酪根基质热演化过程中生氢特征

3.1　干酪根基质升温热解过程中生氢特征

干酪根热解产氢不仅受化学结构的影响，也受孔隙结构的影响。通过对比干酪根结构单元和基质的生氢差异，可分析干酪根孔隙结构对热解生氢的影响。
图10为II-A型干酪根含水基质升温热解中孔隙度与孔隙结构的演化。当温度由500 K升高至800 K，干酪根基质孔隙度由11.78%增加至13.19%，比表面积从1 778.06 m2/g下降至1 523.35 m2/g。这是因为干酪根在低温热解阶段优先形成大分子碎片，形成较大的基质孔隙，从而比表面积减小。当热解温度上升至1 200 K，热解产生的大分子转化为小分子结构，比表面积增大至1 948.72 m2/g，由于小分子填充孔隙空间，形成更多的小孔隙，孔隙度下降至10.37%。当温度继续升高至1 600 K，水分促进干酪根热解，大量热解产物的填充，使得孔隙度大幅度降低至4.97%，比表面积也升高至2 531.29 m2/g。
图10 II-A型干酪根基质升温热解中孔隙度与孔隙结构演化

Fig.10 Evolution of porosity and pore structure during pyrolysis of type II-A kerogen matrix

图11为不同温度下II-A型干酪根单元和基质的热解产氢数量对比。800 K时，干酪根基质与单元均未产生H2。对比1 200 K与1 600 K，1 600 K时在水分的促进作用下，干酪根产氢数量大幅度上升，此时生产大量的小分子热解产物，不仅大幅度增加了模型的比表面积，同时也作为填充物减小了模型的孔隙度，使模型中出现孤立孔隙以及死端孔隙，使得模型的孔径以及连通性均降低，由于H自由基本身具有较高活性，在此基础上被基质俘获重新参与反应的概率进一步提高，自由组合为H2的概率进而减小，表现为基质模型产氢数量相较于单体结构的倍数小于1 200 K时的3倍。
图11 不同温度下单体结构与基质模型的产氢数量

Fig.11 Hydrogen production quantity of monomer structure and matrix model at different temperatures

3.2　不同成熟度干酪根基质的生氢特征

图12为不同成熟度干酪根含水基质的孔隙度与孔隙结构差异。1 200 K下，随热成熟度增加，干酪根孔隙度呈增加趋势，比表面积呈现减小趋势。成熟度较高的干酪根芳香度较高，而1 200 K下芳香结构热解程度较小。因此，相比成熟度更低的干酪根，成熟度高的干酪根孔径、连通性以及孔隙度均更高，比表面积更小。
图12 不同成熟度干酪根基质孔隙度与孔隙结构差异（1 200 K）

Fig.12 Differences in porosity and pore structure of kerogen matrix with different maturity（1 200 K）

图13为不同成熟度干酪根单元和基质的热解产氢数量对比。在1 200 K，相比成熟度较高的干酪根，成熟度较低的干酪根基质和结构单元热解产生的H2数量均更大。然而，相比成熟度较低的干酪根，成熟度较高的干酪根基质模型相比结构单元的产H2数量倍数更大。这是因为高成熟度的干酪根孔隙度更大，孔隙直径更大，连通性更好，比表面积也更小，H自由基组合形成H2的概率高于低成熟度的干酪根。
图13 不同成熟度干酪根单元和基质的热解产氢数量对比（1 200 K）

Fig.13 Comparison of the amount of hydrogen produced by pyrolysis of kerogen units and substrates with different maturity（1 200 K ）

本文首先进行干酪根单体结构的热解模拟，研究不同干酪根单体化学结构对于热解产氢行为的影响机理。随着成熟度的增加，干酪根脂肪链含量降低，芳香度增加，其产氢能力降低但产氢效率增加。随后，将干酪根单体模型扩大至基质模型。基质模型具有孔隙结构，影响了H自由基组合形成H2的概率，因此基质模型中H2数量不能达到单体结构对应的倍数。总得来说随热成熟度增加，干酪根脂肪族单元含量减小，抑制产氢能力，而孔隙度和孔径增加，提高产氢能力。两者共同影响下，高成熟度的干酪根热解产氢能力相比低成熟度的干酪根更低，表明热演化过程中，干酪根化学结构对热解产氢的影响大于孔隙结构。因此，人工改造烃源岩储层时，应更多的选择脂肪族含量较高且孔隙结构更好的储层。

3.3　工程应用启示

本文研究通过分子模拟系统研究了干酪根热演化过程中的热解生氢机理，揭示了干酪根结构、水分作用、热演化进程和热成熟度对热解行为的影响。本节将探讨本研究结果对于寻找天然气藏的启示和为人为改造储层提高产氢效率提供相应指导。
水分子在干酪根热解产氢过程中具有显著的影响。模拟结果显示，水分子通过氢源作用以及促进作用显著提高干酪根产氢能力。因此，在寻找天然氢藏时优先选择被地层水长期饱和的储层。同时，干酪根的化学结构是影响干酪根热解产氢行为的关键因素。模拟结果发现，脂肪族化合物是干酪根热解产氢的主要氢源，且其热解温度相对较低。因此，在勘探天然氢气藏时，应优选含脂肪族化合物丰富的干酪根层。结合地质资料和模拟分析，可以提前筛选出高潜力区域，提高勘探天然氢藏效率。此外，成熟度越高的储层，其孔隙度增加、连通性变好、比表面积降低，H自由基更容易形成H2，且自由孔隙体积增大使其成为天然氢藏的可能性更高。因此，在勘探过程中，应综合评估烃源岩的孔隙特征，优选孔隙结构优良的储层。
同样，由于水分子提高热解产氢能力的显著作用，因此，在储层改造中，应通过注水或保持适宜的含水量，优化地下水分环境，以提高产氢效率。热演化进程是决定干酪根产氢量的关键因素。研究表明，当温度高于1 600 K（实际温度788 K）后，热演化进入高温热解阶段，产氢量大幅度提升。因此，人工改造储层将温度控制在高温热解阶段内可以显著提高产氢能力。
由于模拟工作与现场实际存在一定的差异，因此本研究的适用范围存在一定的局限性。首先，本文通过构建平均化干酪根模型以研究其热解产氢行为，但实际烃源岩储层中干酪根与无机矿物共存，本文构建的模型未考虑无机矿物对干酪根热解产氢行为的潜在影响。目前，关于无机矿物对烃源岩热解产氢行为的影响研究缺乏，部分学者探索了不同类型无机矿物对烃源岩热解产烃行为的影响，发现烃源岩中主要无机矿物类型对干酪根热解产烃的促进作用较为有限［44］。具体而言，绿泥石与伊利石对于本文研究的II型干酪根热解产烃的促进作用较小，高岭石只有在与干酪根的质量比达到60∶1之后才会呈现一定的促进现象。需要说明，在干酪根热解产物组成中，非烃类产物总量远低于烃类产物的总量，而非烃类产物中H2占比非常有限［36］。以上研究表明烃源岩中主要无机矿物类型对干酪根热解产氢的影响较为有限。因此，虽然本研究建立的平均化有机质模型未考虑无机矿物的影响，但该模型在应用于有机质含量较高的烃源岩储层时仍具有一定的适用性。
其次，虽然实际烃源岩储层中水分含量变化幅度较大，但烃源岩干酪根呈混合润湿特征，其内部水分含量普遍较低，因此本研究未考虑高水分含量对干酪根热解过程的影响。此外，ZHAO等［45］在干酪根基质模型中吸附2 088个水分子进行高水分含量下的热解模拟，研究表明高水分含量主要影响干酪根最终产氢量，而对氢气形成途径无明显影响。因此，高水分含量对干酪根热解产氢机理的影响较小，本文研究得到的干酪根在含水条件下的热解产氢机理适用于描述实际烃源岩储层干酪根的热解产氢行为。
最后，虽然本文采用了温度加速技术，一定程度上弥补了实验时间尺度与现场实际改造储层时间尺度的差异，但与地质热演化过程的时间尺度仍具有一定差异，后续将结合跨尺度模拟来进一步弥补这些差异。

4 结论

（1）II⁃A型干酪根升温热演化模拟的低温热解阶段（≤1 200 K）基本不产氢，高温热解阶段（1 600~2 400 K）产氢量大幅提升。不同热成熟阶段中，低成熟阶段干酪根热解产氢能力更强，过成熟阶段干酪根热解产氢效率更高。
（2）干酪根热演化过程产生的氢气主要来自干酪根的脂肪族结构，其在氢源中的占比可超过85%。脂肪族结构上氢原子脱落并互相组合形成氢气是干酪根热解产氢的主要模式。
（3）水分可促进干酪根热解产氢，未成熟干酪根及其在高温热解阶段时，水分促进产氢效果更好。热解中水分通过氢源作用以及催化作用，促进干酪根脂肪族结构热解产氢。
（4）干酪根热演化产氢受化学结构和孔隙结构共同作用。干酪根热演化过程中，脂肪族结构减少可抑制产氢，孔隙度增大可促进产氢，热解产物填充孔隙可抑制产氢，化学结构对生氢能力的影响大于孔隙结构。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
周守为， 朱军龙， 李清平， 等. 科学稳妥实现“双碳”目标， 积极推进能源强国建设［J］. 天然气工业， 2022， 42（12）： 1-11.

ZHOU S W， ZHU J L， LI Q P， et al. Scientifically and prudently achieving the goals of peak carbon emissions and carbon neutrality， actively promoting the construction of an energy power［J］. Natural Gas Industry， 2022， 42（12）： 1-11.

[2]
何东博， 贾成业， 位云生， 等. 世界天然气产业形势与发展趋势［J］. 天然气工业， 2022， 42（11）： 1-12.

HE D B， JIA C Y， WEI Y S， et al. Current situation and development trend of world natural gas industry［J］. Natural Gas Industry， 2022， 42（11）： 1-12.

[3]
韩双彪， 唐致远， 杨春龙， 等. 天然气中氢气成因及能源意义［J］. 天然气地球科学， 2021， 32（9）： 1270-1284.

HAN S B， TANG Z Y， YANG C L， et al. Genesis and energy significance of hydrogen in natural gas［J］. Natural Gas Geoscience， 2021， 32（9）： 1270-1284.

[4]
王强， 党牛， 蒋子龙， 等. 碳中和背景下能源转型研究综述与展望［J］. 地理学报， 2025， 80（3）： 586-604.

WANG Q， DANG N， JIANG Z L， et al. Review and prospects of energy transition research in the context of carbon neutrality［J］. Acta Geoscientica Sinica， 2025， 80（3）： 586-604.

[5]
周莹， 叶嘉乐， 于姗， 等. “双碳”目标下天然气分布式能源发展现状及机遇［J］. 天然气工业， 2024， 44（2）： 23-29.

ZHOU Y， YE J L， YU S， et al. Development status and opportunities of natural gas distributed energy resource under the goal of “dual carbon”［J］. Natural Gas Industry， 2024， 44（2）： 23-29.

[6]
刘金华， 葛政俊， 李晓凤. 苏北盆地CO2气藏成藏规律及富氦成因［J］. 天然气地球科学， 2023， 34（3）： 477-485.

LIU J H，GE Z J，LI X F.The accumulation law and helium rich genesis of CO2 gas reservoir in Subei Basin［J］.Natural Gas Geoscience， 2023， 34（3）： 477-485.

[7]
HANLEY E S， DEANE J P， GALLACHÓIR B P Ó. The role of hydrogen in low carbon energy futures：A review of existing perspectives［J］. Renewable and Sustainable Energy Reviews， 2018， 82： 3027-3045.

[8]
TARKOWSKI R， CZAPOWSKI G. Salt domes in Poland–Potential sites for hydrogen storage in caverns［J］. International Journal of Hydrogen Energy， 2018， 43（46）： 21414-21427.

[9]
TARKOWSKI R. Underground hydrogen storage： Characteristics and prospects［J］. Renewable and Sustainable Energy Reviews， 2019， 105： 86-94.

[10]
WANG L， JING Z J， CHEN X， et al. The origin and occurrence of natural hydrogen［J］. Energies， 2023， 16（5）： 2400.

[11]
张正， 宋凌珺. 电解水制氢技术： 进展、挑战与未来展望［J］. 工程科学学报，2025，47（2）：282-295.

ZHANG Z， SONG L J. Hydrogen production by water electrolysis： Advances， challenges and future prospects［J］. Chinese Journal of Engineering， 2025， 47（2）： 282-295.

[12]
窦立荣， 刘化清， 李博， 等. 全球天然氢气勘探开发利用进展及中国的勘探前景［J］. 岩性油气藏， 2024， 36（2）： 1-14.

DOU L R， LIU H Q， LI B. Global natural hydrogen exploration and development situation and prospects in China［J］.Lithologic Reservoirs， 2024， 36（2）： 1-14.

[13]
宋晗. 低温蛇纹岩化过程中Fe（OH）2相转变产H2的机制研究［D］. 广东：华南理工大学， 2022.

SONG H.Mechanism of Fe（OH）2 Phase Transformation Toproduce H2 During Low Temperature Serpentinization［D］.Guangdong： South China University of Technology， 2022.

[14]
LEFEUVRE N，TRUCHE L，DONZE F V， et al. Native H2 exploration in the western Pyrenean foothills［J］.Geochemistry，Geophysics，Geosystems，2021，22（8）：e2021GC009917.

[15]
LIU Q， WU X， HUANG X， et al. Integrated geochemical identification of natural hydrogen sources［J］. Science Bulletin， 2024， 69（19）： 2993-2996.

[16]
WU Y， JI L， HE C， et al. The effects of pressure and hydrocarbon expulsion on hydrocarbon generation during hydrous pyrolysis of type-I kerogen in source rock［J］. Journal of Natural Gas Science and Engineering， 2016， 34： 1215-1224.

[17]
LEWAN M D， ROY S. Role of water in hydrocarbon generation from Type-I kerogen in Mahogany oil shale of the Green River Formation［J］.Organic Geochemistry，2011，42（1）：31-41.

[18]
LANDAIS P，MICHELS R，Elie M.Are time and temperature the only constraints to the simulation of organic matter maturation［J］.Organic Geochemistry，1994，22（3-5）：617-630.

[19]
FANG C， LI S， MA G， et al. Reaction mechanism and kinetics of pressurized pyrolysis of Chinese oil shale in the presence of water［J］. Petroleum Science， 2012， 9： 532-534.

[20]
LI Q， HAN X， LIU Q， et al. Thermal decomposition of Huadian oil shale. Part 1. Critical organic intermediates［J］. Fuel， 2014， 121： 109-116.

[21]
HUANG Y R， HAN X X， JIANG X M. Comparison of fast pyrolysis characteristics of Huadian oil shales from different mines using Curie-point pyrolysis-GC/MS［J］. Fuel Processing Technology， 2014， 128： 456-460.

[22]
YAN L， BAI Y， ZHAO R， et al. Correlation between coal structure and release of the two organic compounds during pyrolysis［J］. Fuel， 2015， 145： 12-17.

[23]
LIU X， ZHAN J H， LAI D， et al. Initial pyrolysis mechanism of oil shale kerogen with reactive molecular dynamics simulation［J］. Energy & Fuels， 2015， 29（5）： 2987-2997.

[24]
QIAN Y， ZHAN J H， LAI D， et al. Primary understanding of non-isothermal pyrolysis behavior for oil shale kerogen using reactive molecular dynamics simulation［J］.International Journal of Hydrogen Energy， 2016， 41（28）： 12093-12100.

[25]
ZHONG Q，MAO Q，XIAO J，et al.ReaxFF simulations of petroleum coke sulfur removal mechanisms during pyrolysis and combustion［J］.Combustion and Flame，2018，198：146-157.

[26]
ZHANG Z， CHAI J， ZHANG H， et al. Structural model of Longkou oil shale kerogen and the evolution process under steam pyrolysis based on ReaxFF molecular dynamics simulation［J］. Energy Sources， Part A： Recovery， Utilization， and Environmental Effects， 2021， 43（2）： 252-265.

[27]
HAN Q， LI M， LIU X， et al. A maturation scale for molecular simulation of kerogen thermal degradation［J］.Organic Geochemistry， 2023， 175： 104507.

[28]
SHI K， CHEN J， PANG X， et al. Average molecular structure model of shale kerogen： Experimental characterization， structural reconstruction，and pyrolysis analysis［J］.Fuel，2024， 355： 129474.

[29]
QIAN Y， ZHAN J H， LAI D， et al. Primary understanding of non-isothermal pyrolysis behavior for oil shale kerogen using reactive molecular dynamics simulation［J］. International Journal of Hydrogen Energy， 2016， 41（28）： 12093-12100.

[30]
KUANG W， LU M， YEBOAH I， et al. A comprehensive kinetics study on non-isothermal pyrolysis of kerogen from Green River oil shale［J］.Chemical Engineering Journal， 2019，377： 120271-120275.

[31]
ZHANG Z， GUO， ZHANG H， et al. Comparing product distribution and desulfurization during direct pyrolysis and hydropyrolysis of Longkou oil shale kerogen using reactive MD simulations［J］. International Journal of Hydrogen Energy， 2019，44（47）： 25335-25346.

[32]
UNGERER P， COLLELL J， YIANNOURAKOU M. Molecular modeling of the volumetric and thermodynamic properties of kerogen： Influence of organic type and maturity［J］.Energy & Fuels， 2015，29（1）： 91-105.

[33]
CHEN Q，HUANG L，YANG Q，et al.Molecular insights into dual competitive modes of CH4/CO2 in shale nanocomposites：Implications for CO2 sequestration and enhanced gas recovery in deep shale reservoirs［J］.Journal of Molecular Liquids，2024，415： 126359.

[34]
冯鑫霓， 黄亮， 徐侦耀， 等. 不同类型干酪根热解特性及微观机理分子动力学模拟［J/OL］. 成都理工大学学报（自然科学版）， 1-14［2025-03-30］. http：//kns.cnki.net/kcms/detail/51.1634.N.20250311.1059.002.html.

FENG X N， HUANG L， XU Z Y， et al. Molecular dynamics simulation of pyrolysis of different types of kerogen： Characteristics and microscopic mechanisms［J/OL］. Journal of Chengdu University of Technology （Science & Technology Edition）， 1-14［2025-03-30］. http：//kns.cnki.net/kcms/detail/51.1634.N.20250311.1059.002.html.

[35]
黄亮. 基于分子模拟的页岩气多组分竞争吸附机理研究［D］. 北京：中国石油大学（北京）， 2020.

HUANG L.Molecular Simulation Study on Competitive Adsor-ption Mechanism of Multi-components in Shale Gas Reservoir［D］.Beijing：China University of Petroleum（Beijing），2020.

[36]
MA G， CHEN J， YUE C， et al. Reaction molecular dynamics simulation of kerogen hydrogenation catalysis in marine deep shale［J］. Fuel， 2025， 380： 133245.

[37]
CHEN Y， WANG Z， LI B， et al. ReaxFF molecular dynamics study on the microscopic mechanism for kerogen pyrolysis［J］.Langmuir， 2023， 39（50）： 18581-18593.

[38]
LIANG Z， JIANG Z， XUE Z， et al. Evolution of water vapor adsorption behavior on kerogen： Pyrolysis and molecular investigation［J］. Fuel， 2025， 395： 135105.

[39]
PAN S， ZHANG Y， BAI J， et al. Pyrolysis mechanism of kerogen： Model construction and multi-scale molecular simulations［J］.Journal of Analytical and Applied Pyrolysis，2024，183：106837.

[40]
OKIONGBO K S， APLIN A C， LARTER S R. Changes in type II kerogen density as a function of maturity： Evidence from the Kimmeridge Clay Formation［J］.Energy & Fuels，2005，19（6）： 2495-2499.

[41]
LOUCKS R G， REED R M， RUPPEL S C， et al. Morphology， genesis， and distribution of nanometer-scale pores in siliceous mudstones of the Mississippian Barnett Shale［J］.Journal of Sedimentary Research，2009，79（12）：848-861.

[42]
LIU S Q， LIU H， GUO W， et al. Advanced exergy analysis of deep UCG to hydrogen production［J］.Journal of China Coal Society， 2024，49（2）； 1138-1147.

[43]
DE T C， SUAREZ-MARTINEZ I， VALLEJOS-Burgos F， et al. Structural prediction of graphitization and porosity in carbide-derived carbons［J］. Carbon， 2017， 119： 1-9.

[44]
李二庭， 靳军， 高秀伟， 等. 黏土矿物对不同类型烃源岩热解参数影响［J］. 天然气地球科学， 2021， 32（4）： 510-517.

LI E T， JIN J， GAO X W， et al. Effect of clay minerals on pyrolysis parameters of different types of source rocks［J］. Natural Gas Geoscience， 2021， 32（4）： 510-517.

[45]
ZHAO F M， LI B， ZANG L， et al. The mechanism of superheated steam affecting the quality of in-situ pyrolysates of oil shale kerogen： Part A-saturation of pyrolytic organics［J］. Fuel， 2022， 323： 124331.

Options
文章导航

/

〈 〉

陇ICP备2021001824号-7　 甘公网安备 62010202000678号
版权所有 © 《天然气地球科学》编辑部
地址：甘肃省兰州市天水中路8号　
邮编：730000　
电话：(0931)8277790　
E-mail：geogas@lzb.ac.cn
总访问量 　 今日访问 　 在线人数
技术支持：北京玛格泰克科技发展有限公司