楊 陽，吳賢振，劉 浩，周伶杰

(1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州 341000;2.江西省礦業(yè)工程重點(diǎn)實(shí)驗(yàn)室,江西 贛州 341000)

基于歐氏距離的單軸壓縮下粉砂巖熱圖像演化特性研究

楊 陽1，2，吳賢振1，2，劉 浩1，2，周伶杰1，2

(1.江西理工大學(xué)資源與環(huán)境工程學(xué)院，江西 贛州 341000;2.江西省礦業(yè)工程重點(diǎn)實(shí)驗(yàn)室,江西 贛州 341000)

為探尋粉砂巖破裂失穩(wěn)過程紅外異常，引入歐氏距離的方法，通過計(jì)算粉砂巖試樣各相鄰時(shí)刻歐氏距離，對(duì)單軸壓縮條件下粉砂巖紅外輻射溫度場演化特性進(jìn)行研究，結(jié)合熵值、方差兩種指標(biāo)演化趨勢驗(yàn)證歐氏距離方法的可用性。結(jié)果表明：歐氏距離、熵值法和方差均能較好地刻畫粉砂巖紅外輻射溫度場各個(gè)階段；壓密階段至彈性階段前期，歐氏距離存在較多大幅度突升值，表征這一時(shí)段內(nèi)試件自身存在的微裂隙被壓縮，產(chǎn)生熱摩擦效應(yīng)，歐氏距離與熵值、方差刻畫能力相當(dāng)；彈性階段后期至塑性階段，歐氏距離呈較為平穩(wěn)起伏交替，表征裂隙穩(wěn)定發(fā)育，歐氏距離刻畫能力劣于熵值、方差；峰后階段，歐氏距離出現(xiàn)大幅度突升，且幅度、頻率均大于前面的階段，表征粉砂巖試樣即將破裂，歐氏距離刻畫能力優(yōu)于熵值、方差。

破裂失穩(wěn)；紅外熱圖像；歐氏距離；熵值；方差

中國礦產(chǎn)資源儲(chǔ)量豐富，而在礦產(chǎn)資源開發(fā)過程中面臨的災(zāi)害問題正成為制衡資源開采的一大難關(guān)。巖石破裂災(zāi)害是引發(fā)礦山災(zāi)害的主要原因，因此礦山資源開采過程中的巖石破裂監(jiān)測預(yù)警是其防災(zāi)減災(zāi)的關(guān)鍵，成為國內(nèi)外研究的焦點(diǎn)。紅外監(jiān)測技術(shù)具有實(shí)時(shí)、連續(xù)、非接觸等優(yōu)點(diǎn)，在巖石災(zāi)變預(yù)警領(lǐng)域正受到越來越多學(xué)者的青睞。

