高結(jié)旺，黃宜超，鄒杰平

(江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州 314000)

隨著社會經(jīng)濟(jì)發(fā)展和人民物質(zhì)需求的日益增長，國家對于地下資源開采和建設(shè)逐漸加大。而巖土體在自然條件下，風(fēng)化作用導(dǎo)致其內(nèi)部常存在大小、形狀不一的宏微觀缺陷，這種自身存在的缺陷給巖土體安全穩(wěn)定帶來潛在的危害，而在人為擾動、地震等外界因素影響下進(jìn)一步加劇了巖土體的失穩(wěn)破壞，時常給地下施工帶來巨大的生命財產(chǎn)安全威脅[1-3]。因此，對巖土體缺陷結(jié)構(gòu)特征超前感知與識別，可以為地下工程施工提供寶貴的工程經(jīng)驗[4-6]。

近些年，超聲波探傷等技術(shù)被廣泛的運(yùn)用在地下巖土工程中。通過分析超聲波在巖土體傳播過程中的聲學(xué)參數(shù)與巖土體結(jié)構(gòu)特征等之間的關(guān)系，可以識別地下巖體結(jié)構(gòu)、構(gòu)造等基本情況，為減少因地下施工所帶來的安全事故提供前期準(zhǔn)備[7]。M.Kumar等[8]研究發(fā)現(xiàn)當(dāng)介質(zhì)孔洞尺寸遠(yuǎn)小于聲波波長時，孔洞形狀是影響聲波波速的關(guān)鍵。在此基礎(chǔ)上C.Scotellaro等[9]研究指出相比于缺陷其他特征對聲波波速的影響而言，缺陷形狀特征表現(xiàn)的最明顯，巖體剛度與缺陷縱橫比呈正相關(guān)，而波速應(yīng)力敏感性則與缺陷縱橫比呈負(fù)相關(guān)。G.T.Baechle等[10]研究發(fā)現(xiàn)，波速與小孔洞相關(guān)性相比于大孔洞更明顯。Weger等[11]更進(jìn)一步對影響聲學(xué)參數(shù)的缺陷特征進(jìn)行細(xì)化的研究，發(fā)現(xiàn)相比缺陷其他特征而言，聲學(xué)參數(shù)對孔隙結(jié)構(gòu)復(fù)雜性和大孔洞孔徑綜合效應(yīng)更敏感。劉向君等[12]研究了孔洞密度對聲學(xué)參數(shù)之間的表征關(guān)系，指出聲波波速和波幅與孔洞密度呈負(fù)相關(guān)性，并且波幅的變化相對于波速的變化更明顯。梁利喜等[13]基于波動理論和有限差分理論，分別研究了孔洞密度、分布、形狀對聲波衰減系數(shù)影響，發(fā)現(xiàn)孔洞密度對聲波衰減系數(shù)產(chǎn)生最大影響。

大量研究者從室內(nèi)試驗和數(shù)值模型研究了超聲波在介質(zhì)中傳播特性，并指出超聲波衰減系數(shù)、振幅、波速、相位對孔洞結(jié)構(gòu)的變化更敏感[14]。但目前對于超聲波在介質(zhì)中傳播規(guī)律的數(shù)值模擬研究主要集中在二維數(shù)值模型，很難真實的反映超聲波在地下巖土體三維空間中的傳播規(guī)律。為深入研究含孔洞巖樣中超聲波傳播特性，利用COMSOL軟件分別建立了不同孔洞尺寸、角度、形狀三維孔洞壓力聲學(xué)模型，分別從超聲波透射信號、波幅透射系數(shù)、聲壓快照圖、聲壓極大值四個方面分析超聲波傳播規(guī)律。為下一步原位聲波識別復(fù)雜的洞縫組合結(jié)構(gòu)提供理論依據(jù)。

1 數(shù)值模型的建立

