土剪破壞次聲監(jiān)測試驗研究1)

陳 喬 張 闊 陳吉龍 張少杰 呂 君? 周儒夫 韋方強(qiáng),3)

?(礦山地質(zhì)災(zāi)害成災(zāi)機(jī)理與防控重點實驗室,西安 710054)

?(中國科學(xué)院重慶綠色智能技術(shù)研究院,重慶 400714)

??(四川省川建勘察設(shè)計院,成都 610041)

??(中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,成都 610041)

???(中國科學(xué)院聲學(xué)研究所,北京 100190)

???(荊江水文水資源勘測局,湖北荊州 434099)

引言

滑坡是一種全球性的地質(zhì)災(zāi)害,它對世界上很多國家都造成了巨大的危害,促使這些國家的研究人員不斷探索滑坡災(zāi)害監(jiān)測和預(yù)警的技術(shù)手段.追溯滑坡監(jiān)測的研究歷史,學(xué)界比較公認(rèn)的時間起點為日本學(xué)者齋藤迪孝[1]在1963 年提出滑坡預(yù)報經(jīng)驗公式時,之后的幾十年各種高新技術(shù)在滑坡監(jiān)測領(lǐng)域中得到應(yīng)用.在滑坡監(jiān)測技術(shù)的發(fā)展過程中,有不少專家學(xué)者已經(jīng)從不同角度以不同觀點進(jìn)行了系統(tǒng)的總結(jié)[2-3],并對其發(fā)展趨勢做出了展望.而就目前我國的監(jiān)測技術(shù),預(yù)警預(yù)報的成功率仍然較低,如2016 年共發(fā)生地質(zhì)災(zāi)害9710 起,而全國共成功預(yù)報地質(zhì)災(zāi)害只有676 起[4],因此開展滑坡監(jiān)測新技術(shù)的探索具有較高的參考價值.

次聲波是一種在自然界非常常見的聲波,其頻率在0～20 Hz 之間,低頻率和長波長的特點使得次聲波具有傳播距離遠(yuǎn)、衰減小、抗干擾能力強(qiáng)及穿透能力強(qiáng)等優(yōu)點,在傳播過程中大氣等介質(zhì)對它的吸收系數(shù)很小,適用于遠(yuǎn)距離觀測,被廣泛的用于火山噴發(fā)[5-8]、地震[9-11]、泥石流 [12-15] 等自然災(zāi)害監(jiān)測中.前期研究表明巖土體在變形破裂過程中因能量釋放必然會產(chǎn)生彈性波和在空氣中傳播的聲波,滑坡過程會產(chǎn)生明顯的次聲波信號[16],因此次聲監(jiān)測可以作為判斷滑坡的一種技術(shù)手段.

近幾年,國內(nèi)外學(xué)者圍繞基于次聲波的滑坡監(jiān)測技術(shù)展開了大量研究.在次聲滑坡監(jiān)測預(yù)警系統(tǒng)研究方面.Li 等[17]基于譜熵和支持向量機(jī)對次聲信號進(jìn)行分類,可用于滑坡的次聲實時監(jiān)測.任際周等[18-19]從軟硬件設(shè)計的角度研究了基于次聲波技術(shù)的滑坡監(jiān)測預(yù)警系統(tǒng),通過對近年來我國數(shù)十起重大滑坡災(zāi)害臨滑預(yù)警中收集到的數(shù)據(jù)進(jìn)行分析,建立了滑坡預(yù)警模型.由于該模型是利用信號分析和統(tǒng)計的方法建立起來的,缺乏對滑坡次聲波產(chǎn)生機(jī)理的考慮,因此預(yù)警模型的準(zhǔn)確度依賴于大量現(xiàn)場數(shù)據(jù)的統(tǒng)計,而各個區(qū)域滑坡因素的多樣性決定了該模型的現(xiàn)場應(yīng)用具有一定的局限性.