過去的幾十年間，許多學(xué)者對(duì)不同條件下巖石破裂的紅外輻射異常前兆進(jìn)行了廣泛的研究并取得了一定的成就。文獻(xiàn)[1]通過實(shí)驗(yàn)指出，煤巖和砂巖在受壓破裂期間產(chǎn)生了熱圖像和平均紅外溫度(AIRT)異常先兆。文獻(xiàn)[2]通過對(duì)花崗巖受壓破裂紅外熱成像實(shí)驗(yàn)，分別應(yīng)用特征粗糙度、熵、方差作為指標(biāo)來描述紅射溫度場的演化特性。文獻(xiàn)[3]通過實(shí)驗(yàn)研究從升溫、降溫以及幅值、規(guī)模四個(gè)方面對(duì)受壓巖石的紅外輻射溫度場進(jìn)行研究，得出破壞前粉砂巖紅外輻射溫度出現(xiàn)突變異常先兆，其中包含突升異常等四種紅外輻射溫度場異常前兆。文獻(xiàn)[4]通過實(shí)驗(yàn)研究，提出‘紅外溫變場’的概念，指出飽水粉砂巖破壞失穩(wěn)進(jìn)程中紅外溫變場的特征參數(shù)出現(xiàn)陡增的變化異常，同時(shí)提出在延性巖石破壞過程中ITF變化的異常現(xiàn)象具有一定的普遍性。文獻(xiàn)[5]通過實(shí)驗(yàn)研究，將紅外溫度場劃分為包括高溫區(qū)域在內(nèi)的三個(gè)溫度區(qū)域，并以高溫區(qū)域面積比、溫均面積比等形式作為表征形式描述紅外高溫區(qū)域異常前兆。文獻(xiàn)[6]通過實(shí)驗(yàn)，對(duì)潮濕砂巖和干燥砂巖單軸壓縮條件下的紅外輻射特征進(jìn)行了比較；前者平均紅外輻射溫度(AIRT)與應(yīng)力變化步調(diào)一致，曲線波動(dòng)較??；后者平均紅外輻射溫度(AIRT)與應(yīng)力不是簡單的對(duì)應(yīng)關(guān)系，但在破壞瞬間有顯著的增溫現(xiàn)象和熱像特征的明顯變化。文獻(xiàn)[7]通過實(shí)驗(yàn)研究，指出試件加載進(jìn)程中裂隙發(fā)展越劇烈其紅外溫度改變也越大；在試件破壞時(shí)，會(huì)在部分破壞區(qū)域出現(xiàn)高溫輻射條帶的紅外現(xiàn)象。文獻(xiàn)[8-9]結(jié)合實(shí)驗(yàn)研究和數(shù)值模擬，對(duì)試件應(yīng)力狀態(tài)與紅外溫度場變化情況進(jìn)行分析，得出裂紋集中發(fā)育區(qū)域升溫顯著，在試件破裂處呈現(xiàn)出顯著的高溫區(qū)和高溫條帶。

目前針對(duì)巖石加載破裂過程中的紅外溫度場變化情況描述大多基于平均紅外溫度、最大紅外溫度、最小紅外溫度等紅外溫特征參數(shù)變化以及熱圖像來表征巖石內(nèi)部情況，由于巖石結(jié)構(gòu)、物理力學(xué)性質(zhì)的不同，AIRT指標(biāo)只能表示紅外輻射的整體水平，并不能準(zhǔn)確得出巖石破裂的紅外前兆。隨后劉善軍等通過花崗巖受壓破裂紅外溫度觀測實(shí)驗(yàn)提出了用特征粗糙度、熵、方差作為指標(biāo)來分析紅外溫度場的演化特性；吳賢振等通過粉砂巖受壓破裂紅外實(shí)驗(yàn)提出用突升異常、突降異常與突變極值異常、突變點(diǎn)量異常作為巖石破裂的紅外異常前兆特征。

本文根據(jù)裂隙發(fā)育導(dǎo)致巖石溫變這一理論，引入歐氏距離的方法分析粉砂巖在受壓破裂過程中溫度場整體變化情況與應(yīng)力、時(shí)間之間的關(guān)系，并與熵值法、方差進(jìn)行對(duì)比研究，其研究結(jié)果能對(duì)巖石受壓破裂紅外前兆的確定提供一定的指導(dǎo)和推廣作用。

1 基于歐式距離的紅外溫度場表示方法

歐式距離是應(yīng)用范圍最為廣泛的一種距離定義，指在空間坐標(biāo)系中兩個(gè)點(diǎn)或者向量之間的絕對(duì)距離。將巖石受壓破裂過程中同一時(shí)刻試樣表面所有測得的溫度數(shù)據(jù)按照從上到下、從左到右的順序構(gòu)成某一時(shí)刻的溫度向量，由于剪切破裂發(fā)育導(dǎo)致紅外溫度升高、張性破裂發(fā)育導(dǎo)致紅外溫度降低，巖石試樣相鄰時(shí)刻間溫度向量間的歐氏距離表征單位時(shí)間內(nèi)巖石試樣裂隙的發(fā)育情況。歐氏距離由式(1)定義。

(1)

式中：d為相鄰時(shí)刻兩溫度向量的歐氏距離，t為不同時(shí)刻，k為巖石試樣測溫點(diǎn)編號(hào)，xtk為t時(shí)刻k號(hào)測溫點(diǎn)溫度。

