岩体失稳声发射前兆信号提取方法及失稳前兆判别方法技术

本发明专利技术公开一种岩体失稳声发射前兆信号提取方法及失稳前兆判别方法，岩体失稳声发射前兆信号提取方法包括：S1：实时探测并记录当前监测岩体声发射信号的相关参数；S2：将相关参数的时间序列按照预设窗口长度嵌入得到轨迹矩阵；S3：将轨迹矩阵进行SVD分解，得到分解后的轨迹矩阵；S4：计算分解后的轨迹矩阵中特征值的贡献率，并以降序的形式进行排列，得到排列后的特征值贡献率；S5：根据排列后的特征值贡献率的峰度及峰度的均值提取前兆特征值，得到前兆特征值集合；S6：根据前兆特征值集合，重构前兆特征值对应的矩阵，得到前兆特征信号。本发明专利技术能够对岩体状态进行实时跟踪预警预报，且所提取的前兆信号易于识别，判断方式简单。单。单。

【技术实现步骤摘要】
岩体失稳声发射前兆信号提取方法及失稳前兆判别方法


[0001]本专利技术涉及岩土工程监测预警
，具体涉及一种岩体失稳声发射前兆信号提取方法及失稳前兆判别方法。

技术介绍

[0002]随着工程建设不断向西藏等高地应力、构造运动活跃的地区推进，将会面临一系列与岩体稳定性相关的工程问题，如大变形、岩爆、崩塌、滑坡等。因此，实时、准确的对岩体稳定性进行监测、对其失稳破坏进行超前预报、预警，对保障人民生命安全和避免财产损失有着重要意义。
[0003]岩体因自身脆性往往失稳破坏具有突发性、变形量较小等特点，导致将变形作为前兆指标可靠性非常差，前兆不明显。声发射是指材料内部或表面的局部应变能快速释放产生的瞬态弹性波，反映了岩体实时状态的变化，目前大量研究将声发射频率、b值、分形维数，RA/AF、等作为前兆指标，但这些现有技术所采用的前兆指标的波动较大，前兆信号中往往包含较多噪声信号，难以在实际工程中应用。

技术实现思路

