利用颗粒离散元PFC矩张量反演模拟巴西试验中的声发射

5月前浏览16705

本文摘要（由AI生成）：

本文介绍了离散元方法中的颗粒粘结模型（BPM），并通过模拟研究了岩石或类岩石材料在加载过程中的声发射（AE）行为。模拟利用商业软件PFC，考虑了接触粘结和平行粘结两种BPM。文章详细描述了AE模拟和矩张量分解的过程，包括裂纹形成时力的变化和矩张量的计算。模拟结果展示了破坏模式、力-位移曲线和震级分布，并通过后处理得到裂纹类型和AE事件的演化过程。此外，文章还利用AE b值和T-k绘图法评估了裂缝网络的复杂性和震源破裂类型。这些研究为深入理解岩石材料在加载过程中的力学行为和破坏机制提供了重要工具。

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事件存储的单个矩张量是在最大标量力矩的时间步计算的结果。

图1 示例AE事件及矩张量示意图

矩张量可以被认为是源的等效体力，矩张量表示将矩张量矩阵的主值描述为两组向量，其方向和长度分别表示方向和大小（图1（b））。

将公式（1）计算得到的矩张量分解为球张量与偏张量，则声发射破裂类型可根据球张量占比P_iso进行判定

2.3. 模型描述

数值试样的尺寸与实验试样的尺寸相同（如图2所示），直径为50 mm，厚度25 mm。考虑到计算时间，颗粒尺寸遵循0.15 mm至0.25 mm的均匀分布，样品中总共产生14296个颗粒。

图2 数值模型

3. 模拟直接得到结果

模拟直接得出的破坏模式、力-位移曲线和震级分别如图3所示。

(a) 破坏模式

(b) 力-位移曲线

图3 模拟直接得到的结果

模拟得到的dat文件如图4所示。

(a) AE

（b）裂纹文件

（c）矩张量文件

图5 模拟得到的dat文件

3. 后处理得到结果

应用矩张量分析程序对圆盘试样在径向加载作用下的声发射信号进行分析，图6显示了声发射空间分布，利用图5（c）中ios和公式（2）可将破坏模式分为拉伸、剪切和混合。图6中蓝色圆圈表示拉伸型裂纹，红色矩形表示剪切型裂纹，绿色三角形表示混合型裂纹。

图6 不同类型裂纹空间分布图

在模拟过程中，通过矩张量文件中aenum记录了AE事件，进而可以做出应力应变、AE事件数和累计AE事件数之间的关系，如图7所示。

图7 AE事件和AE累积事件的演变过程

4. 通过AE b值评估裂缝网络的复杂性

断裂系统的产生伴随着现场的微震或实验室的AE，这有助于描述断裂分布。同时，我们注意到，Gutenberg和Richter提出的b值不仅可以用于描述AE震级的分形特征，还可以用于描述声发射源的空间分形特征。较大的b值意味着更复杂的裂缝网络，而较小的b值则意味着更简单的裂缝网络。描述地震震级和频率之间关系的古腾堡和里希特公式如下：

其中M为地震震级；N是震级大于M的地震的累积次数；a和b是两个统计常数。根据矩张量文件中的mag列数据可以得出下图，进而求出b值。

图8 AE事件数量和AE震级的关系（直方图显示频率-震级分布，散点显示大于某个震级的AE事件数量）。

5. T-k绘图法

Hudson（1989）提出了矩张量的两个参数T和k，以分析震源破裂类型。参数T用于测量源的偏差部分，其范围从位于−1的纯正CLVD到位于+1的纯负CLVD，与位于原点的纯DC交叉。参数k用于测量各向同性部分，其范围从位于+1的均匀爆炸类型到位于-1的均匀内爆类型。

根据矩张量文件里M11、M12 和M22，利用T-k图MTLAB代码（包含在附件中）可做出图9。由于这是一个二维模型，模型相对均质，因此做出的T-k图呈一条线。三位模型做出的是散点图。

图9 T-k图

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

本文摘要（由AI生成）：本文采用量、矩张量的分解方法研究了微裂纹产生机制及AE事件的震源行为。利用离散元方法中的平行粘结模型进行AE模拟，并通过矩张量分解分析断裂类型。模拟结果显示，破坏模式、力-位移曲线和震级分布均可直接得出。此外，还生成了包含AE事件、裂纹和矩张量信息的dat文件，为进一步分析提供了数据支持。本研究有助于深入理解岩石破坏期间的AE事件震源机制，并为实际工程中的岩石力学问题提供理论依据。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））。力矩张量可以通过在源边界周围积分来计算，对于离散元方法，积分由以下公式计算：其中ΔFi是接触力变化的第i个分量，Rj是接触点和事件形心之间的距离。由于假设力矩张量是对称的，因此力矩张量的非对角分量将被假设为相等，并且将通过等式（1）计算的两个值进行平均来计算。此外，在AE事件持续期间的每个时间步计算力矩张量，为每个AE事件存储的单个力矩张量是在最大标量力矩的时间步计算的结果。矩张量可以被认为是源的等效体力，矩张量表示将矩张量矩阵的主值描述为两组向量，其方向和长度分别表示方向和大小（图1（b））。为了识别断裂机制，提出了几种分解力矩张量的方法。在岩石和类岩石材料的分析中，通常使用两种不同的解释力矩张量的方法。第一种方法是根据特征值分析将其分为拉伸裂纹、剪切裂纹或混合型裂纹。以下分解作为相对比率X、Y和Z获得：中间特征值/最大特征值其中X是剪切分量，Y是偏分量，Z是抗拉静水压平均值。然后，将裂纹视为剪切裂纹（X>60%）、拉伸裂纹（X<40%）和混合模式裂纹（40%<X<60%）。另一种方法是将力矩张量分为各向同性部分和偏分量部分：其中tr(M)是特征值之和，mi*是偏特征值。基于各向同性比，裂缝类型被视为“拉伸”、“内爆”或“剪切”，其计算公式如下：该比率范围从−100%（纯内爆）到100%（纯爆炸），R=0表示纯剪切破坏机制。因此，如果R>30%，则裂纹被视为“拉伸”，如果R<−30%，则裂纹为“内爆”，如果−30%<R<30%，则裂纹视为“剪切”。然而，前者中的混合模式裂纹被称为拉伸和剪切分量的总和，其机制不够明确，因此本研究选择了后者。2.3. 模型描述数值试样的尺寸与实验试样的尺寸相同（如图2所示），高度为100 mm，直径为50 mm。考虑到计算时间，颗粒尺寸遵循0.4 mm至0.6 mm的均匀分布，样品中总共产生8133个颗粒。3. 结果分析模拟直接得出的破坏模式、力-位移曲线和震级分别如图3所示。(a) 破坏模式(b) 力-位移曲线(c) 震级分布图3 模拟直接得到的结果模拟得到的dat文件如图4所示。(a) AE（b）裂纹文件（c）矩张量文件图5 模拟得到的dat文件

有附件