三維微動(dòng)近地表速度結(jié)構(gòu)成像方法在城市地下溶洞探測(cè)中的應(yīng)用

曾 求, 劉 磊, 李成香, 胡云飛, 王向南

(1.資源與生態(tài)環(huán)境地質(zhì)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430034; 2.湖北省地質(zhì)局 地球物理勘探大隊(duì),湖北 武漢 430056;3.中機(jī)三勘巖土工程有限公司,湖北 武漢 430013)

隨著城市建設(shè)的高速發(fā)展,地質(zhì)災(zāi)害成為一個(gè)潛在的安全隱患,例如地下溶洞可能導(dǎo)致地面塌陷、地震和地面沉降等災(zāi)害,對(duì)城市基礎(chǔ)設(shè)施和人員安全構(gòu)成威脅[1-2]。在這種情況下,地震學(xué)三維/四維成像技術(shù)能夠發(fā)揮重要作用,它利用地震背景噪聲互相關(guān)函數(shù)進(jìn)行速度結(jié)構(gòu)成像,能夠準(zhǔn)確地刻畫(huà)地下溶洞的大小和位置,提前發(fā)現(xiàn)潛在的地質(zhì)災(zāi)害風(fēng)險(xiǎn)[3-4],提高城市基礎(chǔ)設(shè)施的安全性和可靠性。

中國(guó)學(xué)者將地震學(xué)方法應(yīng)用在城市中的歷史可以追溯到20世紀(jì)六七十年代[5-6]。此期間,開(kāi)始運(yùn)用地震勘探技術(shù)和地震波形分析方法,研究城市地下結(jié)構(gòu)并開(kāi)展地質(zhì)災(zāi)害風(fēng)險(xiǎn)評(píng)估。隨著科技的進(jìn)步和人類(lèi)需求的增加,這一領(lǐng)域的研究逐漸深入和擴(kuò)展,目前發(fā)展比較成熟的城市淺層地震勘探方法包括主動(dòng)源地震勘探和被動(dòng)源地震勘探方法。其中主動(dòng)源地震勘探方法是一種利用人工地震波進(jìn)行勘探的技術(shù),該方法通過(guò)在地表或地下埋設(shè)震源,人工產(chǎn)生地震波,并利用接收器記錄地震波在地下介質(zhì)中的反射和折射特性來(lái)推斷地下結(jié)構(gòu),進(jìn)而反演重建地下介質(zhì)分布情況,獲取城市地下速度結(jié)構(gòu)、巖層界面和地下構(gòu)造。中國(guó)利用地震學(xué)方法進(jìn)行城市三維地下空間的應(yīng)用取得了顯著成果。例如,北京市應(yīng)用地震臺(tái)陣和記錄儀等設(shè)備,成功獲取了地下速度結(jié)構(gòu)信息,并進(jìn)行了地震災(zāi)害風(fēng)險(xiǎn)評(píng)估[7],同時(shí)廣州、深圳等城市也開(kāi)展了類(lèi)似的研究[8-9]。然而通過(guò)主動(dòng)源地震勘探方法獲取高分辨率的淺層速度結(jié)構(gòu)成本高昂,且存在破壞現(xiàn)場(chǎng)結(jié)構(gòu)的風(fēng)險(xiǎn),因此近年來(lái)被動(dòng)源成像方法開(kāi)始在城市淺層三維速度結(jié)構(gòu)探測(cè)中得到越來(lái)越多的應(yīng)用。

地震背景噪聲成像技術(shù)是一種利用地震背景噪聲信號(hào)進(jìn)行地下結(jié)構(gòu)成像的方法[10]。它基于地球自然震動(dòng)和人類(lèi)活動(dòng),如海浪、風(fēng)、交通、建筑物振動(dòng)等背景噪聲,通過(guò)分析和處理噪聲信號(hào)的波形和頻譜特征,獲取地下結(jié)構(gòu)的相關(guān)信息[11-14]。因這種方法無(wú)需人工產(chǎn)生地震波,對(duì)環(huán)境的影響較小,成本較低,適用于城市地下結(jié)構(gòu)探測(cè)等,得到了廣泛的應(yīng)用。例如,Roux[15]利用地震背景噪聲成像技術(shù)在舊金山市成功探測(cè)到了地下斷層的位置和走向,為地震風(fēng)險(xiǎn)評(píng)估和城市規(guī)劃提供了重要依據(jù);Li et al.[16]利用地震背景噪聲成像技術(shù)在新西蘭某區(qū)域成功探測(cè)到了地下自來(lái)水管線(xiàn)的位置、深度和走向等信息,為城市地下管線(xiàn)的維護(hù)和管理提供了支持。在城市地下溶洞探測(cè)方面,近年來(lái)也有一些相關(guān)研究,例如Chen et al.[17]利用地震背景噪聲數(shù)據(jù)獲取了武漢某地區(qū)地下二維速度結(jié)構(gòu)信息,其結(jié)果與鉆孔數(shù)據(jù)有較好的一致性。