采用φ50 mm×100 mm圓柱體巖樣模型，巖樣密度2 817 kg/m3，超聲波波速為2 500 m/s。COMSOL數(shù)值建模流程如下：1)根據(jù)需要求解問題類型選定恰當(dāng)?shù)奈锢韴鲅芯浚?)建立三維孔洞模型并定義求解區(qū)域的材料屬性，施加聲波邊界條件后對三維孔洞模型進(jìn)行網(wǎng)格剖分；3)邊界條件設(shè)定好后對模型進(jìn)行瞬態(tài)求解計算。通過分析模型計算結(jié)果中超聲波透射信號、波幅透射系數(shù)、聲壓快照圖、聲壓極大值的規(guī)律，反演得到超聲波在孔洞模型中傳播特征規(guī)律。

1.1 模型簡化假設(shè)

為了便于研究超聲波在介質(zhì)中傳播規(guī)律問題，傳統(tǒng)的方法常假設(shè)介質(zhì)為各向同性均質(zhì)材料，及介質(zhì)的力學(xué)和材料性質(zhì)與空間位置無關(guān)[15-16]。

1.2 控制方程

模型不考慮溫度方程的情況下，壓力聲學(xué)瞬態(tài)研究的控制方程為：

(1)

其中，ρ0為介質(zhì)密度；p為聲壓；q為偶極源；c為介質(zhì)中聲速；t為時間；Qm為單極源。

1.3 模型的建立

選用COMSOL數(shù)值軟件壓力聲學(xué)物理場模塊，在瞬態(tài)時域條件下探究超聲波在缺陷巖體中透射規(guī)律試驗。為了達(dá)到數(shù)值試驗?zāi)康模绫?所示，利用帶有規(guī)則孔洞巖樣模型代替缺陷巖體，分別探究孔洞形狀、尺寸、角度三種因素對超聲波透射的影響，進(jìn)而反演得到超聲波在缺陷巖體中傳播特征規(guī)律。

表1 孔洞缺陷模型幾何參數(shù)

1.4 邊界條件設(shè)置

由于設(shè)置合理的初始邊界條件決定著模型最終計算結(jié)果精度，因此含孔洞巖樣模型超聲波透射數(shù)值試驗初始邊界條件設(shè)置主要有如下幾點：

1)設(shè)置孔洞模型表面與內(nèi)部最初聲壓場為零。

2)設(shè)置模型側(cè)面邊界與底部邊界為完美匹配層邊界。

3)模型頂部為聲壓入射邊界，利用高斯脈沖函數(shù)信號作為激勵源調(diào)控聲壓大小。使用脈沖力去描述超聲波的產(chǎn)生是一種簡單常見的方法，如在鋼軌超聲波檢測中常用高斯脈沖函數(shù)描述脈沖信號觸發(fā)[17]。因此，選擇高斯脈沖函數(shù)信號作為超聲波激勵源。高斯脈沖函數(shù)信號的表達(dá)式為[18]：

(2)

其中，ω0=2πf為脈沖信號的中心角頻率；U0為脈沖電壓；Δt為脈沖延時；σ2為控制脈沖寬度。聲波激勵源的發(fā)射頻率選用100 kHz高斯脈沖信號，激勵信號的時域波形圖如圖1所示。

1.5 網(wǎng)格劃分

在利用有限元COMSOL軟件進(jìn)行模型求解時，模型網(wǎng)格劃分需根據(jù)求解問題類型來確定。為了得到超聲波在含孔洞模型傳播過程中解析波這特殊問題，常規(guī)的網(wǎng)格劃分會導(dǎo)致模型計算出現(xiàn)精度低甚至不收斂的現(xiàn)象。為了精確求解含孔洞模型中超聲波傳播過程中解析波，采用奈奎斯特準(zhǔn)則來嚴(yán)格限制模型網(wǎng)格的劃分[19]，如式(3)所示，采用模型網(wǎng)格最大單元大小最大不能超過聲波波長1/10來限制模型網(wǎng)格的劃分。

(3)

其中，L為模型最大單元大?。籿為超聲波波速；fmax為超聲波最大頻率。

2 結(jié)果與分析

