《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》：Contrasting mechanical evolution of bakken shale under hydrous and anhydrous pyrolysis: A nanoindentation-based assessment
编辑推荐：
Bakken组（Bakken Formation）是最大的非常规页岩油储层之一，其中水力压裂（hydraulic fracturing）与水平钻井（horizontal drilling）是提高采收率的关键方法。深入理解其力学性质对于这些工艺的设计至关重要，尤
Bakken组（Bakken Formation）是最大的非常规页岩油储层之一，其中水力压裂（hydraulic fracturing）与水平钻井（horizontal drilling）是提高采收率的关键方法。深入理解其力学性质对于这些工艺的设计至关重要，尤其是因为页岩中的干酪根（kerogen）在热演化过程中会发生化学变化，从而改变其力学行为。该研究主要聚焦于评估富含II型干酪根（Type II kerogen）的Bakken页岩在无水热解（anhydrous pyrolysis, AHP）和含水热解（hydrous pyrolysis, HP）受控热成熟过程中的力学行为如何演化。未成熟样品分别在300、325、350、365、400和450 °C下进行热解。为系统研究页岩样品的力学特征，研究人员在样品表面共进行了775次纳米压痕测试。
结果表明，在无水热解（AHP）过程中，样品杨氏模量（Young's modulus）最初升高至34 GPa以上，随后出现较大幅度下降，并在350–450 °C范围内开始在20–25 GPa之间波动。另一方面，AHP过程中的硬度（hardness）在300 °C升高至接近0.77 GPa，之后稳步下降至365 °C时约0.25 GPa；随后在400 °C回升至0.39 GPa，最终在最后一个温度阶段下降至0.32 GPa。在含水热解（HP）过程中，这些变化趋势明显不同，反映出水存在对干酪根成熟的影响。HP过程中硬度值持续降低，在450 °C时降至约0.15 GPa。相比之下，杨氏模量在300 °C时先降至21.6 GPa，随后在325 °C升至接近25 GPa；之后大幅下降至365 °C时接近11 GPa；随后在400 °C略微回升至12 GPa以上，并在450 °C降至10.5 GPa。
含水热解过程中水的存在增强了热裂解（thermal cracking）现象。这导致生成更多的油和沥青质（bitumen），并使干酪根更加多孔。研究人员认为，这是HP样品力学性质较低的主要原因。反卷积（deconvolution）方法还表明，样品力学性质的分布曲线由三个不同的力学相组合而成。该研究通过明确区分热解作用对干酪根控制力学行为的影响，为页岩在热成熟过程中的微观力学演化提供了新认识，并建立了力学性质与成熟路径之间的直接联系。这为预测富有机质页岩中刚度、脆性及变形响应的变化提供了基础，有助于优化水力压裂策略并改进非常规储层的地质力学建模。
本文发表于《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》，研究对象为Bakken组富含II型干酪根（Type II kerogen）的有机质页岩，核心问题是不同热成熟路径，尤其是含水热解（hydrous pyrolysis, HP）与无水热解（anhydrous pyrolysis, AHP），如何改变页岩的微观力学性质。研究背景在于，Bakken组是重要的非常规页岩油储层，但其低孔低渗特征使得开发高度依赖水平钻井（horizontal drilling）与水力压裂（hydraulic fracturing）。而压裂裂缝的萌生、扩展与支撑剂嵌入风险均与页岩力学性质密切相关，尤其受有机质与无机矿物共存所导致的非均质性（heterogeneity）和各向异性（anisotropy）控制。已有研究表明，干酪根通常较矿物组分更软，且在热演化过程中其化学结构、孔隙发育与相态连接关系持续变化，从而显著影响页岩整体刚度与硬度。然而，关于不同成熟路径下，尤其是有水与无水条件下，干酪根成熟对纳米尺度力学响应的差异性影响，仍缺乏系统认识。因此，这项研究旨在直接针对未经有机质提取的原位页岩样品，揭示热解路径对干酪根主导力学行为的控制作用，并为非常规储层地质力学评价提供更直接依据。
研究人员选取一件Bakken页岩样品开展人工成熟实验。样品成熟前先通过X射线衍射（X-ray diffraction, XRD）表征矿物组成，结果显示样品以石英为主，并含有黏土、黄铁矿及少量铁白云石、钾长石和斜长石。随后，将未成熟样品分别在300、325、350、365、400和450 °C条件下进行72 h热解，构建HP与AHP两类成熟路径。Rock-Eval 6热解用于测定沥青反射率（bitumen reflectance, BRo），证实样品成熟度由初始生油窗（oil window）推进至干气窗（dry-gas window）。在力学测试方面，研究人员对原始样品和12个热解样品共13个样品进行了网格纳米压痕（grid nanoindentation）测试，其中原始样品和AHP样品采用7 × 7网格，HP样品采用9 × 8网格，总计775次压痕。测试结果通过Oliver–Pharr方法计算杨氏模量与硬度，并进一步通过统计反卷积（statistical deconvolution）解析不同力学相的分布特征。研究表明，页岩力学响应不仅随成熟度演化，而且对水参与下的热解路径极为敏感；与AHP相比，HP样品在多数成熟阶段表现出更低的刚度和硬度，这主要与水促进热裂解、增强油和沥青质生成并诱导干酪根孔隙化有关。总体上，该研究建立了成熟路径—干酪根结构变化—微观力学性质演化之间的联系，对压裂优化、储层表征和地质力学建模具有实际意义。
主要技术方法概括：研究人员以Bakken组页岩单一样品为研究材料，首先利用XRD分析矿物组成；随后在1 L Hastelloy C-276反应釜中实施HP与AHP人工热解实验，温度设定为300–450 °C、保温72 h，并用Rock-Eval 6测定BRo以表征成熟度演化；之后采用Hysitron TI 900纳米压痕仪进行网格纳米压痕测试，获取载荷–位移曲线并据Oliver–Pharr模型计算杨氏模量与硬度；最后结合统计反卷积方法，对力学参数分布进行多相分离，识别样品表面的不同力学相及其体积分数特征。
3.1. Nanoindentation curve analysis
该部分展示了不同成熟度页岩的代表性纳米压痕曲线。研究人员从各样品大量曲线中选取最接近平均力学参数的曲线用于比较，结果显示，随着成熟度升高，压痕曲线整体向右偏移，意味着杨氏模量和硬度总体降低。所有样品在加载阶段均表现出弹塑性变形特征。与此同时，部分测试出现异常曲线，包括加载阶段的“pop-in”和卸载或相变相关的“elbow”现象。前者被解释为测试过程中裂纹形成的响应，后者则与局部异质域间刚度突变和机械不稳定性有关。异常曲线约占总测试数的10%，研究人员未将其视为实验误差而剔除，而是认为它们反映了页岩微结构非均质性、矿物—有机质边界应力重分布以及高成熟含水样品中孔隙化干酪根诱发的局部开裂行为。
3.2. Grid nanoindentation
该部分利用网格纳米压痕系统表征不同热解样品的力学空间分布。机械性质图显示，原始样品中高杨氏模量区域通常也具有较高硬度，说明两者存在较直接关联；同时，图谱清晰反映了页岩的复合材料属性，即硬质区域可能对应石英或长石等无机矿物，软质区域则可能与黏土和有机质有关。统计结果表明，未加热样品的平均杨氏模量和硬度分别约为28.5 GPa和0.7 GPa。AHP过程中，杨氏模量先升至34 GPa以上，随后显著下降，并在350–450 °C间维持20–25 GPa波动；硬度则在300 °C达约0.77 GPa后持续下降，在400 °C短暂回升后于450 °C降至0.32 GPa。HP过程中，硬度随温度持续下降，在450 °C约为0.15 GPa；杨氏模量则先降后升，再明显下跌至约11 GPa，并在高温阶段维持低值波动。研究人员据此指出，水的参与显著改变了页岩热成熟期间的力学演化轨迹。
3.3. Deconvolution
该部分使用统计反卷积方法解析页岩表面的力学相组成。对原始样品的杨氏模量和硬度直方图反卷积后，均可分离出三个不同分布峰。第一相代表较软组分，主要归因于有机质或黏土；第三相代表较硬组分，主要与石英、白云石、长石等粉砂级包裹体有关；第二相则被解释为软硬相边界过渡区在随机压痕作用下产生的中间统计相，而非独立矿物学相。对所有样品的分析均显示三相结构持续存在，但不同温度和热解条件下各子峰均值、标准差和体积分数明显变化，说明干酪根成熟确实重塑了页岩力学分相特征。研究同时强调，这种三相划分本质上是统计意义上的力学相划分，并未通过扫描电镜或岩相成像作直接矿物学验证。
3.4. Kerogen maturity and its effect on mechanical properties of shale
这是全文的核心结果部分，重点解释干酪根成熟如何控制页岩力学性质变化。研究人员借助Rock-Eval 6测得BRo，证实样品在HP和AHP过程中随温度升高逐步从初始生油窗演化至干气窗。结果显示，在300 °C这一初始热解温度下，AHP与HP样品表现出显著分异：AHP样品的杨氏模量和硬度均升高，而HP样品则出现下降。研究人员将这种差异归因于HP样品在相同温度下成熟度略高、生成更多沥青质和油，而这些生成物较未成熟干酪根具有更低力学性能；同时，油气排出后干酪根骨架孔隙化增强，也会进一步降低承载能力。随着成熟继续推进，两类样品硬度均下降，且到365 °C即湿气窗末端时差异明显缩小；但杨氏模量方面，HP样品始终低于AHP样品。研究人员指出，水的存在使沥青质保留结构水，从而提供氢源、增强热裂解并促进更多油和沥青质生成，因此HP样品在生油窗阶段的力学软化更为显著。进入干气窗后，两种路径的差异进一步扩大：HP样品力学性质整体维持低位，变化不大；AHP样品则出现杨氏模量和硬度短暂回升后再下降的现象。研究人员将AHP中的这种“峰值”解释为无水条件促进交联反应（cross-linking reactions）和固体热沥青（pyrobitumen）生成，从而带来暂时性结构致密化与刚度增强；而HP中由于水抑制热沥青形成、持续促进裂解反应与高孔隙度保持，因此未出现类似增强。这说明水在高成熟阶段不仅影响成熟速率，更改变干酪根结构重组方式及其对页岩力学框架的贡献。
3.5. Limitations of the study
该部分指出了研究边界。尽管775次压痕构成了较大数据集，但结果仍受限于所选Bakken页岩样品的代表性以及实验室受控热解条件。由于缺少原位储层压力、流体化学环境以及长期成岩作用等因素，研究结论不能直接等同于现场尺度行为。此外，统计反卷积虽然有助于区分不同力学相，但其结果依赖输入假设，可能简化了复杂微结构相互作用。因此，研究人员认为未来需要结合现场数据和多尺度力学测试进一步验证与扩展。
讨论总结：本文的讨论围绕“成熟路径差异如何通过改变干酪根结构与相态分布进而调控页岩微观力学性质”这一主线展开。研究结果表明，Bakken页岩的力学演化不是简单随温度单调变化，而是受热解介质条件显著控制。AHP与HP在低成熟阶段就出现相反响应，说明水不仅改变成熟程度，而且改变了热演化的反应路径。HP中更强的热裂解带来更多油和沥青质生成，并使干酪根孔隙化程度更高，最终导致更低的杨氏模量和硬度；AHP则更易发生交联与热沥青形成，在高成熟阶段表现出短暂增强。反卷积结果进一步说明，这种力学变化并非单一组分均匀改变，而是三类统计力学相的均值和占比同步重塑。研究据此建立了热成熟路径、干酪根结构演化与页岩力学响应之间的对应关系，并强调纳米尺度非均质性对宏观压裂行为和地质力学建模具有基础意义。
研究结论翻译：
该研究系统评估了II型干酪根成熟对Bakken页岩力学性质的影响。未成熟样品分别在有水和无水条件下，于300、325、350、365、400和450 °C进行热解。随后，研究人员采用网格纳米压痕方法测试了共13个样品，包括未加热样品和成熟样品。未加热样品与AHP样品采用7 × 7网格，HP样品采用9 × 8网格，因此总共实施了775次纳米压痕测试。Rock-Eval 6热解结果表明，随着温度升高，样品中的干酪根成熟度不断提高，测试样品覆盖了从初始生油窗到干气窗的成熟范围。反卷积分析表明，所有样品的力学性质均可分解为三个力学相。研究结果可概括如下：II型干酪根成熟显著影响其结构性质并改变其力学特征；II型干酪根的热演化能够显著改变Bakken页岩的力学性质；在整个人工成熟过程中，水的存在对II型干酪根成熟具有关键作用；含水热解过程中，水增强热裂解反应，促进更多油和沥青质生成，并使干酪根更加多孔，这是HP样品力学性质较低的主要原因；结果趋势表明，样品进入生油窗后力学性质开始下降，并在湿气窗末端达到最低值；在干气窗中，HP样品的力学性质基本保持不变，而AHP样品则在该阶段中部出现峰值，这主要归因于无水条件下交联反应更易发生。最后，该研究通过定量揭示不同热成熟路径（含水与无水）如何在受控条件下支配含干酪根体系的纳米尺度弹性性质，阐明了富有机质页岩的力学演化机制。总体而言，流体参与的成熟过程在调控刚度、力学非均质性和相间相互作用方面具有关键作用，并会进一步影响岩石整体力学行为。通过将纳米尺度测量与宏观地质力学意义联系起来，该研究为更准确的非常规储层建模提供了框架，并对储层表征以及水力压裂、水平钻井等现场作业具有直接参考价值。
