11月
JOURNALOFCHINACOALSOCIETY
煤 炭 学 报
Vol.44 No.11 Nov.
2019
张胜,何文超,黎永索,等.基于小波变换时能密度法的隧道空洞充填物识别[J].煤炭学报,2019,44(11):3504ZHANGSheng,HEWenchao,LIYongsuo,etal.Identificationoftunnelcavityfillingsbytime-energydensityanalysis
移动阅读
2019.0921
basedonwavelettransform[J].JournalofChinaCoalSociety,2019,44(11):3504-3514.doi:10.13225/j.cnki.jccs.-3514.doi:10.13225/j.cnki.jccs.2019.0921
基于小波变换时能密度法的隧道空洞充填物识别
张 胜1,2,何文超3,黎永索1,胡 达1,蔡 鑫1
(1.湖南城市学院土木工程学院,湖南益阳 413000;2.休斯敦大学机械工程系,德克萨斯州休斯敦 77204;3.长沙理工大学土木工程学院,湖南长沙 410114)
摘 要:灰岩地区隧道修建过程中,制约安全高效建设的最主要因素是隧道掘进工作面前方的不良地质情况,加之溶洞、断层与破碎带等不良地质具有较强的隐蔽性,如何提高隧道短距离超前地质预报的判释水平显得尤为重要。为便于对地质雷达信号图谱特征进行定量解释,开展了基于时域有限差分法的正演模拟与室内物理模型实验。在小波变换与奇异性检测原理的基础上,构造了与地质雷达发射子波波形相似度高的雷达小波基,添加到小波分析工具箱中,提出了一种新的基于雷达小波基的小波变换时能密度法,将其应用于空洞充填物正演模拟与模型实验地质雷达信号的定量识别,并与波形分析法、Db4小波变换法和雷达小波变换法的识别结果进行比较。研究结果表明,波形分析法虽能有效识别空洞充填物的尺寸大小,但需通过反射系数等先验知识来确定空洞充填物的界面反射位于波峰或波谷;Db4小波变换模极大值法易得到地质雷达信号奇异点,但奇异点的提取位置与不同小波基的时频局部化特征有关,Db4小波变换法识别结果的相对误差约15%;雷达小波变换模极大值法与小波变换时能密度法的识别效果均较好,但小波变换时能密度法的分辨率更高,且不需选择最优尺度。当空洞充填物的尺寸大于电磁波波长时,小波变换时能密度法识别结果的相对误差小于5%。
关键词:隧道工程;空洞充填物;地质雷达;小波变换时能密度法;正演模拟
中图分类号:U45 文献标志码:A 文章编号:0253-9993(2019)11-3504-11
Identificationoftunnelcavityfillingsbytime-energydensityanalysis
basedonwavelettransform
(1.SchoolofCivilEngineering,HunanCityUniversity,Yiyang 413000,China;2.DepartmentofMechanicalEngineering,UniversityofHouston,Houston
ZHANGSheng1,2,HEWenchao3,LIYongsuo1,HUDa1,CAIXin1
77204,USA;3.SchoolofCivilEngineering,ChangshaUniversityofScience&Technology,Changsha 410114,China)
Abstract:Inthetunnelconstructionatlimestonearea,theunfavorablegeologicalconditionsaheadthetunnelfaceisthemostimportantfactortorestrictthesafeandefficientconstructionofthetunnel.Inaddition,theadversegeologicaloftunnelshort-rangeadvancedgeologicalpredictionisparticularlyimportant.Inordertoimprovethequantitativein-terpretationofgroundpenetratingradar(GPR)signalscharacteristics,aforwardmodelingbasedontimedomainfinitebodiessuchaskarsts,faultsandfracturezoneshavestrongconcealment.Thus,howtoimprovetheinterpretationleveldifferencemethodandanexperimentalstudyofindoorphysicalmodelwerecarriedout.Basedontheprincipleof
收稿日期:2019-07-01 修回日期:2019-10-10 责任编辑:郭晓炜
基金项目:国家自然科学基金资助项目(51608183);国家留学基金资助项目(201808430341);湖南省自然科学基金资助项目(2018JJ3022) 作者简介:张 胜(1984—),男,湖南安化人,讲师,博士。E-mail:zhangsheng0403311@163．com 通讯作者:黎永索(1974—),男,湖南汨罗人,教授,博士。E-mail:liyongsuo@126．com
wavelettransformandsingularitydetection,aradarwaveletbasiswithahighsimilaritytoGPRemissionwaveletwas
第11期张 胜等:基于小波变换时能密度法的隧道空洞充填物识别
3505
constructedandaddedtothewaveletanalysistoolbox.Usingtheradarwaveletbasisasafoundation,anewtime-energydensityanalysisbasedonwavelettransformmethod(TEDAWT)wasproposed.TheTEDAWTmethodwasappliedtothequantitativeidentificationofthecavityfillingsintheGPRforwardmodelingandphysicalexperiment,andcom-paredwiththerecognitionresultsofthewaveformanalysismethod,Db4wavelettransformmethodandradarwaveletwhichislocatedinacrestortrough.TheDb4wavelettransformmodulusmaximummethodiseasytoobtainthesingu-characteristicsofdifferentwaveletbases.TherelativeerroroftherecognitionresultsusingDb4wavelettransformmeth-transformmethod.Theresultsshowthatthewaveformanalysismethodcaneffectivelyidentifythesizeofcavityfillings,butreflectioncoefficientandotherpriorknowledgearerequiredtodeterminetheinterfacialreflectionofcavityfillingslarpointsofGPRsignals,buttheextractionpositionofsingularpointsisrelatedtothetime-frequencylocalizationodisabout15%.TherecognitioneffectoftheradarwavelettransformmodulusmaximamethodandtheTEDAWTmethodarebothbetter,whiletheTEDAWTmethodhashigherresolution,anddoesnotneedtoselecttheoptimalscale.Whenthesizeofthecavityfillingsislargerthanthewavelengthoftheelectromagneticwave,therelativeerroroftheTEDAWTmethodislessthan5%.transform;forwardmodeling
Keywords:tunnelengineering;cavityfillings;groundpenetratingradar;time-energydensityanalysisbasedonwavelet
近年来,随着国家基础设施建设的大力发展和高
段。地质雷达信号分析方法普遍采用基于傅里叶理论的仪器自带分析软件,而傅里叶理论没有时频局部化分析能力,无法对隧道掘进工作面前方的地质情况进行精确定位。近年来,随着科学技术的不断进步,小波变换、S变换、复信号分析以及希尔伯特-黄变换等技术在地质雷达信号阈值去噪与特征提取等方面具备了良好的信号解译能力[11-15]。由于小波变换具有许多优良的特性,并且理论比较成熟,因而在工程领域应用最为广泛。然而,采用经典小波对地质雷达信号进行小波时频局部化分析时,选取何种尺度的小波基通常依赖于技术人员的经验,虽然由经验选取的小波基能够对地质雷达信号进行一般化处理,但是小波基天然无法准确选取的事实,阻碍了小波时频局部化分析技术在地质雷达信号处理中的应用。因此,如何利用最新的信号处理方法,构造与地质雷达信号特征相似度高的最优小波基,将是小波理论运用于地质雷达信号处理中亟需突破的难题。
笔者以此为切入点,在小波变换与奇异性检测原理的基础上,构造与地质雷达发射子波波形契合度高的小波基,然后将该小波基添加到MATLAB小波分析工具箱中,提出一种新的基于雷达小波基的小波变换时能密度法,将其应用于隧道空洞充填物地质雷达正演模拟与室内模型实验采集信号的定量识别与分析,以期为隧道短距离超前地质预报定量识别不良地质体提供技术支持。
新技术的不断应用,全国范围内交通工程建设呈现一派繁荣景象。隧道工程因其具有改善线形、克服高程或地形障碍、保护生态环境等优势,已成为交通工程建设中最重要的控制性工程之一。截至2018年底,中国隧道建设包括铁路隧道、公路隧道以及城市轨道交通隧道共计36103km,正在施工的各类隧道我国在隧道工程规模、工程质量与修建技术等方面已成为世界上发展最为迅速的国家。尽管隧道工程建设发展迅速,但在隧道修建过程中仍存在许多难题,其中,行业内公认的隧道建设难点之一是隧道掘进工作面前方的不良地质条件。由于修建隧道所处位置的地质情况较为复杂且地形起伏较大,在修建前期,限于现有的地质勘察技术手段,难以对工程区域内复杂地形地貌进行准确且全面地掌握,且隧道沿线出现的一些不良地质情况(如溶洞、断层与破碎带等)对后期的修建、维护与运营带来了巨大的安全隐患[2-4]。
随着国内外学者在地质探测领域的技术突破,使约20000km,计划修建的隧道约20000km[1]。目前,
得无损探测技术发展迅速且广泛应用于实际工程。隧道施工超前地质预报方法有地质分析法、超前钻探法、地震反射波法、地质雷达法、瞬变电磁法、陆地声纳法、地震波层析成像法等,每种方法都有各自的优、缺点及适用范围[5-8]。地质雷达因其具有分辨率高、图像直观以及对施工干扰少等优点,被广泛应用于隧30m左右)[9-10]。地质雷达信号特征分析是研究地质雷达电磁波在介质中传播规律的前提与基础,也是对地质雷达信号图谱特征进行定量解释的最有效手道短距离超前地质预报中(预测掘进工作面前方
1 识别分析方法
1．1 波形分析法
地质雷达超前地质预报的工作原理是利用地质
3506
煤 炭 学 报
2019年第44卷
雷达发射天线发射高频电磁波以脉冲的形式传播至隧道掘进工作面前方的岩土介质中,遇到溶洞、断层与破碎带等界面时会发生反射,反射波由地质雷达接收天线所接收,最后将在地质雷达主机显示器上获得时-距剖面图像。通过对地质雷达时-距剖面进行处理、分析与解释,达到对隧道掘进工作面前方的地质情况进行短距离超前地质预报的目的[8-10]。
地质雷达回波反射信号可以表示为
小波变换模极大值,具有突出不同时频局部特征信号的能力。因此,采用小波变换模极大值法对信号奇异性特征进行分析时,最重要的前提是选取一个最佳的尺度,在该尺度下才能最大限度显示信号的奇异性与突变成分。
1．3 小波基的构造与添加
小波基函数的选取在小波变换分析中具有不唯一性。采用不同小波基分析同一信号,所得结果也会
x(t)=A(t)cos[2πf0t+θ(t)]
(1)
式中,x(t)为地质雷达回波反射信号;A(t)为振幅函数;t为记录时间;f0为中心频率;θ(t)为相位函数。
2异种不同介质的分界面时地质雷达发射的脉冲电磁波遇到岩土层中存在
,将在界面上产生反射现象,由于相对介电常数的差,其规律符合反射定律,
产生的电磁波能量大小取决于反射系数Γ。
Γ=
εε1-ε式中,ε1+ε2
2(2)
1,ε(2)2分别为分界面上、下介质的相对介电常数。由式可知,隧道围岩与充填物的相对介电常数的差异,直接影响了地质雷达图像的振幅大小、频率高低与相位起伏,这为地质雷达信号特征解译提供了良好的理论基础。1．2 小波变换模极大值法
小波变换模极大值法最突出的特征是选取信号的奇异点。信号奇异点与突变部分往往带有明显的个性化特征,它是一个信号区别于其他信号的重要特征之一。利用小波变换方法分析局部奇异性特征时,小波变换分解系数取决于信号在某一点相邻区域内的特征信号以及小波变换所选用的尺度。将任意信号f(t)的连续小波变换表示为Wf的某一邻域内的任意一点b,(a,b)[16],若对b0 则称(3)
,(aWba)≤有
Wf(0,bWf(a0,b0)
0,0)为f(t)在尺度换模极大值。
at)在尺度a0上的小波变换模极大值点f(0,b0)为f(a0上的小波变如果信号在某点或某区间内可微的次数越高,那么该信号越平滑,奇异性越弱。在数字信号处理中,常用李氏指数来表征信号的局部奇异性特征。信号f(t)在区间(t件为:存在常数1,tk2)>0,上具有一致李氏指数使得对所有的t∈(tα的充要条f(t)的小波变换W1,t2),信号f( 由式(3)~(4)Wa,b)满足如下不等式:
可以得出f(a,b)≤,信号产生突变的特征kaα
(4)
点就是小波变换的模极大值点。由于不同尺度下的
各有差异。构造适合被分析信号特征选用的最优小波基,是目前小波分析中亟需解决的难题。从理论上讲,构造一个新的小波基的前提是只要满足小波基的允许条件。因此,可根据地质雷达发射子波f(t)=t2为中心角频率e-atsin,由发射子波来进行自适应波形匹配ω0t的特点,其中,a为衰减系数;ω,
0构造与地质雷达发射子波波形相似度高的自适应雷达小波基,将该小波基添加到小波分析工具箱中,以期对地质雷达信号进行小波变换时作为小波基进行选用1．4 。
小波变换时能密度法
设信号f(t)连续小波变换系数的幅度平方在时间尺度平面上的加权积分与信号的能量成比例关系-,则满足下列能量守恒关系[1718]C1CψRdaa
2
RW)2db,=即
f(a,bR
f(t)式中,由式ψ为小波基容许条件∫∫;R为实数∫
2
dt(5)
。
(5)可知,小波变换幅度平方的积分同被分
析信号的能量成正比。在分析与处理非平稳信号时,由于海森堡(Heisenberg)测不准原理的限制,在时-频相空间中某一点的瞬时能量密度无法得知。因此,在某一指定时刻瞬时频率的说法是不存在的(a,b)。如果把Wf2/Cψa2认为是时间-尺度平面上的能量密度函数,那么可将Wfa和时间b为中心(、a尺度间隔为,b)2ΔaΔb/ΔCaψa2看作是以尺度、时间间隔为Δb的能量∫
。根据能量的概念,式(5)可以改写成
+¥+¥-¥
f(t)2
dt=
-¥
E(b)db
(6)式中,
E(b)=
小波变换中,尺度因子C1
ψ
∫
∫
+¥
0
Wf(a,b)
2
/a2da
(7)
率。因此,式(7)表示信号所有频带的能量随时间a在某种意义上对应于频b的分布情况,称为时能密度函数。
2 空洞充填物的正演模拟
2．1 时域有限差分法
地质雷达正演模拟是地质雷达图像解译的基础,
第11期张 胜等:基于小波变换时能密度法的隧道空洞充填物识别
3507
解译人员要想对地质雷达信号进行清楚掌握与定量解释,就必须事先了解电磁波在隧道掘进工作面前方的传播规律与图谱特征。时域有限差分法以其具有存储空间小与计算效率高等优势,成为了地质雷达正演模拟的最主要方法[19-21]。表示为
在无源场区域内,Maxwell方程的2个旋度可以
式中,c为电磁波在真空中的传播速度,c=3×2．2 空洞充填物的正演模拟
假设隧道掘进工作面前方为半无限空间连续的均匀介质,电磁波反射与折射均在二维平面内进行。为了研究隧道空洞内充填不同介质的地质雷达正演模拟图像的反射特征,设计了如图1所示的隧道空洞108m/s。
ìÑ×H=ε∂í
∂Et
+σE
î
Ñ×E=-μ∂H
-σmH(8)式中,H为磁场强度,A/m;E∂为电场强度t,V/m;ε为媒质的介电常数;σ为电导率,S/m;t为时间s;μ为相对磁导率,H/m;σm为等效磁导率,w/m。
时域有限差分法利用二阶精度的中心差分形式,
将Maxwell方程中的2个旋度由微分转成差分,电场与磁场在时间顺序上交替抽样,彼此相差半个时间步长。因此,二维电磁波的时域有限差分方程可以表示为
Enx+1
CB
Hny+1/2(i+(1i,/2,j)=CAEnxj)-Hn(y
+1i/,2
(j)
i-+1/2,j)-CB
Hnx+1/2(i,j+1/2)Δx-Hnx+1/2(i,j-1/2)(9)
H
nx
+1/2(i,j+1/2)=CPΔH
ynx
-1/2(i,j+1/2)-
CQ
Enz(i,j+1)Hn+1/2(i+1/2,j)=CPΔHy-Enz(i,j)yn-1/2CQ
Enz(i+x
1,j(i+1/2,j)+Δ)x
-Enz(i,j)(10)
其中,E
nx
+1样,彼此相差半个时间步长,H
nx
+1/2为电场与磁场在时间顺序上交替抽
,n为时间步;Ex方向的电场强度x为坐标轴
CA=2CP=22εε((ii,,jj))-+σσ(;(ii,E,jj)Δz为坐标轴)Δtz方向的电场强度;
t;CB=2ε(i,j)2Δ+σt(i,j)Δt;
和H2μμ((ii,,jj))-+σxσm和m((ii,y,jj)Δ方向上的磁场强度)Δtt;CQ=2μ(i,j)2Δ+σt;Δm(xi,和j)ΔΔt;H
xy分别为y分别为x和y方向上的空间步长;Δt为时间步长;(i,j)为节点坐标。
为了保证离散后时域有限差分方程组解的稳定收敛性,要求时间步长Δt和空间步长Δx与Δy满足以下关系式
Δt≤
æ1
c
2(11)
èΔ1xöø+æèΔ1yö2
ø
充填物的地电模型。图1 隧道空洞充填物的地电模型
隧道掘进工作面前方空洞充填物的地电模型参Fig．1 Geoelectricmodeloftunnelcavityfillings
数设置为:①区域范围:12m×20m,左上角为坐标原点,横坐标为隧道掘进工作面水平距离,纵坐标为探测深度;②目标体:长方形空洞,尺寸为2m×1m,空洞左侧与模拟区域边缘相距5m,埋深为2m;③围
岩灰岩的相对介电常数为6,电导率为0．001S/m,导磁率为1;④空洞充填物分别为空气、冰、干黏土与湿黏土。表1给出了隧道掘进工作面前方空洞内不同充填物的几何与物理参数。
表1 隧道空洞充填物的几何与物理参数Table1 Geometricandphysicalparametersoftunnel
cavityfillings
充填物介电常数
电导率/(S·m-1)
导磁率电磁波速度/(m·ns-1)空气1010．3000冰31．4×10-4干黏土81×10-310．湿黏土
12
1×10-3~~0．0．1111
11
0．17320．10608616
2．3 地质雷达正演模拟特征分析
根据隧道超前地质预报实际检测所使用的参数,模拟时采用的中心频率为100MHz,边界吸收条件为完全匹配层,激励源采用Ricker子波,网格的空间步
3508
煤 炭 学 报
2019年第44卷
长为0．05m,采样步长为0．05m,采样道数为227,总采样时间为200ns。
由于隧道开挖过程中,特别是灰岩地区,隧道掘进工作面前方经常会遇到空洞,且空洞内含有不同的充填物。为了更好反映地质雷达正演模拟的实际情况,本次模拟将在空洞中分别填充空气(ε=1)、冰(ε=3且ε<6)、干黏土(ε=8且ε>6)与湿黏土(ε=12且ε>6)。采用100MHz频率天线对如图拟,得到正演模拟雷达响应图像,并将其减去完整灰岩的地质雷达正演模拟图像,结果如图2所示。处理过程中对每种工况进行标记,如“y1”对应充填物为空气,“y2”对应充填物为冰,“y3”对应充填物
为干黏土,“y4”对应充填物为湿黏土。从图2可以看出:①地质雷达电磁波遇到充填物空气、冰、干黏土与湿黏土等不同介质时,界面处均存在反射,空洞充填物与灰岩的相对介电常数差异较小时,电磁波回波信号的反射不明显;②电磁波由空气进入灰岩,经灰岩再进入充填物,最后进入灰岩,各界面反射波依次被接收天线所接收。当空洞充填物的相对介电常数小于或大于灰岩的相对介电常数时,电磁波会在空洞充填物的第一或第二界面的界面处发生相位反相。而从图2所示的地质雷达时间剖面上无法识别电磁波遇到界面产生反射波的瞬时位置(究竟处于波峰还是波谷)。因此,需采用其他方法作进一步分析。
1所示的地电模型空洞内充填不同介质进行正演模
Fig．2 Time-sectionsofGPRsignalswithdifferenttunnelcavityfillings
图2 隧道空洞充填物的地质雷达信号时间剖面
3 空洞充填物正演模拟的定量识别
3．1 波形分析法的定量识别
地质雷达正演模拟是研究高频电磁波在隧道掘进工作面前方传播规律的有效途径。在本次正演模拟中,当电磁波信号从灰岩传播到空洞充填物时,根据空洞充填物相对介电常数与灰岩的差异,电磁波将在空洞充填物上、下界面处发生反射,由反射系数可以推断反射波相位在界面处的正负,即反射波的波峰或波谷。根据波峰与波谷或波谷与波峰之间的时间差,利用公式H=cΔt/(2ε)可以计算得到空洞充填物的厚度。从图2所示的地质雷达信号时间剖面中抽取最中间的一条信号,由此获得一条反映空洞充填物典型特征的单道信号,结果分别如图3所示。 由图3可知,根据空洞充填物与灰岩相对介电常数的关系,通过计算反射系数来确定反射界面位于反射波波峰还是波谷,从而得到电磁波在不同空洞充填
物上的第一、二反射界面,分别如图3所示的a1,a2,b1,b2,c1,c2和d1,d2。依次计算不同空洞充填物单道信号两特征点的时间差,计算结果分别为11．21,厚度分别为1．68,0．97,1．02和1．02m。采用波形分析法识别不同空洞充填物的相对误差分别为但须事先知道空洞充填物与灰岩相对介电常数的关系,即电磁波遇到异种界面时瞬时反射的位置究竟是反射波波峰还是波谷的位置。否则,难以准确确定反射界面在单道信号上的两特征点。
11．21,19．23和23．59ns。也就是说,空洞充填物的
68%,-3%,2%和2%,说明后三者的识别精度高。
3．2 小波变换模极大值法的定量识别
小波变换模极大值法是一种经典小波去噪方法,信号模极大值点的位置对应于信号奇异点。从理论上讲,小波变换尺度选取得越小,模极大值点与信号奇异点的对应位置越准确,但小尺度下的小波变换系数易受噪声的干扰,产生伪奇异点。
第11期张 胜等:基于小波变换时能密度法的隧道空洞充填物识别
3509
尺度较大时,可使噪声相对平滑,模极大值点位置易于寻找,但会使模极大值点位置产生一定的偏移。因此,在采用小波变换模极大值法识别信号的奇异点时,应选取适宜的尺度以避免小波变换所产生的交迭干扰。
大量试验与研究结果表明,离散小波变换选取合适小波基时,Daubechies系列小波基因其具有良好的
紧支撑性、光滑性与近似对称性等优点被广泛应用于地质雷达信号的处理与分析中[22-25]。经过比较Dau-bechies系列小波变换模极大值法提取信号奇异点的效果后,选取Db4小波作为小波变换模极大值法提取地质雷达信号奇异点的小波基。通过分别对图3所示的地质雷达单道信号进行Db4小波变换模极大值法分析,得到如图4所示结果。
Fig．3 Single-channelGPRsignalsoftunnelcavityfillings
图3 隧道空洞充填物的地质雷达单道信号
由图4可知,空洞不同充填物的第一、二反射界面在地质雷达正演模拟单道信号Db4小波变换模极大值曲线上的两特征点分别为a1,a2,b1,b2,c1,c2
Fig．4 Db4wavelettransformmodulusmaximumcurvesofGPRsignals(a=10)
图4 地质雷达单道信号Db4小波变换模值(a=10)
21．81和26．41ns,即空洞充填物的厚度分别为1．24,1．16,1．16和1．14m,则识别空洞充填物的相对误差分别为24%,16%,16%和14%,相对误差均大于10%。由此说明,小波分析工具箱中已有的小
和d1,d2。两特征点的时间差分别为8．26,13．45,
3510
煤 炭 学 报
2019年第44卷
波基分析地质雷达信号误差较大,亟需构造新的小波3．3 雷达小波变换模极大值法的定量识别
选择合适尺度的小波基是对地质雷达信号进行小波分析与处理的重要环节,其原因是不同小波基分析同一信号会产生不同的结果。在小波分析工具箱中拥有众多小波基可供挑选,但不能保证找基以提高地质雷达探测深度的准确率。
到与地质雷达特征信号相吻合的小波基。因此,需构造与地质雷达特征信号相似度高的自适应雷达小波基,将其添加到小波分析工具箱中,以供进行小波变换时选用。
根据前文1．3节所列方法构造雷达小波基,对空洞不同充填物的地质雷达正演模拟单道信号进行雷达小波变换模极大值法分析,结果如图5所示。
由图5可知,空洞不同充填物的上、下反射界面在地质雷达正演模拟单道信号雷达小波变换模极大值曲线上的两特征点分别为a1,a2,b1,b2,c1,c2
Fig．5 RadarwavelettransformmodulesmaximumcurvesofGPRsignals(a=10)
图5 地质雷达单道信号雷达小波变换模值(a=10)
果的相对误差分别为29%,-4%,2%和2%。
为了比较波形分析法、Db4小波变换模极大值法、雷达小波变换模极大值法和小波变换时能密度法的识别效果,将不同方法得到的测量值与相对误差汇总于表2。
由表2可知,当空洞充填物为空气时,不同方法得到的识别结果相对误差都比较大,依据电磁波在空气中的传播速度0．3m/ns,可以得到电磁波在空气中的波长为3m。由于充填物的厚度仅为1m,可能发生波形混叠,导致无法识别空洞充填物的第二界面;空洞充填物为冰、干黏土与湿黏土时,电磁波在相应介质中的波长分别为1．73,1．06与0．87m。波形分析法、雷达小波变换模极大值法与小波变换时能密度法的识别效果均较好,但波形分析法需预先了解围岩前方的介质属性,通过计算电磁波在界面上反射系数的正负,才能确定空洞充填物第一、二反射界面在电磁波反射信号上的具体位置;通过对图5与图6的比较不难看出,小波变换时能密度法的分辨率要比雷达小波变换模极大值法的高,表明小波变换时能密度
19．22和23．58ns,可以得到空洞不同充填物的厚度分别为1．29,0．96,1．02和1．02m,即识别结果的相对误差分别为29%,-4%,2%和2%。3．4 小波变换时能密度法的定量识别
采用已构造的雷达小波基作为小波变换用的小波基,在MATLAB语言平台工作环境下编写小波变换时能密度分析程序,对空洞不同充填物的地质雷达正演模拟单道信号进行小波变换时能密度法分析,结果如图6所示。
由图6可知,空洞不同充填物的上、下反射界面在地质雷达正演模拟单道信号小波变换时能密度曲
和d1,d2。根据两特征点的时间差8．61,11．09,
线上的两特征点分别为a1,a2,b1,b2,c1,c2和d1,d2。11．09,19．22和23．58ns,即空洞充填物的厚度分别为1．27,0．96,1．02和1．02m,由此可以说明识别结通过计算两特征点的时间间隔,结果依次为8．49,
第11期张 胜等:基于小波变换时能密度法的隧道空洞充填物识别
3511
法具有更好突出地质雷达信号奇异点的位置,从而验证了将地质雷达发射子波作为地质雷达信号小波分析用的小波基是切实可行的,解决了适合地质雷达信
号特征的小波基构造、添加与实现等问题,为分析地质雷达信号小波变换时选择小波基提出了一种新的方法。
Fig．6 Time-energydensityanalysisbasedonwavelettransformcurvesofGPRsignals
表2 不同方法识别空洞充填物的测量值与相对误差
Table2 Measurementvaluesandrelativeerrorsofcavityfillingsunderdifferentmethods
空洞充填物空气冰干黏土湿黏土
波形分析法测量值1．680．971．021．02
相对误差/%
68-322
Db4小波变换模极大值法测量值1．241．161．161．14
相对误差/%
24161614
雷达小波变换模极大值法测量值1．290．961．021．02
相对误差/%
29-422
小波变换时能密度法测量值1．270．961．021．02
相对误差/%
27-422
图6 地质雷达单道信号小波变换时能密度曲线
4 实验研究
4．1 方案设计
在众多研究方法中,室内模型实验是最直观与最可靠的研究手段之一。隧道空洞充填物模型实验中,模型尺寸与材料是影响室内模型实验结果好坏的重要因素。根据已有相关文献与资料,隧道模型尺寸相似比应尽量控制在15~30。采用地质雷达100M天线对隧道掘进工作面前方地质情况进行短距离超前地质预报预测的探测深度一般为30m,本次实验模型箱相似比(深度方向)取为20,即模型箱的高度为1．5m。地质雷达电磁波在空气中的传播速度为0．3m/ns,即波长为3m。根据模型尺寸相似比可将空洞的尺寸设置为10~20cm,本次实验空箱的尺寸取为10cm,除去壁厚与加工误差,空箱的实际尺寸为9．25cm。地质雷达用于识别隧道掘进工作面前
方空洞位置的精确程度主要取决于电磁波的传播速度,即围岩的介电常数。因此,应优先选择与灰岩相对介电常数相同或相近的介质作为模型的实验材料。干砂的相对介电常数为4~6,与灰岩的介电常数较为接近,具有重复使用、试制材料垃圾少以及成本低等优点,本次模型背景材料选择干砂。为便于砂箱物理模型成型,砂箱四周采用混凝土围制而成,长、宽和高分别设置为4,3和1．5m,以尽可能减少混凝土边界对地质雷达采集信号的干扰,具体如图7所示。砂箱中心埋设空箱模拟空洞,空洞充填物分别为空气、冰、干黏土与湿黏土的实物图如图8所示。4．2 地质雷达信号的采集
由于空箱的实际尺寸为9．25cm,在空洞充填物地质雷达检测的物理模型实验中,采用意大利RIS地质雷达,配以1600MHz天线。RIS地质雷达具有轻质便携、天线屏蔽效果好与抗干扰能力强等优点。主
3512
煤 炭 学 报
2019年第44卷
机技术参数如下:扫描速度4761扫/s,脉冲重复频率400kHz,采样点数128~8192,A/D转换16bit,叠加数1~32768,分辨率5psec,动态范围>160dB,信噪比>160dB。按照奈奎斯特采样定理,采样频率2时,雷达信号失真很明显。地质雷达在实际应用过程中,为了达到良好的探测效果,其采样频率应为天线主频的10倍以上,本次使用1600MHz天线所设采样频率为32GHz。根据天线的中心频率,将实验记录时间设为16ns和叠加次数为512。
Fig．7 Physicalexperimentboxoftunnelcavityfillings
图7 隧道空洞充填物的模型箱
至少应为反射波最高频率的2倍,但在频率比仅仅为
依据空洞充填物埋设在模型中的具体位置已知,地质雷达检测时沿一条直线匀速前进。地质雷达数物为空气,“y6”为冰,“y7”为干黏土,“y8”为湿黏土。
Fig．8 Physicaldrawingsofdifferenttunnelcavityfillings
图8 隧道空洞不同充填物的实物
采用地质雷达所携带的分析软件对采集的原始信号进行静校正、去直流点漂移、增益以及带通滤波等处理,得到空洞不同充填物的实测地质雷达信号时间剖面如图9所示。
据采集过程中对每种工况进行标记,如“y5”对应充填
由图9可知,地质雷达电磁波遇到空洞、冰、干黏土与湿黏土等介质时,界面处均存在反射,但不能显示具体的量值,尤其是空洞充填物的第一界面位于电磁波波形的波峰还是波谷无法得知。因此,需采用其他方法作进一步分析。4．3 空洞充填物的定量识别
从如图9所示的隧道空洞充填物实测地质雷达信号时间剖面中抽取最中间的1条信号,获得1条反映空洞充填物典型特征的单道信号。采用前文已构造的雷达小波作为小波变换用的小波基,在MATLAB语言环境下运行小波变换时能密度法分析程序,对空
Fig．9 Time-sectionsofmeasuredGPRsignalswithdifferenttunnelcavityfillings
图9 隧道空洞充填物的实测地质雷达信号时间剖面
洞不同充填物单道信号进行小波变换时能密度法分析,结果如图10所示。
从图10可以看出,空洞不同充填物的上、下反射界面在实测单道信号小波变换时能密度曲线上的两特征点分别为e1,e2,f1,f2,g1,g2与h1,h2。通过计算两特征点的时间间隔(分别0．9084,1．0650,1．9110与
2．4121ns),可以得到空洞充填物的厚度分别为13．63,9．36,9．87与8．92cm,则识别结果的相对误差为47．35%,1．19%,3．35%与-3．57%。尽管采用地质雷达1600MHz高频天线(探测深度浅、精度高),但仍无法识别10cm左右空洞,究其原因是电磁波在空气
第11期张 胜等:基于小波变换时能密度法的隧道空洞充填物识别
3513
中的波长为18．75cm,大于空洞的实际尺寸。采用小波变换时能密度法识别空洞充填物为冰、干黏土与湿黏土的精度高,相对误差均小于5%。由此表明,小波
变换时能密度法能成功识别空洞充填物第一、二界面地质雷达反射信号奇异点的位置,为隧道施工现场超前地质预报图像定量识别提供技术支持。
Fig．10 Time-energydensityanalysisbasedonwavelettransformcurvesofmeasuredGPRsignals
图10 实测单道信号小波变换时能密度曲线
5 结 论
寸大小,但需通过计算反射系数的正或负值等先验知识来确定空洞充填物的界面特征点处于波峰还是波谷,由于空洞充填物的介电常数存在大于或小于围岩介电常数的情况,导致判断过程较繁琐,因此不适用于隧道超前地质预报精细化定量检测。
(2)尽管小波变换模极大值法容易得到信号突(1)波形分析法虽能有效识别空洞充填物的尺
2019,39(4):537-544.
YANJinxiu.Achievementsandchallengesoftunnelingtechnologyin
[2] 钱七虎.隧道工程建设地质预报及信息化技术的主要进展及发
展方向[J].隧道建设,2017,37(3):251-263.
QIANQihu.Maindevelopmentsanddirectionsofgeologicalpredic-Construction,2017,37(3):251-263.
537-544.
Chinaoverpast40years[J].TunnelConstruction,2019,39(4):
tionandinformatizedtechnologyoftunnelconstruction[J].Tunnel[3] 程久龙,李飞,彭苏萍,等.矿井巷道地球物理方法超前探测研
究进展与展望[J].煤炭学报,2014,39(8):1742-1750.
CHENGJiulong,LIFei,PENGSuping,etal.Researchprogressanddevelopmentdirectiononadvanceddetectioninmineroadwayworkingfaceusinggeophysicalmethods[J].JournalofChina[4] ZHANGKai,TANNANTDD,ZHENGWenbo,etal.Prediction
ofkarstfortunnellingusingfuzzyassessmentcombinedwithgeolog-icalinvestigations[J].TunnellingandUndergroundSpaceTechnol-[5] 李术才,刘斌,孙怀凤,等.隧道施工超前地质预报研究现状及
发展趋势[J].岩石力学与工程学报,2014,33(6):1090-1113.LIShucai,LIUBin,SUNHuaifeng,etal.StateofartandtrendsofJournalofRockMechanicsandEngineering,2014,33(6):1090-1113.
ogy,2018,80:64-77.
CoalSociety,2014,39(8):1742-1750.
变的模极大值点,但选取不同小波基,可能会产生截然不同的结果。通过地质雷达正演模拟分析结果得知,小波变换模极大值法识别结果的误差约15%,误差相对较大,其原因是小波基的时频特征与被分析信号时频特征的相似性不高。
(3)在小波变换原理的基础上,从与地质雷达反
射子波波形相匹配出发,成功构造了地质雷达小波基。雷达小波变换模极大值法与小波变换时能密度法的识别效果均较好,但小波变换时能密度法的分辨率更高,压制随机干扰的能力更强,更重要的是不需优选最优尺度。当空洞充填物的尺寸大于电磁波波长时,小波变换时能密度法识别结果的相对误差均小于5%。
参考文献(References):
[1] 严金秀.中国隧道工程技术发展40年[J].隧道建设(中英文),
advancedgeologicalpredictionintunnelconstruction[J].Chinese
[6] LIShucai,LIUBin,XUXinji,etal.Anoverviewofaheadgeological
Technology,2017,63:69-94.
prospectingintunneling[J].TunnellingandUndergroundSpace
[7] 李天斌,孟陆波,朱劲,等.隧道超前地质预报综合分析方法
[J].岩石力学与工程学报,2009,28(12):2429-2436.
3514
煤 炭 学 报
LITianbin,MENGLubo,ZHUJin,etal.Comprehensiveanaly-sismethodforadvancedforecastofgeologyintunnels[J].ChineseJournalofRockMechanicsandEngineering,2009,28(12):2429-2436.
2019年第44卷
iccorrectionmethodbasedonwavelettransformmodulusmaxima360-368.
andthresholds[J].ChineseJournalofGeophysics,2017,60(1):
[8] LAIWWL,DÉROBERTX,ANNANP.Areviewofgroundpene-locatingandtestingtoimaginganddiagnosis[J].NDT&EInterna-tional,2018,96:58-78.
[17] 张胜,凌同华,刘浩然,等.模式自适应小波时能密度法及其在
微差爆破振动信号分析中的应用[J].煤炭学报,2014,39(10):2007-2013.
tratingradarapplicationincivilengineering:A30-yearjourneyfrom
ZHANGSheng,LINGTonghua,LIUHaoran,etal.Patternadaptedwavelettime-energydensitymethodanditsapplicationinmillisec-ondblastvibrationsignalanalysis[J].JournalofChinaCoalSocie-[9] 刘新荣,刘永权,杨忠平,等.基于地质雷达的综合超前预报技
术[J].岩土力学,2015,37(S2):51-56.
LIUvancedXingrong,LIUgeologicalpredictionYongquan,YANGtechnologyZhongping,etfortunnelsal.Syntheticad-[J].Rockand探地雷达探测采煤工作面隐伏钻杆研究SoilMechanics,2015,37(S2):51-56.basedonGPR学报宋劲,2014,39(3):537,王磊.-542.
[J].煤炭
SONGburiedJin,JournaldrillWANGpipebyLei.groundStudypenetratingonthehighradarprecisioninthecoaldetectionface[J].ofZHANGofenhancedWChinaY,HAOCoalSociety,2014,39(3):537-542.
NDTGPRdetectionT,CHANGofRFY,ettaggedal.Time-frequencyburiedplasticpipesanalysis[J].of
刘宗辉and不良地质识别方法,刘毛毛EInternational,2017,(92):88,周东-96.
[J].,等岩土力学.基于探地雷达属性分析的典型岩溶
,2019,40(8):3282-3290.LIUofpenetratingtypicalZonghui,LIUanomaliesMaomao,ZHOUradarattributeinkarsttunnelDong,etconstructional.Recognitionbasedonmethodground2019,40(8):3282-3290.
analysis[J].RockandSoilMechanics,
李尧模拟及其复信号分析,李术才,刘斌,等[J]..钻孔雷达探测地下不良地质体的正演
岩土力学,2017,38(1):300-308.LIsignalYao,LIversegeologicalanalysisShucai,LIUofboreholeBin,etradaral.ForwarddetectionsimulationforMechanics,undergroundandcomplex2017,
ad-38(1):300凌同华-308.
bodies[J].RockandSoilHHT1428.
分析法,张胜[,J].李升冉岩石力学与工程学报.地质雷达隧道超前地质预报检测信号的
,2012,31(7):1422-LINGformsingmethodTonghua,ZHANGSheng,LIShengran.Hilbert-huangtrans-chanicsgroundfordetectionsignaloftunnelgeologicalpredictionu-andEngineering,2012,31(7):1422penetratingradar[J].Chinese-Journal1428.
ofRockMe-周辉林[J].ZHOU电波科学学报,田茂,熊俊志.探地雷达回波信号的特征提取和分类
andneseclassificationHuilin,TIAN,2006,21(4):586ofgroundMao,penetratingXIONG-Junzhi.589.
radarechoFeaturesignals[J].extraction于生宝Journal法的CSAMT,郑建波ofRadio静态校正,高明亮Science,2006,21(4):586[J].,等.-589.Chi-地球物理学报基于小波变换模极大值法和阈值
,2017,60(1):360-368.
YUShengbao,ZHENGJianbo,GAOMingliang,etal.CSAMTstat-
[18] CHENGty,2014,39(10):2007basedJunsheng,YU-onwaveletDejie,YANG2013.
nalysistransform[J].Yu.NDTTime&-energyEInternational,
densitya-[19] 2005,38(7):569邓国文演模拟与超前探测,王齐仁,廖建平-572.[J].,物探与化探等.隧道不良地质现象的探地雷达正
,2015,39(3):651-656.DENGelingnomenaandGuowen,WANGtunnel[advancedJ].GeophysicaldetectionQiren,LIAOofJianping,etal.Forwardmod-&radargeochemicalinadverseexploration,2015,
geologicalphe-[20] 39(3):651刘伟质雷达正演实验研究,周斌,-甘伏平656.
,等[.J].隧道超前预报中不同性质充填溶洞地
现代隧道技术,2014,51(1):153-158.
LIUoffilledWei,ZHOUkarstcavesBin,GANinadvanceFuping,etgeologicalal.GPRpredictionforwardforsimulation[21] [J].钟声ModernTunnelingTechnology,2014,51(1):153-158.
tunnels
征——,—王川婴围岩及充填效应正演分析,吴立新,等.点状不良地质体钻孔雷达响应特
[J].岩土力学,2012,33(4):1191ZHONG-1195.
responseSheng,WANGChuanying,WULixin,etal.Boreholesimulationcharacteristicsofitssurroundingofpointrockunfavorableandfillinggeo-bodies:condition[J].ForwardradarRock[22] andLIUSoilsismethodHaoran,LINGMechanics,2012,33(4):1191forGPRTonghua,LIsignalsbasedDiyuan,et-1195.
optimalal.Abiorthogonalquantitativewavelet
analy-[23] [J].MOHANAJournalofCentralSouthUniversity,2018,25(4):879on-891.
betweenlandmineMA,ABBASandmine-likeAM,GOMAAtargetsusingML,etwaveletsal.Discrimination
andspectral
[24] 黄航如2(1):54analysis[J].-66.
NRIAGJournalofAstronomyandGeophysics,2013,
取的应用,程志平[J].工程地球物理学报,王小龙,等.小波分析在地质雷达反射目标提
,2012,9(1):99-102.HUANGplicationHangru,CHENG[25] Chinese师学明,Journalofwavelet张剑,刘梦花ofEngineeringanalysisZhiping,WANGto,等.小波变换在地质雷达干扰波压制中
Geophysics,2012,9(1):99GPRreflectiveXiaolong,ettargetextract[al.The-102.J].ap-的应用[J].工程地球物理学报,2008,5(3):279-286.
SHIwaveletXueming,ZHANGJian,LIUMenghua,etChinesetransformJournalofinEngineeringde-nosinggroundGeophysics,2008,5(3):279penetratingal.radarApplicationdata[of-286.
J].[10] [11] [12] [13] [14] [15] [16]
因篇幅问题不能全部显示，请点此查看更多更全内容