在滑坡破壞過程中次聲波監(jiān)測機(jī)理研究方面.Chai 等[20]通過單軸壓縮試驗發(fā)現(xiàn)花崗巖破裂過程中也產(chǎn)生次聲信號.朱星等[21-22]以巖體穩(wěn)定性監(jiān)測為目的,對6 種典型的巖石試樣單軸全過程加載試驗過程中的次聲波信號進(jìn)行了實時監(jiān)測,并對數(shù)字次聲波信號進(jìn)行了處理與分析.徐洪等[23]對巖石變形破壞次聲異常的能量特征進(jìn)行分析,研究結(jié)果為巖石破壞前兆預(yù)警提供了重要依據(jù).賈炳等[24-30]以次聲波技術(shù)應(yīng)用于監(jiān)測、預(yù)警礦井煤巖動力災(zāi)害為目的,相繼開展了大量的研究,通過進(jìn)行煤巖試樣的單軸加載試驗,對采集的次聲信號的相對能量和頻率進(jìn)行分析,發(fā)現(xiàn)相對能量和頻率具有明顯的階段性.除了進(jìn)行實驗室物理模型試驗次聲波的研究外,朱星等[31]還對大光包滑坡側(cè)緣拉裂壁巖崩進(jìn)行次聲監(jiān)測,并分析了巖崩次聲的基本特征.上述研究結(jié)果為巖質(zhì)滑坡次聲波監(jiān)測提供了理論參考,但對土質(zhì)滑坡次聲波監(jiān)測的研究鮮有報道.基于此,考慮到滑坡的主要力學(xué)破壞方式為剪切破壞,開展土剪破壞次聲監(jiān)測試驗研究,分析了剪切過程中的微觀聲學(xué)響應(yīng)[32-34],探索土質(zhì)滑坡中次聲產(chǎn)生機(jī)理,為實現(xiàn)土質(zhì)滑坡次聲波監(jiān)測提供理論參考.

1 土體直剪實驗次聲?力?位移聯(lián)合監(jiān)測系統(tǒng)

在實驗中,土體直剪實驗力?聲聯(lián)合監(jiān)測系統(tǒng)主要由變尺寸直剪盒、力學(xué)數(shù)據(jù)采集模塊、聲學(xué)數(shù)據(jù)采集模塊以及位移采集模塊四部分組成,如圖1 所示.

1.1 變尺寸直剪盒設(shè)計

傳統(tǒng)的應(yīng)變控制式直剪儀在進(jìn)行土體剪切實驗時,采用的土樣規(guī)格一般為面積30 cm2高2 cm 的圓柱,可以實時獲取剪切應(yīng)力與剪切位移數(shù)據(jù),但不能采集聲波信號,而在此尺寸下的土樣在剪切過程中所產(chǎn)生的聲波能量又是非常微小,因此本研究需要設(shè)計一個大尺寸且尺寸可變的剪切儀器,在本論文中稱之為變尺寸直剪盒.

如圖2 所示,整個變尺寸剪切盒包括上部剪切盒、上部擋土板、固定上部擋土板的螺紋桿、下部擋土板、導(dǎo)向輪(滾珠花鍵)、下部剪切盒、固定上下剪切盒的螺紋桿、固定上部擋土板的螺紋桿,其中的上部剪切盒、上部擋土板、下部剪切盒、下部擋土板是4 個主要結(jié)構(gòu)體.在設(shè)計變尺寸直剪盒的各個部件和整體的協(xié)調(diào)組合時,根據(jù)研究條件重點考慮了基本的剪切運動、土樣尺寸可變、防銹蝕和防變形的性能.圖3 為加工廠制作后的實際成品.傳統(tǒng)直剪儀的固有結(jié)構(gòu)不利于安裝導(dǎo)波管放大次聲信號,故本次研究所設(shè)計的變尺寸直剪儀更適合土體直剪實驗過程中的次聲?力?位移聯(lián)合監(jiān)測.

圖2 變尺寸剪切盒整體設(shè)計圖Fig.2 Overall design of the variable size shear box

圖3 變尺寸剪切儀成品圖Fig.3 Finished product of variable size shear

1.2 力學(xué)數(shù)據(jù)采集模塊

在實驗中需要對土樣剪切過程中的推力進(jìn)行實時采集,因此需要一套合適的力學(xué)數(shù)據(jù)采集系統(tǒng).本研究中使用的力學(xué)采集儀器主要為平膜盒式力學(xué)傳感器,除此之外,還有配套的數(shù)字轉(zhuǎn)換器以及數(shù)據(jù)處理軟件等,圖4 為整個系統(tǒng)的示意圖.

圖4 力學(xué)數(shù)據(jù)采集系統(tǒng)示意圖Fig.4 Schematic diagram of mechanical data acquisition system

1.3 次聲數(shù)據(jù)采集模塊

