机载探地雷达数值模拟及逆时偏移成像

机载探地雷达可以用于人类无法到达的危险地区、植被严重覆盖的地下目标体探测,然而由于机载探地雷达的特殊性,影响机载探地雷达探测效果的因素包括天线的极化方向、天线的飞行高度以及地表粗糙度等.为了研究这些影响因素与探测效果之间的关系,用三维时间域有限差分模拟电磁波的传播过程,以沙漠地区地下空洞掩体的机载探地雷达探测为实例,分别模拟了不同天线极化方向、天线高度及地表粗糙度情况下的机载探地雷达剖面,分析了各因素对机载探地雷达探测地下空洞目标体的影响.天线极化方向与目标体走向垂直更有利于地下目标体探测;天线距离地表越近,可以获得更高分辨率的雷达剖面;沙漠地表起伏越大,雷达剖面中的散射杂波能量越强,浅部地下目标体信号容易被掩盖.为了消除起伏地形造成的散射杂波,提出用逆时偏移成像技术对共炮集机载探地雷达数据进行偏移成像,成像结果优于基尔霍夫偏移成像结果.     

地质雷达在乌拉山山前断裂探测中的应用

地质雷达是利用电磁波对地下不同电性介质进行探测的地球物理仪器，其探测速率快、分辨率高，可弥补探槽和其他地球物理方法存在探测盲区的缺陷，正在越来越多地应用于活动断层探测领域。本文以乌拉山山前断裂为例开展地质雷达探测工作，使用无人机正射影像技术对测线进行地形校正，获得断层浅部地质雷达图像。研究结果表明，本文研究方法能有效反映探槽揭露的地层单元和断层分布。本次探测中，雷达波形图像特征为:浅地表的土壤层反射波总体较弱；粗粒沉积为主的砾石层反射波总体较强，同相轴连续性好；细粒沉积为主的砂层反射波弱于砾石层，波形以中、高频为主，同相轴具有弱连续性；对于洪冲积地区，地质雷达能分辨具有一定特征的地层单元，这为剖面图像的断层识别提供了标志；通过无人机正射影像技术对地质雷达测线进行地形校正，有利于获得更为准确的探测结果。     

隧道检测中基本形状空洞探地雷达图像特征

空洞是隧道检测的重点探测病害之一,为了进一步提高探测空洞的水平,本文针对隧道衬砌中经常出现空洞的基本形状进行研究.以探地雷达理论为基础,运用时域有限差分法,建立隧道衬砌中圆形、矩形和直角三角形空洞的二维模型,进行探地雷达二维正演模拟.利用PVC管制作圆形空洞,利用纸板制作矩形和直角三角形空洞,依次将物理空洞模型分别埋于砂槽中,应用探地雷达进行探测.并在现场探测各种形状的空洞.数值模拟、实验和实测结果表明:探地雷达可以探测到圆形、矩形和直角三角形空洞的存在,分别可以确定圆心的水平位置,矩形的水平分布范围和直角三角形的水平分布范围.该研究对不同形状空洞的图像解释具有指导和参考作用.     

基于GprMax正演模拟的探地雷达根系探测敏感因素分析

应用探地雷达对植物根系进行探测的有效性已得到证实.但由于根系结构复杂,根围环境异质性强,针对基于探地雷达得到的根系探测数据的解读尚处于经验积累阶段.本研究首先通过对比根系探地雷达实测信号图像和模拟信号图像,证实了利用GprMax模拟探地雷达探测植物根系的有效性.其次通过定义不同根系空间结构和电性参数场景,模拟了不同条件下根目标反射信号的差异,并对影响探地雷达探测植物根系有效性的敏感因素进行了初步分析.模拟结果有助于探地雷达野外根系探测图谱的解译,为探地雷达在植物根系探测中的应用积累经验.     

探地雷达及EH4对钼矿采空区的综合探测

采空区严重威胁着地面建筑物和矿山的安全,对采空区探测方法的研究有着重大的现实意义.本文探讨了探地雷达及EH-4电导率成像法对某露天钼矿采空区探测的方法技术及其应用效果.探地雷达对浅层地段不仅有很高的探测精度和高分辨率,能达到对浅部采空区的精确定位.EH-4电导率成像系统的探测深度大,通过电阻率断面图的局部异常体,可以直观的来判断采空区的位置从而查明深部采空区,结合二者对采空区进行综合探测.     

岩溶勘察中的探地雷达技术及应用