地下溶洞是城市地鐵建設(shè)中的一個(gè)重要問(wèn)題,對(duì)工程的安全運(yùn)行和建設(shè)進(jìn)程都會(huì)產(chǎn)生不利影響。注漿處理是解決這一問(wèn)題常用的方法,不僅可以加固地下溶洞的結(jié)構(gòu),提高地下結(jié)構(gòu)的穩(wěn)定性和承載能力,防止水災(zāi),保護(hù)基礎(chǔ)設(shè)施,也可以提高地下空間的利用率。但在注漿處理前后,需要準(zhǔn)確地確定地下溶洞的大小和位置,獲取溶洞區(qū)域淺層三維精細(xì)速度結(jié)構(gòu)。

目前利用三維微動(dòng)近地表速度結(jié)構(gòu)成像方法(即地震背景噪聲成像方法)獲取城市地下溶洞的大小和位置的研究較少。本研究首先收集地震背景噪聲數(shù)據(jù),并進(jìn)行預(yù)處理;然后計(jì)算地震背景噪聲互相關(guān)函數(shù),并通過(guò)分析互相關(guān)函數(shù)得到速度結(jié)構(gòu)信息;最后利用三維速度結(jié)構(gòu)建模軟件(Paraview)獲取研究區(qū)域近地表三維速度結(jié)構(gòu)并進(jìn)行可視化,通過(guò)與鉆孔數(shù)據(jù)進(jìn)行對(duì)比分析,發(fā)現(xiàn)研究區(qū)地下溶洞的分布特征,可為城市道路建設(shè)中的溶洞注漿提供數(shù)據(jù)指導(dǎo)。

1 數(shù)據(jù)與方法

本研究選取武漢某地鐵線(xiàn)路進(jìn)行實(shí)驗(yàn)(圖1)。研究區(qū)長(zhǎng)約80 m、寬約30 m,鉆孔資料顯示該區(qū)域地層由上至下主要分為4層:人工填土層(～0～2 m)、黏土層(～6～10 m)、黏土夾碎石層(～10～15 m)、灰?guī)r層(～9～15 m),往下為溶洞區(qū)域,溶洞空腔高度未知。該研究區(qū)在前期溶洞探測(cè)中主要采用的方式是鉆孔。圖1中黃色區(qū)域顯示了前期鉆孔施工的主要區(qū)域與推測(cè)的溶洞分布位置,可以看出鉆孔主要集中分布在2線(xiàn)與3線(xiàn)。這是由于1線(xiàn)與4線(xiàn)為居民住宅和商鋪,無(wú)法進(jìn)行鉆孔施工。前期在推測(cè)的溶洞區(qū)域Y1、Y2、Y3區(qū)進(jìn)行過(guò)注漿處理,然而在注漿量達(dá)到推測(cè)溶洞區(qū)域體積上限后仍不能注滿(mǎn)該區(qū)域,推測(cè)研究區(qū)的溶洞可能延伸至1線(xiàn)與4線(xiàn),因居民住宅和商鋪對(duì)鉆孔施工的限制,注漿工程不得不停止。

圖1 研究區(qū)概況與臺(tái)站分布圖

為了獲取研究區(qū)淺層三維速度結(jié)構(gòu),筆者團(tuán)隊(duì)擬采用不對(duì)地面造成破壞的地震背景噪聲成像技術(shù)方法進(jìn)行實(shí)驗(yàn)。在研究區(qū)布設(shè)了131個(gè)地震流動(dòng)臺(tái)站,采用的地震儀是桔燈勘探公司研發(fā)的Atlas智能節(jié)點(diǎn)式地震儀,儀器主頻為4.5 Hz,儀器道間距設(shè)置為2.5 m,采樣率設(shè)置為1 000 Hz。由于路面為土層與硬化路面相間,儀器放置方式采用鋼釘入土與鐵板放置兩種方式混合進(jìn)行,布設(shè)的線(xiàn)路主要沿圖1中的4條線(xiàn)路展開(kāi)。