次聲數(shù)據(jù)采集模塊包括剛性導(dǎo)波管、次聲傳感器、次聲數(shù)據(jù)存儲器以及相應(yīng)的數(shù)據(jù)處理軟件,其中最核心的部件當(dāng)屬次聲傳感器.本次實驗所使用的次聲傳感器來自于中國科學(xué)院聲學(xué)研究所研制的IDS2016 型次聲傳感器,該傳感器具有體積小、靈敏度高、對機(jī)械振動不敏感、頻響好的優(yōu)點.其中3 dB 平穩(wěn)帶寬保持在0.5～200 Hz 范圍內(nèi),靈敏度為50 mV/Pa,在1～100 Hz 頻率范圍內(nèi)靈敏度曲線平直,可以完全覆蓋剪切過程中產(chǎn)生次聲波的整個頻段.

1.4 位移采集模塊

針對傳統(tǒng)的基于傳感器的接觸式測量方法來測量剪切位移,存在儀器安裝困難的現(xiàn)狀,故結(jié)合本研究的實際情況,實驗采用自行研發(fā)的非接觸式光學(xué)測量方法結(jié)合數(shù)字圖像識別技術(shù)來進(jìn)行位移數(shù)據(jù)的采集.

測量的原理是通過計算物體位移變化前后呈現(xiàn)的兩幅圖像之間的相關(guān)性來計算位移值.根據(jù)本實驗的實際情況,利用此測量方法時需要在上部剪切盒表面貼一張長方形白紙(白色與鋼板的深色在二值化后會形成鮮明反差)作為定位標(biāo)記,然后使用數(shù)碼攝像機(jī)錄制整個剪切盒移動過程,得到一串?dāng)?shù)字圖像序列,再利用計算機(jī)程序?qū)γ恳桓眻D像進(jìn)行預(yù)處理,包括圖形矯正、二值化處理、選取特征像素行(特征像素行指,圖像序列中的所有圖像中,作為定位標(biāo)記的像素數(shù)量都相等的那一行像素,之所以直接取行像素作為計算量,是因為在本次研究的實驗中,上部剪切盒只做水平直線運動)等步驟(流程如圖5所示).

圖像預(yù)處理的每一步需使用相同的變換參數(shù).顯然,預(yù)處理后的圖像的特征像素行是一個一維二進(jìn)制行向量(稱為特征行向量).如果將第一張圖像的特征像素行右移(這里的移位為循環(huán)移位,即移位后末位變首位),每次右移后得到新的特征行向量序列,將這一特征行向量序列依次與其他圖像的特征行向量進(jìn)行相關(guān)性計算(見圖6),每次計算后都會得到一組相關(guān)性系數(shù),在一組相關(guān)性系數(shù)中的每一個值都對應(yīng)一個移位的像素數(shù),而移位的像素數(shù)又對應(yīng)一個位移距離,因此,每一組相關(guān)性系數(shù)序列中的最大值對應(yīng)的位移距離就是這幅圖像中被測物體在當(dāng)前時刻的位移值.

由于位移數(shù)據(jù)的采集最初來自數(shù)碼攝像機(jī),所以本測量方法的精度主要依賴于攝像機(jī)的分辨率.本次研究采用的4K 攝像機(jī)橫向分辨率為3840 個像素,在白紙所在的焦距處整個屏幕所囊括的總距離在所有實驗中平均約為1000 mm,由此可以算得測量精度為每像素0.26 mm,在之后的研究中可以采用更高分辨率的攝像機(jī)來提高測量精度.整個流程可以通過Photoshop 圖像處理軟件和Matlab 編程實現(xiàn).

圖5 圖像序列第1 幀和第n 幀的處理流程圖Fig.5 Processing flowchart of the first and nth frames of the image sequence

圖6 相關(guān)性系數(shù)計算示意圖Fig.6 Schematic diagram of correlation coefficient calculation

2 實驗材料

本次實驗選用的土樣均來自三峽庫區(qū)[35-36]重慶市奉節(jié)縣,共采集了5 種土類,包括棕壤(黏土)、黃壤(黏土)、黃棕壤(黏土)、粗骨土(含礫石黏土)、紫色土(含礫石黏土),部分土樣如圖7 所示.其中棕壤取自云霧鄉(xiāng)羅圈村,黃壤取自安坪鄉(xiāng)新浦滑坡后緣,黃棕壤取自長安鄉(xiāng)保安村,粗骨土取自草堂鎮(zhèn)廟灣,紫色土取自康樂鎮(zhèn)平皋村,取樣點位置見圖8,重塑土樣是原狀的黃壤土經(jīng)烘干、碾碎而成,各類土的相關(guān)物理力學(xué)性質(zhì)見表1.