岩溶的空间分布和物性特征为探地雷达技术进行岩溶地球物理探测提供了基本前提条件.在已有的探测实例中应用探地雷达方法基本查明了岩溶的分布范围和充填性质，从而为岩溶区的岩土工程勘察和评价提供了有效手段及依据.     

物探方法在预防路面塌陷中的应用

城市地下空间的开发利用不可避免地在基础设施建设中造成土体流失或土层沉降,形成的地下“病害体”易导致城市道路塌陷的发生。文章以滁州琅琊大道为例,采用地质雷达法和浅层地震勘探法对城市老城区排水管网周边进行探测。结果表明:①在探测深度不超过7.0 m范围内,地质雷达法可以满足探测要求;②在场地条件复杂、干扰较多的情况下,需辅以浅层地震勘探法。利用上述方法,可以更准确地查明地下“病害体”的位置,防止了路面塌陷的发生。     

GPS在探地雷达探测城市道路塌陷隐患中的应用研究

在城市道路塌陷隐患探测中,回溯定位探地雷达剖面中病害异常的现场位置是对其进行复测验证及处治修复的关键.测距轮属于相对定位技术,具有难以消除的累积误差,且探地雷达剖面越长,其累积误差越大,回溯定位病害异常的现场位置越难,而GPS测量定位属于绝对坐标定位技术,没有累积误差,但在复杂的城市环境中,其干扰因素较多,测量误差较大,影响病害现场回溯定位的精度和准确性.联合采用测距轮和GPS对探地雷达进行跟踪测量定位,测距轮按照固定道间距触发探地雷达和GPS同步采集、同步存储,使每一道探地雷达数据与GPS经纬度坐标一一对应、彼此关联、相互追溯;通过绘制GPS路径轨迹,直观形象展示探地雷达实测剖面的测线位置及其工作量完成情况;采用移动平均法对GPS原始数据进行平滑去噪处理,修正测量误差,消除局部突变点,使GPS路径轨迹更平顺,更符合探地雷达检测车实际行驶的路径轨迹,由此计算得到的路径轨迹长度稳定性更高、一致性更好,更符合探地雷达实测剖面长度,同时可以有效提高病害现场回溯定位的精度,显著加快病害异常点加密复测验证工作进度.     

探地雷达时域多分辨法(MRTD)三维正演模拟

应用小波伽略金方法,对Maxwell方程进行离散化,导出了DB2-MRTD算法的探地雷达3D差分公式、数值稳定性条件.在此基础上,开发了探地雷达MRTD(multi-resolution time domain)法正演模拟程序,该程序极大地提高了运算速度,改善了三维探地雷达正演方法,并利用该自制程序,对三角形金属体模型进行了正演模拟,得到了其相应的正演合成三维剖视图及切片图,通过对这些模拟结果进行分析,可以加深对三维雷达反射特征的认识,提高探地雷达探测的可靠性、准确度,同时也说明时域多分辨率法在探地雷达三维正演模拟中的有效性.     

非开挖特深管线的探测技术分析及展望

本文介绍了最近几年对非开挖施工大于3m的深埋管线探测的4种物探探测技术的最新进展,包括地面高精度磁测、井中磁梯度探测、探地雷达、电磁法探测等技术的工作原理、数据处理、应用的成功实例,并总结了各自探测技术的优点、缺点,讨论了物探方法在对非开挖管线的定位技术的发展趋势.     

Evaluation of trenching,ground penetrating radar (GPR) and electrical resistivity tomography (ERT) for sinkhole characterization

This paper explores the suitability and advantages of combining the trenching technique with geophysical surveys [ground penetrating radar (GPR) and electrical resistivity tomography (ERT)] for sinkhole characterization in a mantled karst area. The approach is applied to two active sinkholes concealed by anthropogenic deposits and formed by contrasting subsidence mechanisms; collapse and sagging. The ERT section acquired across the collapse sinkhole images the clayey fill of the depression as an obvious low resistivity area, showing the approximate location of the sinkhole edges. Spatially dense GPR surveys provide information on the position of the boundaries of the concealed subsidence structures and their three‐dimensional (3D) internal geometry, revealing the dominant subsidence mechanism. We illustrate the impact of several factors on the quality of the GPR data such as sinkhole size, nominal frequency of the antennas, antenna shielding, and the presence of backfilled excavations and above‐surface objects. Trenches provided detailed information on the subsurface structure of the sinkhole, subsidence magnitude, partitioning of the strain, and the position of the sinkhole edges, especially when they are deep enough and excavated across the central sector and perpendicular to the boundaries. The stratigraphic and structural relationships observed in the trench were then used to infer the spatial evolution of the sinkholes (e.g. enlargement), their kinematic behavior (episodic versus progressive), and to differentiate discrete subsidence events and their associated magnitude. Numerical dates were used to estimate average subsidence rates and the recurrence of subsidence events. Such integrated data sets may be used as an objective basis to forecast the future behavior of potentially damaging sinkholes and to assess the associated hazard and risk. Copyright © 2013 John Wiley & Sons, Ltd.     

