利用窗函数进行爆破振动功率谱的平滑化.pdfVIP

2008年 12月 利用窗函数进行爆破振动功率谱的平滑化 李洪涛等 利用窗函数进行爆破振动功率谱的平滑化’ 李洪涛 杨兴国 四川大学水利水电学院(四川成都，610065) [摘 要] 文章针对爆破振动频谱分析中谱峰值的准确位置不易确定，频带能量分布特征难以直观显现的问题， 利用信号处理中的窗函数，采用滑动平均法对爆破振动功率谱进行平滑化，实现了爆破振动频带范围和频带能量 分布的直观显现。进一步研究表明，在对功率谱进行平滑化的过程中，可以采用关窗法来进行窗函数带宽的选择， 保证既不丢掉功率谱原有的趋势特点，又能很好地将重要的峰点凸现出来。 [关键词] 爆破振动 功率谱 平滑 窗函数 Matlab [分类号] TIY235．1 1 引言 就等于频率 [f—df／2，厂+df／2】范围内的均方值， 爆破地震效应及其对建筑物影响问题是工程爆 ( 即表示随机信号 丁)在每单位频带上分量的 破的重要课题之一。建(构)筑物的爆破振动破坏 均方值。在一般意义上，功和能量与振幅的平方成 与爆破地震波的频率密切相关，不同频率成分的振 正比例，所以谱密度 ( 可以看作随机信号在单 动对于建(构)筑物的影响是不同的。基于此，很多 位频带内的谐波分量的能量按频率厂分布的度量， 学者和研究人员开展 了爆破振动频率特性 的研 将其称为功率谱密度函数 。 究 ¨ ；《爆破安全规程》(GB6722-2003)中关于地 目前爆破振动测试普遍采用了数字式记录，从 面建筑爆破振动安全控制标准，也考虑 了爆破振动 爆破振动记录文件中提取出的记录数据，一般是具 频率的影响6‘J。对于爆破振动频率特性研究，最基 有一定时间间隔的离散信号，采用离散傅立叶变换 本的方法就是对爆破振动波形进行频谱分析，求得 可以对爆破振动进行频谱分析。 傅立叶或功率谱。但是由于爆破振动是一种典型的 在Matlab中提供 了快速傅立叶变换的工具箱 非稳态随机过程 ，采用频谱分析方法求得的谱分布 函数，只需通过比较简单的编程，即可以实现爆破振 往往呈现锯齿状，在进行频率特性研究时，不仅谱峰 动的功率谱分析。本文后面的功率谱平滑化过程 ， 值的准确位置不易确定，而且对于频带分布范围，以 也是利用Matlab语言编制计算程序实现的。 及各频带内的能量分布，更是难以直观显现。本文 下面将以一个爆破振动实测资料来说明爆破振 利用信号处理中的窗函数，采用滑动平均方法对爆 动频谱分析过程，以及爆破振动功率谱的特点。 破振动功率谱进行平滑化，取得了良好的效果。 图1为丹江口大坝加高右岸土建施工某次爆破 2 爆破振动功率谱分析 过程中，实测竖直 向速度历程曲线，及利用Matlab 在爆破振动频谱分析中，用得比较多的是功率 编程计算得到的功率谱图。爆破振动信号采样频率 谱。功率谱函数可分为 自功率谱和互功率谱，在爆 为5000Hz，开挖爆破采用手风钻钻孔，孔径 42 破振动谱分析中，主要应用的是 自功率谱函数。 m131，孔深2．0m，炸药使用032mm乳化炸药药卷。 自相关函数间接地给出了包含在随机过程中的 从图1中的功率谱图可以看出，爆破振动的频 频率信息，根据自相关函数的傅立叶变换及其逆变 率成分是非常丰富的，其功率谱的突峰比较多，在很 换，并依据 自相关函数的基本定义，可以得到 】： 小的频带范围内，功率谱密度也存在很大的突变性， E ()]=』 ∽ df