潘紀(jì)順, 高東攀, 冷元寶, 李傳金 , 呂芝全, 朱四新, 朱 瑞, 趙 詳

(1.華北水利水電大學(xué),河南 鄭州 450045; 2.黃河水利科學(xué)研究院,河南 鄭州 450003)

堤壩隱患的天然源面波成像試驗(yàn)研究及應(yīng)用

潘紀(jì)順1, 高東攀1, 冷元寶2, 李傳金1, 呂芝全1, 朱四新1, 朱 瑞1, 趙 詳1

(1.華北水利水電大學(xué),河南 鄭州 450045; 2.黃河水利科學(xué)研究院,河南 鄭州 450003)

諸如空洞、軟弱層、滲漏、裂縫等病害嚴(yán)重威脅堤壩的安全,因此,堤壩隱患的探測具有重要的現(xiàn)實(shí)意義。當(dāng)前應(yīng)用到堤壩隱患探測的物理方法主要有直流電法、電磁法、彈性波法、放射性法和流場法等,而這些方法都是有源的,或者向地下供電,或者發(fā)射電磁波,或者產(chǎn)生地震波。而天然源面波探測是無源的,它利用地球表面無時不刻存在的微弱震動,從里邊提取面波信息,通過對面波頻散曲線進(jìn)一步反演,可以獲得地下介質(zhì)速度結(jié)構(gòu)。通過黃河大堤隱患原型試驗(yàn)和長江大堤的實(shí)際探測,表明該技術(shù)對探測諸如軟弱層等隱患有一定的效果和實(shí)際應(yīng)用價值。

天然源;面波;堤壩隱患;黃河原型試驗(yàn);長江大堤

古今中外,歷朝歷代,水利建設(shè)都會得到足夠的重視。時至今日,中國大中型水庫已經(jīng)超過8萬余座,水利堤防工程長度已經(jīng)超過26萬km[1-2],長江三峽水電樞紐工程、黃河小浪底水利水電樞紐工程和南水北調(diào)等水利工程舉世矚目。自2011年起,中央財政每年向水利行業(yè)投放4 000個億元資金,且連續(xù)投放10年,由此可以預(yù)見,水利水電事業(yè)將處于一個快速發(fā)展時期,數(shù)量眾多的水利水電工程將應(yīng)運(yùn)而生,這必將對中國經(jīng)濟(jì)發(fā)展和社會安定產(chǎn)生巨大影響。然而,許多堤防工程久經(jīng)風(fēng)吹雨打、河流沖刷、年久失修。有的水利工程水患后雖經(jīng)反復(fù)修繕、加固,亦或由于存在施工質(zhì)量問題而導(dǎo)致存在裂縫、洞穴、松散、軟弱層、滲漏層等各類隱患,這些隱患在洪水期間極易形成滲水、管涌、漏洞、散浸、跌窩等險情,嚴(yán)重威脅堤壩的安全。因此,如何有效、快捷地探明這些堤壩的隱患,成為關(guān)乎水利堤壩安全的一個科學(xué)問題。

