岩石破坏时频分析及瞬时频率前兆探析.docVIP

岩石破坏时频分析及瞬时频率前兆探析 波形分析是指以声发射信号的时域波形为对象，选用一定的信号处理方法来获取声发射源信息。声发射波形信号中蕴含着大量声发射源的信息。因此，近年来，波形信号分析成为了声发射技术研究的新热点。声发射波形分析方法主要有基于傅里叶变换的频谱分析和时频分析两种方法。在频谱分析方面最突出的研究成果，是在室内实验[13-14]与现场监测[15-16]中发现了不同尺度岩石破坏时出现的声发射信号由高频向低频移动的“频移”现象。傅里叶变换的不足之处是信号的“整体变换”，即不能揭示频率分量随时间变化的情况。因此，只适用于平稳信号的分析[17]。时频分析(time-frequencyanalysis)的基本思想是设计时间和频率的联合函数，用于描述信号在不同时刻的谱分量[18]。时频分析将一维时域信号变换为二维时频图像，能够展现出更为丰富的信号内涵。时频分析包括短时傅里叶变换和小波变换两类线性时频变换，和以维格纳变换为代表的双线性时频变换[18]。小波分析是声发信号多尺度分析的有力工具[19-20]。然而，小波变换给出的是时间ndash;尺度(或分辨率)之间的关系[19]，对于要求频率信息的场合，时频变换具有不可替代的作用。时频变换在岩石破坏前兆信号分析方面有诸多应用：例如，对岩石试件多轴加载的声发射信号进行分析，发现低应力水平时声发射信号为高频小幅值、岩爆或地震前兆为低频大幅值。 与数字信号处理的其他领域相比，岩石破坏声发射信号时频分析的研究还很不充分，现有的岩石破坏频域前兆包含的信息仅有时间、频率和幅度，缺少动力学信息的描述；二次型时频分布的“交叉项”识别以及频率的精确定位缺少综合性的方法，因而与实际应用的要求还有不小的距离，其主要问题是：(1)大数据：即声发射原始波形的采样为海量数据，必须分割后才能进行处理；对大数据的分割的研究少有报道。(2)对声发射信号滤波及时频分析算法的研究不够充分。(3)对时频图像的进一步处理与分析的研究，尚不能达到细致描述岩石破坏动态过程的要求。本文以花岗岩单轴压缩声发射数据为实例，开展大数据样本的分割理论、时频分析窗函数滤波、时频分析算法、和以时ndash;频联合分布、谱分量分布模式表征动态过程为内涵的，多维度瞬时频率前兆信息的研究，以期为岩石破坏声发射信号前兆信息的识别提供原理上的支持。 实验系统及原理 岩石试样为山东莱州花岗岩，呈浅肉红色，粗粒块状结构，肉眼观察试件致密、无裂纹(见图1(b))。花岗岩是典型的硬岩，由于工程开挖，在高应力作用下，会产生突然猛烈的破坏。将花岗岩块首先用岩石切割机进行切割，再用双面磨石机磨平，制成规格为150mmtimes;60mmtimes;30mm(高times;长times;宽)的长方形试样，试样两侧安装有2个声发射传感器，直接固定在试样表面(见图1(a))。矿物成分的X射线衍射(X-raydiffraction，XRD)分析表明，花岗岩中晶体矿物主要为石英(27.0%)、钾长石(37.0%)、钭长石(31.0%)和黏土矿物(5.0%)。扫描电镜(scanningelectronicmicroscope，SEM)观察分析表明：花岗岩是非均匀、多孔性介质；石英、长石、云母接触边界明显可见；石英、长石中存在着大量的微米级孔洞。岩石的非均匀性源自于其中存在着不同硬度的成分、存在着微观孔洞与不同成分的边界。当岩石加载时，这些微观缺陷会扩展、传播、最终汇合成宏观裂纹并导致岩石的完全破坏。实验是在中国矿业大学(北京)煤炭资源与安全开采国家重点实验室进行，实验系统如图2所示。花岗岩试件的单轴加载实验是在EHFndash;UG500全数字液压伺服试验机进行，最大静态载荷为759kN，荷载精度小于0.5%，采用轴向位移控制，加载速率为0.002mm/s，自动记录实验过程中的力和位移数据。在对花岗岩样的加载过程中，同步进行声发射探测与红外探测。本文主要研究声发射信号的处理。声发射探测采用中国矿业大学(北京)深部岩土与地下工程国家重点实验室的PXWAE声发射监测系统。系统由传感器、放大器、信号采集处理和记录与显示(采集卡和微机)组成(见图3[2])，最高采样速率为20MHz，采样精度12位。采集卡与计算机的数据通讯采用PCI总线，可实现速率达132MB/s的波形信号传输。传感器接收从试件内发出的波动信号，使声发射传感器压电陶瓷变形，电压发生变化。这种变化的电压信号通过前置放大器放大，再经过信号调制变成稳定的声发射信号存储至微机中并实时显示波形等参数。实验采用的声发射传感器的共振频率为150kHz，前置放大器40dB，增益设置为10倍。实验中，系统的采样率设置为1MHz。 声发射信号分割 1.信号分割理论 图4给出了一个拉长了声发射信号随时间t