Clay content evaluation in soils through GPR signal processing

The mechanical behavior of soils is partly affected by their clay content, which arises some important issues in many fields of employment, such as civil and environmental engineering, geology, and agriculture. This work focuses on pavement engineering, although the method applies to other fields of interest. Clay content in bearing courses of road pavement frequently causes damages and defects (e.g., cracks, deformations, and ruts). Therefore, the road safety and operability decreases, directly affecting the increase of expected accidents. In this study, different ground-penetrating radar (GPR) methods and techniques were used to non-destructively investigate the clay content in sub-asphalt compacted soils. Experimental layout provided the use of typical road materials, employed for road bearing courses construction. Three types of soils classified by the American Association of State Highway and Transportation Officials (AASHTO) as A1, A2, and A3 were used and adequately compacted in electrically and hydraulically isolated test boxes. Percentages of bentonite clay were gradually added, ranging from 2% to 25% by weight. Analyses were carried out for each clay content using two different GPR instruments. A pulse radar with ground-coupled antennae at 500 MHz centre frequency and a vector network analyzer spanning the 1–3 GHz frequency range were used. Signals were processed in both time and frequency domains, and the consistency of results was validated by the Rayleigh scattering method, the full-waveform inversion, and the signal picking techniques. Promising results were obtained for the detection of clay content affecting the bearing capacity of sub-asphalt layers.     

核电站防波堤探地雷达检测的适用性研究

在分析核电站防波堤常见安全隐患及其产生机理的基础上,采用目前较为先进的探地雷达检测技术,对实际的堤防进行模拟检测,探究探地雷达检测防波堤的适用性和准确性.分析表明,探地雷达作为一种无损检测方法对防波堤的检测是可行的,且检测准确度满足实际工程的要求,可以用于核电站海堤的日常安全检测和维修加固后的成果检验.     

The application of GPR and ERI in combination with exposure logging and retrodeformation analysis to characterize sinkholes and reconstruct their impact on fluvial sedimentation

This work illustrates the practicality of investigating sinkholes integrating data gathered by ground penetrating radar (GPR), electrical resistivity imaging (ERI) and trenching or direct logging of the subsidence‐affected sediments in combination with retrodeformation analysis. This mutidisciplinary approach has been tested in a large paleosinkhole developed during the deposition of a Quaternary terrace on salt‐bearing evaporites. The subsidence structure, exposed in an artificial excavation, is located next to Puilatos, a village that was abandoned in the 1970s due to severe subsidence damage. Detailed logging of the exposure revealed that the subsidence structure corresponds to an asymmetric sagging and collapse paleosinkhole with no clear evidence of recent activity. The sedimentological and structural relationships together with the retrodeformation analysis indicate that synsedimentary subsidence controlled channel location, the development of a palustrine environment and local changes in the channel pattern. GPR profiles were acquired using an array of systems with different antenna frequencies, including some recently developed shielded antennas with improved vertical resolution and penetration depth. Although radargrams imaged the faulted sagging structure and provided valuable data on fault throw, they did not satisfactorily image the complex architecture of the fluvial deposit. ERI showed lower resolution but higher penetration depth when compared to GPR, roughly capturing the subsidence structure and yielding information on the thickness of the high‐resistivity alluvium and the nature of the underlying low‐resistivity karstic residue developed on top of the halite‐bearing evaporitic bedrock. Data comparison allows the assessment of the advantages and limitations of these complementary techniques, highly useful for site‐specific sinkhole risk management. Copyright © 2016 John Wiley & Sons, Ltd.     

机载探地雷达的进展以及数值模拟

机载探地雷达可能解决危险环境或广域条件下的近地表探测问题,用于解决环境、生态或军事方面的问题.然而由于种种原因,该技术的发展却显得比较慢.为了推进该技术的发展,本文介绍了目前世界范围内机载探地雷达的进展,并利用时间域有限差分法对一些典型模型进行数值模拟,并用特定的偏移成像方法对模拟结果进行成像.目前存在的机载探地雷达主要有三种类型:第一种为将常规探地雷达天线悬挂在直升飞机上,第二种为针对机载探地雷达开发的雷达系统,第三种为具有探地能力的合成孔径雷达.数值模拟结果表明,不管是水平地面的情况下,还是起伏界面的情况下,机载探地雷达都能清楚探测一定深度范围内的地下目标.可见,机载探地雷达是存满希望的一种方法.     