圖7 土樣類別圖Fig.7 Types of soil samples

3 實驗方案

3.1 實驗流程

首先,將野外取回的所有土樣分批放入烘箱中,以105?C 的溫度持續(xù)24 h 鼓風(fēng)干燥,并將烘干后的土樣研磨并過2 mm 篩備用; 接下來,按具體實驗?zāi)康?將變尺寸剪切儀的上部擋土板和下部擋土板固定在規(guī)定的位置,采用分層壓實的方法將處理好的土樣堆填于剪切盒中;然后,在土樣頂部插入剛性導(dǎo)波管,在導(dǎo)波管頂部安裝次聲傳感器探頭,將力學(xué)傳感器放置在千斤頂與上部剪切盒之間,將記錄位移用的攝像頭擺放在合適位置;最后,記錄實驗開始時間,并打開數(shù)據(jù)采集軟件,在記錄30 s 背景噪音的基礎(chǔ)上,開始按壓千斤頂為上部剪切盒加載直至土樣剪切破壞.

圖8 土樣取樣點位置圖Fig.8 Location of soil sampling points

表1 各類土的相關(guān)物理力學(xué)性質(zhì)Table 1 Related physical and mechanical properties of various soils

3.2 分析方法研究

3.2.1 背景噪聲分析與處理

在室內(nèi)試驗過程中,會出現(xiàn)很多外在條件干擾的情況.例如開、關(guān)門和抽風(fēng)機(jī)運轉(zhuǎn)時會導(dǎo)致空氣不穩(wěn)定流通,實驗用水所發(fā)出的超低頻聲等.針對上述可控的低頻聲源,本試驗時間均選定在凌晨,關(guān)閉與試驗無直接關(guān)系的電器設(shè)備,以減少外界帶來的干擾.在預(yù)處理的基礎(chǔ)上,考慮到本次研究中關(guān)注的重點是聲波信號的次聲部分,所以需要對聲波信號進(jìn)行了低通濾波處理.在文中利用巴特沃斯濾波器對信號進(jìn)行處理,設(shè)置通帶截止頻率為20 Hz,阻帶截止頻率為25 Hz,通帶內(nèi)最大衰減為0.5 dB,阻帶衰減為40 dB.同時,利用背景噪聲強(qiáng)度較低的現(xiàn)狀,設(shè)置聲壓強(qiáng)度閾值為0.000 5 Pa,最后,對閾值處理后小波信號進(jìn)行小波逆變換.

3.2.2 次聲事件拾取

為了提取次聲事件,將整個時域按次聲事件分段,本文提出了一種按次聲事件由計算機(jī)程序自動劃分時域的方法.首先根據(jù)每個采樣點的信噪比初步判斷該采樣點是否位于一個次聲事件的時段內(nèi),再歸并比較密集的有效采樣點形成一個次聲事件,具體步驟如下.

第一步,初步判斷一個采樣點是否位于一個次聲事件的時段內(nèi).

根據(jù)離散聲波信號的能量公式

其中,E表示信號能量;s表示聲波通過的面積,可以認(rèn)為是常量; ρ 表示聲波傳播介質(zhì)的密度,這里即指空氣的密度;v表示音速;I表示采樣點總數(shù);x(i) 表示第i個采樣點的聲壓值,可以得到單個采樣點所代表的能量

根據(jù)功率的定義(單位時間內(nèi)的能量),本文定義離散聲波信號中第i個采樣點的功率為

本文對“均方功率”的定義由式(1)右側(cè)去掉s,ρ 和v再開方得到,因為式中的s,ρ 和v為常量,所以去掉這3 個量不影響信號點間的相對取值且避免了對它們的測量,再開方的目的是為了使最后得到的值與聲壓值量綱相同.

由式(3)定義一個判據(jù)式

其中,Pj表示第j個采樣點的功率,Pmax表示前30 s背景噪聲的最大功率,L表示用于計算信噪比的時窗長度.

