脆性巖石應力
——應變全過程聲學特征演化規律

張國凱，李海波，王明洋，鄧樹新，王振

(1.南京理工大學機械工程學院，江蘇南京，210094；2.中國科學院武漢巖土力學研究所巖土力學與工程國家重點實驗室，湖北武漢，430071；3.陸軍工程大學爆炸沖擊防災減災國家重點實驗室，江蘇南京，210007)

巖石破裂過程一直是巖石力學研究的熱點。在大規模地下工程建設過程中，伴隨著巖體的開挖卸荷，圍巖內部缺陷的孕育萌生、擴展、聚合等導致巖體力學特性劣化，甚至引起巖體工程的失穩破壞，因此，研究巖石在加(卸)載下裂紋擴展規律及對應的宏觀及細觀表征具有重要意義。聲發射(AE)是材料斷裂導致能量以應力波形式釋放的現象，可用于研究脆性材料破裂及損傷過程[1-2]。目前，AE信號后處理主要包括AE特征參數分析法和波形分析法[3]。最常用的是AE特征參數分析，如AE能量[1-2]、參數b(與震源活動特征相關的參數)[4]、振鈴計數和峰值頻率[5]、AE事件[6]等。紀洪廣[5]通過高頻和低頻傳感器同步測試花崗巖AE 演化，發現低頻與高頻AE 變化趨勢基本一致，但主頻帶在破裂過程中發生改變。EBERHARDT等[7]研究了Lac du Bonnet花崗巖破壞過程中AE信號平均持時、振鈴下降計數、能量等特征參數變化規律，發現AE活動頻繁時對應新生裂紋的起裂。GANNE等[8]通過AE測試，發現巖石破裂過程中AE能量累積的4個階段分別對應裂紋的產生、擴展、累積以及最終破壞的發生，且在材料破壞前出現明顯的能量加速釋放。LOCKNER[9]通過AE 測試發現，巖石破壞前AE事件率與非彈性應變密切相關，AE 震源定位很好地表征了微破裂的萌生、匯聚發展過程。EBERHARDT等[10]采用AE 事件量化了巖石的損傷演化規律。KIM等[11]對比了宏觀變形和AE 表征的花崗巖損傷，發現AE 能量對裂紋起裂擴展更加敏感。然而，相比AE特征參數，波形分析能夠更直觀、全面地反映巖石破裂特征[12]。何滿潮等[3]對巖爆過程中關鍵拐點處AE時域波形及功率譜圖進行了分析，發現AE頻譜經歷了由低頻單峰向高頻多峰轉變，最終又回歸低頻單峰的演化過程。同時，AE 頻率特性為研究巖石微觀斷裂及分析微裂紋的破裂機制提供了重要依據[13]。聲波測試是一種受環境影響小、操作快速便捷的無損檢測手段，由于其波形中攜帶了受傳播介質物理力學狀態及微結構(如微裂紋、孔隙等)影響的相關信息，已被廣泛應用于巖土工程測試[14]。我國“水工建筑物巖石基礎開挖工程施工技術規范”以開挖前后波速變化率超過15%作為巖石損傷破壞的閾值[14]。目前，關于巖石室內加載下的波速變化規律已有較多研究。李新平等[15]基于模型試驗研究了在不同圍壓下裂隙巖體在峰值強度前的波速變化規律。GRAHAM 等[16]研究了三軸壓縮下花崗巖波速變化規律，發現初始加載軸向波速增幅遠大于側向加載軸向波速增幅；隨著裂紋擴展，側向波速率先下降，且降幅逐漸增大，在應力峰值附近側向波速下降幅度超過0.17。而花崗巖在單軸壓縮下峰值前波速降幅小于15%[17]。李浩然等[18]研究了花崗巖單軸壓縮下軸向縱橫波速與聲發射演化特征，并分析了二者變化的相關性。PETRU?áLEK等[19]研究了單軸壓縮下巖石聲波及聲發射變化，通過實測波速反演AE震源定位，發現波速及波速模型的選取直接影響AE震源的分布規律。由于花崗巖是一種典型的強度高、脆性突出的巖石，在傳統加載方式下，巖石在峰后往往發生快速劇烈破壞，由于聲波發射頻率和巖石劇烈破壞時產生的強擾動，導致花崗巖應力-應變全過程聲波測試變得困難。目前，聲波研究主要針對波速開展，而關于巖石破裂過程中聲波頻譜特征、應力-應變全過程多參量AE 特征及不同參量間對應關系的研究相對較少。本文通過MTS 加載系統獲得花崗巖單軸壓縮全應力-應變曲線，并同步測試全過程AE 及聲波變化規律。對于AE 測試，分別采用特征參數和波形分析方法，研究巖石破裂全過程AE 事件、能量分布、波形及頻譜特征等變化規律；其次，分析峰前和峰后卸載段不同方向波速及波速各向異性變化規律，研究破裂不同階段聲波頻譜特征變化規律；同時，基于時變波速優化AE震源時空演化分布，結合宏觀變形、AE 及聲波特征綜合分析花崗巖應力-應變曲線全過程的破裂特征。