探地雷达小波域三维波动方程偏移

阐述了矩阵多分辨分析理论中的标准形式与非标准形式,并以Hilbert算子为例,说明了算子多分辨表示的压缩效果,为小波域偏移算法奠定了理论基础.从三维雷达波动方程出发,利用爆炸反射原理和浮动坐标变换,推导出三维探地雷达波动方程差分格式,并通过方程分裂算法及小波多分辨算法,在小波域求解波场外推矩阵,进而得到探地雷达小波域三维波动方程偏移算法,在此基础上,开发了探地雷达小波域偏移处理程序,并把该程序应用于三个球体空洞的3-D正演结果及实际的雷达数据中,通过对比偏移处理前后的雷达资料,得知该三维偏移算法能使3-D正演剖面中的反射波归位、绕射波收敛,极大地提高了雷达剖面的分辨率,有利于探地雷达资料的地质解释.     

Quaternary deformation mapping with ground penetrating radar

Ground penetrating radar (GPR) was used in a programme of geological and hydrogeological investigations into the Quaternary of west Cumbria, UK. The investigations were part of an extensive programme to determine the suitability of the area for a deep radioactive waste repository. The hydrogeological characteristics of the drift deposits are important since they affect both recharge and discharge. The glacially derived Quaternary sediments include a variety of deformation structures related to their mode of deposition and subsequent modification by glacial and periglacial processes. These deformation structures range from variable scale thrusts (centimeter to tens of meter displacement due to proglacial thrusting) to small to medium scale faults, folds and collapse structures which are early or syn-sedimentary adjustments. From the GPR data a number of dipping reflectors, having dip angles in the range of 8–23°, were interpreted as thrust planes and related faults and folds have been clearly observed. Inwardly dipping reflectors, having a maximum dip angle of 36°, creating a characteristic `V' shaped anomaly, were also observed and interpreted to be linear collapse structures where glacitectonically emplaced slices of ice melted out causing sediments to slump. Lateral continuations, up to 90 m, of GPR-interpreted dislocations were mapped between survey lines, suggesting that such features could be laterally continuous on the scale of hundreds of metres.     

Detection of Hazardous Cavities Below a Road Using Combined Geophysical Methods

Assessment of the risk arising from near-surface natural hazard is a crucial step in safeguarding the security of the roads in karst areas. It helps authorities and other related parties to apply suitable procedures for ground treatment, mitigate potential natural hazards and minimize human and economic losses. Karstic terrains in the Salento Peninsula (Apulia region—South Italy) is a major challenge to engineering constructions and roads due to extensive occurrence of cavities and/or sinkholes that cause ground subsidence and both roads and building collapse. Cavities are air/sediment-filled underground voids, commonly developed in calcarenite sedimentary rocks by the infiltration of rainwater into the ground, opening up, over a long period of time, holes and tunnels. Mitigation of natural hazards can best be achieved through careful geoscientific studies. Traditionally, engineers use destructive probing techniques for the detection of cavities across regular grids or random distances. Such probing is insufficient on its own to provide confidence that cavities will not be encountered. Frequency of probing and depth of investigation may become more expensive. Besides, probing is intrusive, non-continuous, slow, expensive and cannot provide a complete lateral picture of the subsurface geology. Near-surface cavities usually can be easily detected by surface geophysical methods. Traditional and recently developed measuring techniques in seismic, geoelectrics and georadar are suitable for economical investigation of hazardous, potentially collapsing cavities. The presented research focused on an integrated geophysical survey that was carried out in a near-coast road located at Porto Cesareo, a small village a few kilometers south west of Lecce (south Italy). The roads in this area are intensively affected by dangerous surface cracks that cause structural instability. The survey aimed to image the shallow subsurface structures, including karstic features, and evaluate their extent, as they may cause rock instability and lead to cracking of the road. Seismic refraction tomography and ground-penetrating radar surveys were carried out along several parallel traverses extending about 100 m on the cracked road. The acquired data were processed and interpreted integrally to elucidate the shallow structural setting of the site. Integrated interpretation led to the delineation of hazard zones rich with karstic features in the area. Most of these karstic features are associated with vertical and subvertical linear features and cavities. These features are the main reason of the rock instability that resulted in potentially dangerous cracking of road.