一种声发射/微震传感器校准系数量化方法

技术特征：
1.一种声发射/微震传感器校准系数量化方法，其特征在于，包括如下步骤：步骤1、落球实验及其三维数值模型的建立：利用声发射/微震传感器采集实验中落球在岩石中激发的应力波，建立落球实验对应的三维实体模型，划分网格并对落球及岩石的网格单元进行物理力学参数的赋值，建成三维数值模型；步骤2、落球激发应力波传播过程的数值模拟：对三维数值模型施加边界条件，并模拟球体下落、撞击岩石并在岩石中激发应力波的传播过程；步骤3、比较落球实验及数值模拟结果，确定各声发射/微震传感器的校准系数。2.如权利要求1所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤1具体为：步骤1.1、落球实验：在制作好的直径为5cm、高为10cm的圆柱体的标准岩样的表面标记点处涂抹凡士林并贴置声发射/微震传感器；将直径为0.8cm的钢球置于岩样上表面中心点正上方10cm处；准备就绪后，让钢球做自由落体撞击岩样上表面中心点并实时采集落球激发的应力波波形信息，重复3次上述实验以减小实验误差；步骤1.2、三维实体建模：采用数值模拟软件建立以球体、圆柱体为基本构件，球体位于圆柱体正上方10cm处为相对空间位置的三维实体模型，其中球体、圆柱体尺寸与声发射实验选取的钢球和标准岩样的物理模型尺寸相一致；步骤1.3、网格划分：三维实体模型结构简单，尺寸较小，网格尺度选择极细化处理以保证模型计算精度；步骤1.4、物理力学参数的赋值：通过数值模拟软件分别对钢球和标准岩样进行物理力学参数赋值，包括杨氏模量、泊松比、密度，由于多数数值模拟软件无法直接对波速赋值，通过赋予材料弹性模量及密度来间接赋值波速，波速与弹性模量、密度间的关系如下式表示：式中，e为标准岩样的弹性模量，单位pa；ρ为标准岩样的密度，单位kg/m3。3.如权利要求1所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤2中对三维数值模型施加边界条件具体为：首先，钢球撞击标准岩样上表面，钢球及标准岩样表面区域独立且不会发生表面穿透，因此将钢球与标准岩样上表面设置为接触对，钢球设置为源，标准岩样上表面设置为目标；其次，为简化模拟过程，使钢球以匀速状态撞击标准岩样上表面中心点，设置钢球初始速度为v
p
：式中：v
p
为钢球的初始速度，单位为m/s；d为钢球下表面到标准岩样上表面的初值距离，单位为m；g为重力加速度；最后，需解算标准岩样表面的固体力学速度场波形数据，只对岩石的底面设置固定约束。4.如权利要求1所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤3具体包括如下步骤：步骤3.1：提取落球实验中的应力波初动振幅；
步骤3.2：提取数值模拟结果中的应力波初动振幅；步骤3.3：分别解算应力波在岩石中传播的衰减率和传感器、岩壁耦合界面对应力波的衰减系数，对声发射/微震传感器精确校准。5.如权利要求4所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤3.1具体包括如下步骤：步骤3.1.1：对声发射/微震传感器采集的波形数据进行预处理，将电信号波形转换为速度波形(传感器电压-速度转换系数v/(m/s)：501.19)，去掉尾部四分之一波形数据，并截取第一个非0点至最后一个非0点之间的波形数据作为最终波形数据；步骤3.1.2：采用赤池信息准则aic方法依次计算初始到时：选取不同的波形节点，截取不同长度的波形数据，计算不同截取方式对应的aic值，提取aic最小值对应的波形节点作为初始到时：aic＝k
×
lgs
12
k2×
lgs
12
式中：k为截取的波形长度，为截取的波形数据的方差，k2为剩余波形长度，剩余波形数据的方差；步骤3.1.3：根据下式提取速度波形初动振幅：au'
n
＝max||bo
n
||
∞
式中：au'
n
为落球实验中第n个声发射/微震传感器采集的应力波的初动振幅，bo
n
为第n个声发射/微震传感器采集到的波形数据中从初始到时至下一个非0数据之间的波速矩阵，||bo
n
||
∞
表示对波速矩阵bo
n
中的元素取绝对值，max||bo
n
||
∞
表示选取波速矩阵的元素绝对值的最大值。6.如权利要求5所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤3.2具体包括如下步骤：步骤3.2.1：根据数值模拟结果，提取数值模拟中各个模拟传感器监测点处的速度分量矩阵，并依次计算出沿模拟传感器监测点方向的速度波形矩阵：v
m
＝[v
xm
,v
ym
,v
zm
]v
′
m
＝i
m
×vm
式中，v
m
为数值模拟中第m个模拟传感器监测点处的速度分量矩阵，v
xm
,v
ym
,v
zm
分别为x、y、z方向的速度分量，单位m/s，v
′
m
为垂直于第m个模拟传感器监测点方向的速度矩阵，i
m
为各模拟传感器监测点的方向向量；步骤3.2.2：根据下式计算数值模拟中各波形数据的初动振幅：au
m
＝max||bo
m
||
∞
式中，au
m
为数值模拟中第m个模拟传感器监测点采集的应力波的初动振幅；m为模拟传感器监测点的序号，bo
m
为v
′
m
中截取从初始到时到下一个非0数据点之间的波速矩阵，||bo
m
||
∞
表示对波速矩阵bo
m
中的元素取绝对值，max||bo
m
||
∞
表示选取波速矩阵的元素绝对值的最大值，初始到时的计算方法与步骤3.1.2相同。7.如权利要求1所述的声发射/微震传感器校准系数量化方法，其特征在于，所述步骤3.3具体为：步骤3.3.1：声发射/微震传感器的校准系数计算：在模拟计算中岩石材料往往设置为均值材料，而真实的岩石材料内部构造十分复杂，应力波在岩石内部传播过程中衰减更加
明显，设置应力波在该岩石材料中传播单位长度的幅值衰减系数为：式中：n表示落球实验中声发射/微震传感器的数量和数值模拟中模拟传感器监测点数量；d为应力波在岩石中传播的衰减率；l
i
为声源到第i个声发射/微震传感器的距离，l
j
为声源到第j个声发射/微震传感器的距离；au
i
为数值模拟中第i个模拟传感器监测点采集的应力波的初动振幅，au
j
为数值模拟中第j个模拟传感器监测点采集的应力波的初动振幅；au
′
i
为落球实验中第i个声发射/微震传感器采集的应力波的初动振幅，au'
j
为落球实验中第j个声发射/微震传感器采集的应力波的初动振幅；步骤3.3.2：在实验中因传感器与岩壁耦合效果差原因导致提取的初动振幅降低，其振幅的衰减系数为t，即传感器、岩壁耦合界面对应力波的衰减率，因此定义数值模拟计算应力波初动振幅au与实验提取的应力波初动振幅au'关系为：au'＝t
×
(au-d
×
l)式中：t为传感器、岩壁耦合界面对应力波的衰减率；l
n
为声源到第n个声发射/微震传感器n的距离；分别解算应力波在岩石中传播的衰减率和传感器、岩壁耦合界面对应力波的衰减系数，对声发射/微震传感器精确校准。步骤3.3.3：声发射/微震传感器的校准计算公式为：其中，用d和t表示声发射/微震传感器的校准系数校准系数。

技术总结
本发明的一种声发射/微震传感器校准系数量化方法，包括如下步骤：进行落球试验及建立落球试验三维数值模型；对建立的三维数值模型施加边界条件，并模拟球体下落、撞击岩石并在岩石中激发应力波的传播过程；分别提取落球试验以及数值模拟结果中声发射/微震传感器对应位置处的应力波初动振幅确定各声发射/微震传感器的校准系数。本发明结合物理试验及数值模拟，同时解算出应力波在岩石中传播的衰减程度及各传感器的校准系数，消除了因传感器与岩壁耦合效果不同带来的误差，有利于提高声发射/微震波的采集精度，为反演岩石内部裂纹大小与演化规律提供更加精确的基础参数。演化规律提供更加精确的基础参数。演化规律提供更加精确的基础参数。


技术研发人员：张鹏海 马庆山 侯俊旭 刘溪鸽 邓文学 杨天鸿
受保护的技术使用者：东北大学
技术研发日：2022.05.30
技术公布日：2022/8/16