由于整個(gè)巖石加載破裂過程中紅外輻射溫度場出現(xiàn)的溫變現(xiàn)象，主要是由于剪切破裂與張性破裂發(fā)育導(dǎo)致的，相鄰時(shí)刻溫度向量的歐氏距離越大表示裂隙發(fā)育越劇烈，反之則是裂隙發(fā)育越輕微，所以可以利用相鄰時(shí)刻溫度向量的歐氏距離表征紅外溫度場的溫變現(xiàn)象。

2 方法的驗(yàn)證

為證實(shí)歐氏距離方法在表征巖石破裂過程中紅外溫度場溫變現(xiàn)象的可行性，進(jìn)行了粉砂巖加載紅外熱成像試驗(yàn)。

2.1 試驗(yàn)設(shè)計(jì)

采用粉砂巖作為試驗(yàn)試樣，將其加工成50mm×50mm×100mm規(guī)格的標(biāo)準(zhǔn)柱體試件3塊，分別編號(hào)為S1、S2、S3，打磨上下兩端平面，使上下兩端面平整度符合要求。分別將3個(gè)試件放置在50℃的溫水中浸泡半個(gè)小時(shí)，取出待試件表面干燥后座密封處理。試驗(yàn)采用單軸壓縮，在RLW-3000型伺服試驗(yàn)系統(tǒng)上完成，使用位移控制方式加載，先預(yù)加壓到1.5kN,然后以1.5mm/min恒速加載，直到試件破壞。使用美國制造的SC3000熱成像儀對(duì)試件進(jìn)行正面實(shí)時(shí)監(jiān)測，得到試件加載直至破壞的熱圖像。

2.2 試驗(yàn)數(shù)據(jù)處理

根據(jù)需要，對(duì)試驗(yàn)中取得的熱圖像進(jìn)行1幀/s的重采樣，并按照時(shí)間升序排成紅外溫度矩陣序列；然后應(yīng)用matlab軟件將每個(gè)時(shí)刻紅外溫度矩陣按照從上到下、從左到右的順序轉(zhuǎn)化為紅外溫度向量，通過編程分別求解出三個(gè)試件時(shí)間升序紅外溫度向量間的歐氏距離、熵值、方差，并對(duì)得到的數(shù)據(jù)進(jìn)行定量分析。

2.3 粉砂巖紅外輻射溫度場歐氏距離

歐氏距離表示的是相鄰時(shí)刻試件紅外溫度場變化情況，也就是裂隙的發(fā)育程度，求出每個(gè)相鄰時(shí)刻試件紅外溫度場的歐氏距離，見圖1。

圖1 粉砂巖歐氏距離和應(yīng)力隨時(shí)間變化曲線

根據(jù)圖1，相鄰時(shí)刻試件紅外溫度場歐氏距離呈上下起伏交替，當(dāng)試件完全破壞后歐氏距離出現(xiàn)一個(gè)峰值并迅速回落至較低程度。加載前期，粉砂巖位于壓密-彈性前期階段，歐氏距離在該階段存在較多突升值，突升值變化越大表示相鄰時(shí)刻試件紅外溫度場變化程度越劇烈，主要是由于巖石內(nèi)部裂隙被壓密，在壓密過程中產(chǎn)生熱摩擦效應(yīng)，導(dǎo)致壓密位置溫度上升，從而使得歐氏距離出現(xiàn)較多的突升值。A到B點(diǎn)期間，粉砂巖處于彈性后期-塑性變形階段，歐氏距離在該階段表現(xiàn)為平穩(wěn)的起伏交替，表示該期間相鄰時(shí)刻試件紅外溫度場變化較為平穩(wěn)，剪切破裂和張性破裂發(fā)育比較穩(wěn)定。當(dāng)施加應(yīng)力到達(dá)峰值點(diǎn)C之后，歐氏距離突升值再次增多，隨著繼續(xù)加壓，粉砂巖內(nèi)部剪切破裂和張性破裂大量發(fā)育，歐氏距離的突升值幅度和頻率大量增加，在試件S1上壓力峰值點(diǎn)C之后一段時(shí)間歐氏距離數(shù)值產(chǎn)生劇烈波動(dòng)，并且其幅值遠(yuǎn)遠(yuǎn)大于前面的時(shí)刻，在S2、S3試件上表現(xiàn)的不是很明顯。D點(diǎn)以后階段，由于裂隙相互交叉且相互聯(lián)合形成宏觀斷裂面，歐氏距離最后時(shí)刻出現(xiàn)的峰值點(diǎn)說明粉砂巖試件已經(jīng)破壞，峰值點(diǎn)出現(xiàn)主要是由于粉砂巖試件破壞后內(nèi)部剪切裂隙和張性裂隙基本不再發(fā)育，試件通過與空氣進(jìn)行熱交換使得溫度迅速下降，導(dǎo)致歐氏距離出現(xiàn)峰值點(diǎn)。

