李新寶，陳 霞，柴緒令，方玉滿，郭民龍

（1.吉林省水利水電勘測(cè)設(shè)計(jì)研究院，吉林 長春 130012；2.遼寧省冶金地質(zhì)勘查局地質(zhì)勘查研究院，遼寧 鞍山 114038）

0 引言

當(dāng)巖體存在破碎或局部不連續(xù)時(shí)，破壞了巖體的均一性，降低了強(qiáng)度，工程中需要灌漿以達(dá)到增強(qiáng)巖體完整性，提高巖體強(qiáng)度的目的。灌漿是否達(dá)到預(yù)定的效果，需要采用物探方法進(jìn)行檢測(cè)，彈性波CT技術(shù)就是有效的檢測(cè)方法之一。

彈性波CT方法，又稱彈性波層析成像技術(shù)，是通過在孔中發(fā)射和接收彈性波，通過一定的觀測(cè)采集方式和數(shù)據(jù)處理方法，得到彈性波在巖體內(nèi)傳播速度信息的一種技術(shù)方法。由于彈性波傳播速度的大小與巖體的堅(jiān)硬和完整程度相關(guān)，同種巖性下，堅(jiān)硬且完整的石具有較高的波速，因此，利用波速分布信息的成像結(jié)果，可以分析和確定出低速區(qū)域，評(píng)價(jià)灌漿效果，或確定灌漿位置。

1 基本原理

彈性波CT成像，首先通過扇形測(cè)試獲得大量的首波走時(shí)數(shù)據(jù)，然后通過求解大型矩陣方程獲取兩孔之間速度剖面圖像。設(shè)在成像剖面內(nèi)測(cè)有N條射線，首先根據(jù)測(cè)試精度將剖面分為M個(gè)網(wǎng)格，以射線理論為基礎(chǔ)的成像方法歸結(jié)為求解如下方程：

式中l(wèi)ij是第i條射線在第j個(gè)單元內(nèi)的路徑長度；Sj是第j個(gè)單元的慢度（1/V）值；ti是第i條射線的走時(shí)值。

用SIRT解上述方程步驟如下：

1）用直射線方法計(jì)算各單元平均慢度值S(0)。設(shè)某一單元內(nèi)有n條射線通過，lm是其中第m條通過單元的射線長度，且其射線總長為Lm，走時(shí)為Tm，則通過單元所用的時(shí)間分配tm=Tmlm/Lm，n條射線通過單元的總時(shí)間分配為，總長度分配為，則該單元的平均慢度S=tn/ln。

2）用聯(lián)合迭代（SIRT）算法校正各單元慢度值（S）。設(shè)某一單元內(nèi)共有n條射線通過，Tm′，Tm分別是通過該單元的第m條射線的計(jì)算走時(shí)和實(shí)測(cè)走時(shí)，lm（可用直射線或者彎曲射線追蹤算法求出）是射線通過該單元內(nèi)的長度。則分配給該單元的走時(shí)誤差 εm=（Tm-Tm′）lm/Lm，n 條射線通過單元內(nèi)總走時(shí)誤差為，總射線長度為，單元慢度S用下式校正：

計(jì)算時(shí)可根據(jù)測(cè)區(qū)性質(zhì)加一限制條件，使（1/Vmax）＜S（k+1）＜（1/Vmin）提高計(jì)算精度。

2 應(yīng)用實(shí)例

遼寧某基礎(chǔ)處理工程，測(cè)區(qū)基巖為大理巖及變質(zhì)巖，由于長期的地質(zhì)運(yùn)動(dòng)及人類采礦等活動(dòng)造成了基巖巖體破碎、裂隙發(fā)育、巖體強(qiáng)度降低等問題。為了提高巖體的完整性、加大巖體整體強(qiáng)度、增強(qiáng)巖體的穩(wěn)定性，需要對(duì)測(cè)區(qū)巖體進(jìn)行灌漿處理。為了解灌漿效果，在灌漿前后進(jìn)行了地震CT觀測(cè)。

2.1 工作方法與技術(shù)

根據(jù)測(cè)試要求及實(shí)際情況，本次彈性波CT成像測(cè)試工作采用在一個(gè)鉆孔中用電火花震源激震，另一鉆孔中用12個(gè)檢波器的檢波器串進(jìn)行接收的孔間透射方式，激發(fā)點(diǎn)距1 m、接收點(diǎn)距1 m。每一個(gè)激發(fā)及接收點(diǎn)均形成一個(gè)扇形區(qū)域。再根據(jù)不同路徑的彈性波走時(shí)應(yīng)用WYS2003軟件進(jìn)行反演計(jì)算并繪制孔間速度變化的色譜圖像。

使用的儀器為 SE2404EE型地震儀和電火花震源。

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

在外業(yè)實(shí)測(cè)的彈性波CT原始記錄上，讀取每一炮點(diǎn)的每一接收道的直達(dá)波初至?xí)r間。應(yīng)用“WYS2003彈性波CT層析成像軟件”進(jìn)行輸入處理，形成觀測(cè)系統(tǒng)圖、數(shù)據(jù)曲線圖、成像剖面圖。

