孫 蒙，侯興民，張一林

(煙臺(tái)大學(xué) 土木工程學(xué)院， 山東 煙臺(tái) 264005)

巖塊波速不僅能用于求得巖體完整性系數(shù)進(jìn)而評(píng)價(jià)巖體完整程度，還能綜合反映巖石的結(jié)構(gòu)和物理性質(zhì)?！豆こ處r體試驗(yàn)方法標(biāo)準(zhǔn)》[1](GB/T 50266—2013)(以下簡(jiǎn)稱“標(biāo)準(zhǔn)”)中提到的巖體波速測(cè)試方法很多，如地震波法、超聲波法等，而規(guī)定巖塊的波速測(cè)試只有超聲波法，對(duì)于振動(dòng)法測(cè)巖塊波速，一直缺乏相關(guān)研究。

地震波和超聲波兩種彈性波在巖石介質(zhì)中傳播時(shí)，其波速存在一定的差異，超聲波速度高于地震波，主要是由于兩種彈性波自身頻率不同造成的[2-3]。工程巖體波速測(cè)試采用地震波法，若巖塊波速測(cè)試仍采用聲波測(cè)試，所得到的巖體完整性系數(shù)具有一定誤差，不能準(zhǔn)確地判定巖體完整程度，根據(jù)同一種彈性波激發(fā)所測(cè)出的波速得到的巖體完整性系數(shù)才較為精確[4]。一直以來，巖塊波速都是采用超聲波法測(cè)試，國(guó)內(nèi)外學(xué)者對(duì)其研究較多，也取得了很多成果。

Mckenzie等[5]于1982年研究了超聲波測(cè)試技術(shù)的可行性和儀器規(guī)格，利用從測(cè)試或鉆孔中獲得的聲波速度和衰減參數(shù)Q(通過測(cè)量傳播過程中脈沖形狀的變化來估算)來表征巖體完整性。李成波等[6]研究了巖石加卸載時(shí)超聲波速與各彈性參數(shù)之間的變化規(guī)律，通過實(shí)測(cè)與理論推導(dǎo)的對(duì)比以及對(duì)該實(shí)驗(yàn)的回歸分析，得出波速與泊松比等彈性參數(shù)有密切的關(guān)系。陳旭等[7]利用超聲波法測(cè)定了三巖樣的縱波波速及波形的時(shí)域特征，從而反映了巖石內(nèi)部的情況。邵長(zhǎng)云等[8]也是利用超聲波法對(duì)巖樣進(jìn)行室內(nèi)波速測(cè)試，利用測(cè)定的波速比來劃分巖石的風(fēng)化程度。由此看出巖塊縱波波速在表征巖石損傷程度、綜合評(píng)估巖石脆裂性、巖塊物理力學(xué)性質(zhì)判斷和巖石質(zhì)量指標(biāo)評(píng)價(jià)等工作中也得到了廣泛應(yīng)用，因此巖塊縱波波速測(cè)試的準(zhǔn)確測(cè)量具有重要的意義。許強(qiáng)[9]則用地震波法和超聲波法對(duì)巖體波速進(jìn)行測(cè)試，對(duì)采用兩種方法測(cè)得的數(shù)據(jù)對(duì)比研究差異規(guī)律，提出了將聲波速度非線性折減為地震波速度計(jì)算巖體完整性系數(shù)。吳振祥等[10]對(duì)砂巖進(jìn)行循環(huán)荷載下的疲勞加載，同時(shí)進(jìn)行超聲波速測(cè)量，并對(duì)提取的超聲波信號(hào)進(jìn)行小波變換，得到了波譜參數(shù)的敏感變化。李廣場(chǎng)[11]通過多組巖塊聲波測(cè)試以及地震波測(cè)試巖體波速，地震波明顯小于聲波波速，且地震波速度變化幅度高于聲波。