研究區(qū)位于鬧市區(qū)的主干道,為成像工作提供了充足的噪聲信號(hào)。為了降低車(chē)輛行人對(duì)設(shè)備采集工作的影響,所采用的連續(xù)記錄時(shí)間段(8:00—20:00)避開(kāi)了車(chē)流量峰值,選在了車(chē)流量稍小的時(shí)間段(20:00—24:00)進(jìn)行數(shù)據(jù)的連續(xù)采集。圖2中紅色和黑色波形分別給出了第21道節(jié)點(diǎn)式地震儀記錄到的前2 h(20:00—22:00)和后2 h(22:00—24:00)噪聲頻譜,可以看出在數(shù)據(jù)記錄的4 h內(nèi),噪聲的主頻主要都集中在3～40 Hz。根據(jù)Lak et al.[18]的研究可知,在研究區(qū)域內(nèi)頻率為10～40 Hz的噪聲主要來(lái)自于交通噪聲,10 Hz以下的噪聲主要來(lái)自于人類(lèi)活動(dòng)。同時(shí)也觀(guān)察到數(shù)據(jù)記錄前2 h的噪聲水平明顯大于后2 h,推測(cè)可能與晚上22:00后的車(chē)流量以及人類(lèi)活動(dòng)大幅下降有關(guān)。由于噪聲互相關(guān)波形包含的面波信號(hào)能量一定程度上依賴(lài)持續(xù)穩(wěn)定高能的背景噪聲[19],因此本研究選擇噪聲水平相對(duì)較高的時(shí)間段20:00—22:00進(jìn)行數(shù)據(jù)的連續(xù)采集。

圖2 第21道節(jié)點(diǎn)式地震儀噪聲頻譜圖

本研究所使用的地震背景噪聲互相關(guān)函數(shù)計(jì)算處理流程與Bensen et al.[20]類(lèi)似,這一數(shù)據(jù)處理流程在淺層速度結(jié)構(gòu)成像的研究中已經(jīng)得到了廣泛應(yīng)用[21-23]。由于主要針對(duì)100 m以淺的地震波速度結(jié)構(gòu)進(jìn)行研究,根據(jù)曾求等[24]對(duì)于高頻互相關(guān)面波提取的實(shí)驗(yàn)結(jié)果,本次選擇將數(shù)據(jù)裁剪到300 s并且降采樣到100 Hz以提高計(jì)算效率,同時(shí)進(jìn)行時(shí)間域的歸一化以及頻率域的譜白化處理,然后計(jì)算地震背景噪聲互相關(guān)函數(shù),最后進(jìn)行疊加以提高互相關(guān)函數(shù)的面波信號(hào)噪聲比。由于在計(jì)算互相關(guān)函數(shù)前選擇窄帶濾波能夠在一定程度上提高互相關(guān)函數(shù)的面波信噪比,因此根據(jù)噪聲水平(圖2),筆者團(tuán)隊(duì)選擇在歸一化前對(duì)原始數(shù)據(jù)進(jìn)行3～40 Hz的帶通濾波。對(duì)于5 Hz以上的面波互相關(guān)函數(shù)的提取,采取通過(guò)長(zhǎng)時(shí)間的疊加來(lái)提高互相關(guān)函數(shù)信噪比,在疊加時(shí)間為0～1 h時(shí)信噪比有較大幅度的提升,在疊加時(shí)間為2 h時(shí)信噪比已經(jīng)趨于穩(wěn)定[24],因此本研究所使用的疊加時(shí)間為2 h。