數(shù)據(jù)處理時(shí)，根據(jù)實(shí)際情況選擇各種參數(shù)，一般采用直射線聯(lián)合迭代法，迭代次數(shù)一般小于50次，計(jì)算殘差小于3.0。

2.3 資料分析

1）灌漿前勘察資料分析

圖1 ZK9—ZK10孔彈性波CT成像剖面圖

圖1為灌漿前勘察獲得的ZK9-ZK10孔之間的CT剖面圖，從圖中可以看到，測(cè)量范圍內(nèi)巖體波速均勻性很差，巖體波速集中在2500～4500 m/s的區(qū)域，占72.6%，低于2500 m/s部分占16.9%、高于4500 m/s的部分占10.5%。分析認(rèn)為：灌漿前的巖體整體波速較低，完整性較差。

2）灌漿后測(cè)試資料分析

從灌漿后檢測(cè)鉆孔j8-j13的彈性波CT觀測(cè)系統(tǒng)圖（由于施工后場地變化j8，j13未在原勘察孔位置，但距ZK9，ZK10孔不超過6 m），可以看出，激發(fā)和接收點(diǎn)距均為1 m，根據(jù)彈性波CT分辨率應(yīng)滿足：設(shè)d為異常體尺寸,Δx為激發(fā)點(diǎn)距離和接收點(diǎn)距離,觀測(cè)系統(tǒng)激發(fā)點(diǎn)和接收點(diǎn)布置應(yīng)該盡可能滿足Δx≤d/3[1，5]。 次觀測(cè)系統(tǒng)安排合理，能夠滿足測(cè)試要求。

從j8～j13彈性波CT數(shù)據(jù)曲線中可以看到，在第24，31，37組數(shù)據(jù)中各有1個(gè)離散值，其他數(shù)據(jù)正常，剔除3個(gè)離散值后進(jìn)行反演計(jì)算，得到j(luò)8～j13孔彈性波CT成像剖面圖。

灌漿后的檢測(cè)鉆孔為J8～J13，孔間距16.0 m，激發(fā)和接受點(diǎn)距均為1 m。在獲得的觀測(cè)數(shù)據(jù)中，其第一組數(shù)據(jù)的擬合曲線，首波走時(shí)離散度較小，其他各組數(shù)據(jù)較均勻。

圖2 j13-j8孔彈性波CT成像剖面圖

圖2為J8-J13地震CT成像結(jié)果，由圖可見，兩孔間巖體波速較均勻，局部波速小于4650 m/s占 1.8%；大部在 4650～5500 m/s占 93.7%；部分巖體波速大于5500 m/s占4.3%。

由此可見，經(jīng)灌漿處理后，巖體的波速大大提高，巖體的整體性得到了增強(qiáng)，巖體更加均勻，其強(qiáng)度等物性指標(biāo)也得到了提高，表明灌漿工作取得了較好的效果。

3 注意的問題

1）在工作開始前，應(yīng)根據(jù)任務(wù)要求、現(xiàn)場實(shí)際情況等設(shè)計(jì)好觀測(cè)系統(tǒng)。在彈性波CT工作中觀測(cè)系統(tǒng)十分重要，因?yàn)椴煌纳渚€密度和射線角度直接影響著成像質(zhì)量和能夠發(fā)現(xiàn)的最小異常。

2）在數(shù)據(jù)的采集過程中，要特別注意：震源起震與儀器記錄起跳時(shí)間的一致性。

3）在用WYS2003軟件繼續(xù)數(shù)據(jù)反演時(shí)，各項(xiàng)參數(shù)應(yīng)根據(jù)實(shí)際情況選擇，初始單元列數(shù)不宜過大或太小，否則將引起失真，圖形應(yīng)用插值圓滑處理。

4 結(jié)語

應(yīng)用實(shí)例表明，彈性波CT技術(shù)具有設(shè)備比較簡單，施工方便，信息量大，對(duì)空間位置的反應(yīng)精細(xì)。在彈性波CT剖面圖上可以對(duì)巖體進(jìn)行比較全面的判斷，可以比較直觀的看到巖體內(nèi)各波速帶的位置、形態(tài)等特征，是一種能夠?qū)Ρ容^大的巖體進(jìn)行波速測(cè)試的技術(shù)，具有一定的先進(jìn)性和實(shí)用性。

[1]王運(yùn)生，顧漢明，王家映彈性波CT關(guān)鍵技術(shù)與應(yīng)用實(shí)例[J].工程勘察，2005，（3）.

[2]王運(yùn)生.最佳路徑算法在計(jì)算波路中的應(yīng)用[J].物探化探計(jì)算技術(shù),1992，（1）.

[3]董清華，朱介壽.井間電阻率層析成像及其應(yīng)用[J].計(jì)算物理，1999（5）.

[4]董清華.井間地球物理成像與工程勘察[J].工程勘察，1999，（4）.

[5]陳國金.井間觀測(cè)系統(tǒng)的討論[J].石油物探,1996（3）.