1 試驗巖樣及技術

試驗巖樣為強度高、脆性突出、質地均勻的中粒徑花崗巖，通過不同方向聲波測試選取波速各向異性較小、直徑×高度為50 mm×100 mm 的標準圓柱形試樣，如圖1所示。花崗巖主要成分為長石、石英和云母，巖樣的平均密度為2 650 kg/m3，縱波波速為3 895～4 015 m/s。

圖1 試驗巖樣Fig.1 Rock samples for testing

試驗采用美國MTS-815 電液伺服實驗機系統進行單軸加載，該系統可同時記錄荷載、應力和應變等。具體加載方式如下：首先采用軸向荷載控制，加載速率為0.4 kN/s；為了獲得脆性花崗巖破壞的全過程應力-應變曲線，當加載接近峰值應力時采用環向變形控制，環向變形加載速率為0.5 μm/s。

聲學測試采用中國科學院武漢巖土力學研究所研發的聲波、聲發射一體化裝置，該系統可實現信號的連續波形采集記錄。本文共設置12 個P 波傳感器，如圖2(a)和(b)所示。圖2中，傳感器10，11 和12 作為超聲激發探頭激發超聲脈沖，傳感器1～9作為聲波和聲發射接收傳感器。聲波測試采用傳統透射法，通過儀器內置電壓信號發射卡向激發通道施加電壓脈沖，依靠激發與接收探頭測試不同方向波速的變化。巖樣側面從傳感器1，5和9所在軸線展開，聲波測試路徑見圖2(c)(其中α為波傳播方向與徑向的夾角)。聲波和聲發射傳感器采用Nano30 探頭，該傳感器對應的頻帶響應寬度為125～750 kHz，通道采集門檻值設為40 dB，前置放大器為40 dB，通道采樣速率為10 MHz。試驗開始之前，檢查傳感器與巖樣的耦合情況，一切正常后才可進行試驗。

2 試驗結果及分析

2.1 應力--應變過程分析

通過測試巖樣的軸向應變ε1及側向應變ε3，計算得到巖樣體積應變εv為

巖樣的變形包含巖石基質的變形及微裂紋引起的變形，在單軸壓縮下，巖石內部裂紋演化引起的變形為

式中：E為彈性模量，υ為泊松比；σ1為單軸壓縮下的軸向應力。

圖2 測試傳感器布置Fig.2 Arrangement of testing sensors

圖3所示為試驗花崗巖A1 的全應力-應變曲線，根據特征應力可劃分為幾個不同階段。首先是賦存在巖石中的天然裂紋及孔隙的閉合，該階段裂紋體應變隨著應力增大而逐漸減小，宏觀表現為應力-應變呈非線性變化，即階段Ⅰ，且該非線性變化程度與原始微裂紋及孔隙的分布密切相關。由于花崗巖的初始壓密主要發生在軸向，該階段橫向變形遠小于軸向變形，軸向應變主導該變形過程。隨著原始裂紋的逐漸壓密閉合，裂紋體應變逐漸減小為零，此時對應裂紋閉合應力σcc。之后進入彈性階段Ⅱ，應力-應變曲線基本保持線性增長，且裂紋體應變保持為零。當應力超過斷裂強度后，微裂紋開始萌生、擴展，裂紋體應變由零開始逐漸增加，對應巖石擴容的起始，即裂紋起裂應力σci(為95 MPa，即0.56σc，σc為峰值應力)。此后，微裂紋開始擴展并進入穩定擴展階段Ⅲ，此階段微裂紋主要沿軸向擴展，橫向變形速率逐漸增大，同時橫向變形增量與軸向變形增量的比值也不斷增加。當應力達到體應變的反彎點(最大值)時，花崗巖由體積壓縮轉變為體積膨脹，此時對應裂紋的損傷應力σcd[6-7](127 MPa，0.74σc)。之后進入裂紋不穩定擴展階段Ⅳ，巖樣擴容速度逐漸增大，且擴容以橫向變形為主導，橫向變形增加速率迅速增大，并逐漸占據主導優勢。當達到峰值強度時，軸向應力-應變曲線快速跌落，呈Ⅱ型破壞，峰后沿軸向擴展的張拉裂紋快速發展，導致軸向應變增量相對較小，而橫向變形則快速增大，進而引起花崗巖體應變的快速膨脹及裂紋體應變的非線性增大。