根據(jù)以上參數(shù)設(shè)置,本文計(jì)算了研究區(qū)20:00—22:00 2 h的地震流動(dòng)臺(tái)站垂直—垂直分量互相關(guān)函數(shù)。圖3為以S01道為虛擬源、S02-S33道為接收臺(tái)站的互相關(guān)函數(shù)剖面圖,濾波頻帶為3～40 Hz。從圖3中可以看出有明顯面波頻散現(xiàn)象,其中較高頻的面波信號(hào)視速度約200 m/s,較低頻的面波信號(hào)視速度約300 m/s。使用基于波形時(shí)頻分析(Frequency-Time Analysis,FTAN)的do_mft程序[25-26]對(duì)互相關(guān)波形進(jìn)行面波頻散提取,通過(guò)將互相關(guān)波形轉(zhuǎn)換到時(shí)間域,根據(jù)波形能量給出頻譜能量分布圖,然后參考能量峰值挑選對(duì)應(yīng)的頻散點(diǎn)作為互相關(guān)函數(shù)的面波群速度頻散曲線(xiàn)參與淺層速度結(jié)構(gòu)的反演。圖4為3～40 Hz的S21-S26道頻散測(cè)量示意圖,其中群速度(U)主要集中在100～350 m/s,并且表現(xiàn)出明顯的正頻散特征,其中能夠獲取穩(wěn)定頻散點(diǎn)的頻帶范圍約3～25 Hz。

圖3 S01道與其他道互相關(guān)函數(shù)剖面圖

圖4 S21-S26道頻散測(cè)量示意圖

2 速度結(jié)構(gòu)反演與三維模型可視化

本文采用移動(dòng)虛擬源和接收排列獲得的相速度頻散曲線(xiàn),反演得到排列中點(diǎn)下方分層速度結(jié)構(gòu)[27]。再結(jié)合頻散提取結(jié)果,根據(jù)瑞雷波(Rayleigh Wave)反演的最大深度約0.5～0.67倍波長(zhǎng)[28],能夠獲取最低頻頻散點(diǎn)約3 Hz、視速度約380 m/s、對(duì)應(yīng)的敏感深度約85 m,最高頻頻散點(diǎn)約25 Hz、視速度約150 m/s、對(duì)應(yīng)的敏感深度約3 m,因此在進(jìn)行速度結(jié)構(gòu)反演時(shí),構(gòu)建一維反演模型的深度范圍為0～100 m,以覆蓋本研究頻散點(diǎn)對(duì)應(yīng)的所有敏感深度范圍。為了得到合適的一維S波速度反演的初始模型,本文首先對(duì)區(qū)域的平均頻散曲線(xiàn)(圖5)進(jìn)行反演得到平均一維S波速度模型,這里采用Hermann和Ammon的線(xiàn)性反演程序(surf96)[26],并且給予不同層不同的反演權(quán)重[29],初始速度模型由當(dāng)?shù)氐貙臃謱忧闆r結(jié)合群速度給出。

圖5 平均一維S波速度反演實(shí)例

基于平均一維S波速度反演模型(圖5),筆者團(tuán)隊(duì)對(duì)研究區(qū)所有測(cè)量點(diǎn)的一維S波速度結(jié)構(gòu)進(jìn)行了反演,獲取了淺層深度0～100 m的三維速度結(jié)構(gòu)結(jié)果。圖6為深度0～60 m的常規(guī)二維速度結(jié)構(gòu)剖面成像結(jié)果。整體而言,淺部人工填土層(～0～2 m)以及黏土層(～6～10 m)的S波速度較小,均<250 m/s,厚度約10～20 m;相較而言,在深度>20 m時(shí)S波速度開(kāi)始逐漸增大,對(duì)應(yīng)下層黏土夾碎石層(～10～15 m)和灰?guī)r層(～9～15 m),S波速度約250～400 m/s。從不同剖面之間的對(duì)比來(lái)看,0～50 m的S波速度結(jié)構(gòu)表現(xiàn)出明顯的橫向不均勻性,這種橫向不均勻性在深度20～50 m間表現(xiàn)得尤為明顯,每個(gè)速度剖面在該深度區(qū)間都存在較大的低速異常區(qū)域。

圖6 0～60 m垂向二維速度結(jié)構(gòu)剖面圖