2.1 孔洞尺寸對超聲波傳播特性的影響

圖2為巖樣孔洞形狀為圓形、孔洞尺寸大小分別為0 mm，5 mm，10 mm，15 mm，20 mm的模型透射信號曲線。分析得到超聲波在介質(zhì)傳播過程中聲壓波幅與孔洞尺寸呈現(xiàn)明顯的負(fù)相關(guān)，即模型孔洞尺寸越大，聲壓波幅越小。如當(dāng)模型為完整時，聲壓波幅為390 Pa；孔洞大小為10 mm時，聲壓波幅為337 Pa；孔洞大小為20 mm時，聲壓波幅為235 Pa。出現(xiàn)該現(xiàn)象原因在于介質(zhì)中孔洞對超聲波傳播產(chǎn)生了散射、繞射、放射等干擾現(xiàn)象，導(dǎo)致超聲波能量在介質(zhì)傳播過程中產(chǎn)生損失，致使聲壓波幅出現(xiàn)大幅下降現(xiàn)象。

當(dāng)模型中存在孔洞時，由于超聲波在傳播過程中介質(zhì)屬性發(fā)生改變，聲波能量必定發(fā)生不同量級的損失。因此引入波幅透射系數(shù)α來定量的表征超聲波衰減，如式(4)所示：

(4)

其中，α為波幅透射系數(shù)；Ai,A0分別為有孔洞和無孔洞模型接收信號幅值。波幅透射系數(shù)越接近1，則超聲波能量在缺陷巖體內(nèi)衰減度越小，即超聲波透射性能越好。

圖3為透射系數(shù)α與孔洞尺寸的擬合關(guān)系曲線，分析發(fā)現(xiàn)模型孔洞大小與波幅透射系數(shù)α存在線性關(guān)系，線性擬合相關(guān)系數(shù)R2達(dá)到0.998。當(dāng)完整模型波幅透射系數(shù)為1，說明超聲波透射過程中沒有發(fā)生衰減現(xiàn)象；而當(dāng)孔洞尺寸為5 mm時，透射系數(shù)衰減到0.987較無孔洞下降了1.3%；孔洞尺寸增大為20 mm時，透射系數(shù)進(jìn)一步衰減到0.607較無孔洞下降了39.7%。這說明超聲波在模型傳播過程中，超聲波透射能力極易受孔洞大小的影響，即模型孔洞尺寸越大，超聲波透射系數(shù)越低，超聲波能量衰減越大，相應(yīng)超聲波透射能力越低。

通過對透射邊界上超聲波透射信號和透射系數(shù)研究，驗證了孔洞尺寸大小會干擾超聲波的傳播并且產(chǎn)生嚴(yán)重的能量損失。為了深入分析超聲波在孔洞模型傳播過程中聲壓分布規(guī)律，選取10 mm孔洞模型中聲壓快照圖(見圖4)。發(fā)現(xiàn)超聲波在孔洞模型傳播過程中聲壓分布隨時間的變化主要分為兩個階段：第一階段，聲波從入射邊界傳播到孔洞之前時，聲波聲壓分布呈現(xiàn)有規(guī)律的直線傳播，在此期間沒有出現(xiàn)繞射、反射現(xiàn)象；第二階段當(dāng)聲波傳播穿透過孔洞到達(dá)接收邊界時，由于傳播介質(zhì)的變化，致使聲波在介質(zhì)中的傳播發(fā)生明顯的繞射、反射現(xiàn)象。

