胡振華，李云鵬，夏潤(rùn)賢，何震洲，朱虹

（浙江省錢(qián)塘江流域中心，浙江 杭州 310016）

1 問(wèn)題的提出

錢(qián)塘江是浙江人民的母親河，錢(qián)塘江海塘則是沿江地區(qū)防洪御潮的重要屏障。位于海寧市、海鹽縣以及杭州市蕭山區(qū)境內(nèi)的古海塘有400多年歷史，已列入國(guó)家級(jí)或省級(jí)保護(hù)文物。海塘多年來(lái)受到外界作用（涌潮、強(qiáng)降水）易存在發(fā)育隱患，這些隱患往往具有隱蔽性、漸變性，需要定期對(duì)海塘進(jìn)行常態(tài)化監(jiān)測(cè)以維持其長(zhǎng)期運(yùn)行。傳統(tǒng)人工監(jiān)測(cè)多采用巡查、觀測(cè)等方式，能夠?qū)扔谢蛘谘葑兊碾[患起到監(jiān)控作用，但很難滿(mǎn)足大規(guī)模、高效率、高精度以及無(wú)損養(yǎng)護(hù)監(jiān)測(cè)的需求。地質(zhì)雷達(dá)方法高效、無(wú)損等優(yōu)點(diǎn)恰恰彌補(bǔ)了人工監(jiān)測(cè)方法的不足，尤其是由二維單通道地質(zhì)雷達(dá)發(fā)展而來(lái)的三維車(chē)載地質(zhì)雷達(dá)的應(yīng)用，使常態(tài)化監(jiān)測(cè)成為可能，有利于保障海塘安全。

地質(zhì)雷達(dá)儀器設(shè)備伴隨著電磁理論的發(fā)展而出現(xiàn)，近三十年來(lái)發(fā)展迅速，由于應(yīng)用場(chǎng)景不同，產(chǎn)生多天線組合的陣列式及步進(jìn)變頻三維地質(zhì)雷達(dá)。美國(guó)GSSI公司1990年設(shè)計(jì)出世界上第一臺(tái)四通道數(shù)字化地質(zhì)雷達(dá)，這是陣列式三維地質(zhì)雷達(dá)的雛形；隨后，1994年最先推出高速公路測(cè)量雷達(dá)SIR-10系列，并在我國(guó)首先使用。2001年，挪威3D-Radar公司推出GeoScope步進(jìn)變頻式三維地質(zhì)雷達(dá)，通過(guò)發(fā)射步進(jìn)頻率的電磁波兼顧分辨率和勘探深度。2016年，中公高科養(yǎng)護(hù)科技股份有限公司和大連中?？萍及l(fā)展有限公司聯(lián)合開(kāi)發(fā)52通道三維陣列地質(zhì)雷達(dá)，我國(guó)地質(zhì)雷達(dá)實(shí)現(xiàn)由二維單個(gè)天線到三維陣列式天線的跨越。

理論研究方面，李世念等[1]利用基于三維時(shí)域有限差分法原理的模擬軟件GprMax 3D，開(kāi)展道路路基空洞模擬和電磁波三維探地正演模擬，對(duì)獲取的空洞模型三維數(shù)據(jù)體進(jìn)行分析，研究2種不同充填條件下空洞模型的三維地質(zhì)雷達(dá)圖像特征，有助于三維地質(zhì)雷達(dá)實(shí)測(cè)圖像解譯，提高判定監(jiān)測(cè)結(jié)果的可靠性，為開(kāi)展三維地質(zhì)雷達(dá)在道路空洞探測(cè)中的研究積累經(jīng)驗(yàn)。在實(shí)際應(yīng)用領(lǐng)域，唐嘉明[2]于2020年基于三維地質(zhì)雷達(dá)對(duì)瀝青路面施工質(zhì)量評(píng)價(jià)與控制進(jìn)行數(shù)值模擬及實(shí)測(cè)數(shù)據(jù)分析研究。王洪華等[3]為全面了解上海某道路下方塌陷空洞的三維空間分布，采用意大利 IDS 公司生產(chǎn)的RIS-K2型地質(zhì)雷達(dá)系統(tǒng)實(shí)現(xiàn)道路塌陷隱患三維探測(cè)的可視化和定量化，使塌陷隱患的反映更直觀、解譯更準(zhǔn)確，大大克服了二維雷達(dá)探測(cè)的局限性。趙鐠[4]利用三維地質(zhì)雷達(dá)對(duì)北京某實(shí)訓(xùn)場(chǎng)實(shí)施道路結(jié)構(gòu)及管線探測(cè)，取得理想的應(yīng)用效果。

目前，相關(guān)儀器設(shè)備研發(fā)及理論研究已具備一定基礎(chǔ)，但國(guó)內(nèi)研究機(jī)構(gòu)尚未在海塘堤壩監(jiān)測(cè)中開(kāi)展相關(guān)工作。本文采用數(shù)值模擬和現(xiàn)場(chǎng)實(shí)測(cè)相結(jié)合的方法對(duì)錢(qián)塘江海塘進(jìn)行相關(guān)研究，研究成果有利于三維地質(zhì)雷達(dá)技術(shù)的推廣應(yīng)用。

2 方法原理

二維地質(zhì)雷達(dá)利用高頻電磁波信號(hào)對(duì)地下目標(biāo)體進(jìn)行掃描獲得探測(cè)影像，電磁波信號(hào)遇到具有波陰抗差異的界面會(huì)形成反射，且界面間波陰抗差異越大，反射信號(hào)越強(qiáng)。與發(fā)射天線同步移動(dòng)的接收天線接收到反射的回波信號(hào)，通過(guò)數(shù)字信號(hào)處理，形成探測(cè)剖面圖。對(duì)雷達(dá)圖像進(jìn)行分析、判斷，可推斷出地下介質(zhì)分布特征。工作原理見(jiàn)圖1。

圖1 二維單通道地質(zhì)雷達(dá)工作原理圖

三維地質(zhì)雷達(dá)采用陣列天線激發(fā)、接收技術(shù)，通過(guò)對(duì)某區(qū)域探測(cè)，能夠形成高密度三維立體電磁波數(shù)據(jù)。該技術(shù)在二維地質(zhì)雷達(dá)的基礎(chǔ)上發(fā)展起來(lái)，通過(guò)發(fā)射天線和接收天線進(jìn)行交叉發(fā)射-接收，實(shí)現(xiàn)一發(fā)雙收的目的（見(jiàn)圖2）。組成的天線陣列可同時(shí)采集多個(gè)通道的數(shù)據(jù)，達(dá)到三維勘探效果，屬于三維探測(cè)技術(shù)。

圖2 三維地質(zhì)雷達(dá)的天線陣列工作原理示意圖

與二維地質(zhì)雷達(dá)相比，三維地質(zhì)雷達(dá)有以下優(yōu)勢(shì)：

（1）采用三維陣列天線技術(shù)，可采集到高密度、無(wú)縫拼接的海量數(shù)據(jù)，不會(huì)造成地下信息缺失；

（2）高密度采集獲得的縱、橫向數(shù)據(jù)間距接近天線中心波長(zhǎng)的1/4，滿(mǎn)足高分辨率要求；

（3）通過(guò)三維陣列式天線實(shí)現(xiàn)真三維采集，為將來(lái)全三維處理、解釋提供硬件實(shí)現(xiàn)的可能；

（4）三維陣列式天線發(fā)射脈沖頻率一般可選擇200 MHz~1 300 MHz，高頻天線可做到超淺層的高分辨，低頻天線可保證一定的探測(cè)深度；

（5）通過(guò)帶基站的動(dòng)態(tài)GNSS對(duì)天線陣列進(jìn)行高精度定位，保證雷達(dá)數(shù)據(jù)的精確歸位，直觀解釋圖像。

3 海塘隱患模型的正演計(jì)算

正演計(jì)算是確定方法有效性的常用手段。常見(jiàn)海塘隱患有層間脫空、局部發(fā)育空洞等，該類(lèi)病害多位于海塘表層以下，故一般無(wú)大量積水分布，具體表現(xiàn)為上、下兩層結(jié)構(gòu)之間存在由空氣充填的空隙。本次正演計(jì)算以常見(jiàn)的海塘結(jié)構(gòu)層間脫空為對(duì)象，建立簡(jiǎn)化地電結(jié)構(gòu)模型（見(jiàn)圖3），設(shè)置網(wǎng)絡(luò)空間大小3.00 m×1.50 m，網(wǎng)格間距0.03 m×0.03 m，電磁波中心頻率200 MHz，各物性參數(shù)按表1賦值，采用GPRMax3.0正演模擬軟件進(jìn)行正演計(jì)算，結(jié)果見(jiàn)圖4。

表1 模型介質(zhì)物性參數(shù)表

圖3 海塘結(jié)構(gòu)斷面圖與正演模型圖 單位：cm

圖4 正演計(jì)算結(jié)果圖 單位：cm

由圖4可見(jiàn)，脫空模型位置在地質(zhì)雷達(dá)剖面上呈現(xiàn)弧形繞射異常。越接近模型中心異常強(qiáng)度越大，趨向兩側(cè)則異常漸弱，且形態(tài)趨于寬緩。計(jì)算結(jié)果表明，三維地質(zhì)雷達(dá)用于海塘隱患探測(cè)具有可行性。

4 異常特征及驗(yàn)證

基于正演計(jì)算結(jié)果，在錢(qián)塘江海塘里程91K+400 m段，針對(duì)漿砌條石下方0.20~1.00 m深度范圍的脫空發(fā)育情況開(kāi)展試驗(yàn)。

根據(jù)設(shè)計(jì)資料，古海塘試驗(yàn)段塘面下方從上至下分別為漿砌條石、碎石回填土、水泥摻粉煤灰灌漿。其中條石厚約0.25 m；碎石回填層厚約2.00 m，呈倒梯形，自西向東厚度不斷增大，碎石由反濾層包裹，反濾層西側(cè)為回填粉土；水泥摻粉煤灰灌漿厚約3.50 m。

車(chē)載三維地質(zhì)雷達(dá)沿海塘道路進(jìn)行覆蓋探測(cè)，測(cè)線沿道路方向布置，探測(cè)天線頻率200 MHz，點(diǎn)距4 cm，測(cè)線線距4 cm，剖面見(jiàn)圖5。

圖5 三維地質(zhì)雷達(dá)剖面圖

由探測(cè)結(jié)果可知：在三維地質(zhì)深度切片圖中，海塘隱患出現(xiàn)不規(guī)則強(qiáng)振幅反射異常，相位局部存在倒轉(zhuǎn)；在正交的2個(gè)二維時(shí)間剖面中，異常區(qū)振幅增強(qiáng)，有多次波發(fā)育，與周邊存在明顯差異。根據(jù)以上異常特征，圈定該處為疑似脫空異常隱患位置，采用800 MHz天線布置5條二維地質(zhì)雷達(dá)對(duì)異常點(diǎn)開(kāi)展復(fù)測(cè)，工作布置見(jiàn)圖6。

圖6 工作布置圖

以L5、L7測(cè)線為例（見(jiàn)圖7~8），兩線分別相交于2.7/L5、14.8/L7位置。L5測(cè)線剖面中同相軸零星錯(cuò)斷，0.5~9.5 m處（虛線框區(qū)域）同相軸局部上翹；L7測(cè)線處相位局部倒轉(zhuǎn)，9.0~20.0 m處（虛線框區(qū)域）同相軸也存在上翹情況。2條正交測(cè)線的異常特征表明：該處異常呈現(xiàn)板狀分布特征，中部存在脫空現(xiàn)象，導(dǎo)致電磁波傳播速度加快，同相軸形成弧形上翹的特征。

圖7 L5測(cè)線剖面圖

圖8 L7測(cè)線剖面圖

為進(jìn)一步確認(rèn)隱患情況，分別采用釬探和開(kāi)挖的形式進(jìn)行驗(yàn)證（見(jiàn)圖9）。釬探結(jié)果顯示脫空高約25 cm，隨后圍繞釬探點(diǎn)開(kāi)挖，以進(jìn)一步確認(rèn)病害異常。根據(jù)開(kāi)挖情況，脫空隱患由深部基礎(chǔ)物質(zhì)流失造成。后續(xù)對(duì)該隱患點(diǎn)進(jìn)行回填恢復(fù)處理。

圖9 開(kāi)挖驗(yàn)證情況圖

5 結(jié) 語(yǔ)

地質(zhì)雷達(dá)法是監(jiān)測(cè)海塘隱患常用的方法，三維地質(zhì)雷達(dá)法由于實(shí)現(xiàn)了三維勘探，能夠更加高精度、高保真、高速度地對(duì)海塘隱患進(jìn)行探測(cè)定位。本文應(yīng)用三維地質(zhì)雷達(dá)探測(cè)海塘堤壩隱患，并通過(guò)正演模擬和實(shí)際探測(cè)論證該方法的異常特征及有效性。

結(jié)合實(shí)際探測(cè)情況，由于三維地質(zhì)雷達(dá)多采用車(chē)載方式進(jìn)行作業(yè)，發(fā)射-接收天線存在空氣耦合問(wèn)題，影響勘探深度，所以常用車(chē)載三維地質(zhì)雷達(dá)對(duì)海塘隱患進(jìn)行普查，發(fā)現(xiàn)定位異常后采用二維地質(zhì)雷達(dá)或其他高精度物理探測(cè)方法對(duì)異常情況進(jìn)行詳細(xì)解剖，這樣的作業(yè)模式更能提高探測(cè)結(jié)果的準(zhǔn)確度和探測(cè)效率。