圖3 單軸壓縮花崗巖全應力-應變曲線Fig.3 Complete stress-strain curves of granite under uniaxial compression

2.2 聲發射特征演化規律

花崗巖破裂過程中AE 事件率及AE 波形對應的能量分布如圖4所示。由圖4可知：在初始加載階段，原始微裂隙的閉合、摩擦導致少量AE事件的產生，此階段AE 事件率極小且呈不穩定變化，AE 能量基本保持在10×10-18J 以下。當應力達到75.5 MPa(0.44σc)時，此時AE事件率顯著增大；同時，能量超過10×10-18J 的AE 信號出現小幅凸起，對應微裂紋的初始萌生擴展，即起裂應力σci；隨著應力的逐漸增加，AE 事件率逐漸增多，大能量AE 數量也逐漸增多，當應力達到135 MPa時，事件率發生陡增，同時伴隨著AE 信號能量的劇烈凸起(對應圖4中標記a)，AE攜帶能量最大值由低于1×10-16J激增到1×10-15J左右。在峰值應力處(對應圖4中標記b)，由于裂紋的快速擴展引起應力快速下降，花崗巖內存儲的應變能迅速釋放，導致AE事件率快速增加及能量激增，但由于采用環向變形控制，花崗巖并未出現劇烈破壞。在經歷快速能量釋放后，花崗巖內部回到短暫平衡狀態，應力也出現小幅增加，并在較短時間后再次出現快速下降，同時引起AE 事件率和AE 能量分布凸起(對應圖4中標記c)。之后應力呈緩慢下降，而AE事件累積速率也逐漸趨緩，此時，裂紋擴展引起的花崗巖損傷發展相對較慢。

圖4 AE事件率及其能量分布Fig.4 Distribution of AE event rate and energy

波形分析是AE信號后處理的另一種方法，波形分析反映的巖石破裂特征更全面、直觀[12]。不同應力下AE 信號的時頻變化如圖5所示。圖5(a)中，初始壓密階段AE 信號主頻約為115 kHz。在σci附近，AE頻譜振幅和頻帶寬度均增大，主頻增加到230 kHz，次主頻則達到310 kHz，頻譜表現出明顯雙峰值特征。圖5(c)所示頻譜對應裂紋穩定擴展階段，此時AE 主頻增加到300 kHz，而對應的次主頻則轉變為230 kHz，頻譜振幅繼續增大。隨著外荷載不斷增加，裂紋擴展速率和AE事件率逐漸增大，在峰值應力附近，AE頻帶寬度迅速增大，頻譜中額外增加了100～200 kHz 的低頻成分和400～500 kHz 的高頻成分，呈現出多峰值特征，且頻譜振幅達到最大。

隨著花崗巖的快速擴容進入峰后段，環向變形控制下花崗巖峰后呈漸進性破壞，應力呈階梯狀下降，圖4中c處(0.76σc)AE 頻譜表現為雙峰值，主頻下降到290 kHz，次主頻減小為240 kHz，對應的頻譜振幅也迅速下降。花崗巖峰后段AE主頻下降的原因如下：一方面，由于裂紋擴展聚集，導致裂紋尺度增大；另一方面，隨著裂紋密度和擴展尺度增加，AE信號在損傷巖石中的衰減逐漸增強，且高頻信號的衰減要比低頻信號的衰減大，兩方面共同作用導致AE主頻不斷減小。在峰后段(0.43σc)附近，AE主頻下降到230 kHz，高頻成分振幅逐漸衰減，頻帶寬度也逐漸變窄。

圖5 不同應力水平下AE波形時頻譜圖Fig.5 Time-frequency spectra of AE waveforms recorded under different stress levels