鑒于此，本文提出了利用振動(dòng)測(cè)試的方法來對(duì)巖塊縱波波速進(jìn)行測(cè)試。振動(dòng)測(cè)試法克服了因超聲波頻率不同所引起的巖塊波速差異以及超聲波受干擾和試驗(yàn)條件限制的缺陷，將振動(dòng)傳感器直接粘到巖塊上，再配合上測(cè)試設(shè)備測(cè)出波速，通過對(duì)比超聲波測(cè)試結(jié)果，表明本文方法可以精確地測(cè)得巖塊縱波波速。

1 巖塊縱波波速測(cè)試原理及數(shù)據(jù)處理方法

目前“標(biāo)準(zhǔn)”中給出的巖塊波速測(cè)試的方法只有超聲波法，但是該方法因波動(dòng)走時(shí)難以準(zhǔn)確獲取，所得到的波速會(huì)有一定誤差。巖塊是由各礦物構(gòu)成的含節(jié)理、裂隙的非均勻介質(zhì)，由于頻散效應(yīng)，彈性波在巖塊中的波速大小與其頻率相關(guān)。對(duì)于同一巖塊，激發(fā)不同頻率的彈性波，所測(cè)得的縱波波速有一定差異。超聲波測(cè)試得到的波速略高，一般高于地震波5%～15%左右，主要由彈性波的頻散特性所致[2-3]。

而大多數(shù)實(shí)驗(yàn)室內(nèi)進(jìn)行測(cè)試的巖塊尺寸較小，相比于巖體，巖塊中所含有的節(jié)理微裂隙等結(jié)構(gòu)面發(fā)育程度較低，可近似將其看作均勻彈性介質(zhì)看來考慮。振動(dòng)波與超聲波都屬于彈性波，其本質(zhì)都是機(jī)械振動(dòng)在介質(zhì)中的傳播過程，根據(jù)彈性波理論，均勻介質(zhì)的縱波波速可由該表達(dá)式求得：

(1)

式中:E為介質(zhì)的彈性模量；υ為泊松比，ρ為介質(zhì)密度。除了上述參數(shù)外，巖塊的巖性、結(jié)構(gòu)面發(fā)育程度、孔隙度、密度、軟硬程度等不同，所得的波速就不同。

無論是用超聲波法還是本文提出的振動(dòng)法測(cè)量巖塊縱波波速，都是激發(fā)巖塊內(nèi)的直達(dá)波來測(cè)得縱波波速。使用超聲波法測(cè)量巖塊縱波波速時(shí)，由于其頻率較高，對(duì)波速的影響不能忽略。黃世強(qiáng)等通過實(shí)驗(yàn)室測(cè)試發(fā)現(xiàn)了聲波頻率從25 kHz增加到1 000 kHz時(shí)測(cè)得的玄武巖巖塊縱波波速由4 700 m/s增大到了6 000 m/s，呈明顯的頻散特征[12]。除此之外，換能器的耦合、波動(dòng)走時(shí)的準(zhǔn)確獲取、防干擾措施不足等都會(huì)影響測(cè)試結(jié)果。巖塊尺寸較小，波程較短，由于存在以上原因，超聲波法在測(cè)量中很容易出現(xiàn)較大誤差。對(duì)此，采用本文提出的振動(dòng)法測(cè)試巖塊縱波波速，各種實(shí)驗(yàn)條件布置得合理精細(xì)，在一定程度上使測(cè)試結(jié)果更加精確。

1.1 振動(dòng)法測(cè)試的基本原理

本文提出的振動(dòng)測(cè)試方法與超聲波法不同的是，通過敲擊巖塊，激發(fā)振動(dòng)波，在巖塊兩端橫截面中心布設(shè)的高通頻帶振動(dòng)加速度傳感器，采集高頻振動(dòng)信號(hào)，再根據(jù)敲擊點(diǎn)與傳感器之間的距離關(guān)系對(duì)所求得的波速進(jìn)行修正，確保所得的巖塊縱波波速準(zhǔn)確。振動(dòng)法測(cè)試示意圖如圖1所示。

圖1 振動(dòng)法測(cè)試示意圖