式(4)左側(cè)的意義為以第i個采樣點為中心的L時長范圍內(nèi)的信號能量,右側(cè)的意義為在長度同樣為L的時段內(nèi)噪聲的最大能量,整個式子表示將局部的信號能量與相同長度噪聲能量進(jìn)行比較.針對新信號中的每一個點代入式(5) 中判斷,結(jié)果為真時,把該點標(biāo)記為1,否則標(biāo)記為0,這樣就形成一個由0 和1 組成的序列,稱為第一輪標(biāo)記序列,其長度等于原始信號.于是整個時域按信噪比被初步劃分出來.此時的劃分是不夠的,因為第一輪標(biāo)記序列中的1 值點過于零碎,需要將間隔很短的信號段再一次歸并.

第二步,對第一輪標(biāo)記序列進(jìn)行歸并整合.這一步會用到一個所有值均為1 的行向量r作為參考向量,把這個參考向量作為時窗在第一輪標(biāo)記序列的時域中平移,并與和時窗長度相同的局部按下式做計算

其中,ci表示判別該采樣點(第i個點)是否需要重新標(biāo)記為1 或0 的值,J表示移動時窗的采樣點序列總數(shù),sj表示第一輪標(biāo)記序列中的第j個點,??表示向上取整,??表示向下取整.

由式(6)計算得到結(jié)果ci后,再判別是否需要重置采樣點的聲壓值,當(dāng)ci> 0 時,說明第一輪標(biāo)記序列中與時窗長度相同的局部向量與作為參考樣本的向量相互匹配,該點可以標(biāo)記為1,否則當(dāng)ci< 0 時應(yīng)該標(biāo)記為0.最后再次得到一個基于時域的由若干0 和1 分別連續(xù)分布的新序列,稱為次聲事件標(biāo)記序列,其中連續(xù)的1 指明該時段內(nèi)存在能量明顯大于背景噪音的次聲信號,而連續(xù)的0 指明該時段內(nèi)不存在有效的次聲信號.

式(5)中的時窗長度L和式(6)中時窗長度J是兩個對最終劃分結(jié)果影響較大的計算因子.其中L應(yīng)設(shè)置一個較小值,若設(shè)置值過大,一些能量較小的次聲事件將不會被識別出.而時窗長度J影響最終次聲事件劃分出來的數(shù)量,若J取值過小,次聲事件會劃分的過密,無疑會增加后期分析的工作量,而J若取值過大,對次聲事件的劃分又太過于簡單,一些間隔很遠(yuǎn)的有效信號也被歸并為一個次聲事件.經(jīng)過反復(fù)實驗,本文認(rèn)為L=5,J=100 是合適的.最后根據(jù)式(5)和式(6),使用MATLAB 編制了對信號按次聲事件自動劃分時域的程序.

4 試驗結(jié)果分析

4.1 應(yīng)變軟化類黏性土的漸進(jìn)性破壞過程

本次共完成了不同尺寸、不同土性和不同飽和度共13 組土剪切實驗(見表2).圖9 展示了所有實驗的“推力?時間”曲線和聲波信號時域圖土樣的剪切力與位移原始關(guān)系.從推力變化的整體趨勢來看,每組推力都隨時間先陡然升高,超過峰值后逐漸降低.通過對推力進(jìn)行規(guī)格化處理,統(tǒng)一在0～7000 N 范圍內(nèi),可以進(jìn)一步反映推力在時域上的變化趨勢.以具有代表性的sy3 實驗為例(如圖10),圖中實線為實測推力值,推力曲線呈現(xiàn)鋸齒狀是因為本次實驗加載方式采用的是千斤頂人工加載,所以導(dǎo)致加載不連續(xù),反映出了該類土的應(yīng)力松弛特點.

人工加載對剪切速度的影響可以認(rèn)為是微乎其微的,在本實驗中推力與剪切帶受到的剪力幾乎相等,由牛頓第二定律

其中,f是土體受到的剪力,單位kN;F是千斤頂?shù)耐屏?單位kN;m是土體上半部分質(zhì)量,實測sy3 的土體總質(zhì)量49.24 kg;a為位移加速度,單位為m/s2.經(jīng)計算ma的值在?0.046 6 到0.049 9 之間,遠(yuǎn)遠(yuǎn)小于推力大小,因此可大致認(rèn)為在千斤頂加載的過程中上部剪切盒保持勻速運動.

表2 實驗組別一覽Table 2 List of experimental groups