在地下溶洞結(jié)構(gòu)探測(cè)中,準(zhǔn)確判斷異常區(qū)域之間的連通性顯得尤為重要,然而傳統(tǒng)二維速度結(jié)構(gòu)剖面只能提供一個(gè)平面上的速度變化信息,無(wú)法直觀(guān)地展示另一維度方向上的變化,因此很難準(zhǔn)確判斷異常區(qū)域之間的連通性。三維速度結(jié)構(gòu)圖可以提供地下速度結(jié)構(gòu)在三個(gè)維度上的變化信息,包括水平和垂直方向,這樣可以更準(zhǔn)確地確定地下異常區(qū)域的位置、形狀和大小,并能直觀(guān)地觀(guān)察到地下異常區(qū)域的連通性,因此本研究根據(jù)軟件VTK(Visualization Tool Kit)模型規(guī)范,基于三維點(diǎn)云數(shù)據(jù)構(gòu)建了點(diǎn)云數(shù)據(jù)到VTK三維地震速度結(jié)構(gòu)模型的建立方法。通過(guò)輸入獲取的淺層S波三維速度結(jié)構(gòu)點(diǎn)云數(shù)據(jù)直接生成VTK模型,并通過(guò)Paraview軟件實(shí)現(xiàn)三維可視化(圖7)。VTK模型可以輕松實(shí)現(xiàn)模型的平移、剖切、數(shù)值篩選、透明化等操作,適合對(duì)淺層地震速度異常區(qū)域的分析。

圖7 三維速度結(jié)構(gòu)模型可視化示意圖

3 分析與討論

在反演研究區(qū)三維速度結(jié)構(gòu)后,需要驗(yàn)證該結(jié)果的準(zhǔn)確性,因此將獲取的3線(xiàn)二維速度結(jié)構(gòu)剖面與該區(qū)域的鉆孔資料進(jìn)行對(duì)比驗(yàn)證(圖8)。圖8-a中黑色和白色透明長(zhǎng)條分別展示了未挖到溶洞以及挖到溶洞的鉆孔深度,圖8-b展示了不同鉆孔資料與橫向低速異常區(qū)域的深度差異,同時(shí)藍(lán)色誤差棒展示了對(duì)應(yīng)鉆孔的差異大小,綜合計(jì)算鉆孔驗(yàn)證得到的深度誤差約±1.846 62 m,滿(mǎn)足一般地鐵注漿回填的深度誤差要求。

圖8 3線(xiàn)二維速度結(jié)構(gòu)剖面與鉆孔資料對(duì)比圖

針對(duì)研究區(qū)溶洞區(qū)域的分析與劃分,一般而言由于體積效應(yīng),在地震速度結(jié)構(gòu)圖上溶洞通常表現(xiàn)為低速異常,然而縱剖面只能展示速度在深度方向的變化情況,無(wú)法很好地與研究區(qū)前期通過(guò)鉆孔推測(cè)的溶洞區(qū)域進(jìn)行對(duì)比。筆者團(tuán)隊(duì)繪制了研究區(qū)從5～45 m每間隔5 m的水平二維速度結(jié)構(gòu)在深度上的切面分布情況(圖9),圖9中白色透明區(qū)域推測(cè)為溶洞分布區(qū)域。

圖9 水平二維速度結(jié)構(gòu)剖面圖

從縱向上看,研究區(qū)20 m以?xún)?nèi)的速度結(jié)構(gòu)整體偏小(<250 m/s),這與人工填土層與黏土層的參考S波速度(150～250 m/s)相吻合;當(dāng)深度為20 m以下時(shí),研究區(qū)內(nèi)的東西部S波速度開(kāi)始上升至約400 m/s,伴隨著在推測(cè)溶洞區(qū)域出現(xiàn)相對(duì)低速區(qū)域,但在該深度上低速區(qū)域大面積成片出現(xiàn),推測(cè)可能是由于淺層低速層的沉積作用;當(dāng)深度為25～40 m時(shí),在推測(cè)溶洞區(qū)域Y1、Y2、Y3區(qū)均出現(xiàn)明顯的低速異常,其中Y2、Y3區(qū)深度較Y1淺,Y2、Y3的頂部在25 m左右,而Y1區(qū)的頂部在35 m左右,這與該區(qū)域的鉆孔數(shù)據(jù)也比較吻合;當(dāng)深度>45 m后,Y1、Y2、Y3區(qū)均恢復(fù)正常速度,推測(cè)在這3個(gè)區(qū)域內(nèi)的溶洞底部約45 m。橫向上看,Y1的低速區(qū)在南部與3線(xiàn)的居民區(qū)相連,Y2、Y3區(qū)在南部可能與3線(xiàn)的居民區(qū)相連、東北部與1線(xiàn)的居民區(qū)相連。

