基于颗粒离散元PFC的岩石单轴压缩声发射矩张量模拟

1年前浏览7163

1. 引言

应力扰动，无论是自然的还是由人类活动引起的，都会引起地下地层的地震活动。当受到诱导应力时，完整的岩石或预先存在的裂缝可能会突然破坏。岩石灾难性破坏涉及以声发射（AE）、微震或诱发地震的形式发出不同震级地震活动。当进行实验室实验以分析断裂定位时，通过微裂纹的形成和相互作用可以产生声发射。

AE数据为研究岩石破坏的机理提供了宝贵的信息。AE事件的监测已用于研究导致岩石破坏的微裂纹聚结过程。由于完整岩石的破裂性质和断裂带的凹凸不平，因此对声发射数据进行了分析，以研究更大范围内的机制。在岩石力学试验（如三轴压缩试验）期间，可以使用岩石试样周围的一系列传感器记录实验室中的AE事件。AE震源的空间-时间分布能够表征完整岩石或具有预先存在裂缝的岩石的破坏。跟踪微裂纹位置和相互作用可以绘制剪切裂缝的成核、扩展和定向图。此外，可以通过在不同围压下进行三轴压缩试验来研究不同带的形成（例如，压实或剪切）。

除了找到AE事件的位置外，对AE事件源机制的分析还提供了信息，以确定导致岩石破坏的破坏模式。近年来，大量研究工作者进行了不同的实验工作，以研究AE事件的震源机制。矩张量分析也被用作解释AE事件源机制的定量方法。

微尺度数值技术也可用于推断基本的岩石破坏机制。通过使用离散单元法（DEM），我们可以模拟完整和断裂岩石在压缩下的破坏，以研究断裂过程中涉及的主要机制。不同参数的影响，例如，颗粒形状和尺寸、晶粒互锁通过单个裂纹的能量释放诱导更多裂纹，颗粒的不规则性和角度对通过DEM模拟岩石力学行为和破坏模式进行了广泛研究。通过DEM研究晶粒位移和微裂纹发展，还可以模拟AE事件，并将其定位在微尺度数值样本中。

在这项工作中，我们通过确定岩石破坏产生的AE事件的矩张量、矩张量的分解来研究微裂纹产生的机制，以及解释每个AE事件的不同震源机制，重点对AE事件的震源行为进行定量分析。由于试验设计的复杂性，研究震源机制的实验工作可能限于特定场景（例如，岩石类型和围压）。使用机械地震耦合数值模型可以实现更广泛的场景。在这种情况下，我们可以更好地理解岩石破坏期间AE事件的震源机制。通过模拟颗粒和颗粒之间的接触，我们可以研究不同类型的岩石破坏。

2. 建模方法

AE模拟算法在PFC（particle Flow Code）中的平行粘结模型（BPM）中实现，PFC是Itasca开发的显式离散元方法（DEM）的一个分支。这里，假设粒子是刚性的（不可变形的），并通过模拟粒子之间相互作用行为的接触模型结合。假设每一个粘结断裂都是一个微裂纹，并且每个微裂纹都会在测试期间导致AE活动。通过在测试期间监测所有粘结断裂，可以收集AE信息，并计算力矩张量（MT），以分析断裂类型和破坏机制。

2.1. 平行粘结模型

平行粘结模型（BPM）作为一种离散元方法，最初由Cundall引入，并由Cundal和Strack开发。Potyondy等人详细描述了BPM的原理。该模型模拟了由二维或三维粘结的非均匀尺寸圆形或球形颗粒表示的固体材料的运动和相互作用。模拟材料的特性通常由颗粒和粘结键的刚度和强度微参数决定。商业软件PFC中提供了两种类型的BPM，即接触粘结模型和平行粘结模型。接触粘结作用于接触点，只能传递法向力和剪切力。平行粘结作用于两个颗粒之间的横截面，可以传递力和力矩，这是岩石或类岩石材料的更现实的粘结模型。

2.2. AE模拟和矩张量分解

当粘结键断裂时，源极粒子将移动，相邻触点将发生一些变形。因此，由于裂纹的形成，周围接触处的力将发生变化。力矩张量的分量可以通过将接触点处的力变化乘以接触点到裂纹位置的距离来计算。值得注意的是，AE事件可以由多个微裂纹组成。在空间和时间上紧密发生的微裂纹被视为一个AE事件，事件质心被假设为事件的几何中心（图1（a））。力矩张量可以通过在源边界周围积分来计算，对于离散元方法，积分由以下公式计算：

其中ΔF_i是接触力变化的第i个分量，R_j是接触点和事件形心之间的距离。由于假设力矩张量是对称的，因此力矩张量的非对角分量将被假设为相等，并且将通过等式（1）计算的两个值进行平均来计算。此外，在AE事件持续期间的每个时间步计算力矩张量，为每个AE事件存储的单个力矩张量是在最大标量力矩的时间步计算的结果。