技术领域
本发明涉及岩石物理力学性质测试技术领域，特别涉及一种基于优势声发射震源能量特征的岩石长期强度确定方法。
背景技术
岩石蠕变表征的是在某一恒定的荷载作用下岩石变形随时间增大的现象。工程实践表明：金属矿山地下开采过程中岩体的破坏与岩石蠕变性质密切相关，岩石蠕变性质是诱发金属矿山采场冒顶、采空区垮塌等典型灾害的主要因素之一。岩石的强度是随着荷载作用时间的增大而降低。岩石长期强度定义为荷载作用时间趋于无穷大时对应的强度最小值。因此，当蠕变应力小于长期强度σ
目前确定岩石长期强度σ
岩石蠕变破坏其实质是微裂纹持续演化的结果。当蠕变应力小于岩石长期强度σ
现有技术存在的问题如下：
(1)直接法存在蠕变历时长，时效性差的局限性。
(2)间接法未充分考虑岩石轴向与横向蠕变特征具有相对明显的差异，单方向蠕变特征并不能代表岩石整体蠕变特征。同时，在选取相关曲线“拐点”时也存在一定的主观性与随意性。
发明内容
本发明提供一种基于优势声发射震源能量特征的岩石长期强度确定方法，本发明尝试通过建立岩石蠕变破坏过程中微裂纹(声发射震源)演化特征与长期强度σ
本发明提供了一种基于优势声发射震源能量特征的岩石长期强度确定方法，包括以下步骤：
对标准岩石试件进行不少于6次纵波波速测试，得到岩石纵波波速的平均值；
在标准岩石试件上布置不少于6个声发射传感器，并记录每个声发射传感器的位置坐标；
将声发射传感器的位置坐标及纵波波速的平均值作为实验输入参数，对标准岩石试件进行分级蠕变破坏声发射震源定位试验，得到声发射震源的坐标；
分别提取声发射震源对应的波形文件，计算得到各个声发射声发射震源对应的P波到时与初动振幅，在此基础上，对各个声发射震源对应的矩张量进行求解，得到各矩张量对应的特征值；
根据特征值对声发射震源进行分类；
通过声发射震源的类型确定优势声发射震源；
对分级蠕变破坏声发射震源定位试验中的各级优势声发射震源的绝对能量进行拾取；
根据绝对能量计算各级蠕变过程中优势声发射震源的绝对能量的平均值；
根据各级蠕变过程中优势声发射震源的绝对能量平均值，绘制绝对能量平均值-蠕变应力曲线；
将绝对能量平均值-蠕变应力曲线上最低点与最低点之后第一个数据点对应的蠕变应力平均值作为岩石长期强度σ
可选的，分级蠕变破坏声发射震源定位试验加载方式采用分级加载，加载速率控制在0.1MPa/s～0.5MPa/s之间，各级蠕变应力取βσ
可选的，矩张量的求解具体为通过提取六个有效信号的初动振幅A
式中：A(x)—为初动振幅；
R,r—声源与传感器之的距离，及其方向余弦，由传感器坐标与震源坐标计算得到；
C
Ref(t,r)—反应系数，可通过断铅试验进行标定，取值为2，
由于矩张量
可选的，声发射震源进行分类具体为：
(1)由矩张量
(2)对3个特征值进行归一化处理，得到X,Y,Z，其中X＝λ
(3)求解方程组
(4)根据X值，按式(3)对震源类型进行辨识
可选的，蠕变破坏过程优势声发射震源类型辨识的方法包括：
(1)统计各级蠕变过程中，剪切震源、张拉震源与混合震源数量，记S
(2)统计整个蠕变过程中，剪切震源、张拉震源与混合震源数量，记S
式中：m为分级次数；
(3)比较S
可选的，优势声发射震源的绝对能量进行拾取包括：
对于某一个确定的优势声发射震源其绝对能量的拾取，以接收到声发射信号P波到时最小值min(T
可选的，分级蠕变破坏声发射震源定位试验中，记N
可选的，岩石长期强度σ
(1)取蠕变应力-优势声发射震源平均绝对能量曲线最小值为A点，其坐标表示为
(2)取蠕变应力-优势声发射震源平均绝对能量曲线A点之后第一个数据点为B，其坐标可表示为
(3)取A点X轴坐标值σ
与现有技术相比，本发明的有益效果在于：本发明提供的基于优势声发射震源能量特征的岩石长期强度确定方法，该发明以声发射震源定位为基础，以矩张量分析方法为手段，将岩石蠕变过程中声发射震源定性分为剪切、张拉与混合型震源。在此基础上，以优势声发射震源类型对应的绝对能量变化特征来确定岩石长期强度。该方法充分考虑了不同种类岩石在蠕变过程中，声发射震源类型与岩石宏观破坏模式之间的关系，搭建了岩石长期强度与岩石破坏过程中微裂纹之间的联系，具有物理意义明确、易操作的特点，是岩石长期强度确定方法的有效补充。
附图说明
图1为本发明实施例1提供的标准圆柱体试件；
图2为本发明实施例1提供的标准圆柱体试件相对坐标编制；
图3为本发明实施例1提供的传感器布置图；
图4为本发明实施例1提供的加载方式与加载速率示意图；
图5为本发明实施例1提供的震源空间坐标计算结果；
图6为本发明实施例1提供的P波到达时间t
图7为本发明实施例1提供的优势声发射震源绝对能量拾取示意图；
图8为本发明实施例1提供的蠕变应力-优势声发射震源平均绝对能量曲线；
图9为本发明实施例1提供的长期强度σ
图10为本发明实施例2提供的标准圆柱体岩石试件；
图11为本发明实施例2提供的主要仪器设备示意图，其中，图(a)为三轴流变仪，图(b)为声发射传感器；
图12为本发明实施例2提供的标准圆柱体试件相对坐标编制；
图13为本发明实施例2提供的Nano30型声发射传感器布置图；
图14为本发明实施例2提供的岩石纵波波速测试的结构示意图；
图15为本发明实施例2提供的震源空间坐标计算结果；
图16为本发明实施例2提供的剪切震源辨识结果；
图17为本发明实施例2提供的张拉震源辨识结果；
图18为本发明实施例2提供的混合震源辨识结果；
图19为本发明实施例2提供的蠕变应力-优势声发射震源平均绝对能量曲线；
图20为本发明实施例2提供的长期强度σ
图21为本发明实施例2提供的稳态蠕变速率法长期强度σ
具体实施方式
下面结合附图，对本发明的一个具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明的技术方案和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
目前确定岩石长期强度σ
岩石蠕变破坏其实质是微裂纹持续演化的结果。当蠕变应力小于岩石长期强度σ
现有技术存在的问题如下：
(1)直接法存在蠕变历时长，时效性差的局限性。
(2)间接法未充分考虑岩石轴向与横向蠕变特征具有相对明显的差异，单方向蠕变特征并不能代表岩石整体蠕变特征。同时，在选取相关曲线“拐点”时也存在一定的主观性与随意性。
基于上述问题，本发明提供一种基于优势声发射震源能量特征的岩石长期强度确定方法，尝试通过建立岩石蠕变破坏过程中微裂纹(声发射震源)演化特征与长期强度σ
本发明提供了一种基于优势声发射震源能量特征的岩石长期强度确定方法，包括以下步骤：
对标准岩石试件进行不少于6次纵波波速测试，得到岩石纵波波速的平均值；
在标准岩石试件上布置不少于6个声发射传感器，并记录每个声发射传感器的位置坐标；
将声发射传感器的位置坐标及纵波波速的平均值作为实验输入参数，对标准岩石试件进行分级蠕变破坏声发射震源定位试验，得到声发射震源的坐标；
分别提取声发射震源对应的波形文件，计算得到各个声发射声发射震源对应的P波到时与初动振幅，在此基础上，对各个声发射震源对应的矩张量进行求解，得到各矩张量对应的特征值；
根据特征值对声发射震源进行分类；
通过声发射震源的类型确定优势声发射震源；
对分级蠕变破坏声发射震源定位试验中的各级优势声发射震源的绝对能量进行拾取；
根据绝对能量计算各级蠕变过程中优势声发射震源的绝对能量的平均值；
根据各级蠕变过程中优势声发射震源的绝对能量平均值，绘制绝对能量平均值-蠕变应力曲线；
将绝对能量平均值-蠕变应力曲线上最低点与最低点之后第一个数据点对应的蠕变应力平均值作为岩石长期强度σ
本发明提供的基于优势声发射震源能量特征的岩石长期强度确定方法，该发明以声发射震源定位为基础，以矩张量分析方法为手段，将岩石蠕变过程中声发射震源定性分为剪切、张拉与混合型震源。在此基础上，以优势声发射震源类型对应的绝对能量变化特征来确定岩石长期强度。该方法充分考虑了不同种类岩石在蠕变过程中，声发射震源类型与岩石宏观破坏模式之间的关系，搭建了岩石长期强度与岩石破坏过程中微裂纹之间的联系，具有物理意义明确、易操作的特点，是岩石长期强度确定方法的有效补充。
本发明适用于单轴压缩条件下岩石长期强度σ
实施例1：
1.1标准岩石试件的制备
按《工程岩体试验方法标准》(GB/T50266-2013)进行试件制备。试件可由钻孔岩芯或岩块制备。以标准直径50mm，高度100mm的圆柱体为例。试件制备所需的仪器设备包括：取芯机、切割机、磨石机等。其中，取芯机用于钻取直径为50mm，高度大于100mm的岩芯。切割机用于对经钻取后的岩芯进行切割，以得到直径为50mm，高度略大于100mm的圆柱体岩芯。磨石机用对切割后圆柱体岩芯进行端面进行打磨处理，并最终制备得到直径50mm，高度100mm的圆柱体。为保证试验精确，试件两端面不平行度小于0.05mm，沿试件高度方向，直径的误差小于0.3mm，试件端面应垂直于试件轴线，偏差小于0.25°。标准圆柱试件见图1。
1.2岩石蠕变过程声发射震源定位试验
1.2.1仪器与设备
所需的仪器与设备包括：具有蠕变功能的刚性压力机、声发射仪、声发射传感器及相关固定装置、声波仪、应变片等。其中，刚性压力机应具备自动采集轴向应力(力)与轴向应变(位移)的功能。如轴向应变(位移)无法采集，可在试件上粘贴应变片替代。声发射仪应具备三维定位功能，建议采用PAC声发射仪备PCI-2系统。声发射传感器建议采用Nano30传感器。声波仪用于测试岩石试件的纵波波速，建议采用RSM-RCT(B)声波仪。
1.2.2标准圆柱体相对坐标编制
为便于声发射传感器安装定位，以及震源空间坐标的计算，需进行标准圆柱体相对坐标编制。步骤如下：
(1)以试件底面中心为坐标原点，以试件下端面为X-Y平面，以试件轴向为Z轴，建立空间坐标系。
(2)以X轴为起点，并标记为0°刻度线，逆时针沿试件表面每隔30°作Z轴平行线。
(3)以X-Y平面为起点，沿Z轴方向每隔10mm作圆周线。试件相对坐标编制，见图2。
1.2.3声发射传感器布置
共需对称布置8个声发射传感器，布置步骤如下：
(1)用铅笔沿试件表面轴线方向作直线L1。采用游标卡尺沿直线方向精确量取距上、下端面10mm的点，并将此两点(2号点、6号点)作为声发射传感器布置位置的中心点。
(2)以试件轴线作L1的对称线L2，做2号点、6号点的对称点，作为4号、8号传感器布置中心点。
(3)同理，分别沿试件表面作直线L3、L4，使L3、L4构成的平面垂直于L1、L2构成的平面，并确定剩余的1号与5号，3号与7号传感器布置中心点。
(4)最后，以试件底面中心为坐标原点，分别计算、记录以上8个传感器布置中心点坐标。震源定位传感器布置示意，见图3。
1.2.4岩石纵波波速测试
采用RSM-RCT(B)声波仪对岩石纵波波速进行测试，测试次数不少于5次，并取平均值作为震源定位的计算参数。
1.2.5岩石蠕变破坏声发射震源定位试验实施
(1)加载方式与加载速率
蠕变声发射震源定位试验加载方式采用分级加载，加载速率建议控制在0.1MPa/s～0.5MPa/s之间。各级蠕变应力建议取βσ
(2)声发射采集参数设置
声发射参数的设置与岩石种类与测试环境有关。不同的岩石在不同测试环境下声发射参数不尽相同。以红砂岩声发射测试为例，声发射参数设置，见表1。
表1声发射参数设置
(3)完成试件纵波波速测试，将试件置于具有蠕变功能的刚性压力机，完成加载速率与声发射采集参数等设置，进行岩石破坏声发射震源定位试验。其中，需保证刚性压力机与声发射采集时间同步进行，即两者开始采集的时间相同。
(4)震源空间坐标计算
声发射震源坐标可通过最二乘法、单纯形算法等计算得到。同时，也可由声发射仪直接得到。
作为示例，采用PAC声发射仪定位结果，震源定位结果如图5。
1.3矩张量求解
声发射震源类型辨识包括：波形文件的选取、P波到达时间t
1.3.1波形文件的选取
在声发射震源定位中，对于一个声发射震源其产生时间与对应空间坐标的计算，最少需要4个声发射传感器所采集的波形文件。但基于矩张量分析的声发射震源类型辨识，最少需要6声发射传感器所采集的波形文件才能实现声发射震源类型辨识。因此，在理论上用于声发射震源类型所对应的震源数小于声发射震源定位所对应的震源数。
根据实际声发射震源定位数量，通常由两种方法提取单个声发射震源类型辨识所对应的波形文件。第一种，对于震源定位数量少，可直接通过人工方式提取。第二种，对于震源定位数量多，可根据声发射仪波形文件存储规则，通过编制相关程序批量进行提取。
1.3.2P波到达时间t
以6个声发射传感器为例，即该震源相关信息“同时”被6个声发射传感器采集。
(1)根据赤池信息理准则(AIC)，采用F
AIC(k
式中：R
var、n
(2)将AIC函数最小点作为P波到达时间t
(3)以t
因此，对于6个声发射传感器，可求出6组P波到达时间t
以某一个传感为例：P波到达时间t
1.3.3矩张量求解
通过提取六个有效信号的初动振幅A
式中：A(x)—为初动振幅；
R,r—声源与传感器之的距离，及其方向余弦，可由传感器坐标与震源坐标计算得到；
C
Ref(t,r)—反应系数，可通过断铅试验进行标定，通常取2。
由于矩张量
1.4震源类型辨识
(1)由矩张量
(2)对3个特征值进行归一化处理，得到X,Y,Z。其中X＝λ
(3)求解方程组(2)
(4)根据X值，按式(3)对震源类型进行辨识
1.5长期强度σ
1.5.1蠕变破坏过程优势声发射震源类型辨识
已有的试验表明：在形状、尺寸与加载方式相同条件下，某一类型的声发射震源将“贯穿”整个岩石破坏过程，并且其数量最多。本发明定义该类震源为优势震源，其辨识方法如下：
(1)统计各级蠕变过程中，剪切震源、张拉震源与混合震源数量。记S
(2)统计整个蠕变过程中，剪切震源、张拉震源与混合震源数量。记S
式中：m—分级次数。
(3)比较S
1.5.2优势声发射震源绝对能量拾取
由震源类型辨识方法可知：对于某一个确定的优势声发射震源，其最少对应于6个声发射波形文件。即对于某一个确定的优势声发射震源，其最少存在6个绝对能量值。由于单个声发射震源空间坐标常具有一定的“随机性”。因此，在通常情况下单个震源到各接收到信号的传感器的距离L
1.5.3各级蠕变过程中优势声发射震源平均绝对能量计算
统计各级蠕变过程中，优势声发射震源数量与绝对能量。记N
1.5.4各级蠕变应力条件下优势声发射震源平均绝对能量曲线绘制
以各级蠕变应力σ
1.5.4长期强度σ
(1)取蠕变应力-优势声发射震源平均绝对能量曲线最小值为A点，其坐标可表示为
(2)取蠕变应力-优势声发射震源平均绝对能量曲线A点之后第一个数据点为B，其坐标可表示为
(3)取A点X轴坐标值σ
长期强度σ
实施例2
2.1标准岩石试件的制备
按《工程岩体试验方法标准》(GB/T50266-2013)制备直径50mm，高度100mm的圆柱体岩石试件，见图10。
2.2岩石蠕变过程声发射震源定位试验
2.2.1仪器与设备
试验加载仪器即刚性压力同采用GDS VIS 400kN HPTAS三轴流变仪，声发射定位由PCI-2声发射系统完成，同时配Nano30型声发射传感器，岩石纵波波速测试采用RSM-RCT(B)声波仪。主要仪器设备见图11。
2.2.2标准圆柱体相对坐标编制
为便于声发射传感器安装定位，以及震源空间坐标的计算，需进行标准圆柱体相对坐标编制。步骤如下：
(1)以试件底面中心为坐标原点，以试件下端面为X-Y平面，以试件轴向为Z轴，建立空间坐标系。
(2)以X轴为起点，并标记为0°刻度线，逆时针沿试件表面每隔30°作Z轴平行线。
(3)以X-Y平面为起点，沿Z轴方向每隔10mm作圆周线。
试件相对坐标编制，见图12。
2.2.3声发射传感器布置
共需对称布置8个声发射传感器，布置步骤如下：
(1)用铅笔沿试件表面轴线方向作直线L1。采用游标卡尺沿直线方向精确量取距上、下端面10mm的点，并将此两点(2号点、6号点)作为声发射传感器布置位置的中心点。
(2)以试件轴线作L1的对称线L2，做2号点、6号点的对称点，作为4号、8号传感器布置中心点。
(3)同理，分别沿试件表面作直线L3、L4，使L3、L4构成的平面垂直于L1、L2构成的平面，并确定剩余的1号与5号，3号与7号传感器布置中心点。
(4)最后，以试件底面中心为坐标原点，分别计算、记录以上8个传感器布置中心点坐标。震源定位传感器布置示意，见图13。
2.2.4岩石纵波波速测试
采用RSM-RCT(B)声波仪对岩石岩石纵波波速进行测试，测试次数不少于5次，并取平均值作为震源定位的计算参数。测试过程，见图14。
2.2.5岩石蠕变过程声发射震源定位试验实施
声发射震源定位由GDS VIS 400kN HPTAS三轴流变仪与PCI-2声发射系统完成。蠕变声发射震源定位试验加载方式采用分级加载，加载速率为0.1MPa/s。以红砂岩为例，由单轴压缩试验测得岩石平均单轴抗压强度为约60.35MPa，起裂应力约为30.00MPa。试验中，分级蠕变应力分别取23MPa、27MPa、31MPa、35MPa、39MPa、42MPa与45MPa，其数值大小分别约为起裂应力的0.80、0.90、1.03、1.16、1.30、1.40与1.50倍。声发射参数设置，参照表1。震源空间坐标计算由PCI-2声发射系统完成。图15给出了某一蠕变应力条件下震源定位结果。
2.3矩张量求解
在声发射震源定位的基础上，筛选出震源对应的6个及以的波形文件。赤池信息理准则(AIC)，将AIC函数最小点作为P波到达时间t
以某个震源对应6个波形文件为例，计算结果见表2。其中，各传感器实际接收信号时间T
表2 P波到时与初动振幅计算结果
按式(6)求解矩张量M六个独立分量：
式中：A(x)—为初动振幅；
R,r—声源与传感器之的距离，及其方向余弦，可由传感器坐标与震源坐标计算得到；
C
Ref(t,r)—反应系数，可通过断铅试验进行标定，通常取2。
求解结果见式(7)：
2.4震源类型辨识
式(7)经归一化后的特征值为[-0.3510.0141.597]。按式(8)对该震源类型进行辨识。计算结果表明：其中X的相对比例为22.865％，Y的相对比例为50.833％，Z的相对比例为26.302％。因此，该震源为张拉型震源。
对应于图15作为示例，图16～18给出了某一蠕变应力条件下剪切、张拉与混合震源辨识结果。
2.5长期强度σ
2.5.1蠕变破坏过程优势声发射震源类型辨识
(1)统计各级蠕变过程中，剪切震源、张拉震源与混合震源数量。记S
(2)统计整个蠕变过程中，剪切震源、张拉震源与混合震源数量。记S
式中：m—分级次数。
(3)比较S
表3给出了各级蠕变应力条件下各类型声发射震源数量，从表中可以看出剪切型震源最多。因此，剪切开型震源为优势震源。
表3各级蠕变应力条件下各类型声发射震源数量
2.5.2优势声发射震源绝对能量拾取
表2给出了某震源信号被6个声发射传感所采集，相应的生成了6个波形文件。从表2可以看出：经对6个传感器所生成的6个波形文件进行P波到时t
2.5.3各级蠕变过程中优势声发射震源平均绝对能量计算
统计各级蠕变过程中，优势声发射震源数量与绝对能量。记N
2.5.4各级蠕变应力条件下优势声发射震源平均绝对能量曲线绘制
以各级蠕变应力σ
2.5.4长期强度σ
图19显示，当蠕变应力小于等于39MPa时，随着蠕变应力的增大优势声发射震源平均绝对能量逐步减小；当蠕变应力大于等于39MPa时，随着蠕变应力的增大优势声发射震源平均绝对能量逐步增大；蠕变等于39MPa时优势声发射震源平均绝对能量最小，其数值约为39aJ。
(1)取蠕变应力-优势声发射震源平均绝对能量曲线最小值为A点。
(2)取蠕变应力-优势声发射震源平均绝对能量曲线A点之后第一个数据点为B。
(3)取A点X轴坐标值与B点X轴坐标值的平均值作为岩石长期强度σ
作为对比，图21给出了稳态蠕变速率法长期强度σ
以上公开的仅为本发明的几个具体实施例，但是，本发明实施例并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。
机译： 三轴压缩过程中岩石长期强度的确定方法
机译： 三轴压缩过程中岩石长期强度的确定方法
机译： 岩石长期强度指标的确定方法
