地震層析成像研究及在巖溶發(fā)育區(qū)應(yīng)用(全文)

精品文檔-下載后可編輯地震層析成像研究及在巖溶發(fā)育區(qū)應(yīng)用(全文)首先通過扇形觀測(cè)系統(tǒng)獲取鉆孔或平硐的首波走時(shí)數(shù)據(jù)(ti),然后通過求解大型矩陣方程來獲取兩孔之間速度剖面圖像,根據(jù)速度剖面圖像可以直觀準(zhǔn)確地判斷異常體大小分布。層析成像中的反演方法可分為線性方法和非線性方法兩種,非線性方法主要有:遺傳算法、模擬退火法和神經(jīng)網(wǎng)絡(luò)法等。在體波層析成像中,使用線性反演方法的較多,如奇異值分解法(SVD)、共軛梯度法(CG)和最小二乘法(LSQR)等(見圖2)。典型的聯(lián)合迭代圖像重建法(SIRT)不利于同一迭代的第j條射線的修改值來計(jì)算第j+1條射線的修改值,而是把第k輪迭代中由所有射線得到的修改值保存起來,在本輪對(duì)射線迭代結(jié)束時(shí)求某種平均$f,然后由fk+1=fk+$f對(duì)每個(gè)象元的慢度作修改,并留作下一輪迭代使用。最小二乘法(LSQR)求解的是ATAx=ATb,具有計(jì)算量小的優(yōu)點(diǎn),并且能很容易地利用矩陣的稀疏性簡(jiǎn)化計(jì)算,因而適合求解大型稀疏問題。

數(shù)值模擬

本文為了比較SIRT算法和LSQR算法求解大型稀疏矩陣方程的能力,建立如圖3所示模型進(jìn)行速度層析成像反演。左邊設(shè)置4個(gè)激發(fā)點(diǎn),右邊12個(gè)檢波點(diǎn)接收。模型背景速度為3000m/s,異常區(qū)域速度分別為5000m/s和2000m/s,其反演結(jié)果見圖4。從圖4中分析得出,在SIRT算法中,獲得的數(shù)據(jù)量較少,而未知像素較多,得到的系數(shù)矩陣通常是大型的稀疏矩陣,因而方程組都是欠定的,甚至是嚴(yán)重欠定的。所以采用SIRT的反演方法,波速異常區(qū)域雖然可以大致反映出來,但與原始模型有一定誤差,如右圖中較小的速度變大,較大的速度變小,此誤差源于SIRT方法的局限性;LSQR的反演方法,具有數(shù)值穩(wěn)定、能充分利用矩陣稀疏性減少計(jì)算量等特點(diǎn),所以采用LSQR法取得了較好的結(jié)果,波速異常區(qū)域在圖4(a)中明顯反映出來,與原始模型吻合,且波速誤差較小。通過兩種算法模型的對(duì)比分析得出,LSQR算法在精度、效率和穩(wěn)定性三個(gè)方面都能達(dá)到較好的平衡,運(yùn)用LSQR算法進(jìn)行層析處理能獲得較為滿意的圖像。

工程實(shí)例

1.工區(qū)概況及勘探目的

壩址區(qū)地層主要出露因民組(Pyt2)和落雪組(Plt2)兩組地層,為一套淺變質(zhì)的碳酸鹽巖。壩址發(fā)育的地表巖溶類型有:溶縫、溶溝、溶槽、溶洞、溶蝕洼地。本次地震CT探測(cè)主要目的是探查壩址區(qū)右岸K25巖溶區(qū)發(fā)育分布情況。巖體中彈性波的傳播特征主要取決于巖體的物理力學(xué)性質(zhì),除本身強(qiáng)度外,在很大程度上受巖體結(jié)構(gòu)面性質(zhì)及其組合關(guān)系的控制。受風(fēng)化卸荷影響的巖體,其物理力學(xué)性質(zhì)及強(qiáng)度與新鮮基巖存在較大差異,且風(fēng)化與卸荷程度不同的巖體結(jié)構(gòu)和強(qiáng)度也會(huì)不同,這些差異和變化都可以通過CT成像中的波速分布情況來判斷。該水電站平硐、鉆孔聲波測(cè)試成果表明,堅(jiān)硬完整的碳酸鹽巖聲波速度較高,一般在5000m/s左右,受溶洞的影響,則呈現(xiàn)低波速,一般在3500m/s以下,這種明顯的波速差異為平硐地震CT法的應(yīng)用提供了地球物理前提。

2.工作方法及成果分析

本次地震CT采用一硐發(fā)射、另一硐接收的觀測(cè)方式。具體步驟是觀測(cè)射線網(wǎng)絡(luò);根據(jù)射線的疏密程度及成像精度劃分規(guī)則的成像單元;運(yùn)用射線追蹤理論橢圓約束法;采用的最小二乘法形成被測(cè)區(qū)域的波速圖像;根據(jù)圖像中的波速分布情況來重建探測(cè)區(qū)域內(nèi)巖體空間分布。信號(hào)采集器為R24地震儀,以及與之配套的高靈敏地震檢波器。激發(fā)方式采用炸藥震源,采樣間隔100us,道間距為310m,主頻30kHz檢波器接收做全硐觀測(cè)。圖5為右岸850m高程地震CT成果圖,測(cè)試區(qū)域跨PD56、PD58和LPD24平硐,由剖面的波速圖像可看到,完整性較好的巖體波速以4000~5000m/s為主,而位于PD58平硐硐深22m左右的上游硐壁附近,存在波速為2000~2250m/s低速異常區(qū)(k25)。后通過對(duì)此低速異常區(qū)進(jìn)行開挖驗(yàn)證,近似橢圓型巖溶區(qū),其內(nèi)溶蝕角礫發(fā)育、有泥質(zhì)填充,與周圍巖體速度差異大,異常明顯。利用硐間層析成像技術(shù),直觀劃分了巖體速度分布情況,準(zhǔn)確的探明了巖溶發(fā)育區(qū)的空間分布及走向情況。

3.與平硐波速測(cè)試對(duì)比

圖6為PD58平硐內(nèi)進(jìn)行的地震波速、單孔聲波及跨孔聲波的巖體波速測(cè)試。從圖中可以看出,三種方法得到的巖體波速大致呈相同的趨勢(shì),位于PD58平硐硐深20m左右的位置都出現(xiàn)明顯速度降低,與地震CT探測(cè)成果相吻合。

結(jié)語

本文利用硐間層析成像技術(shù),直觀準(zhǔn)確地判斷了巖溶發(fā)育區(qū)某水電站右岸壩址區(qū)巖體質(zhì)量,探明了喀斯特洞穴的空間分布及走向情況,為水電站壩址設(shè)計(jì)提供了必要的地球物理參數(shù)。層析成像成果與該地區(qū)單孔聲波及跨孔聲波所測(cè)的巖體波速相吻合,表明了地震層析成像技術(shù)是一種高分辨率探測(cè)巖溶發(fā)育區(qū)的有效手段。(1)CT圖像重建與許多因素有關(guān),在反演方法上,進(jìn)行了模型試算,得出最小二乘正交分解法(LSQR)的效果較聯(lián)合迭代重建法(SIRT)有一定的提高,LSQR算法在精度、效率和穩(wěn)定性三個(gè)方面都能達(dá)到較好的平衡。(2)在巖溶地區(qū)的工程勘察中,地震層析成像技術(shù)在劃分巖體速度分布情況及探查巖溶或者巖體破碎帶等方面具有良好的效果,并與地震波CT和單孔聲波測(cè)試、跨孔聲波測(cè)試等物探方法