AE 信號的頻率特性尤其是主頻與微破裂密切相關，加載全過程AE波形峰值頻率及幅值隨時間的變化如圖6所示，其中圓圈直徑與信號幅值呈正比。在加載初始階段，AE 幅值極小，峰值頻率主要集中在87～98 kHz 頻帶。隨著應力增大，低頻信號逐漸較少，相反，頻率大于250 kHz 的高頻信號逐漸產生，起裂應力附近高頻信號數量和比例均快速增加。根據文獻[20]中裂紋類型與AE 特征的聯系，此時高頻(>273 kHz)小幅值AE 對應于沿晶或穿晶微破裂。在損傷應力(見圖6中a點)前，低頻(87～98 kHz)小幅值信號數量逐漸增多，同時，高頻大幅值AE數量也逐漸增多，并在損傷應力附近急劇增加，部分AE幅值高達2 V，而高頻大幅值信號多對應中尺度裂紋，也說明裂紋擴展尺度增大。之后，高頻大幅值AE數量急劇增加，中等尺度裂紋的快速擴展導致裂紋尺度繼續增加，導致峰值應力前出現一定數量的低頻大幅值AE，對應更大程度破裂的發生。

圖6 AE波形峰值頻率及幅值隨時間的變化Fig.6 Change of AE waveform peak frequency and amplitude with time

在峰值應力附近，AE 頻帶寬度也快速增大，增加了頻率小于87 kHz，100～250 kHz(中等頻率)及大于310 kHz 的AE 信號。中頻小幅值信號主要是破裂面的摩擦造成的，因此，伴隨著峰后快速應力降，從小尺度穿晶裂紋到大尺度宏觀裂紋及破裂面的摩擦均有發生，裂紋破裂類型較復雜，該現象從峰值應力一直持續到c點，且花崗巖破裂的大幅值AE 信號主要發生在該階段。c點之后應力緩慢下降，宏觀破裂面已基本形成，AE 事件率快速下降，AE 頻帶寬度減小，其中以中等頻率信號最明顯。

由上述分析可知，峰值頻率在87～97 kHz和273～292 kHz 范圍內的AE 信號貫穿花崗巖破裂全過程，是AE 信號的優勢頻率，且大幅度AE 信號主要集中在該頻帶。由于頻域特征往往具有本征性、唯一性，說明該頻帶與巖石礦物成分相關，是花崗巖的特征頻率[3]。另外，花崗巖破裂全過程大幅值AE 信號主要發生在峰后2次大幅應力降之間，快速應力降對應裂紋的快速擴展，同時導致AE 頻帶寬度及大幅度AE數量快速增多。因此，AE 信號頻率成分的復雜程度與花崗巖破壞劇烈程度和破裂模式密切相關。

2.3 聲波特征變化規律

根據圖2中的傳感器布置，以巖樣A1 為例，研究不同方向波速與應力-應變的關系，如圖7所示。由圖7可知：在加載初期，花崗巖各方向波速差異較小，隨著外荷載增大，波速逐漸呈非線性增加，同時，不同方向波速差異也逐漸增大，波傳播方向與徑向夾角越大，波速增速越大，這是原始微缺陷在偏應力加載下的導向性閉合所致[17]。當應力達到82.5 MPa(圖7中虛線l)，側向波速v0達到峰值，相比初始波速增加了5.8%，而v31(與徑向夾角為31°)、v58(與徑向夾角為58°)分別增加13%和18.8%。此后，v0在較小范圍內保持穩定，之后開始逐漸下降，而v31和v58則繼續增加，并在體應變反彎點附近達到峰值。單軸壓縮下裂紋的導向性起裂、擴展是導致波速變化差異的主要原因[17]。

圖7 不同方向波速與應力-應變關系Fig.7 Relationship between velocities in different directions and stress-strain

每間隔一定時間對聲波進行1次測試，在靠近峰值應力(159.4 MPa，即93%σc，對應圖7中虛線m)前，相比波速峰值，此時v0，v31和v58下降幅度分別為4.5%，2.4%和2.0%，其中側向波速下降幅度最大。在峰值應力附近，由于采用環向變形控制加載，花崗巖并沒有出現快速破壞，巖石損傷累積也繼續呈漸進式發展；在峰后階段，應力緩慢下降，軸向卸載下波速以較小速率持續減小，其變化趨勢與應力變化趨勢基本一致；且峰后段v0下降速率逐漸增大，相比波速峰值，加載結束附近v0，v31和v58下降幅度分別達到17.4%，14.9%和16%。