2.4 粉砂巖紅外溫度場熵值

熵值最先被應(yīng)用于信息論中，被學(xué)者廣泛用來判斷一個(gè)指標(biāo)的離散水平。對(duì)粉砂巖加載至破壞過程中每張熱圖像在時(shí)間序列下的熵值，熵值越大熱圖像上的溫度分布起伏越大，反之則越小。求出試件紅外溫度場熵值隨時(shí)間的變化圖，作出圖2。

圖2 粉砂巖熵值和應(yīng)力隨時(shí)間變化曲線

根據(jù)圖2，粉砂巖試樣熵值隨著時(shí)間序列呈上下起伏交替，當(dāng)試件完全破壞后熵值產(chǎn)生一個(gè)大幅度降低。加載初期，粉砂巖試樣處于壓密階段，熵值迅速增加，主要是由于剪切裂隙和張性裂隙發(fā)育位置溫度快速突變，導(dǎo)致熵值增大。當(dāng)粉砂巖試件處于彈性后期-塑性階段時(shí)，熵值緩慢增加，說明此時(shí)裂隙發(fā)育較為平穩(wěn)，粉砂巖試樣溫度突變速度放緩。當(dāng)粉砂巖試樣處于應(yīng)力峰值點(diǎn)C后階段，熵值劇烈增加且增加速度遠(yuǎn)遠(yuǎn)大于前面的階段，主要是裂隙大量發(fā)育導(dǎo)致試樣部分區(qū)域溫度劇烈升高，熵值增長速度加快。在D點(diǎn)附近粉砂巖試樣被破壞，內(nèi)部裂隙發(fā)育迅速減小，試樣間的熱傳遞使得熵值迅速降低。

2.5 粉砂巖紅外輻射溫度場方差

1918年，羅納德·費(fèi)雪提出了方差的概念，用來在概率論中度量隨機(jī)變量和其數(shù)學(xué)期望之間的偏離程度。對(duì)粉砂巖加載至破壞過程中每張熱圖像在時(shí)間序列下的方差，方差越大熱圖像上的溫度與其平均溫度的偏離程度越大，反之則越小。求出試件紅外溫度場方差隨時(shí)間變化圖，作出圖3。

圖3 粉砂巖方差和應(yīng)力隨時(shí)間變化曲線

根據(jù)圖3，粉砂巖試樣處于壓密階段，方差迅速增加，說明圖像上的溫度與其平均溫度的偏離程度迅速增大。當(dāng)粉砂巖試樣進(jìn)入彈性后期-塑性階段區(qū)間時(shí)，方差開始緩慢降低，主要是由于裂隙發(fā)育較為緩慢導(dǎo)致。當(dāng)壓力繼續(xù)增加，進(jìn)入塑性-峰后階段方差迅速降低，此時(shí)因?yàn)榱严断嗷ソ徊媲蚁嗷ヂ?lián)合形成宏觀斷裂面，促使方差急劇下降。在D點(diǎn)之后，方差發(fā)生一次較大幅度的突變，試件1、試件3方差升高，試件2方差降低。

3 討論