在聲波傳播過程中模型孔洞尺寸不僅會對聲壓分布產(chǎn)生影響，也會對聲壓大小產(chǎn)生一定的影響。為了研究孔洞尺寸對超聲波傳播過程中聲壓大小的影響，提取了五種孔洞尺寸模型下超聲波聲壓極大值與時間關(guān)系曲線(見圖5)。發(fā)現(xiàn)可以把不同孔洞尺寸模型聲壓極大值的時變規(guī)律分成三階段：傳播時間0 μs～20 μs期間為第一階段，由于在此階段五種不同孔洞尺寸模型超聲波傳播并未透過孔洞，所以各模型聲壓極大值重合成一條曲線且當(dāng)傳播時間達(dá)到20 μs時均達(dá)到1 095 Pa。原因是此階段超聲波在各孔洞模型傳播過程中，聲壓極大值由于介質(zhì)振動數(shù)量和速度增加而增加；傳播時間20 μs～30 μs期間為第二階段，在此階段超聲波完全透過孔洞，各孔洞尺寸模型出現(xiàn)第一次聲壓極大值峰值時間點在28 μs左右。出現(xiàn)該現(xiàn)象是因為當(dāng)超聲波透過孔洞到達(dá)孔洞下邊界時，由于介質(zhì)的變化導(dǎo)致孔洞下邊界上質(zhì)點振動速度增大進(jìn)而聲壓增大。進(jìn)一步分析發(fā)現(xiàn)聲壓極大值峰值時間點以及大小都與孔洞尺寸呈現(xiàn)正相關(guān)性，原因是孔洞尺寸越大導(dǎo)致超聲波到達(dá)孔洞上邊界的傳播時間越短，相應(yīng)的邊界質(zhì)點振動時間提前，出現(xiàn)峰值時間越早峰值點的值也越大；傳播時間30 μs之后為第三階段，在此期間聲壓極大值出現(xiàn)一次增加后再減少，原因是聲波在到達(dá)接收邊界上時，導(dǎo)致邊界上受振動質(zhì)子的速度和個數(shù)增加而后聲波逐漸被接收邊界吸收。

2.2 孔洞角度對超聲波傳播特性的影響

為了探究孔洞角度與超聲波傳播特性之間的關(guān)系，選取了0°，30°，45°，60°，90°和完整(無孔洞)模型透射邊界上聲壓波幅大小進(jìn)行對比分析。由圖6可知，模型超聲波透射幅值與孔洞角度的變化呈現(xiàn)負(fù)相關(guān)性。完整模型的聲壓波幅最大為384 Pa，而0°模型較完整模型波幅下降了10.7%，為343 Pa，當(dāng)模型孔洞角度為90°時，較完整模型波幅下降了44.5%，為213 Pa。原因是孔洞角度決定著聲波垂直向孔洞傳播方向面積的投影，即孔洞角度增大時投影面積增大，造成對聲波透射干擾增強(qiáng)，聲波傳播過程中的能量損失增大。

圖7為5種不同孔洞角度與波幅透射系數(shù)的關(guān)系，分析發(fā)現(xiàn)孔洞角度與透射系數(shù)存在線性關(guān)系，擬合相關(guān)系數(shù)R2為0.97，這進(jìn)一步說明隨著孔洞角度增加即垂直聲波傳播方向投影面積增大，聲波在模型傳播過程中能量衰減越大，導(dǎo)致波幅透射系數(shù)降低，超聲波在孔洞模型中透射能力越弱。

由于超聲波在介質(zhì)傳播過程中缺陷(孔洞、裂縫等)會對聲波聲場傳播產(chǎn)生一定的干擾，如圖8所示，選取了孔洞角度為45°條件下超聲波傳播過程中聲壓分布情況。分析發(fā)現(xiàn)隨著聲波在孔洞模型中傳播聲壓值和分布出現(xiàn)了明顯差異。在聲波從入射邊界進(jìn)入但并未完全透過孔洞0 μs～20 μs期間，聲波沿直線傳播沒有出現(xiàn)干擾現(xiàn)象；當(dāng)聲波傳播穿過孔洞20 μs～30 μs期間，由于傳播介質(zhì)的改變干擾了聲波的傳播，此階段聲波出現(xiàn)明顯的反射現(xiàn)象；聲波傳播到被透射邊界完全吸收的30 μs～40 μs期間，此階段聲波出現(xiàn)明顯的散射和繞射現(xiàn)象，直至聲波傳播到透射邊界后被完全吸收，至此結(jié)束了超聲波的傳播過程。