新中國建國之初,基本上沿襲清代河務(wù)機(jī)構(gòu)的“簽堤”方法探測堤壩隱患,即憑借“鐵簽”入土的快慢難易程度,根據(jù)人的直覺來判斷。“簽堤”方法探測深度淺,對規(guī)模小尺度的隱患容易遺漏。20世紀(jì)50年代末—80年代初,“錐探”方法在堤壩隱患探測中起主導(dǎo)作用,但是該方法依舊存在速度慢、效率低、易漏掉小尺寸隱患的缺點(diǎn)。然而,伴隨中國地球物理探測技術(shù)的發(fā)展,堤壩隱患探測技術(shù)得到了較大發(fā)展,并演化成“地球物理探測為主、鉆探為輔”的綜合探測體系。地球物理探測技術(shù)在堤壩隱患探測的應(yīng)用追索到20世紀(jì)50—60年代的山東大學(xué)和山東省河務(wù)局,兩者合作嘗試?yán)免?0放射性探測黃河的堤壩隱患。20世紀(jì)70年代,鞍山電子研究所研制了YB-1型暗縫探測儀,取得一定探測效果。80年代,山東省水利科學(xué)研究所成功研制了ZDT-1型堤壩隱患探測儀,并于次年獲得水利部科技進(jìn)步二等獎[3]。在90年代中國加大了資金的投入,地球物理探測技術(shù)在堤壩隱患探測領(lǐng)域呈現(xiàn)蓬勃發(fā)展的趨勢。1990年鄧習(xí)珠等人利用其研制的TTY-1型便攜式智能堤壩探測儀成功探測了堤壩中的蟻穴隱患[4]。1992年,“堤防隱患探測技術(shù)研究”列入國家“八五”重點(diǎn)科技項(xiàng)目。1993年葛建國等人將淺層反射地震的方法應(yīng)用到堤壩隱患探測中[5]。1993年劉康和等人將K剖面法成功應(yīng)用到堤壩隱患探測中[6]。1994年陳紹求等人利用雙頻激電儀探測堤壩隱患[7]。1998年吳相安等人對地質(zhì)雷達(dá)在堤壩隱患探測中的應(yīng)用效果進(jìn)行了研究[8]。1998年“三江洪水”爆發(fā),當(dāng)年水利部將“堤防隱患和險情探測儀器開發(fā)重大科技攻關(guān)項(xiàng)目”列入“988”科研計劃中。中國水利水電科學(xué)研究院研制了SDC-2型堤壩滲漏探測儀器,隨后,水利部長江水利委員會對地質(zhì)雷達(dá)的觀測系統(tǒng)進(jìn)行了革新,研制出了“雙頻多普勒相控陣探地雷達(dá)三維掃描成像系統(tǒng)”,水利部黃河水利委員會研制了ZDT-1型智能堤壩隱患探測儀,這些儀器在中國水利工程堤壩隱患探測中都有上乘的表現(xiàn)。2000年,國家自然基金資助了“洪水特性與減災(zāi)方法研究”項(xiàng)目。2000年何繼善院士發(fā)明了利用水流場流向和相對流速來判斷堤壩隱患存在與否的“流場法”,該方法能夠有效地探測到滲漏、管涌類的隱患[9]。2005年陳建生等利用溫度場方法來研究堤壩滲漏通道,取得了一定的效果[10]。所有這些大大促進(jìn)了中國堤壩隱患探測技術(shù)的發(fā)展,據(jù)統(tǒng)計,目前已經(jīng)應(yīng)用地球物理方法探測堤防工程2 000多千米,水庫100多座。

綜上所述,越來越多的地球物理探測技術(shù)被應(yīng)用到堤壩隱患的探測中,每種方法都有其優(yōu)缺點(diǎn)。但是,歸納起來發(fā)現(xiàn),當(dāng)前所有應(yīng)用到隱患探測的地球物理方法都需要場源,為主動源探測。本文試圖另辟蹊徑,作為一種嘗試,將天然源面波成像技術(shù)應(yīng)用到堤壩隱患探測中,試圖給出一種新的探測方法,給出了其在黃河隱患原型實(shí)驗(yàn)和長江實(shí)際堤壩的探測結(jié)果。

1 SPAC法天然源面波成像技術(shù)

天然源面波成像技術(shù)是從地面觀測到的隨機(jī)噪聲,利用相關(guān)技術(shù)從中提取面波信息,繼而對面波頻散曲線進(jìn)行反演得到地下介質(zhì)的S波速度結(jié)構(gòu)。該方法技術(shù)的理論是由著名地震學(xué)家Aki于1957年提出的[11],日本學(xué)者Okada H.大大發(fā)展了這項(xiàng)技術(shù)[12],中國許多學(xué)者也做了有益的探索[13-19]。當(dāng)前占主流的SPAC方法要求布設(shè)2個同心圓,每個圓周上均勻布設(shè)3臺觀測儀器,圓心上布設(shè)1臺觀測儀器,采用極坐標(biāo)系,微動信號是位置(r,θ)和時間t的函數(shù),在地球表面存在A(0,0)和B(r,θ)兩點(diǎn),A點(diǎn)為圓心,B點(diǎn)為圓周上一點(diǎn),其天然源觀測信號表述為:

(1)

(2)

A、B兩點(diǎn)的空間自相關(guān)函數(shù)可以定義為:

(3)

圓周上B點(diǎn)和中心點(diǎn)A點(diǎn)觀測信號的標(biāo)準(zhǔn)化自相關(guān)函數(shù)的方位平均值可表示為:

(4)

式中:x=2πfr/c(f),J0(x)是第1類0階貝賽爾函數(shù);c(f)或c(ω)為波的傳播速度;θ是波的入射角度;S(ω,r,θ)是圓心處的觀測信號和圓周上某點(diǎn)的觀測信號的交叉譜;S0(ω,0,θ)和Sr(ω,r,θ)是圓心處和圓周上某點(diǎn)的信號功率譜;ρ(ω,r)稱為空間自相關(guān)系數(shù)。