根據不同方向波速差異，采用測試波速的最大值vmax和最小值vmin來計算巖石波速各向異性系數K[19]：

圖8所示為巖樣A2和A3不同方向波速變化。由圖8可見：巖樣A2 和A3 與A1 波速變化的趨勢基本一致；在初始加載階段，隨著波傳播導向傾角的增加，波速增幅基本呈逐漸增大趨勢；同時，波速各向異性系數K也不斷增大，其中巖樣A3增加速率較大，這與巖樣原始微裂隙的幾何分布有關。

圖8 巖樣A2和A3不同方向的波速變化Fig.8 Change of velocities in different directions for samples A2 and A3

隨著應力增加，側向波速v0逐漸達到峰值，然后開始減小，而K的增加速率逐漸增大。在峰值應力前，各方向波速均開始下降；在峰值應力處(圖8中σcA2和σcA3)，巖樣A3的側向波速下降幅度最大，導致K快速增加。對于巖樣A2，峰值應力后各方向波速下降速率均逐漸增大，尤其是v0和v45(與徑向的夾角為45°)；但當加載到1 400 s 時，v45和v58下降速率比v0的大，說明裂紋此時局部擴展加劇。當加載結束時，巖樣A2各方向波速下降幅度為0.30～0.44，而A3波速下降幅度為0.13～0.36。

聲波信號包含了大量與巖石物理力學特征相關的信息。為了分離反射脈沖，選取首個到達的脈沖來分析加載過程中的衰減。根據已有研究，采用單一振幅的半余弦函數對接收的波形進行處理[21-22]，如圖9(a)所示。由于微裂紋的萌生擴展主要沿軸向擴展，這里僅對側向透射波進行處理，如圖9(b)所示(巖樣A1)。

圖9 不同應力下側向透射波頻譜特征變化Fig.9 Amplitude spectra of lateral transmitted waves at different stress levels

當應力約為2 MPa 時，透射波主頻為213 kHz，伴隨著應力的增加，原始微裂隙及微孔洞逐漸壓密閉合，透射波主頻呈小幅增加。當應力達到62 MPa時，主頻增加到244 kHz，同時，頻譜幅值從0.46 V 增加到0.58 V，增幅更大，此時，巖石內部微裂紋基本閉合(見圖3)。當應力達到82 MPa 時，主頻基本不變，頻譜幅值出現小幅增加。由于AE事件率逐漸呈穩定增加(見圖4)，微裂紋開始萌生擴展，進而導致透射波頻譜幅值快速下降，在峰值應力前(160 MPa)，頻譜幅值下降到0.45 V。進入峰后階段，由于裂紋繼續擴展，透射波頻譜幅值繼續下降，同時伴隨的是主頻小幅減小。當峰后應力減小到56 MPa 時，主頻下降到182 kHz，頻譜幅值則減小為0.29 V。

與A1相比，巖樣A3透射波頻譜特征在不同階段差異更大。同樣，在裂紋初始壓密階段，透射波主頻和幅值逐漸增大。在閉合應力附近達到峰值，之后由于裂紋的孕育、萌生擴展，透射波主頻和幅值快速下降。當應力達到101 MPa時，頻譜幅值由0.93 V快速下降到0.57 V，此后一直持續到峰值應力，并在峰前(148 MPa)減小到0.24 V，峰前透射波的頻譜特征下降幅度遠比巖樣A1的下降幅度大。進入到峰后階段，由于宏觀裂紋形成，導致側向波速快速下降(見圖8)，同時，透射波頻譜幅值下降到較小值(0.09 V)。經對比發現，透射波頻譜幅值變化幅度明顯大于波速和頻譜主頻的變化幅度。

2.4 AE震源時空演化規律

AE 測試技術的另一個優勢是能夠有效追蹤裂紋擴展區域[4]。這里采用蓋格爾時差定位法[23]分析花崗巖破裂全過程AE震源演化分布。基于實時測試的波速變化，將AE 震源演化分為3 個階段，不同應力水平下AE震源分布如圖10所示，其中圓圈直徑與震源幅值呈正比。