如圖1所示，試驗(yàn)設(shè)備主要包含了夾緊系統(tǒng)、信號(hào)采集系統(tǒng)、數(shù)據(jù)處理系統(tǒng)和激振系統(tǒng)。振動(dòng)加速度傳感器質(zhì)量輕、體積小，便于與巖塊試樣粘貼，同時(shí)其靈敏度為100 mv/g，通頻帶為0.5 kHz～10 kHz，能準(zhǔn)確清晰地采集巖塊高頻振動(dòng)信號(hào)。夾緊系統(tǒng)能有效地將巖塊固定住，且與巖塊試樣的接觸面積較小，能減少對(duì)測(cè)試結(jié)果的干擾。激發(fā)系統(tǒng)為BY-L200型力錘，用力錘在巖塊橫截面處沿垂直巖塊橫截面方向敲擊激振。由于巖塊尺寸較小，若采樣頻率設(shè)置不合理，采集的信號(hào)不能準(zhǔn)確反映原始信號(hào)，使后續(xù)數(shù)據(jù)處理中識(shí)別的到時(shí)差精度不夠，影響波速計(jì)算。因此在進(jìn)行信號(hào)采集時(shí)，采樣頻率要達(dá)到10 MHz。

1.2 數(shù)據(jù)處理

精確計(jì)算巖塊縱波波速的關(guān)鍵是準(zhǔn)確識(shí)別兩個(gè)傳感器接收振動(dòng)信號(hào)之間的到時(shí)差Δt，主要采用互相關(guān)分析法?；ハ嚓P(guān)函數(shù)是描述輸入信號(hào)與輸出信號(hào)在任意兩個(gè)不同時(shí)刻的相關(guān)程度。用互相關(guān)法進(jìn)行時(shí)延估計(jì)，則是通過計(jì)算輸入信號(hào)與輸出信號(hào)的互相關(guān)函數(shù)的峰值點(diǎn)的位置來確定兩函數(shù)的偏移量，從而得到到時(shí)差[13]。

設(shè)兩個(gè)信號(hào)x(t)和y(t)由被測(cè)信號(hào)s(t)和噪聲n(t)組成，即

x(t)=s(t)+n(t)

(2)

y(t)=s(t+τ)+n(t)

(3)

式中:τ為延遲時(shí)間。根據(jù)前面的相關(guān)性的定義可知，隨機(jī)信號(hào)x(t)和y(t)的互相關(guān)函數(shù)可定義為：

(4)

式中:T是隨機(jī)信號(hào)x(t)和y(t)的采樣時(shí)間?；ハ嚓P(guān)函數(shù)Rxy(τ)是τ的函數(shù)，它完整地描述了兩個(gè)信號(hào)之間的相關(guān)情況或取值依賴關(guān)系。

對(duì)振動(dòng)信號(hào)x(t)和y(t)進(jìn)行互相關(guān)分析，即對(duì)振動(dòng)信號(hào)x(t)和延遲后的信號(hào)y(t+τ)進(jìn)行卷積運(yùn)算，得到的結(jié)果就是振動(dòng)信號(hào)x(t)和y(t)的互相關(guān)函數(shù)?；ハ嚓P(guān)函數(shù)最大值對(duì)應(yīng)的延時(shí)即為到時(shí)差Δt，進(jìn)而根據(jù)式(5)求出巖塊縱波波速。

(5)

式中:Vp為巖塊縱波波速;L為巖塊長(zhǎng)度;Δt為兩個(gè)傳感器接收到振動(dòng)信號(hào)的到時(shí)差。

2 巖樣制備及試驗(yàn)方法

2.1 巖樣特征

選用大連液流電池儲(chǔ)能調(diào)峰電站取到的輝綠巖，青島金能洞庫取的花崗巖、泥質(zhì)灰?guī)r三種巖性8個(gè)巖塊作為測(cè)試巖樣，將所要測(cè)量的巖塊加工成圓柱巖樣，兩端面加工光滑。