2 微動觀測儀器設(shè)備

本次試驗(yàn)的設(shè)備系統(tǒng)由中國地震局地震研究所(武漢)研制,共有8套微動儀,每臺儀器由外接電源、GPS時鐘、數(shù)據(jù)采集器和拾振器四部分組成,如圖1所示。

(2) 基于層次分析法的主觀權(quán)重的確定。根據(jù)表4的專家打分，最終得出各指標(biāo)的權(quán)重Wc={0.338，0.063，0.210，0.063，0.124，0.081，0.121}。

圖1 微動觀測儀器組成Fig.1 Composition of micro observation instruments

3 天然源面波成像技術(shù)在黃河隱患原型堤壩上的應(yīng)用

在距離黃河花園口景區(qū)以西約3.7 km處的崗李水庫建設(shè)兩道平行的實(shí)驗(yàn)堤進(jìn)行病害設(shè)置,主要服務(wù)于病害探測實(shí)驗(yàn)研究,監(jiān)測實(shí)驗(yàn)研究,堤防除險加固、潰決搶險實(shí)驗(yàn)研究。

病害的種類設(shè)置主要有:①松軟層2處,埋深分別是1 m和3 m,每處長10 m,寬2 m,厚0.5 m,用干沙、石子、秸稈等混合雜物填充。②孔洞,管件規(guī)格為外徑0.5 m,長2 m,壁厚2 cm以內(nèi),在不同位置埋置深度1、2、3、4、5 m。③裂縫,采用高2 m、厚0.05 m的2塊絕緣泡沫板模擬,其中與實(shí)驗(yàn)堤軸線方向垂直的稱為垂直裂縫,與實(shí)驗(yàn)堤軸線方向成30°、45°交角的稱為傾斜裂縫。

由于堤頂狹窄,無法采用常用的同心圓臺陣觀測臺陣,實(shí)際觀測采用由兩個等邊三角形共邊組合而成的菱形臺陣觀測臺陣,如圖2所示,圓點(diǎn)表示臺站,一共有7個臺站。圓點(diǎn)上方為臺站的編號,分別為S1,S2,…,S7。S5位于觀測系統(tǒng)的中心,代表測點(diǎn)的位置。7個臺站所處位置的幾何關(guān)系為:臺站S1,S4,S6組成一個等邊三角形,其中心處放置觀測臺站S3;臺站S2,S4,S6組成另一個等邊三角形,其中心處放置觀測臺站S7;臺站S4和S6連線的中間處放置臺站S5;臺站S1,S3,S5,S7,S2在一條直線上。由于各測點(diǎn)的場地不同,S5-S4的距離有兩種,分別為1 m和1.5 m。

圖2 菱形觀測臺陣Fig.2 Diamond observatory array

圖3是7個臺站同時記錄60 s得到的波形圖。圖4顯示了7 個臺站同時記錄到的微動信號的功率譜。

3.1 對軟弱層的探測試驗(yàn)

模型在樁號65～75 m段距壩頂3 m深度埋設(shè)厚度0.5 m的軟弱層,在85～95 m段距壩頂1 m深度埋設(shè)厚度0.5 m的軟弱層。觀測試驗(yàn)中心點(diǎn)的位置分別為:66、68、72、74、76、78、80、82、84、86、91、93、95、97、99,也就是說,測點(diǎn)66、68、72、74位于3 m埋深軟弱層的上方,測點(diǎn)86、91、93、95位于1 m埋深軟弱層的上方,86、91、93測點(diǎn)位于兩軟弱層的中間位置。從反演圖像(圖5)可以看出,1 m、3 m埋深的軟弱層都能夠得到較好的反映。同時從圖像反映出來,6 m處有一個界面,為新筑堤壩的底界面,其下11 m附近同樣有一個界面,這表明天然源面波成像技術(shù)在分層結(jié)構(gòu)的探測應(yīng)用中有較好的效果。

圖3 觀測時間為60 s的一段波形圖Fig.3 Oscillogram of 60 seconds observation time

圖4 7個臺站同時記錄到的微動信號的功率譜Fig.4 Power spectrum of micro signal by seven stations

圖5 軟弱層的探測試驗(yàn)反演圖像Fig.5 Inversion image of probe test in weak layer

3.2 對孔洞的探測試驗(yàn)