由圖10可見：在起裂應力出現之前，出現少量小幅值AE 震源，并呈隨機分布。在起裂應力階段后，裂紋率先在花崗巖上端部萌生，震源釋放速率逐漸增大，同時大幅值AE震源數量逐漸增多。在峰值應力處，大振幅震源在花崗巖中上部聚集成核，宏觀表現為裂紋擴展的局部化，初步沿軸向向下形成1條破裂帶。花崗巖的局部化破裂成核導致應力下降，即進入峰后階段。首先，震源繼續沿巖石中上部擴展分布，同時逐漸向巖石下端部擴展，由于峰后卸載下花崗巖屬于漸進性破壞，AE 震源的幅值呈逐漸下降趨勢，大幅值震源的數量明顯比峰值應力附近的大幅值震源數量少。

2.5 不同參量特征對比分析

裂紋體應變是巖石破裂過程內部微裂紋演化的宏觀表現，而聲波和聲發射特征與巖石斷裂及損傷演化密切相關[1,19]。波速與AE 事件率隨時間的變化如圖11所示。由圖11可見：在加載初期，原始微裂紋的閉合及沿裂紋壁面的滑移引起AE事件率間歇性地增加，而波速則逐漸增加；當應力達到75.5 MPa 時，AE 事件率顯著增大，對應微裂紋的起裂，但波速并未開始下降，而是繼續緩慢增加，直到101 MPa 附近v0才開始逐漸減小。在標記點a處(見圖4和圖6)，由于裂紋加速擴展，AE 事件率快速增加，同時，波速下降速率增大，尤其在峰值應力處，裂紋的快速擴展引起AE事件率增大和波速加速下降。峰后段的2次顯著應力降(峰值應力處和標記點d)均伴隨著波速快速下降，但相比AE事件率的變化幅度，波速下降幅度較小，變化更平緩。第2 次顯著應力降(標記點d)之后，AE 事件率迅速減小并基本呈穩定發展，v31和v58下降速率也逐漸平緩，但v0繼續以較大速率減小，此時，v0對裂紋擴展的響應比AE 事件對裂紋擴展的響應大。

圖10 不同階段AE震源演化分布Fig.10 Spatiotemporal distribution of AE hypocenters at different stages

圖11 AE事件率與波速隨時間的變化Fig.11 Change of AE event rate and wave velocities with time

圖12所示為宏觀變形與累計AE事件隨時間的變化。由圖12可見：累計AE事件增加的起點出現在花崗巖宏觀線彈性段(即裂紋體應變為零)，明顯早于裂紋體應變增加的起點(擴容起始點)，這表明花崗巖應力-應變曲線的宏觀線彈性段已經出現微破裂的孕育、萌生。擴容起始點后，AE 事件和裂紋體應變均呈非線性增加，在體應變的反彎點(體應變最大值)附近，AE 事件率迅速增大，這與裂紋體應變的變化趨勢基本一致。在峰后階段，由于花崗巖快速擴容，體應變發生陡降并迅速歸零，裂紋體應變也出現跳躍式增加，但AE事件變化速率并未出現大幅增大。值得關注的是，第2 次大幅應力降之后，AE 事件率繼續增大，但增大速率明顯減小，而體應變和裂紋體應變繼續呈大速率增加，二者對裂紋擴展的響應明顯比AE事件對裂紋擴展的響應大。

圖12 裂紋體應變與累計AE事件隨時間的變化Fig.12 Change of crack volumetric strain and accumulated AE event with time

3 討論

本次試驗加載結束時，所有巖樣波速下降幅度范圍為0.18～0.44，且波速降低主要發生在峰后破裂階段。在軸向位移控制加載條件下，強度高、脆性突出的實驗花崗巖多在峰值強度后很短時間內就發生劇烈破壞，巖石支撐結構的損傷也主要發生在該階段。由于聲波測試受到發射頻率和強擾動等的影響，導致花崗巖劇烈破壞階段聲波測試存在一定困難。在環向變形控制加載條件下，脆性花崗巖在破壞前的波速下降幅度能達到較大值，巖樣加載結束時最大波速降幅達0.44。波速降η與損傷D的量化關系為[24]

加載結束時，基于波速量化的花崗巖損傷達到0.69。另外，值得說明的是，由于室內巖石試件聲波傳感器布置存在數量限制，只能測試幾條指定方向的波速變化，而脆性巖石的破壞多是由微破裂的彌散型分布轉變為聚集、局部化的分布[9]，因此，在此基礎上對巖石波速各向異性的分析也只是一種簡化分析。

