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

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

【技术实现步骤摘要】
一种声发射/微震传感器校准系数量化方法


[0001]专利技术属于声发射/微震传感器数据校准
，具体涉及一种声发射/微震传感器校准系数量化方法。

技术介绍

[0002]随着我国大型土建工程规模的不断增大，如矿山工程、核电工程、隧道工程等，岩石力学理论及相应更加精细的岩体内部破裂机理的研究愈发重要。声发射是岩石材料在外载荷作用下由于自身变形断裂导致储存的应变能以应力波的形式快速释放的一种现象，声发射现象反映了岩石微裂纹萌生扩展及贯通破坏过程，是探究岩石破裂机理的一种重要监测手段，利用声发射波形信息对深入研究和理解岩石的破坏特征、损伤演化等具有重要意义。
[0003]随着岩石力学的不断发展，对岩石内部破裂机制的研究逐渐深入，近年来越来越多的学者将矩张量理论引入到岩石内部破裂的震源机制研究中，以期得到岩石裂纹萌生、扩展、贯通及相互作用的内部损伤演化过程和规律及其与应力之间的对应关系，从而应用于工程中指导生产。而波形信息采集的准确性是矩张量反演方法的合理性的必要条件，初动振幅幅值作为矩张量反演的初始参数其准确性直接影响反演结果。
[0004]然而，实验采集声发射应力波初动振幅时，由于传感器与岩壁耦合效果差等因素，导致传感器接收到的应力波初动振幅与真实应力波初动振幅大小存在差异，因此，亟需一种传感器校准系数量化方法，提高初动振幅提取精度，为真实解算岩石内部裂纹大小与演化规律奠定基础。

技术实现思路

[0005]为了解决上述问题，本专利技术提供一种声发射/微震传感器校准系数量化方法。
[0006]本专利技术提供一种声发射/微震传感器校准系数量化方法，包括如下步骤：
[0007]步骤1、落球实验及其三维数值模型的建立：利用声发射/微震传感器采集实验中落球在岩石中激发的应力波，建立落球实验对应的三维实体模型，划分网格并对落球及岩石的网格单元进行物理力学参数的赋值，建成三维数值模型；
[0008]步骤2、落球激发应力波传播过程的数值模拟：对三维数值模型施加边界条件，并模拟球体下落、撞击岩石并在岩石中激发应力波的传播过程；
[0009]步骤3、比较落球实验及数值模拟结果，确定各声发射/微震传感器的校准系数。
[0010]在本专利技术的声发射/微震传感器校准系数量化方法中，所述步骤1 具体为：
[0011]步骤1.1、落球实验：在制作好的直径为5cm、高为10cm的圆柱体的标准岩样的表面标记点处涂抹凡士林并贴置声发射/微震传感器；将直径为0.8cm的钢球置于岩样上表面中心点正上方10cm处；准备就绪后，让钢球做自由落体撞击岩样上表面中心点并实时采集落球激发的应力波波形信息，重复3次上述实验以减小实验误差；
[0012]步骤1.2、三维实体建模：采用数值模拟软件建立以球体、圆柱体为基本构件，球体
位于圆柱体正上方10cm处为相对空间位置的三维实体模型，其中球体、圆柱体尺寸与声发射实验选取的钢球和标准岩样的物理模型尺寸相一致；
[0013]步骤1.3、网格划分：三维实体模型结构简单，尺寸较小，网格尺度选择极细化处理以保证模型计算精度；
[0014]步骤1.4、物理力学参数的赋值：通过数值模拟软件分别对钢球和标准岩样进行物理力学参数赋值，包括杨氏模量、泊松比、密度，由于多数数值模拟软件无法直接对波速赋值，通过赋予材料弹性模量及密度来间接赋值波速，波速与弹性模量、密度间的关系如下式表示：
[0015][0016]式中，E为标准岩样的弹性模量，单位Pa；ρ为标准岩样的密度，单位kg/m3。
[0017]在本专利技术的声发射/微震传感器校准系数量化方法中，所述步骤2 中对三维数值模型施加边界条件具体为：
[0018]首先，钢球撞击标准岩样上表面，钢球及标准岩样表面区域独立且不会发生表面穿透，因此将钢球与标准岩样上表面设置为接触对，钢球设置为源，标准岩样上表面设置为目标；
[0019]其次，为简化模拟过程，使钢球以匀速状态撞击标准岩样上表面中心点，设置钢球初始速度为v
P
：
[0020][0021]式中：v
P
为钢球的初始速度，单位为m/s；d为钢球下表面到标准岩样上表面的初值距离，单位为m；g为重力加速度；
[0022]最后，需解算标准岩样表面的固体力学速度场波形数据，只对岩石的底面设置固定约束。
[0023]在本专利技术的声发射/微震传感器校准系数量化方法中，所述步骤3 具体包括如下步骤：
[0024]步骤3.1：提取落球实验中的应力波初动振幅；
[0025]步骤3.2：提取数值模拟结果中的应力波初动振幅；
[0026]步骤3.3：分别解算应力波在岩石中传播的衰减率和传感器、岩壁耦合界面对应力波的衰减系数，对声发射/微震传感器精确校准。
[0027]在本专利技术的声发射/微震传感器校准系数量化方法中，所述步骤 3.1具体包括如下步骤：
[0028]步骤3.1.1：对声发射/微震传感器采集的波形数据进行预处理，将电信号波形转换为速度波形(传感器电压
‑
速度转换系数V/(m/s)：501.19)，去掉尾部四分之一波形数据，并截取第一个非0点至最后一个非0点之间的波形数据作为最终波形数据；
[0029]步骤3.1.2：采用赤池信息准则AIC方法依次计算初始到时：选取不同的波形节点，截取不同长度的波形数据，计算不同截取方式对应的AIC值，提取AIC最小值对应的波形节点作为初始到时：
[0030]AIC＝k
×
lgs
12
+k2×
lgs
12
[0031]式中：k为截取的波形长度，为截取的波形数据的方差，k2为剩余波形长度，剩余波形数据的方差；
[0032]步骤3.1.3：根据下式提取速度波形初动振幅：
[0033]Au'
n
＝max||BO
n
||
∞
[0034]式中：Au'
n
为落球实验中第n个声发射/微震传感器采集的应力波的初动振幅，BO
n
为第n个声发射/微震传感器采集到的波形数据中从初始到时至下一个非0数据之间的波速矩阵，||BO
n
||
∞
表示对波速矩阵BO
n
中的元素取绝对值，max||BO
n
||
∞
表示选取波速矩阵的元素绝对值的最大值。
[0035]在本专利技术的声发射/微震传感器校准系数量化方法中，所述步骤 3.2具体包括如下步骤：
[0036]步骤3.2.1：根据数值模拟结果，提取数值模拟中各个模拟传感器监测点处的速度分量矩阵，并依次计算出沿模拟传感器监测点方向的速度波形矩阵：
[0037]V
m
＝[V
xm
,V
ym
,V
zm
][0038]V
′
m
...

【技术保护点】

【技术特征摘要】
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数据之...

【专利技术属性】
技术研发人员：张鹏海，马庆山，侯俊旭，刘溪鸽，邓文学，杨天鸿，
申请(专利权)人：东北大学，
类型：发明
国别省市：