孔洞的探測試驗(yàn)主要針對3 m埋深和1 m埋深的2種類型試驗(yàn),其中,樁號12 m的測點(diǎn)正好位于3 m 埋深管子的正上方,針對3 m埋深管子,做了觀測點(diǎn)中心在9.5 m、11.5 m、12 m、12.5 m的觀測試驗(yàn),其圖形見圖6-(a)。盡管分辨不好,但是可以明顯看出3 m埋深管子形成的高橫波速度的團(tuán)狀反映。而對于埋深1 m的管子,其中心點(diǎn)位于觀測點(diǎn)5 m的正下方,做了觀測點(diǎn)中心在4.5 m、5 m、6.5 m的觀測試驗(yàn),結(jié)果見圖6-(b)。1 m埋深管子的形態(tài)有稍許反映,但是,已經(jīng)遠(yuǎn)沒有3 m埋深的探測效果好。

圖6 孔洞的探測試驗(yàn)圖像Fig.6 Probe test image of holes(a).埋深3 m;(b).埋深1 m。

對于垂直裂縫,其位于樁號60 m的正下方,同時在樁號59 m、60 m、60.5 m、61.5 m做了觀測,反演后構(gòu)制的成果圖見圖7,從圖上并未發(fā)現(xiàn)裂縫的存在,微動探測分辨不了裂縫隱患。

圖7 裂縫的探測試驗(yàn)反演圖像Fig.7 Inversion image of probe test of fissure

4 天然源面波成像技術(shù)在張家港長江堤壩上的應(yīng)用

張家港長江大堤朝東圩港段現(xiàn)有堤壩歷史上經(jīng)歷多次改造,在做防洪規(guī)劃設(shè)計之前需對當(dāng)前堤壩的情況做一個詳細(xì)地了解,勘察單位布設(shè)了鉆孔,應(yīng)用地質(zhì)雷達(dá)、高密度電法等物探手段進(jìn)行了探測,也進(jìn)行了天然面波探測試驗(yàn)。

實(shí)驗(yàn)在晚上進(jìn)行,天氣晴朗,微風(fēng)?；旧蠜]有干擾,采集的數(shù)據(jù)質(zhì)量也比較高。實(shí)驗(yàn)臺陣布設(shè)如圖8所示,臺陣半徑R=3.2 m、r=1.6 m,拾振器使用1檔,共進(jìn)行了8次觀測,每次觀測30 min左右。

圖8 朝東圩港堤段臺陣布設(shè)Fig.8 Array layout of Chaodong section

實(shí)驗(yàn)時數(shù)據(jù)采集的條件比較好,各個點(diǎn)都得到了結(jié)果。對每一個臺陣觀測到的頻散曲線進(jìn)行反演,得到了8條一維的剪切波速度結(jié)構(gòu),使用Kriging法進(jìn)行插值,得到了臺陣下方的剪切波速度剖面,如圖9所示。

由圖9可以看出,長江大堤朝東圩港段在6 m深度存在一明顯界面,分析為新舊堤的一個分界面,在淺部淺色部分顯示堤壩壩體淺表存在局部不均勻。

4 結(jié)論

通過原型堤壩試驗(yàn)和在張家港長江大堤的實(shí)際檢測,可以得到以下結(jié)論:

(1) 微動探測技術(shù)在堤壩隱患探測中可以作為一種新的手段來研究,其最大的優(yōu)點(diǎn)是無需提供物理源,并且探測對堤壩是無損的。

圖9 長江大堤探測成果圖Fig.9 Detection results figure of Changjiang grand dyke

(2) 原型堤壩隱患的微動探測試驗(yàn)表明,微動技術(shù)能夠較好地解決軟弱層的探測問題,對孔洞的探測也有一定的反映,但是對裂縫的探測卻沒有效果。

(3) 微動探測技術(shù)能夠較好地進(jìn)行堤壩的分層,尤其是對不同時期堤壩的區(qū)分有現(xiàn)實(shí)意義。

[1] 冷元寶,朱文仲,何劍,等.我國堤壩隱患及滲漏探測技術(shù)現(xiàn)狀及展望[J].水利水電科技進(jìn)展,2002,22(2):59-62.

[2] 冷元寶,任建東,王銳,等.堤防隱患與監(jiān)測技術(shù)展望[J].工程地球物理學(xué)報,2004,1(1):74-77.

[3] 山東省水利科學(xué)研究所.電法探測堤防隱患技術(shù)[J].人民黃河,1992(10):59.

[4] 鄧習(xí)珠,肖明,曹思文.堤壩隱患電法探測技術(shù)的應(yīng)用[J].江西水利科技,1990(2):35-38.