[0004]本专利技术的目的在于提供一种岩体失稳声发射前兆信号提取方法及失稳前兆判别方法，以能够对岩体状态进行实时跟踪预警预报，且所提取的前兆信号易于识别，判断方式简单。
[0005]本专利技术解决上述技术问题的技术方案如下：
[0006]本专利技术提供一种岩体失稳声发射前兆信号提取方法，所述岩体失稳声发射前兆信号提取方法包括：
[0007]S1：实时探测并记录当前监测岩体声发射信号的相关参数；
[0008]S2：将所述相关参数的时间序列按照预设窗口长度嵌入得到轨迹矩阵；
[0009]S3：将所述轨迹矩阵进行SVD分解，得到分解后的轨迹矩阵；
[0010]S4：计算所述分解后的轨迹矩阵中特征值的贡献率，并以降序的形式进行排列，得到排列后的特征值贡献率；
[0011]S5：根据排列后的特征值贡献率的峰度及峰度的均值提取前兆特征值，得到前兆特征值集合；
[0012]S6：根据所述前兆特征值集合，重构所述前兆特征值对应的矩阵，得到前兆特征信号。
[0013]可选择地，所述步骤S1中，所述相关参数包括声发射上升时间/幅值参数和平均频率参数。
[0014]可选择地，所述步骤S2中，所述预设窗口长度L满足：2≤L≤N
‑
1，其中，N表示为时间序列数据长度。
[0015]可选择地，所述步骤S2中，所述轨迹矩阵X为：
[0016][0017]其中，K为轨迹矩阵的列数，且K＝N
‑
L+1，N表示为时间序列数据长度，L为窗口长度参数。
[0018]可选择地，所述步骤S3包括：
[0019]S31：计算XX
T
的初始特征值并将所述初始特征值按递减的顺序排列，得到特征值；
[0020]S32：获取轨迹矩阵X的左特征向量和右特征向量；
[0021]S33：根据所述轨迹矩阵X的左特征向量和右特征向量，以及所述特征值，得到所述分解后的轨迹矩阵。
[0022]可选择地，所述分解后的轨迹矩阵为：
[0023][0024]其中，d为非零特征值的个数，且d＝max{i,λ
i
>0}＝rank X，为奇异值，λ
i
(i＝1,2,
…
,L)为矩阵S＝XX
T
按降序排列后的特征值，U
i
和V
i
分别为轨迹矩阵X的左特征向量和右特征向量。
[0025]可选择地，所述步骤S4中，所述特征值的贡献率η
i
为：
[0026][0027]其中，λ
i
(i＝1,2,
…
,L)为矩阵S＝XX
T
按降序排列后的特征值，d为非零特征值的个数，且d＝max{i,λ
i
>0}＝rank X，X为轨迹矩阵。
[0028]可选择地，所述步骤S4中，所述特征值贡献率的峰度Ku为：
[0029][0030]其中，λ
i
(i＝1,2,
…
,L)为矩阵S＝XX
T
按降序排列后的特征值，d为非零特征值的个数，且d＝max{i,λ
i
>0}＝rank X，X为轨迹矩阵，为特征值贡献率的均值。
[0031]可选择地，所述步骤S6即为将每一个前兆特征值重构为长度为N的新时间序列。
[0032]本专利技术还提供一种岩体失稳前兆判别方法，所述岩体失稳前兆判别方法包括上述的岩体失稳声发射前兆信号提取方法，还包括：
[0033]利用以下失稳前兆条件之一对所述前兆特征信号进行判别，得到判别结果：
[0034]当平均频率参数产生突降和/或当声发射上升时间/幅值参数产生突增时，输出所述前兆特征信号有效，否则，输出所述前兆特征信号无效。
[0035]本专利技术具有以下有益效果：
[0036]声发射波形可以由声发射监测设备直接测出，而声发射上升时间参数可以从声发
射波形图中直接获取，获取非常方便快捷，难度低；声发射上升时间的变异系数、偏度和峰度的计算过程以及后续判断非常简单，在常规的处理器中嵌入计算程序即可实现，无需构建复杂的算法模型，岩体失稳破坏的前兆点易于识别，并且计算误差小，稳定性高，一般的工程施工单位均具备应用条件。
[0037]此外，本专利技术的前兆指标多元、且与岩土破裂类型紧密相关，可提高预警的准确性。
[0038]综上，本专利技术可实时表征岩体破坏类型(裂纹类型)的发展趋势，可在隧道工程、边坡工程、采矿工程、水利水电工程等领域中的岩体预警预报中广泛使用。
附图说明
[0039]图1为本专利技术岩体失稳声发射前兆信号提取方法的流程图；
[0040]图2为单轴荷载作用下含预制裂缝岩体失稳破坏过程中采集到的声发射AF曲线图；
[0041]图3为单轴荷载作用下含预制裂缝岩体失稳破坏过程中采集到的声发射RA曲线图；
[0042]图4为声发射前兆特征值提取示意图；
[0043]图5为所提取的声发射AF、RA前兆信号及前兆点判别图。
具体实施方式
[0044]以下结合附图对本专利技术的原理和特征进行描述，所举实例只用于解释本专利技术，并非用于限定本专利技术的范围。
[0045]本专利技术提供一种岩体失稳声发射前兆信号提取方法，参考图1所示，所述岩体失稳声发射前兆信号提取方法包括：
[0046]S1：实时探测并记录当前监测岩体声发射信号的相关参数；
[0047]本专利技术对当前监测岩体声发射信号的相关参数的实时探测设备不做具体限制，作为一种实施方式，在本专利技术中，采用声发射监测设备获取相关参数，其中，相关参数包括声发射上升时间/幅值参数RF和平均频率参数AF。安装固定好声发射探头后即可监测岩体的声发射信号，并获取RA、AF参数，操作相对便捷。
[0048]S2：将所述相关参数的时间本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种岩体失稳声发射前兆信号提取方法，其特征在于，所述岩体失稳声发射前兆信号提取方法包括：S1：实时探测并记录当前监测岩体声发射信号的相关参数；S2：将所述相关参数的时间序列按照预设窗口长度嵌入得到轨迹矩阵；S3：将所述轨迹矩阵进行SVD分解，得到分解后的轨迹矩阵；S4：计算所述分解后的轨迹矩阵中特征值的贡献率，并以降序的形式进行排列，得到排列后的特征值贡献率；S5：根据排列后的特征值贡献率的峰度及峰度的均值提取前兆特征值，得到前兆特征值集合；S6：根据所述前兆特征值集合，重构所述前兆特征值对应的矩阵，得到前兆特征信号。2.根据权利要求1所述的岩体失稳声发射前兆信号提取方法，其特征在于，所述步骤S1中，所述相关参数包括声发射上升时间/幅值参数和平均频率参数。3.根据权利要求1所述的岩体失稳声发射前兆信号提取方法，其特征在于，所述步骤S2中，所述预设窗口长度L满足：2≤L≤N
‑
1，其中，N表示为时间序列数据长度。4.根据权利要求1所述的岩体失稳声发射前兆信号提取方法，其特征在于，所述步骤S2中，所述轨迹矩阵X为：其中，K为轨迹矩阵的列数，且K＝N
‑
L+1，N表示为时间序列数据长度，L为窗口长度参数。5.根据权利要求4所述的岩体失稳声发射前兆信号提取方法，其特征在于，所述步骤S3包括：S31：计算XX
T
的初始特征值并将所述初始特征值按递减的顺序排列，得到特征值；S32：获取轨迹矩阵X的左特征向量和右特征向量；S33：根据所述轨迹矩阵X的左特征向量和右特征向量，以及所述特征值，得到所述分解后的轨迹矩阵。6.根据权利要求1或5所述的岩体失稳声发射前兆信号提取方法，其特征在于，所述分解后的轨迹矩阵为：其中，d为非零特征值个数，且d＝max{i,λ
i
>0}＝rank X...

【专利技术属性】
技术研发人员：许强，朱星，唐垚，樊柱军，段毅，杨生彬，高晓，贺鸣，冯军伟，
申请(专利权)人：中国电力工程顾问集团西北电力设计院有限公司，
类型：发明
国别省市：