页岩的力学性能对于页岩气和致密油的开采至关重要。1,2,3随着全球能源需求的增加，水力压裂被广泛用于提高页岩层的渗透率。4然而，页岩层的各向异性和异质性使得压裂设计和井筒稳定性的优化变得复杂。5早期研究主要关注页岩层的宏观性能，因此忽略了其微观结构的复杂性。6,7这种现象限制了水力压裂中页岩应力响应的准确性，而微观结构在其中起着关键作用。
页岩的微观结构作为功能性的“基因组”，决定了其力学性能，并且是评估页岩气生产可压裂性的关键。8微观结构由矿物颗粒、有机质、孔隙-裂缝网络和层理组成；这些成分的分布、丰度和耦合共同决定了页岩的力学性能和传输能力。9,10许多研究探讨了页岩微观结构与其力学性能之间的关系。10, 11, 12脆性矿物的含量决定了岩石的硬度，高石英含量的海相页岩表现出更高的硬度和破裂压力。13, 14, 15随着热成熟度的增加，有机质发生芳香化，导致页岩的整体模量增加。16,17富含粘土的层理会削弱岩石强度并产生明显的力学各向异性。3,18页岩中的微裂缝倾向于沿层理界面传播，从而改变破坏模式。19,20在矿物颗粒界面和颗粒边界方面也取得了显著进展，研究从宏观观察发展到微观力学分析。21, 22, 23页岩的颗粒边界结构（也称为界面）由三部分组成：未风化的致密界面、风化矿物和微裂缝。24, 25, 26颗粒边界通过Hall–Petch强化、控制塑性及断裂响应的位错-边界相互作用以及提供微裂缝成核点的边界应力集中来影响力学性能。
传统的宏观技术（如单轴和三轴压缩试验以及超声波速度测量）已被用于研究页岩的力学性能。27, 28这些技术依赖于岩心样本，只能提供复合岩石的平均力学性能。12到目前为止，纳米压痕试验和电子背散射衍射（EBSD）分析已将研究方向转向了页岩的微观力学性能。纳米压痕试验可以测量颗粒尺度上的弹性模量和硬度，揭示不同矿物组成、颗粒形状和颗粒边界的微观力学性能。1,29这些研究表明，力学性能也会因颗粒类型而异；即使在同一颗粒内的不同位置，也表现出不同的力学性能。30,31脆性矿物（如黄铁矿和长石）比粘土等较软的矿物更坚固。同时，EBSD研究可以描绘颗粒边界的特征参数，例如错位角、距离和长度。EBSD通常用于单相材料（如金属和陶瓷）。然而，页岩由多种矿物复合体组成，具有复杂的矿物晶体结构和极细的尺寸。因此，易受损的颗粒会导致EBSD信号质量较差，Kikuchi图案模糊或消失，给测量和解释带来挑战。先前的研究表明，低角度颗粒边界（小于15°）提供更高的强度，而高角度边界则倾向于削弱材料并促进裂纹形成。32页岩的微观力学性能反映了其颗粒微观结构和状态；然而，纳米压痕和EBSD等技术在精确表征单个颗粒内的特征方面存在局限性，特别是在单一矿物颗粒内部的变化方面。33
已经探索了多尺度和统计建模方法，以将微观力学特征放大到更大尺度行为。33,34然而，从纳米尺度到宏观尺度的预测仍是一个挑战。35,36一些研究建议将孔隙和裂缝网络纳入模型，而其他研究则强调矿物学和层理面的作用。37尽管取得了进展，但仍然缺乏对颗粒边界、接触类型和颗粒形态如何相互作用以影响页岩宏观力学性能的完整理解。以往的研究通常关注矿物组成或边界取向等单个特征，而没有考虑这些元素如何共同塑造页岩的整体行为。38,39此外，大多数研究尚未将这些微观结构因素纳入水力压裂中页岩宏观行为的模型中。36,38,40本研究的创新之处在于其多方法结合，将颗粒边界错位、颗粒形状和接触类型等微观结构特征与页岩的宏观力学性能联系起来。此外，使用Mori-Tanaka模型来放大微观结构性能，提供了更全面的理解，即微观结构如何影响宏观行为。40
微观力学放大方法旨在根据已知的组成特征和微观结构特征预测材料的有效介观性能。广泛认可的有四种主要方法：混合规则法、自洽法、Mori-Tanaka法和微分法。22,41混合规则法基于Voigt（上限）和Reuss（下限）界限，提供了最直接的方法。22它假设每种成分分别承受均匀应变或均匀应力。21自洽法模拟嵌入在未知有效性能的均匀介质中的单个夹杂物，这些性能通过自洽方式确定。这种方法特别适用于多晶体建模，其中单个颗粒可以被视为夹杂物，也适用于夹杂物体积分数较高的复合材料。24然而，其局限性在于控制方程是隐式的，需要迭代数值解。因此，在夹杂物体积分数非常高时，其预测的可靠性可能会降低。Mori-Tanaka方法可以被视为工程中最成功和最广泛使用的方法之一，它是Eshelby解的优雅扩展。21该方法通过用有限夹杂物浓度下的平均基体相场替换远场应力或应变，扩展了Eshelby解，从而考虑了夹杂物之间的相互作用。2,30微分方案有效地捕捉了增强作用，尽管需要求解控制微分方程，计算成本较高。Mori-Tanaka方法即使在夹杂物体积分数达到40%或更高时也能提供合理且稳定的预测。相比之下，自洽方法在高分数（超过50%）时可能会产生不稳定甚至不可靠的结果，因为它忽略了连续的基体相。尽管微分方案在较高夹杂物含量下也表现良好，但其计算复杂性更高且效率较低。32在所有方法中，Mori-Tanaka方法被广泛用于估算页岩的等效力学性能。42, 43, 44它符合页岩典型的软基体嵌入硬矿物颗粒的结构，能够捕捉中等至高矿物含量下的夹杂物相互作用。因此，Mori-Tanaka方法通常提供更高的准确性。
在本研究中，我们旨在通过研究颗粒边界特征、接触类型和颗粒形态如何影响页岩的力学性能来填补这一空白。我们结合了EBSD、FESEM和纳米压痕技术来分析页岩的这些微观结构特征。EBSD测试用于研究颗粒边界错位和位错分布，从而了解局部变形机制。纳米压痕测量颗粒尺度上的微观弹性模量和硬度，而FESEM有助于可视化页岩的颗粒形态。Mori-Tanaka方法用于将微观结构性能放大到宏观层面，建立微观结构与宏观行为之间的关系。通过将这些微观结构特征整合到预测模型中，工程师可以优化水力压裂设计并改善页岩储层开发。这项研究还通过提供一种新的方式来理解页岩等异质材料中微观结构与力学性能之间的复杂关系，从而推动了岩石力学的发展。