為了深入探究不同孔洞角度下超聲波傳播過程中聲壓大小時變規(guī)律，如圖9所示提取了傳播全過程中聲壓極大值與時間關(guān)系曲線進(jìn)行分析。

可知在聲壓極大值方面，不同孔洞角度對超聲波傳播的響應(yīng)主要有兩方面：

1)聲壓極大值出現(xiàn)多峰值。當(dāng)聲波從入射邊界傳播到孔洞邊界之前時，各孔洞模型聲壓極大值均隨時間等值增大；當(dāng)聲波透過邊界進(jìn)入孔洞29 μs左右時，由于聲波傳播到孔洞邊界導(dǎo)致邊界上受振動質(zhì)點的數(shù)量和速度增大，各孔洞模型相繼出現(xiàn)第一次也是全部傳播時間段內(nèi)最大聲壓極大值峰值點；當(dāng)聲波透出孔洞到最后被模型底部完美匹配層吸收時，聲壓極大值出現(xiàn)二次峰值點，但聲壓值總體趨勢呈遞減趨勢發(fā)展。

2)各孔洞模型中出現(xiàn)聲壓極大值峰值大小與時間點和孔洞角度相關(guān)聯(lián)。如圖9所示，當(dāng)各模型出現(xiàn)第一次也是最大一次聲壓極大值峰值點時，孔洞模型角度越大出現(xiàn)聲壓極大值越大，而相應(yīng)的時間點越晚，比如孔洞角度為0°，45°，90°時，出現(xiàn)峰值點的大小和時間分別為27.1 μs和1 785 Pa，28.9 μs和2 093 Pa，29.7 μs和2 234 Pa。

2.3 孔洞形狀對超聲波傳播特性的影響

為了探究孔洞形狀對超聲波傳播特性的影響，如圖10所示，選取了5種不同形狀孔洞模型中超聲波聲壓透射信號進(jìn)行分析，發(fā)現(xiàn)單從孔洞形狀方面分析聲波波幅信號沒有明顯的規(guī)律性。而從孔洞垂直聲波傳播方向的尺寸分析發(fā)現(xiàn)，垂直聲波傳播方向尺寸分別為5 mm，10 mm，20 mm的長方形A、正方形、等邊三角形B模型中聲壓波幅為364 Pa，291 Pa，233 Pa。說明當(dāng)模型孔隙率相同但孔洞形狀不同時，孔洞垂直聲波傳播方向的尺寸是影響波幅大小的關(guān)鍵，即孔洞垂直聲波傳播方向的尺寸越大，聲波波幅越小。

基于孔洞形狀與波幅信號沒有明顯的規(guī)律性，而與垂直聲波傳播方向的孔洞尺寸存在明顯的負(fù)相關(guān)，因此提取圖11波幅透射系數(shù)與孔洞垂直聲波傳播方向上投影面積關(guān)系深入分析。利用線性擬合的方法將五種不同孔洞形狀垂直聲波傳播方向的投影面積與透射系數(shù)進(jìn)行分析，線性擬合系數(shù)R2為0.777擬合效果一般，但這也在一定程度上反映了在孔隙率相同而孔洞形狀不同條件下，孔洞垂直聲波傳播方向的投影面積越大，波幅透射系數(shù)越小，聲波透射能力越弱。

為了深入探究超聲波傳播過程中聲場分布與孔洞形狀之間的關(guān)系，如圖12所示，分別選取了正方形和三角形A兩種典型聲波傳播過程中的聲壓快照圖進(jìn)行研究。由圖12發(fā)現(xiàn)，聲波傳播到孔洞尺寸、邊界、棱角都會出現(xiàn)聲壓場紊亂的現(xiàn)象，說明聲波傳播受到了孔洞干擾，但孔洞形狀對聲壓場干擾的形式不同。如正方形孔洞(如圖12(a)所示)超聲波在孔洞上邊界主要發(fā)生反射現(xiàn)象，而繞射現(xiàn)象主要發(fā)生在孔洞左右邊界。而三角形A孔洞(如圖12(b)所示)在聲場中主要發(fā)生散射和繞射現(xiàn)象。

