3.3 計(jì)算機(jī)仿真實(shí)例
為驗(yàn)證小波變換方法在巖石聲波信號(hào)分析中的
有效性,確定最優(yōu)小波基與分析尺度,利用本次巖
芯測(cè)試探頭對(duì)零激發(fā)信號(hào)(圖1),采用不同類(lèi)型小波
基對(duì)其進(jìn)行不同尺度下的小波變換與重構(gòu),進(jìn)而將
重構(gòu)信號(hào)與原始信號(hào)進(jìn)行誤差對(duì)比分析,最終確定
誤差最小的db4 小波為最優(yōu)小波基,進(jìn)行8 尺度小
波分析(詳細(xì)過(guò)程見(jiàn)作者[3]的研究),分析結(jié)果如圖3
所示,由于屬高頻干擾,故圖中未列出第1,2 層小
波分量。
由圖3 可見(jiàn),聲波信號(hào)通過(guò)小波變換后被分解
為不同頻帶通道的信號(hào),彌補(bǔ)了傳統(tǒng)的Fourier 變換
時(shí)窗口大小不隨頻率變化而改變的缺陷,從而可以
分別對(duì)每個(gè)頻帶通道的信號(hào)進(jìn)行獨(dú)立分析。
(a) D3(時(shí)域) (b) D3(頻域)
(c) D4(時(shí)域) (d) D4(頻域)
A
A/104
A
A/104
第 26 卷增1 俞縉,等. 基于小波變換的巖芯卸荷擾動(dòng)聲學(xué)反應(yīng)分析 . 3561 .
A/10 4
AA
A
(e) D5( 時(shí)域) (f) D5( 頻域)
A
A/10 4
(g) D6( 時(shí)域) (h) D6( 頻域)
A/10 4
(i) D7( 時(shí)域) (j) D7( 頻域)
A/10 4
600 200 400
f/kHz
(k) D8( 時(shí)域) (l) D8( 頻域)
圖3 超聲激發(fā)信號(hào)的小波變換
Fig.3 Wavelet transformation of acoustic exciting wave signal of transducer
表2 XND16 孔巖芯特征和測(cè)試結(jié)果
Table 2 Characters of rock cores and testing results in
borehole XND16
巖芯卸荷擾動(dòng)聲學(xué)反應(yīng)分析
本次研究對(duì)巖芯進(jìn)行了大量的聲波信號(hào)小波分
析,現(xiàn)以XND16 孔(角礫巖)為例具體分析。該孔巖
芯特征與測(cè)試結(jié)果見(jiàn)表2。
運(yùn)用db4 小波基,對(duì)各巖芯聲波信號(hào)進(jìn)行8 尺
度小波分析。因結(jié)果圖數(shù)量很大,本文僅列出第5
層分析結(jié)果,如圖4 所示。
由表2、圖4 可以看出,XND16–4,XND16–7
樣品強(qiáng)度較低,波速與巖體測(cè)井波速相比有明顯降
低。其聲波信號(hào)的第5 層小波分量波動(dòng)幅度明顯低
于其余巖芯,頻譜上只在200 kHz 上有微弱起伏,350 kHz 處之后的波動(dòng)已經(jīng)巖芯濾波作用衰減干凈。
觀察巖芯,XND16–4 樣品含兩組裂隙,夾角10°~
50°,微閉~微張,局部發(fā)育小溶孔;XND16–7
樣品含一組緩傾裂隙,泥質(zhì)膠結(jié)。其余巖芯強(qiáng)度、
- 巖石聲波參數(shù)測(cè)試儀
HS-YS4A
- 全文觀看
基于小波變換的巖芯卸荷擾動(dòng)聲學(xué)反應(yīng)分析 另存為下載 文章應(yīng)用本所DB系列
|