許東溝露天礦爆破震動(dòng)波頻譜分析

楊 號(hào) 常來(lái)山 王鑫瑀2

（1.遼寧科技大學(xué)礦業(yè)工程學(xué)院，遼寧鞍山114051；2.河北鋼鐵集團(tuán)礦業(yè)有限公司，河北唐山063000）

爆破震動(dòng)信號(hào)是一種非周期性的瞬態(tài)隨機(jī)波形，從其波形中可以直觀地分析出振幅的大小、爆破震動(dòng)持續(xù)時(shí)間的長(zhǎng)短等物理量［1-2］，但是頻率的高低和分布情況卻很難得到。分析爆破震動(dòng)波的頻率特性，主要原因是較小的爆破震動(dòng)可能產(chǎn)生較強(qiáng)的共振效應(yīng)［3-5］。頻率的結(jié)構(gòu)和分布情況，對(duì)合理選擇爆破藥量、爆破方法和確定爆破安全距離提供了科學(xué)依據(jù)［6-8］，因此在爆破震動(dòng)分析中，獲得爆破震動(dòng)波的主要頻率成分是十分必要的［9］。

根據(jù)許東溝露天礦現(xiàn)場(chǎng)爆破震動(dòng)實(shí)測(cè)數(shù)據(jù)，基于頻譜分析FFT原理，利用Matlab軟件中的頻譜分析功能，繪出震動(dòng)信號(hào)的傅里葉頻譜圖和功率譜圖，對(duì)許東溝露天礦爆破震動(dòng)信號(hào)的頻譜、功率譜分布特性進(jìn)行了分析研究。

1 震動(dòng)信號(hào)頻譜分析

1.1 傅里葉快速變換

傅里葉變換可將爆破震動(dòng)信號(hào)由時(shí)域轉(zhuǎn)換為頻域，以分析頻率相關(guān)的信息。函數(shù)a（t）在R2上的連續(xù)傅里葉變換［10］定義為

爆破震動(dòng)信號(hào)為離散信號(hào)，需要對(duì)等時(shí)間間隔的信號(hào)進(jìn)行積分以獲得離散傅里葉變換（DFT）［11］：

逆變換公式為

FFT是利用計(jì)算機(jī)將DFT快速實(shí)現(xiàn)的方法，它巧妙地將震動(dòng)信號(hào)分解成若干頻率正弦波的疊加。

1.2 爆破震動(dòng)信號(hào)的譜密度分析

爆破震動(dòng)信號(hào)有很大的隨機(jī)性，不適合作幅值譜、相位譜及離散譜分析，往往只做功率譜分析。功率譜密度函數(shù)可定義為

功率譜密度函數(shù)在實(shí)際工程中用單邊功率譜表示，即

2 爆破震動(dòng)現(xiàn)場(chǎng)測(cè)試分析

2.1 工程概況與測(cè)點(diǎn)布置

許東溝采場(chǎng)位于遼寧省鞍山市東南15～20 km。爆區(qū)位于采場(chǎng)南端，巖性為磁鐵石英巖，爆區(qū)呈塊狀分布，寬32.79 m，長(zhǎng)74.61 m，鉆孔采用45R牙輪鉆機(jī)，孔徑250 mm，孔深15 m，孔距7 m，排距7 m，采用的炸藥為銨油乳化炸藥，連續(xù)柱狀裝藥，逐孔微差起爆。炮孔布置與起爆網(wǎng)絡(luò)如圖1?，F(xiàn)場(chǎng)測(cè)試采用IDTS-3850型爆破震動(dòng)儀，爆區(qū)與測(cè)點(diǎn)布置如圖2所示。

2.2 爆破震動(dòng)實(shí)測(cè)數(shù)據(jù)

根據(jù)現(xiàn)場(chǎng)爆破試驗(yàn)，用爆破震動(dòng)儀記錄許東溝采場(chǎng)爆區(qū)的數(shù)據(jù)匯總?cè)绫?所示。

從表1可以看出，爆破震動(dòng)峰值速度總體呈衰減趨勢(shì)，垂向相對(duì)于徑向衰減較快，隨著爆心距的增加，其垂直方向的爆破震動(dòng)速度將小于徑向爆破震動(dòng)的數(shù)值，主振方向有向徑向轉(zhuǎn)化的趨勢(shì)；爆破震動(dòng)信號(hào)的徑向主頻分量大于垂直方向，其主頻具有向低頻發(fā)展的趨勢(shì)，且徑向較為明顯。

2.3 頻譜分析

利用傅里葉變換對(duì)采集的爆破震動(dòng)信號(hào)進(jìn)行頻譜分析時(shí)，在Matlab中編制相應(yīng)的信號(hào)處理和分析程序，得到該爆破震動(dòng)實(shí)測(cè)信號(hào)的頻譜圖和功率譜密度圖?；贛atlab的FFT設(shè)計(jì)流程圖如圖3所示，圖4為爆破震動(dòng)信號(hào)各測(cè)點(diǎn)的垂向震動(dòng)頻譜圖，圖5為各測(cè)點(diǎn)的徑向震動(dòng)頻譜圖（測(cè)點(diǎn)1徑向測(cè)振失?。?/p>

爆破震動(dòng)信號(hào)波形曲線的復(fù)雜程度可根據(jù)它的頻率成分體現(xiàn)出來(lái)（如圖4、圖5），礦山爆破震動(dòng)信號(hào)具有較寬的頻帶分布，常常存在多峰值現(xiàn)象，主頻并不突出。垂向震動(dòng)的主頻主要分布在低中頻區(qū)域，范圍在8.5～54.9 Hz，主要集中在12～33 Hz；徑向震動(dòng)的主頻主要分布在低頻區(qū)域，范圍在6.1～23.2 Hz，主要集中在12～20 Hz。垂向震動(dòng)主頻相對(duì)于徑向主頻頻帶較寬，但其幅值相對(duì)較小。

2.4 功率譜分析

由Matlab頻譜分析功能繪制的爆區(qū)各測(cè)點(diǎn)的垂向震動(dòng)功率譜如圖6所示，圖7為各測(cè)點(diǎn)的徑向震動(dòng)功率譜圖。

礦山爆破震動(dòng)具有典型的非平穩(wěn)隨機(jī)特征。由圖6和圖7的功率譜圖可以看出，爆破震動(dòng)能量在頻域上雖然展布比較廣泛，但大部分能量都集中在低中頻范圍內(nèi)，低中頻震動(dòng)信號(hào)在徑向方向上的能量高于垂直分量的能量。傳播時(shí)能量衰減較快，且高頻快于低頻。

3 結(jié)論

（1）許東溝露天礦爆破震動(dòng)峰值速度總體呈衰減趨勢(shì)，隨著爆心距的增加，高頻信號(hào)衰減較快，震動(dòng)主頻有向低頻區(qū)域發(fā)展的趨勢(shì)。

（2）在爆破震動(dòng)信號(hào)的頻譜分析中，爆破震動(dòng)信號(hào)具有一定的頻帶寬度，有時(shí)主頻并不突出，包含多個(gè)子頻；垂向震動(dòng)主頻相對(duì)于徑向主頻頻帶較寬，垂向震動(dòng)主頻主要集中在12～33 Hz，徑向震動(dòng)主頻主要集中在12～20 Hz。