基于不同孔洞形狀會導(dǎo)致模型中超聲波聲壓場出現(xiàn)不同的紊亂現(xiàn)象，為了探究不同孔洞形狀中聲波傳播過程中聲壓極大值與孔洞形狀的變化關(guān)系，選取了如圖13所示五種不同孔洞形狀聲壓極大值時間曲線圖。發(fā)現(xiàn)不同孔洞形狀聲壓極大值出現(xiàn)明顯的多峰值現(xiàn)象，但出現(xiàn)聲壓極大值峰值時間點和大小與孔洞形狀沒有明顯的關(guān)系，而在孔隙率、透射下邊界尺寸相同的不同孔洞形狀模型中，出現(xiàn)聲波聲壓極大值最大峰值點的大小和時間與孔洞平行聲波傳播方向的尺寸呈負(fù)相關(guān)，即孔洞平行聲波傳播方向的尺寸越大，聲壓極大值最大峰值點的值越小，相應(yīng)出現(xiàn)峰值點的時間點越晚。如圖13所示，正方形孔洞和三角形A孔洞模型孔隙率、透射下邊界尺寸相同，出現(xiàn)聲壓極大值最大峰值點的時間和大小分別是29.1 μs和2 417 Pa，33.4 μs和1 814 Pa，在相同條件下長方形B和三角形B模型中出現(xiàn)聲壓極大值最大峰值點的時間和大小分別是30.3 μs和2 421 Pa，31.9 μs和2 293 Pa。所以在孔隙率、透射下邊界尺寸相同的不同孔洞形狀模型中，平行聲波傳播方向的尺寸是造成出現(xiàn)聲壓極大值最大峰值點的時間和大小存在規(guī)律性差異的關(guān)鍵，也從側(cè)面反映了孔洞形狀會造成超聲波在傳播過程中存在差異。

3 結(jié)論

為了探究含孔洞巖樣中超聲波傳播特性規(guī)律，采用COMSOL建立了不同孔洞尺寸、角度、形狀三維孔洞模型，利用壓力聲學(xué)模塊進(jìn)行求解計算，從聲壓透射信號、波幅透射系數(shù)、聲壓快照圖、聲壓極大值四個聲學(xué)參數(shù)方面分別探究了不同孔洞模型超聲波傳播特性規(guī)律。主要結(jié)論有如下幾點：1)從聲壓透射信號、波幅透射系數(shù)、聲壓快照圖、聲壓極大值四個聲學(xué)參數(shù)與不同孔洞模型整體分析得到，三種孔洞模型對超聲波傳播影響關(guān)系為孔洞尺寸>孔洞角度>孔洞形狀。2)超聲波在三種孔洞模型傳播中，超聲波透射信號和波幅透射系數(shù)與垂直于聲波傳播方向孔洞尺寸或面積存在負(fù)相關(guān)性，說明垂直于聲波傳播方向孔洞尺寸或面積是造成聲波能量衰減的主要原因。因此，從超聲波透射信號和波幅透射系數(shù)方面可以定量的評價地下巖土體孔洞的大小。3)巖體存在缺陷時，超聲波在傳播過程中會產(chǎn)生反射、繞射、散射影響，導(dǎo)致聲壓場出現(xiàn)不同程度的紊亂現(xiàn)象。在孔洞尺寸方面，聲壓極大值峰值時間點和峰值大小與孔洞尺寸呈正相關(guān)；孔洞角度方面，聲壓極大值峰值時間點與孔洞角度呈負(fù)相關(guān)，與峰值大小呈正相關(guān)；在孔洞形狀方面，聲壓極大值最大峰值時間點和峰值大小與孔洞平行聲波傳播方向的尺寸呈正相關(guān)。因此，分析聲壓極大值最大峰值時間點和峰值大小關(guān)系可以定性的評價地下巖體孔洞形態(tài)。