[5] 葛建國,陶建基.淺層地震反射波探測法在堤壩上的應(yīng)用[J].水利工程管理技術(shù),1993(1):30-32.

[6] 劉康和.應(yīng)用K剖面法探測堤壩隱患[J].地質(zhì)與勘探,1993,29(2):40-43.

[7] 陳紹求.雙頻激電法在堤壩隱患探測中的應(yīng)用[J].中南礦業(yè)學(xué)院學(xué)報,1994,25(增刊):25-28.

[8] 吳相安,徐興新,吳晉,等.水利隱患GPR探測方法研究[J].地質(zhì)與勘探,1998,34(3):47-52.

[9] 何繼善.堤防滲漏管涌“流場法”探測技術(shù)[J].銅業(yè)工程,2000(1):5-8.

[10] 陳建生,董海洲,吳慶林.虛擬熱源法研究壩基裂隙巖體滲漏通道[J].巖石力學(xué)與工程學(xué)報,2005,24(22):4019-4024.

[11] Aki K.Space and time spectra of stationary stochastic waves,with special reference to microtremors[J].Bull.Earthq.Res.Inst.,1957,35:415-456.

[12] Okada H.Theory of efficient array observations of microtremors with special reference to the SPAC method[J].Explor.Geophys.,2006,37:73-85.

[13] 王振東.微動的空間自相關(guān)法及其實(shí)用技術(shù)[J].物探與化探,1986,10(2):123-133.

[14] 馮少孔.微動勘探技術(shù)及其在土木工程中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報,2003,22(6):1029-1036.

[15] 陶夏新.工程場地條件評定中的地脈動研究[J].地震工程與工程振動,2001,21(4):18-23.

[16] 葉太蘭.微動臺陣探測技術(shù)及其應(yīng)用研究[J].中國地震,2004,20(1):47-52.

[17] 廖成旺,鄧濤,丁煒,等.微動觀測臺站定位方式研究[J].大地測量與地球動力學(xué),2007,27(6):61-64.

[18] 何正勤,丁志峰,賈輝,等.用微動中的面波信息探測地殼淺部的速度結(jié)構(gòu)[J].地球物理學(xué)報,2007,50(2):492-498.

[19] 徐佩芬,李傳金,凌甦群,等.利用微動勘察方法探測煤礦陷落柱[J].地球物理學(xué)報,2009,52(7):1923-1930.

(責(zé)任編輯：于繼紅)

Natural Source Surface Wave Imaging on the Hidden Trouble of Dikes and Dams

PAN Jishun1, GAO Dongpan1, LENG Yuanbao2, LI Zhuanjin1, LV Zhiquan1, ZHU Sixin1, ZHU Rui1, ZHAO Xiang1

(1.NorthChinaUniversityofWaterConservancyandHydroelectricPower,Zhengzhou,Henan450045; 2.YellowRiverInstituteofHydraulicResearch,Zhengzhou,Henan450003)

The troubles such as caves,weak layer,leakage and cracks threaten the safety of dikes and dams. So it is important to detect the troubles. Nowadays the geophysical method used in detection of the troubles including DC electrical method,electromagnetic method,elastic wave method,radioactive method and flow fluid method. All these methods are active,some need supply electric current,some need eject electromagnetic wave,the other need create seismic wave. However passive surface is passive. The passive surface method gets the surface wave from the microtremor which exists on the earth surface all the time. The subsurface velocity structure will be obtained by inverting the wave dispersive curve. Both Yellow River dam trouble prototype test and the practical detections in the Yangtze River dam prove the passive surface method is effective to detect the weak layer trouble and it has some practical application value.

natural source; surface wave; dikes and dams troubles,Yellow River prototype test; the Yangtze River dam

2016-05-16;改回日期:2016-05-30

國家自然科學(xué)基金(41174052);河南省基礎(chǔ)與前沿技術(shù)研究計劃項(xiàng)目(112300410252)。

潘紀(jì)順(1968-),男,教授,博士,地球物理專業(yè),從事工程地球物理勘探、地震數(shù)據(jù)處理及地球物理反演等的教學(xué)與研究。E-mail:jspan123@126.com

TV871

A

1671-1211(2016)03-0306-05

10.16536/j.cnki.issn.1671-1211.2016.03.014

數(shù)字出版網(wǎng)址:http://www.cnki.net/kcms/detail/42.1736.X.20160530.0937.018.html 數(shù)字出版日期:2016-05-30 09:37