巖石加載至破壞過程中，一般經(jīng)過彈性階段、塑性階段、塑性—峰后階段。在巖石加載過程中會(huì)產(chǎn)生微裂隙的發(fā)育，而微裂隙發(fā)育會(huì)導(dǎo)致紅外溫變，其中剪切裂隙的發(fā)育會(huì)導(dǎo)致紅外溫度升高、張性裂隙的發(fā)育會(huì)導(dǎo)致紅外溫度降低。因此隨著巖石加載至破壞，巖石試樣紅外輻射溫度場會(huì)產(chǎn)生溫度突變和離散現(xiàn)象，從而導(dǎo)致相鄰時(shí)刻紅外溫度場的歐氏距離增大。本文的粉砂巖加載過程紅外熱像試驗(yàn)表明，歐氏距離能較好的反應(yīng)粉砂巖紅外溫度場時(shí)間序列下的溫變情況，同時(shí)也能反應(yīng)出裂隙實(shí)時(shí)發(fā)育情況，定量的描述粉砂巖紅外溫度場的變化情況。

針對(duì)歐氏距離與熵值法、方差兩種指標(biāo)對(duì)紅外輻射溫度場刻畫的情況可以得出，歐氏距離對(duì)于粉砂巖破裂紅外前兆的識(shí)別能力更強(qiáng)，在粉砂巖臨近破裂時(shí)歐氏距離會(huì)出現(xiàn)頻繁的較大值，特別在S1試件的圖像上可以看出臨破裂時(shí)的歐氏距離遠(yuǎn)遠(yuǎn)大于前面各個(gè)階段。而熵值法在反映粉砂巖紅外溫度場的階段性變化特征效果更好，因此可以考慮結(jié)合熵值法與歐氏距離兩種方法對(duì)巖石破裂紅外溫度場各階段特性進(jìn)行綜合描述。

4 結(jié)論

本文提出歐氏距離的方法對(duì)巖石加載至破壞各階段紅外輻射溫度場變化情況進(jìn)行研究，同時(shí)通過與熵值法、方差兩種刻畫紅外輻射溫度場的指標(biāo)進(jìn)行比較，得出以下結(jié)論。

1)整體來看，歐氏距離、熵值法和方差均能較好地粉砂巖紅外輻射溫度場各個(gè)階段的性質(zhì)，歐氏距離對(duì)于時(shí)間序列下裂隙的發(fā)育情況描述更為具體，同時(shí)對(duì)粉砂巖紅外溫變前兆描述更為清晰。

2)壓密階段至彈性階段前期，歐氏距離數(shù)值存在較多大幅度突升值，表征這一時(shí)段內(nèi)試件自身存在的微裂隙被壓縮，產(chǎn)生熱摩擦效應(yīng)，導(dǎo)致紅外輻射溫度場歐氏距離大幅度突升。熵值、方差兩種指標(biāo)在該階段存在明顯上升，與歐氏距離的刻畫趨勢是一致的，表明歐氏距離可以較好的刻畫該階段紅外輻射溫度場的性質(zhì)。

3)在彈性階段后期至塑性階段，歐氏距離數(shù)值呈現(xiàn)較為平穩(wěn)的起伏交替，表征這一時(shí)段內(nèi)裂隙發(fā)育比較穩(wěn)定，紅外輻射溫度場歐氏距離不會(huì)產(chǎn)生較大的變化。在該階段熵值緩慢上升、方差緩慢下降，與歐氏距離的刻畫趨勢基本一致，但刻畫能力稍好于歐氏距離。

4)在峰后階段，歐氏距離呈現(xiàn)大幅度突升，且突升程度和頻率均大于前面的階段，此時(shí)說明裂隙劇烈發(fā)育，促使紅外輻射溫度場歐氏距離劇烈增大，其后隨著繼續(xù)加壓裂隙相互交叉且相互聯(lián)合形成宏觀斷裂面，粉砂巖試樣完全破壞。在該階段熵值加速上升、方差加速下降，與歐氏距離刻畫趨勢是一致的，但刻畫效果沒有歐氏距離顯著。