巖塊“標(biāo)準(zhǔn)”要求進(jìn)行制備,圓柱體試件直徑取48 mm～54 mm，高徑比為2.0～2.5。巖塊兩端面不平行度誤差不大于0.05 mm，高度、直徑和邊長(zhǎng)的誤差不大于0.3 mm，端面垂直試件軸線且最大偏差不超過0.25°[1],同時(shí)也選取了幾個(gè)不同直徑及高徑比的巖塊進(jìn)行對(duì)比。對(duì)于巖心長(zhǎng)度和直徑的測(cè)量工作，用游標(biāo)卡尺直接測(cè)得,通過對(duì)多次測(cè)量結(jié)果取平均值來減小測(cè)量誤差。本次采集加工的巖塊樣品的巖性、尺寸及風(fēng)化程度如表1所示。

表1 巖塊樣品的巖性、尺寸及風(fēng)化程度

2.2 振動(dòng)法測(cè)試

實(shí)驗(yàn)前先將COINV數(shù)據(jù)采集儀、電腦、振動(dòng)加速度傳感器(頻率響應(yīng)0.5 Hz～10 kHz)、力錘等設(shè)備連接好，用桌虎鉗將巖塊固定。為了保證振動(dòng)加速度傳感器與巖塊試樣有較好的接觸，用砂紙對(duì)巖塊表面進(jìn)行打磨，使巖塊兩端面清潔、無雜質(zhì)，然后分別將兩個(gè)傳感器通過雙面膠粘到巖塊兩端橫截面中心，于測(cè)點(diǎn)1粘貼傳感器1采集輸入信號(hào)，在測(cè)點(diǎn)2粘貼傳感器2采集輸出信號(hào)。

調(diào)試設(shè)備及測(cè)試參數(shù)，測(cè)試時(shí)采樣頻率設(shè)為10 MHz，采樣時(shí)間為2 s，用力錘在巖塊橫截面處沿垂直巖塊橫截面方向敲擊激振，配合上COINV數(shù)據(jù)采集儀、電腦記錄顯示振動(dòng)信號(hào)和波形，測(cè)得巖塊縱波波速。

2.3 超聲波測(cè)試

為了驗(yàn)證本文提出的振動(dòng)測(cè)試方法的準(zhǔn)確性，采用北京智博聯(lián)的ZBL-U510非金屬聲波測(cè)試儀測(cè)量自然狀態(tài)下各取樣巖塊的縱波波速，換能器工作頻率為50 kHz，換能器與巖塊之間采用TM-100型超聲耦合劑耦合。每個(gè)巖塊重復(fù)測(cè)量3次，取平均值作為其縱波波速，記錄巖塊縱波波速。

3 試驗(yàn)結(jié)果分析

本次振動(dòng)法和超聲波法測(cè)試的結(jié)果如表2所示。圖2、圖3為振動(dòng)法實(shí)測(cè)的試樣H-1—試樣L-2八個(gè)巖塊試樣輸入與輸出端振動(dòng)加速度傳感器采集下的水平互相關(guān)分析圖以及加速度幅值時(shí)程曲線圖。

表2 巖塊縱波波速測(cè)試結(jié)果

圖2 各試樣水平互相關(guān)分析圖

圖3 各試樣水平加速度幅值時(shí)程曲線圖

由表2可以看出，對(duì)于試驗(yàn)采用的三種不同巖性的巖塊，測(cè)得的縱波波速有明顯差異，但都在一定范圍內(nèi)波動(dòng)?；◢弾r的縱波波速一般為4 614.7 m/s～5 512.3 m/s，輝綠巖孔隙致密、堅(jiān)硬，其縱波波速最高，一般為4 915.7 m/s～5 532.1 m/s，泥質(zhì)灰?guī)r的縱波波速大致為3 281.3 m/s～3 442.7 m/s。室內(nèi)所測(cè)得的輝綠巖縱波波速高于花崗巖縱波波速，泥質(zhì)灰?guī)r由于其所含方解石等礦物較多為碎屑狀，松軟易碎，波在巖塊內(nèi)傳播時(shí)間相對(duì)較長(zhǎng)，其縱波波速在三種巖性的巖塊中最低。因此，巖性是影響巖塊縱波波速最主要的因素。由于巖性不同，所加工出的巖樣其固有頻率、阻尼比、振型等固有特性不同，而巖樣的固有頻率決定了傳感器的選擇，只有當(dāng)傳感器的通頻帶高于巖樣的固有頻率時(shí)才能準(zhǔn)確地采集下振動(dòng)信號(hào)。因?yàn)閹r性是影響巖塊縱波波速最主要的因素，所以傳感器的選擇至關(guān)重要。