圖10通過(guò)Paraview軟件展示了研究區(qū)的三維速度低速異常分布圖,利用速度分割工具對(duì)三維模型進(jìn)行剖切并表示出了在25、30、35、40 m的低速異常區(qū)。圖10-b中紅色虛線(xiàn)區(qū)域表示基于三維模型推測(cè)出的溶洞范圍。通過(guò)與前期鉆孔推測(cè)出的溶洞范圍(黃色透明區(qū)域)Y1、Y2、Y3區(qū)進(jìn)行對(duì)比可以發(fā)現(xiàn),在35 m處兩者結(jié)果較吻合;此外還可以發(fā)現(xiàn)Y1、Y2、Y3區(qū)并非資料推測(cè)的三塊獨(dú)立區(qū)域,并且Y2、Y3區(qū)的溶洞腔比Y1區(qū)更淺,這與鉆孔數(shù)據(jù)對(duì)應(yīng)較好。圖10-a中的紅圈區(qū)域是一塊細(xì)長(zhǎng)的通道,連接了Y1與Y3區(qū)并向南部4線(xiàn)居民區(qū)延展,推測(cè)Y1和Y3處部分注漿往南部匯聚并擴(kuò)散到4線(xiàn)居民區(qū)下的溶洞處;同時(shí)綠色區(qū)域?qū)?yīng)東北部的1線(xiàn)居民處,推測(cè)該區(qū)域可能是Y2與Y3區(qū)向北部的延伸;在研究區(qū)東北部可能存在另一較大的溶洞,該推測(cè)結(jié)果在后續(xù)地鐵對(duì)1線(xiàn)和3線(xiàn)圍擋鉆孔施工工作中得到了驗(yàn)證,并且水平向的溶洞推測(cè)區(qū)域誤差約2 m。

圖10 三維速度低速異常分布圖

4 結(jié)論

本文利用節(jié)點(diǎn)式地震儀在武漢某地鐵施工區(qū)域進(jìn)行地震背景噪聲觀(guān)測(cè)實(shí)驗(yàn),并通過(guò)地震背景噪聲成像與自動(dòng)化三維模型建立程序獲得以下結(jié)論:

(1) 通過(guò)對(duì)連續(xù)地震觀(guān)測(cè)記錄進(jìn)行地震背景噪聲互相關(guān)函數(shù)計(jì)算,得到了研究區(qū)3～25 Hz的瑞雷面波頻散數(shù)據(jù);隨后通過(guò)頻散數(shù)據(jù)反演成像,得到了淺層地震三維散點(diǎn)結(jié)果;接著利用自動(dòng)化三維模型建立程序?qū)υ摻Y(jié)果進(jìn)行三維成像,得到了研究區(qū)0～50 m的三維微動(dòng)近地表速度結(jié)構(gòu)成像結(jié)果。

(2) 通過(guò)與鉆孔資料對(duì)比,發(fā)現(xiàn)垂向二維速度結(jié)構(gòu)剖面低速異常區(qū)與鉆孔資料較為一致。水平向二維速度結(jié)構(gòu)剖面結(jié)果顯示,研究區(qū)淺層25 m以下存在低速異常區(qū),該區(qū)域的位置分布與前期鉆孔資料推測(cè)的溶洞位置水平分布吻合。

(3) 三維微動(dòng)近地表速度結(jié)構(gòu)成像結(jié)果顯示,前期鉆孔資料推測(cè)的溶洞區(qū)域往南和東北方向延伸,推測(cè)1線(xiàn)和2線(xiàn)、3線(xiàn)和4線(xiàn)之間存在溶洞通道,導(dǎo)致注漿時(shí)向南北方向發(fā)散。后續(xù)在1線(xiàn)和4線(xiàn)的鉆孔資料中證實(shí)了這一推測(cè)。

本文成功地揭示了地鐵施工區(qū)的溶洞分布特征,這些溶洞的大小和位置對(duì)注漿處理具有重要的指導(dǎo)作用。研究結(jié)果表明,利用地震背景噪聲互相關(guān)函數(shù)可以有效地成像地鐵建設(shè)區(qū)域的速度結(jié)構(gòu),為城市道路建設(shè)中的地下溶洞注漿處理提供了一種新的方法。