5)綜合上述各階段三種方法刻畫效果，歐氏距離在試樣受壓破裂過程中裂隙發(fā)育和紅外前兆判定方面具有顯著的優(yōu)勢，研究中可以考慮將歐氏距離與熵值方法結(jié)合起來應(yīng)用。

[1] 吳立新,王金莊. 煤巖受壓紅外熱象與輻射溫度特征實(shí)驗(yàn)[J]. 中國科學(xué)(D輯:地球科學(xué)),1998(1):41-46.

[2] 劉善軍,魏嘉磊,黃建偉,等. 巖石加載過程中紅外輻射溫度場演化的定量分析方法[J]. 巖石力學(xué)與工程學(xué)報(bào),2015(S1):2968-2976.

[3] 吳賢振,高祥,劉祥鑫,等. 飽水粉砂巖破裂過程中紅外溫度突變異常[J]. 煤炭學(xué)報(bào),2015(S2):328-336.

[4] 吳賢振,高祥,趙奎,等. 巖石破裂過程中紅外溫度場瞬時(shí)變化異常探究[J]. 巖石力學(xué)與工程學(xué)報(bào),2016(8):1578-1594.

[5] 高祥,吳賢振,徐金輝,等. 粉砂巖破裂失穩(wěn)過程的紅外高溫區(qū)域異常前兆[J]. 中國礦業(yè),2015,24(12):172-176.

[6] 劉善軍,吳立新,張艷博,等. 潮濕巖石受力過程紅外輻射的變化特征[J]. 東北大學(xué)學(xué)報(bào)：自然科學(xué)版,2010(2):265-268.

[7] 董玉芬,王來貴,劉向峰,等. 巖石變形過程中紅外輻射的實(shí)驗(yàn)研究[J]. 巖土力學(xué),2001(2):134-137.

[8] 李國良. 加載巖石紅外輻射實(shí)驗(yàn)與數(shù)值模擬對(duì)比分析[J]. 河北理工大學(xué)學(xué)報(bào)：自然科學(xué)版,2008(3):5-9.

[9] 譚志宏,閆江偉,韓正林. 裂紋擴(kuò)展過程紅外熱像實(shí)驗(yàn)及數(shù)值模擬研究[J]. 河南理工大學(xué)學(xué)報(bào)：自然科學(xué)版,2011(4):401-405.

Study on thermal image evolution characteristics of silt sandstone under uniaxial compression based on euclidean distance

YANG Yang1，2，WU Xianzhen1，2，LIU Hao1，2，ZHOU Lingjie1，2

(1.College of Resources and Environmental Engineering，Jiangxi University of Science and Technology，Ganzhou 341000,China；2.Jiangxi Provincial Key Laboratory of Mining Engineering,Ganzhou 341000,China)

In order to explore the infrared anomaly in the process of siltstone failure, the method of Euclidean distance was introduced. By calculating the Euclidean distance of siltstone samples, the evolution characteristics of siltstone infrared radiation temperature field under uniaxial compression were studied. Entropy and variance are used to verify the availability of Euclidean distance method. The results show that the Euclidean distance, the entropy method and the variance can describe the different stages of the infrared radiation temperature field of siltstone well. From the compaction stage to the early stage of the elastic stage, the Euclidean distance has a large sudden rise value, The Euclidean distance is similar to the entropy value and the variance characterization. In the late stage of elasticity to the plastic stage, the Euclidean distance is relatively smooth and alternated alternately, which indicates that the fracture is stable and stable. And the range and frequency are larger than those in the previous stage, which indicate that the siltstone sample is about to burst, and the Euclidean distance characterization is better than the entropy value,variance.

fracture instability;infrared thermal image;euclidean distance;entropy;variance

2016-10-12

江西省教育廳科技計(jì)劃項(xiàng)目資助(編號(hào)：GJJ12336);江西省研究生創(chuàng)新專項(xiàng)資金項(xiàng)目資助(編號(hào)：YC2016-S303)

楊陽(1989-)，男，碩士研究生，主要從事巖石力學(xué)，采礦理論與技術(shù)研究，E-mail：625406641@qq.com。

TD315

A

1004-4051(2017)03-0132-04