AR模型在爆破震动信号频谱分析中的应用.pdf
第 33 卷 第 1 期 2016 年 3 月 爆 破 BLASTING Vol 33 No 1 Mar 2016 doi 10 3963/ j issn 1001 -487X 2016 01 008 AR 模型在爆破震动信号频谱分析中的应用* 孙新建 1, 2, 曾亚平1, 苏振妍1 (1 青海大学, 西宁 810016; 2 清华大学 水沙科学与水利水电工程国家重点实验室, 北京 100084) 摘 要 爆破震动信号的频谱分析以傅氏变换为基础的经典谱估计方法为主, 该方法存在方差性能差、 易 产生虚假峰值、 谱曲线起伏大等缺点。AR 模型谱估计方法计算精确度高、 且简单, 并能弥补经典谱估计在信 号分析中的不足。采用经典、 AR 模型谱估计方法对不同类型的工程实测爆破震动信号进行频谱分析, 其中 以 AR 模型谱估计方法得到爆破震动信号的频谱曲线平滑、 方差小、 主频明确, 且主频附近无虚假峰值, 表明 AR 模型现代谱估计比经典谱估计在爆破震动信号的频谱分析中更具优越性。 关键词 爆破震动信号;频谱分析;AR 模型;现代谱估计 中图分类号 TD235 3 文献标识码 A 文章编号 1001 -487X (2016) 01 -0040 -05 Application of AR Model in Spectrum Analysis of Blasting Vibration Signals SUN Xin-jian1, 2, ZENG Ya-Ping1, SU Zhen-yan1 (1 Qinghai University, Xining 810016, China; 2 State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China) Abstract The of spectrum analysis of blasting vibration signal is a classical spectral estimation mainly based on Fourier trans, which has some disadvantages, including poor perance of variance, easier to produce false peak and high curve fluctuation While AR model spectrum estimation has the advantages of high accuracy and simplicity, and makes up for the shortcomings of the classical spectrum estimation in signal analy- sis The spectrum analysis of the vibration signals of different types of projects was carried out through classical and AR model spectral estimation Results show that spectrum curve of blasting vibration signal with AR model spectrum estimation has the following characteristics such as smoothing, small variance, clear dominant fre- quency and no false peak near the dominant frequency Therefore, the AR model modern spectrum estimation is more superior to the classical spectral estimation on Spectrum analysis of blasting vibration signals Key words blasting vibration signals;frequency spectrum analysis;AR model;modern spectrum estimation 收稿日期 2015 -11 -14 作者简介 孙新建 (1976 - ) , 男, 山东鄄城, 博士后、 副教授, 主要从 事爆破震动安全研究,(E-mail) sxj000918 sina com。 基金项目 国家重点基础研究发展计划 (973) 项目 (2015CB057904) 爆破震动信号是一种典型强非线性、 非平稳的 随机信号, 爆破震动信号频谱分析是研究爆破震动 特性的重要手段之一。功率谱是信号频谱分析的主 要手段, 其反映了信号的功率能量在频域上的分布 情况, 进而揭示了信号的功率谱峰及隐含的周期性 等大量信息, 在诸多研究领域都发挥了重要作 用 [1-5]。经典谱估计与现代谱估计统称为功率谱估 计, 其中非参数化谱估计为经典谱估计, 主要包括周 期图法与自相关法两种, 频率分辨率低以及方差性 能较差是其主要缺陷。现代谱估计又称为参数化谱 估计, 主要包括 AR 模型法、 MA 模型法、 ARMA 模型 法等, 现代谱估计弥补了经典谱估计的分辨率低及 方差性能差等缺陷 [6], 现代谱估计对信号频谱分析 具有重大意义。因 AR 模型参数法简单, 目前广泛 应用在导航系统故障检测、 望远镜目标跟踪误差分 析、 结构损伤预警等多个研究领域 [7-11]。目前 AR 模型现代谱估计在爆破震动领域应用较少, 采用 AR 模型谱估计方法在频域分析爆破震动特性具有 重大意义。 1 经典谱估计 目前, 自相关法、 周期图法等经典谱估计都是以 傅氏变换为基础的无参数模型来实现信号的频谱特 征的, 由于傅氏变换不考虑随机信号的特性, 使得经 典谱估计自身存在方差性能较差、 分辨率低、 谱曲线 起伏大、 易产生虚假峰值等缺点。 1 1 自相关法 自相关法是序列 x (n) 利用维纳-辛钦定理估出 自相关函数 R (n) , 并傅氏变换后, 得到 x (n) 的功率 谱估计。 R ∧ x [n] 1 N ∑ N-1-n k 0 x [k] x [k n] (1) P ∧ x (Ω) ∑ L N -LR ∧ x [n] e-jnΩ,L ≤ N - 1 (2) 1 2 周期图法 周期图法是把对随机信号 x (n) 监测的 N 个数 据看作一个有限的能量序列, x (n) 经离散傅氏变换 后, 得 x (n), 并对其幅值进行平方, 再除上 N, 可得 x (n) 的功率谱估计。 xN[k] x [k], k 0, 1, L, N - 1 0, ﹛ 其它 (3) XN(e jΩ ) DTFT {xN[k] }∑ N-1 k 0 xN[k] e-jΩk(4) P ∧ x (Ω) 1 N XN(e jΩ )2 (5) 2 AR 模型法 现代谱估计针对经典谱估计在信号分析中的不 足, 近几年提出了 AR 模型、 MA 模型、 ARMA 模型的 现代谱估计。由于 AR 模型计算精确度高, 且计算 简单,目前 AR 模型得到广泛应用。 AR 模型可假定信号序列 x (n) 被输入序列 u (n) 所激励的线性系统的输出, 用 x (n) 或其自相 关函数 rx(m) 估计系统转移函数 H (z) 的参数, 再由 H (z) 的参数估计 x (n) 的功率谱等步骤。 x (n) 若为随机序列, 则 u (n) 为白噪声序列, 设 σ2为 u (n) 的方差为, 则 u (n) 和 x (n) 之间的关系为 x (n)∑ ∞ k 0 h (k) u (n - k)(6) 式中, h (k) 为信号系统单位抽样响应, 用差分方程 表示 x (n) 为 x (n) -∑ P k 1 akx (n - k) u (n)(7) 式中 ak为系数; P 为阶次, 对式 (6) 、 式 (7) Z 变换 后得 H (z) 1 A (z) (8) 式中 A (z) 1 ∑ P k 1 akz-k; H (z)∑ ∞ k 0 h (k) z-k。 由随机信号通过线性系统的理论可得信号 x (n) 的功率谱 Px(e jω ) Pu(e jω )H (z)2 σ2 1 ∑ake-jωk 2 (9) 式中 Px(e jω ) 、 Pu(e jω ) 分别为 x (n) 、 u (n) 的功率谱; 式中的未知系数可由 Yule-Walker 方程求得 rx(0)rx(1)rx(2)rx(P) rx(1)rx(0)rx(1)rx(P - 1) rx(2)rx(1)rx(0)rx(P - 2) rx(P) rx(P - 1) rx(P - 2) rx(0 ) 1 a1 a2 a P σ2 0 0 0 (10) AR 模型的阶次 P 对 AR 模型谱估计精确度影 响非常大, 可由 Akaike 信息理论准确求出合适的 P 值, 计算公式如下 AIC (P) ln σ2 2P/ N (11) 式中, N 为 x (n) 的长度, 当 P 从 1 开始逐次增加, 当 AIC 出现最小值时的 P 值为 AR 模型的合适阶次。 将 P 代入式 (10) 求解, 便可得到随机信号的 AR 模 型谱估计。 3 爆破震动信号频谱分析比较 爆破震动信号为典型强非线性、 非平稳的随机 信号, 为比较经典谱估计与现代谱估计在爆破震动 信号频谱分析中的性能, 用经典谱估计和 AR 模型 谱估计两种方法分别对某水电工程实测左坝肩开 挖、 右岸坝肩开挖、 地下主变室开挖、 岩锚梁开挖、 尾 水洞开挖、 地下主厂房开挖爆破震动信号 A、 B、 C、 D、 E、 F 进行频谱分析, 其中爆破震动速度信号 A F 的时程曲线如图 1 图 6 所示。 14第 33 卷 第 1 期 孙新建, 曾亚平, 苏振妍 AR 模型在爆破震动信号频谱分析中的应用 图 1 A 信号时程曲线 Fig 1 The time history curve of A signal 图 2 B 信号时程曲线 Fig 2 The time history curve of B signal 图 3 C 信号时程曲线 Fig 3 The time history curve of C signal 图 4 D 信号时程曲线 Fig 4 The time history curve of D signal 图 5 E 信号时程曲线 Fig 5 The time history curve of E signal 图 6 F 信号时程曲线 Fig 6 The time history curve of F signal 由式 (11) 对爆破震动信号 A、 B、 C、 D、 E、 F 进行 计算可得到相应的 AIC 值和 AR 模型阶次 P, 其中 P 由 1 阶取到 2000 阶。经计算 AR 模型合适阶次 P 如表 1 所示, 经式 (10) 、 式 (11) 联合求解可分别得 到各爆破震动信号的 AR 模型谱估计。另外, 由相 关法求得各爆破震动信号的经典谱估计, 为对比分 析两种谱估计, 对每个爆破震动信号采用自相关法 和 AR 模型法分别求得的功率谱绘制在一个图中, 如图 7 图 12 所示。 表 1 A F 爆破震动信号 AR 模型合适阶次 p Table 1 AR model appropriate order p of A F blasting vibration signal 信号ABCDEF P60036118857350 对图 7 图 12 进行分析研究, 可知两种谱估计 方法求得的破震动信号 A、 B、 C、 D、 E、 F 优势频率相 同, 功率谱密度变化趋势基本一致。但经典谱曲线 起伏剧烈, 方差大, 且多个虚假峰值存在于主频附 近; 而 AR 模型现代谱曲线比较平滑, 方差也较小, 且主频附近没有虚假峰值, 主频率明确。 24爆 破 2016 年 3 月 图 7 A 信号的功率谱密度 Fig 7 The power spectrum density of A signal 图 8 B 信号的功率谱密度 Fig 8 The power spectrum density of B signal 图 9 C 信号的功率谱密度 Fig 9 The power spectrum density of C signal 图 10 D 信号的功率谱密度 Fig 10 The power spectrum density of D signal 图 11 E 信号的功率谱密度 Fig 11 The power spectrum density of E signal 图 12 F 信号的功率谱密度 Fig 12 The power spectrum density of F signal 4 结语 通过对比分析经典谱估计和 AR 模型的现代谱 估计两种方法在实际工程中的应用效果, 证明 AR 模型谱估计在爆破震动信号频谱分析中有很强的适 应性, 能够弥补了经典谱估计自身存在的虚假峰值、 谱曲线起伏大等缺陷。由此, 表明 AR 模型谱估计 在爆破震动信号频谱特性分析中更具优越性。 参考文献 (References) [1] 肖 立, 米 彦, 陈 攀, 等 现代功率谱估计在绝缘 子污秽诊断中的应用 [j] 高电压技术, 2006 (2) 32- 33, 56 [1] XIAO Li, MI Yan, CHEN Pan, et al Application of mod- ern spectral estimation in diagnosis of discharge of pollu- ted insulator [j] High Voltage Engineering, 2006 (2) 32- 33, 56 (in Chinese) [2] 李文超 现代功率谱估计在设备故障诊断应用 [j] 传 感器世界, 2012 (2) 25-29 [2] LI Wen-chao Application of modern power spectral esti- mation on malfunction diagnosis of equipment [j] Signal Process and System, 2012 (2) 25-29 (in Chinese) [3] 许萍萍, 韩雪山, 何渊井 运输过程振动信号功率谱估 计的实现 [j] 技术与方法, 2013 (7) 189-191 34第 33 卷 第 1 期 孙新建, 曾亚平, 苏振妍 AR 模型在爆破震动信号频谱分析中的应用 [3] XU Ping-ping, HAN Xue-shan, HE Yuan-jing Realization of PSD estimation of vibration signals in transportation process [ j] Techniques and s, 2013 (7) 189- 191 (in Chinese) [4] 伊 鑫, 曲爱华 基于 Welch 算法的经典功率谱估计 的 Matlab 分析 [j] 现代电子技术, 2010 (3) 7-9 [4] YI Xin, QU Ai-hua Matlab simulation analysis of power specrum estimation based on welch [j] Modern Electronics Technique, 2010 (3) 7-9 (in Chinese) [5] 傅广操, 樊明捷 Matlab 在现代功率谱估计中的应用 [j] 电脑学习, 2003 (6) 6-7 [5] FU Guang-cao, FAN Ming-jie The application of matlab in modern spectral estimation[ j] Computer Study, 2003 (6) 6-7 (in Chinese) [6] 冯 磊 经典功率谱估计与现代功率谱估计的对比 [j] 商业文化 (学术版) , 2009 (5) 239, 223 [6] FENG Lei The classic power spectrum estimation and com- parison of modern power spectral estimation [j] Business Culture(Academic Version) , 2009 (5) 239, 223(in Chi- nese) [7] 金 天, 丁东方, 丛 丽 基于 AR 模型参数的导航系 统故障检测 [j] 电子设计工程, 2014 (15) 168-170 [7] jIN Tian, DING Dong-fang, CONG Li Navigation system fail- ure detection based on ar model parameters [j] Electronic Design Engineering, 2014 (15) 168-170 (in Chinese) [8] 廖宏宇, 马晓燠, 郭友明, 等 基于 AR 模型搜索迭代算 法的望远镜跟踪误差分析 [j] 物理学报, 2014 (17) 444-452 [8] LIAO Hong-yuan, MA Xiao-yu, GUO You-ming, et al A- nalysis of tracking error of telescope based on ar-search- iteration algorithm [ j] Acta Phys Sin, 2014, 17 444- 452 (in Chinese) [9] 刁延松, 任 红 基于 AR 模型和因子分析的结构损伤 预警研究 [j] 振动与冲击, 2014 (18) 115-119, 144 [9] DIAO Yan-song, REN Hong Structural damage early warn- ing based on ar model and factor analysis [j] journal of Vibration and Shock, 2014 (18) 115-119, 144(in Chi- nese) [10] 尚小晶, 吴忠伟, 徐成波 改进 AR 模型特征提取及分 类 [j] 吉林建筑工程学院学报, 2014 (3) 74-77 [10] SHANG Xiao-jing, WU Zhong-wei, XU Cheng-bo The improved AR model feature extraction and classification [j] journal of jilin Institute of Architecture and Civil Engineering, 2014 (3) 74-77 (in Chinese) [11] 郑凤霞 基于 GM (1, 1) 和 AR 模型的四川高等教育 在校学生数预测 [j] 西昌学院学报 (自然科学版) , 2014 (3) 36-37, 41 [11] ZHENG Feng-xia Research on the number of enrollment students of colleges and universities in Sichuan based on GM (1, 1)model and AR model [j] journal of Xichang College Natural Science Edition, 2014 (3) 36-37, 41 (in Chinese) (上接第 33 页) [5] LU Wen-bo, YANG jian-hua, CHEN Ming, et al Mecha- nism and equivalent numerical simulation of transient re- lease of excavation load for deep tunnel [j] Chinese jour- nal of Rock Mechanics and Engineering, 2011, 30 (6) 1089-1096 (in Chinese) [6] XIAO jian-Qing, DING De-Xin, jIANG Fu-Liang, et al Fatigue damage variable and evolution of rock subjected to cyclic loading [j] International journal of Rock Mechan- ics and Mining Sciences, 2010, 47 (3) 461-468 [7] 闫长斌, 徐国元, 杨 飞 爆破动荷载作用下围岩累积 损伤效应声波测试研究 [j] 岩土工程学报, 2007, 29 (1) 88-93 [7] YAN Chang-bin, XUE Guo-yuan, YANG Fei Measure- ment of sound waves to study cumulative damage effect on surrounding rock under blasting load [j] jouranl of Rock and Soil Engineering, 2007, 29 (1) 88-93 (in Chinese) [8] 孙金山, 周传波, 卢文波, 等 大冶金矿爆破开采巷道 围岩累积损伤规律研究 [j] 爆破, 2013, 30 (3) 10-14 [8] SUN jin-shan, ZHOU Chuan-bo, LU Wen-bo, et al Study on cumulative damage law of rock mass under blasting loading in Daye iron mine [j] Blasting, 2013, 30 (3) 10- 14 (in Chinese) [9] 徐 颖, 庄新炉, 贾 虎 井巷掘进爆破对邻近围岩的 损伤及锚杆锚固力影响的数值模拟研究 [j] 爆破, 2007, 24 (S1) 1-5 [9] XU Ying, ZHUANG Xin-lu, jIA Hu Numerical simulation of influence of shaft and adjacent surrounding rock and bolt anchorage force [j] Blasting, 2007, 24(S1) 1-5 (in Chinese) [10] 马文伟 岩石动态力学特性及钻爆掘进巷道破岩过 程分析 [D] 淮南 安徽理工大学, 2014 [10] MA Wen-wei The analysis of rock dynamic mechanical properties and rock breaking process in tunnel blasting excavation [D] Huainan Anhui University of Science and Technology, 2014 (in Chinese) [11] HOLMQUIST T j, jOHNSON G R A computational con- stitutive model for concrete subjected to larger strains, high strain rates and high pressure [C] ∥jackson N, Dickert S 14thInternational Symposium Ballistics USA American De- fense Preparedness′ Association, 1995 591 600 [12] 何满潮, 钱七虎 深部岩体力学基础 [M] 北京 科学 出版社, 2010 160-162 44爆 破 2016 年 3 月