圖9 “推力?時間”曲線和聲波信號時域圖Fig.9 “Thrust-time”curve and time domain diagram of acoustic signal

圖10 sy3 推力?位移曲線Fig.10 Thrust-displacement curve of sy3

同時,推力的趨勢線隨位移變化呈現(xiàn)升高后降低的規(guī)律,當(dāng)sy3 發(fā)生剪切破壞時,首先在某一點處超過抗剪強(qiáng)度而發(fā)生破裂,然后從此破裂點開始向整個剪切面蔓延擴(kuò)展,隨著剪切面的蔓延擴(kuò)展,整個土體的抗剪強(qiáng)度也在降低,應(yīng)變局部化現(xiàn)象增強(qiáng),直至整個剪切面貫通,整體表現(xiàn)出明顯的應(yīng)變軟化特性.

4.2 剪切過程的次聲響應(yīng)

根據(jù)漸進(jìn)破壞理論和推力趨勢線的形態(tài),本文將其劃分為了彈性變形階段、彈塑性變形階段以及峰值后的應(yīng)變軟化階段3 個階段.為了探究與聲學(xué)信號的關(guān)系,將推力?位移曲線在時域上重新描述,繪制出的“推力?時間”曲線圖,以sy3 為例(見圖11),推力趨勢線在經(jīng)過前30 s 的背景噪音錄制階段后隨時間逐漸增大,加載后先是呈直線狀態(tài)增大,處于彈性變形階段;在74 s 后曲線斜率開始減小,曲線本身呈向下彎曲狀,處于彈塑性變形階段; 在118 s 處推力趨勢線達(dá)到峰值6.1 kN 后開始減少,隨后隨著時間推移逐漸減小,對應(yīng)著應(yīng)變軟化階段.

sy3 在次聲事件拾取的基礎(chǔ)上獲取了26 個次聲事件(見圖11),將每個事件對應(yīng)的次聲信號均方功率峰值用虛線連接起來形成均方功率的包絡(luò)線,隨著時間推移包絡(luò)線開始逐漸上升,在105 s 處均方功率達(dá)到0.047 9 Pa 后開始下降,之后偶有起伏但都沒有超過105 s 處的峰值,并且很明顯的,該包絡(luò)線的時域形態(tài)與推力趨勢線較為相似.

圖11 sy3 均方功率圖Fig.11 sy3 mean square power diagram

從圖11 來看,剪切試驗過程中的時域信號在105 s 處能量達(dá)到了峰值,同時應(yīng)力曲線也達(dá)到峰值,實際上從現(xiàn)場結(jié)果來看試樣在達(dá)到應(yīng)力峰值時已經(jīng)破壞,但此時應(yīng)力并沒有大幅度下降,表明土樣還有一定殘余強(qiáng)度.破壞過程中的次聲信號強(qiáng)度比破裂后剪切面滑移所產(chǎn)生的次聲信號強(qiáng)度大,說明土樣內(nèi)部裂隙發(fā)展所產(chǎn)生的次聲信號與土樣顆粒摩擦所產(chǎn)生的次聲信號有一定的區(qū)別.

在此基礎(chǔ)上,將所有組實驗的均方峰值與推力峰值之間的時間差進(jìn)行計算,從結(jié)果來看(見表3 和圖12),次聲信號的均方功率峰值都在推力峰值之前,提前的時間分布在2.30～47.11 s,平均值為22.95 s.該結(jié)論可為土質(zhì)滑坡次聲預(yù)警模型研究提供理論參考,具有重大意義和工程應(yīng)用價值.

表3 各組實驗次聲事件次數(shù)與時間差計算結(jié)果Table 1 Calculation results of the number and time difference of infrasound events in each group of experiments

圖12 所有實驗均方功率圖Fig.12 Mean square power plots for all experiments

在“推力?位移” 曲線的3 個階段中,彈性變形階段推力趨勢線呈近直線型,斜率幾乎不變,彈塑性變形階段趨勢線的增量逐漸降低,線形向下凹,應(yīng)變軟化階段趨勢線經(jīng)過峰值后逐漸降低,但較彈性變形階段下降趨勢更為平緩.根據(jù)應(yīng)變軟化類黏性土在剪切破壞中滿足漸進(jìn)性破壞理論,剪切面破壞面的形成呈漸進(jìn)性,另外參考一些對大型直剪實驗剪切帶形態(tài)的相關(guān)研究[37],本文認(rèn)為彈性變形階段對應(yīng)于圖13 中的剪切前到初始剪切之間的狀態(tài),彈塑性變形階段和應(yīng)變軟化階段對應(yīng)初始剪切到剪切完成之間的狀態(tài).