AE 事件是空間分布的多個傳感器對破裂源的反演定位，事件率的顯著增大對應微裂紋的萌生擴展，而裂紋的擴展則會引起波傳播衰減，但測試結果顯示AE增加早于波速下降，且在峰值應力前，AE事件對局部裂紋快速擴展的響應比波速的大。其原因是：1)巖石試件聲波傳感器布置存在限制，只能測試有限幾條路徑區域的波速，而巖石的破裂多呈局部化，指定聲波測試路徑不一定恰好穿過裂紋萌生起裂及高度發育區；而聲發射則是對巖石全局破裂進行探測；2)聲波測試主頻約為250 kHz，對應的波長約為15 mm，裂紋初始起裂長度遠小于該波長，根據波傳播理論，遠小于波長的微裂紋擴展對波傳播影響極小，而聲發射對破裂初始的小尺度微裂紋的敏感性極強。另外，對于花崗巖峰后卸載破壞段，AE 事件率下降較明顯，側向波速v0則繼續以較大速率下降，可見AE 事件表征的花崗巖漸進性破壞過程反而更明顯；而宏觀體應變和裂紋體應變受峰后裂紋擴展影響后的變化幅度更大，巖石損傷累積更明顯。這是因為峰后卸載段宏觀破裂面已經形成(見圖9)，環向變形控制下破裂面的擴展張開引起的破裂強度并非十分突出(AE 事件增加速率相對較小)，而側向波速和體應變尤其是橫向變形對破裂面的張開十分敏感，導致宏觀變形和波速的變化反而更加劇烈。因此，不同特參量表征的花崗巖破裂過程及活動性與花崗巖破裂模式和破壞劇烈程度有關。

通過對巖石破裂過程中聲波透射波進行處理獲得不同應力狀態下的頻譜曲線，結果表明透射波頻譜特征，尤其是頻譜幅值受巖石裂紋演化的影響較大。在花崗巖的初始加載階段，微裂紋的壓密閉合導致頻譜主頻和幅值增加，而隨著應力增大，新生裂紋的孕育、萌生擴展則引起特征值快速下降。同時，對比波速和頻譜特征值變化發現，頻譜幅值受裂紋演化的變化幅度要比波速的變化幅度大。

4 結論

1)花崗巖破裂全過程AE頻譜表現出雙峰值或多峰值，AE 主頻及幅值先增后減，主頻由初始加載段的115 kHz 增加到穩定破裂段的300 kHz，峰后快速下降到230 kHz。AE事件率的峰值與能量分布的凸起較吻合，對應裂紋加速擴展，應力突降伴隨著AE事件率的劇烈凸起和能量快速釋放。

2)在加載初期，AE 峰值頻率集中在87～97 kHz，起裂應力σci附近，高頻小幅值AE快速增加；損傷應力σcd后高頻大幅值AE迅速增多，對應中尺度裂紋擴展；峰值應力σc附近AE頻譜帶寬快速增大，低頻大幅值AE急劇釋放，對應大尺度破裂的發生；峰后應力快速降低，從小尺度裂紋到大尺度宏觀主裂紋相互伴生；峰值頻率為87～97 kHz 與273～292 kHz 的AE貫穿破裂全過程，且大幅值AE 主要集中在該頻帶內，是花崗巖的特征頻率，AE 頻率成分的復雜程度與破裂劇烈程度和模式密切相關。

3)單軸壓縮下波速先增加后減小，橫向波速率先下降且降幅最大，破壞前波速降幅高達0.44，裂紋擴展的導向性引起波速各向異性逐漸增大；AE 震源演化較好地描繪了裂紋擴展分布區域，峰值應力附近局部大幅值震源聚集成核，是花崗巖破壞的誘因。

4)超聲透射波頻譜幅值和主頻隨著微裂紋的壓密閉合而逐漸增大，之后，由于裂紋的萌生擴展而逐漸減小，且頻譜幅值受裂紋演化的變化幅度遠比波速和頻譜主頻的大。

5)AE 事件在花崗巖的宏觀線彈性段開始穩定增加，此時微裂紋開始孕育萌生，AE 事件增加早于波速下降，在峰前對裂紋擴展的響應及敏感度比波速的大；但在峰后卸載段，宏觀變形和波速對裂紋擴展的響應反而比AE事件對裂紋擴展的響應大，這與花崗巖峰后的破裂模式直接相關。