而對(duì)于同一巖性的巖塊，由于受到節(jié)理、裂隙、孔隙度、軟硬程度等因素的影響，其波速也有一定程度的差異[14-15]。當(dāng)激發(fā)的縱波在傳播路徑上遇到上述影響因素時(shí)，會(huì)發(fā)生一定程度的繞射、反射、能量衰減等現(xiàn)象，增大了縱波的傳播路徑，使得測(cè)試得到的波速偏低，節(jié)理、裂隙越多，會(huì)大大降低波速值。相同尺寸新鮮未風(fēng)化的花崗巖與微風(fēng)化的花崗巖縱波波速采用振動(dòng)法測(cè)試的最大相對(duì)誤差為8.3%，采用超聲波法測(cè)試的最大相對(duì)誤差為 6.7%，可見風(fēng)化程度影響了整體測(cè)試結(jié)果[16]。所以在評(píng)價(jià)巖體完整性工作中應(yīng)選取新鮮未風(fēng)化的巖塊進(jìn)行波速測(cè)試。

巖塊試樣的高徑比也是影響巖塊縱波波速的一個(gè)重要因素。由表2可以看出，同一巖性、風(fēng)化程度的巖塊試樣，隨著巖塊高徑比增大，振動(dòng)法和超聲波法測(cè)試結(jié)果有所降低，但相對(duì)于振動(dòng)法，超聲波法測(cè)試結(jié)果的變化幅度較大。大多數(shù)情況下，較大的高徑比的巖塊試樣中存在孔隙、節(jié)理、裂隙相對(duì)較多，會(huì)導(dǎo)致波速衰減，而標(biāo)準(zhǔn)試件或高徑比較小的巖塊，其所含的孔隙、節(jié)理、裂隙較少，能減小其對(duì)縱波波速測(cè)量結(jié)果的影響。

測(cè)試方法對(duì)巖塊縱波波速也有一定的影響。在受到上述因素影響時(shí)，雖然振動(dòng)法測(cè)試結(jié)果與超聲波法測(cè)試結(jié)果有相同的變化趨勢(shì)，但由于頻散效應(yīng)，在巖塊中激發(fā)的彈性波頻率不同，所測(cè)得的速度有一定差異，振動(dòng)法測(cè)得的巖塊縱波波速與超聲波法測(cè)試結(jié)果對(duì)比，兩者相對(duì)誤差達(dá)到-8.57%。

4 結(jié) 論

(1) 本文利用振動(dòng)法測(cè)試巖塊縱波波速，與超聲波法測(cè)試結(jié)果對(duì)比，兩者相對(duì)誤差達(dá)到-8.57%。采用振動(dòng)法測(cè)試巖塊縱波波速能有效減少頻散效應(yīng)對(duì)波速的影響，提高測(cè)試精度。

(2) 本文通過實(shí)測(cè)分析了各因素對(duì)巖塊縱波波速的影響規(guī)律，同時(shí)發(fā)現(xiàn)在受到巖性、節(jié)理裂隙、風(fēng)化程度等因素影響時(shí)，振動(dòng)法與超聲波法測(cè)試結(jié)果有著相同的變化趨勢(shì)。

(3) 以往無論巖體波速測(cè)試采用振動(dòng)法還是超聲波法，巖塊波速測(cè)試都是采用超聲波法，據(jù)此所得到的巖體完整性系數(shù)Kv不準(zhǔn)確。采用本文提出的振動(dòng)法測(cè)試巖塊波速能夠有效克服因彈性波頻率不同所引起的波速差異，從而精確地求得Kv值，更客觀的評(píng)價(jià)巖體完整性。