圖13 剪切過程示意圖(陰影區(qū)域表示壓密區(qū),1#、2#和3#三根粗直線表示鋁絲,用于表征土樣平行于剪切方向的截面的形變)Fig.13 Schematic diagram of the cutting process

在剪切過程初期,土體整體變形以上部試樣中壓密區(qū)土體的彈性變形為主,剪切帶中的剪切變形為輔,因此推力趨勢線呈近直線狀,在此階段,土顆粒主要發(fā)生彈性變形,顆粒間接觸面相互擠壓產(chǎn)生的聲波能量較小,但隨著剪切過程的發(fā)展,土體整體變形進(jìn)入到以剪切帶塑性變形為主的階段,此階段中主要以黏土膠體團(tuán)粒本身的變形為主,這種膠體團(tuán)粒的變形比彈性變形所產(chǎn)生的聲波能量要大,從以土顆粒彈性變形為主到黏土團(tuán)粒本身變形為主,后者產(chǎn)生的聲波信號成分也逐漸超過前者,所以在宏觀上表現(xiàn)為彈塑性階段次聲信號幅度值比彈性階段更大(見圖14),均方功率峰值一般也發(fā)生在此階段.在推力峰值時,土體剪切面開始漸進(jìn)性破壞,并進(jìn)入到應(yīng)變軟化階段,在此階段剪切面從某一點發(fā)生破裂并逐步蔓延到整個剪切面,此時剪切帶的運動形式在微觀上發(fā)展成以黏土團(tuán)粒拉張破裂后土顆粒間的摩擦為主,黏土團(tuán)粒的破裂為輔,此時的信號能量也逐漸減弱,偶爾局部團(tuán)粒破裂活動較為集中時,在時域信號上還會產(chǎn)生相應(yīng)的脈沖,但包含此脈沖的次聲事件的幅值一般不會再超過彈性階段的幅值.

圖14 不同階段聲波幅值變化情況Fig.14 Changes in sound wave amplitude at different stages

5 結(jié)論

(1)針對常規(guī)土體應(yīng)變式直剪儀用于聲發(fā)射實驗時聲學(xué)信號微弱的缺陷,研發(fā)了專門用于探測次聲信號的變尺寸直剪儀,設(shè)計了非接觸式位移數(shù)據(jù)的自動采集模塊,搭建了一套大尺寸土體直剪實驗力?聲?位移聯(lián)合采集系統(tǒng),為土質(zhì)滑坡次聲監(jiān)測技術(shù)室內(nèi)實驗研究奠定了基礎(chǔ).

(2) 將聲學(xué)信號結(jié)合力學(xué)和位移數(shù)據(jù)進(jìn)行分析,發(fā)現(xiàn)次聲信號在時域上的脈沖形態(tài)與直剪實驗中的“推力?時間” 曲線具有明顯的相關(guān)性,并在應(yīng)變軟化類黏性土直剪過程中剪切帶的漸進(jìn)性破壞理論的基礎(chǔ)上,分析了剪切過程中彈性階段、彈塑性階段和應(yīng)變軟化階段的微觀聲學(xué)響應(yīng),認(rèn)為黏性土剪切破壞時發(fā)出的次聲信號主要來源于彈性階段中黏土顆粒的相互擠壓、彈塑性階段和應(yīng)變軟化階段中由黏土顆粒構(gòu)成的膠體團(tuán)粒本身的變形以及黏土膠體團(tuán)粒拉張破裂后包含其中的土顆粒間的相互摩擦三種形式.

(3)在土剪不同階段,次聲信號的幅值大小呈現(xiàn)出應(yīng)變軟化階段< 彈性階段< 彈塑性階段的規(guī)律,同時,次聲信號的均方功率峰值都在推力峰值之前,平均提前為22.95 s,這說明將次聲技術(shù)作為土質(zhì)滑坡預(yù)警手段的具有理論可行性.

本研究對于土質(zhì)邊坡的動態(tài)監(jiān)測和穩(wěn)定性預(yù)測都具有重要的指導(dǎo)意義,可為進(jìn)一步基于次聲的土質(zhì)滑坡預(yù)警模型研究提供一定的理論參考.