探地雷达测量土壤水含量的进展

包括探地雷达(Ground Penetrating Radar)在内的许多电磁方法正越来越多的应用于近地表水文地质领域。探地雷达技术基于介质的电磁特性,通过介质介电常数和速度的变化来反映土壤的体积水含量。水的相对介电常数为81,而周围介质的相对介电常数约为2~30,这样大的介电常数差异使探地雷达成为监测和分析地下水的有效工具。介绍4种应用探地雷达测量土壤水含量的方法,包括:探地雷达反射波法、探地雷达地面波法、钻孔探地雷达方法、以及探地雷达地表反射系数方法4种方法各有不同的适用条件,其中反射波方法最为常用。     

Imaging low-velocity anomalies with the aid of seismic tomography

Theoretical considerations (Snell's law) suggest that low-velocity fanomalies are undersampled and therefore should be poorly resolved by inversion schemes based on ray-tracing methods. A synthetic study considering a 40×20 m low-velocity anomaly (300 m/s) placed at the center of a 400×50 m block with gradient background velocity model (from 3000 m/s at the surface to 4000 m/s at the base) indicates that the low ray density in ray-tracing coverage diagrams of tomographic inversion studies can be used as evidence for the existence of low-velocity anomalies. Combined normal incidence seismic reflection images and the velocity models obtained by tomographic inversions of first-arrival travel times form an efficient scheme to resolve low-velocity anomalies such as fracture zones. Furthermore, the velocity models derived from tomographic inversions are used in a wave equation datuming algorithm to account for statics caused by a strongly laterally variable shallow surface (weathering) layer and provide seismic reflection images of fracture zones (low-velocity anomaly) within a granitic pluton.     

GPR studies over the tsunami affected Karaikal beach,Tamil Nadu,south India

In this study, results of GPR profiling related to mapping of subsurface sedimentary layers at tsunami affected Karaikal beach are presented . A 400 MHz antenna was used for profiling along 262 m stretch of transect from beach to backshore areas with penetration of about 2.0 m depth (50 ns two-way travel time). The velocity analysis was carried out to estimate the depth information along the GPR profile. Based on the significant changes in the reflection amplitude, three different zones are marked and the upper zone is noticed with less moisture compared to other two (saturated) zones. The water table is noticed to vary from 0.5 to 0.75 m depth (12–15 ns) as moving away from the coastline. Buried erosional surface is observed at 1.5 m depth (40–42 ns), which represents the limit up to which the extreme event acted upon. In other words, it is the depth to which the tsunami sediments have been piled up to about 1.5 m thickness. Three field test pits were made along the transect and sedimentary sequences were recorded. The sand layers, especially, heavy mineral layers, recorded in the test pits indicate a positive correlation with the amplitude and velocity changes in the GPR profile. Such interpretation seems to be difficult in the middle zone due to its water saturation condition. But it is fairly clear in the lower zone located just below the erosional surface where the strata is comparatively more compact. The inferences from the GPR profile thus provide a lucid insight to the subsurface sediment sequences of the tsunami sediments in the Karaikal beach.     

三维介质中探地雷达(GPR)波传播逐段迭代射线追踪方法研究和应用

基于F erm at原理和Sne ll定理对三维任意界面情况下的两点间GPR波传播射线追踪问题进行了研究,采用了一个适用于任意界面情况下计算反(折)射点的一阶近似公式,从任意给定的初始路径出发进行逐段迭代计算,给出了三维逐段迭代射线追踪算法的计算格式,并进行了三维介质下射线追踪计算。结果表明计算速度相当快,且其计算精度可以根据需要达到任意要求。在此基础上,根据雷达波传播规律,考虑每一段的能量损失。根据地质雷达的发射子波,最终获得三维介质中雷达波传播和接收记录。并将结果应用于实际剖面拟合中,提高了解释的精度。     

Ground‐penetrating radar and sedimentological analysis of Holocene floodplains: Insight from the Tuross valley,New South Wales

Ground‐penetrating radar (GPR) has been used on an array of floodplain types on the lower Tuross River, in southeastern New South Wales, as part of an investigation into controls on channel‐floodplain relationships. Ground‐penetrating radar transects from two floodplains are presented, along with sedimentological detail from trenches dug along the profiles at key locations. Sedimentological investigations showed that 100 MHz antenna gave an approximation of overall bedding trends in the upper 3 m when automatic gain control processing was used. Spreading and exponential compensation processing provided insight into textural changes associated with increased silt content distal of the levee crest. One trench showed that thinning beds were responsible for onlapping reflectors. Signal attenuation at ～4 m depth below the raised floodplain surface resulted from a >50 cm‐thick bed of sandy clay. The close integration of GPR and sedimentological data produced an excellent dataset, that enabled form‐process associations and floodplain evolution to be established for these sandy floodplains. However, accurate subsurface assessment and interpretation must stem from carefully combined GPR and sedimentological datasets.     

Ground penetrating radar and 2D electric resistivity studies for tracing hydrocarbon leakage site,close to Abha City: a case study

This work presents the application of the ground-penetrating radar (GPR) method and electric resistivity tomography (ERT) technique in outlining a zone of contamination due to the light non-aqueous phase liquid (LNAPL) plume underground in the area of an impacted fuel station, close to Abha City. The GPR has been performed using SIR3000 unit with the 100 and 400 MHz antennas. The main objective of the GPR survey was to evaluate the lateral extension of contamination. The complex GPR signature of the plume was well characterized. Low reflectivity zone corresponds to hydrocarbon vapor phase in the vadose zone. Enhanced reflections are associated with free and residual products in the fractured saturated zone directly above the water table. An electric resistivity tomography (ERT) survey was performed on four profiles within the site to investigate the vertical and horizontal extent of the contamination plume and to define the bottom of the landfill. The 2D electric profiles show the presence of low-resistivity (4O to 37 Ω m) anomalies that refers to the presence of accumulated hydrocarbons. From the interpretation of the GPR and ERT profile, it was possible to locate the top and bottom of the contamination plume of the waste disposal site. The radar signal penetrated deep enough and enabled the identification of a second reflector at approximately 10-m deep, interpreted as the hard basement surface which causes the strong amplitude reflection in the GPR profile. The results of GPR and ERT showed good agreement.     

Geotechnical investigation of the El-Elb dam site,northwest Riyadh,Saudi Arabia,using 2D resistivity and ground-penetrating radar techniques

In the present study, the existence of cavities, voids, and fractures was verified at the site of the El-Elb Dam, which is located to the northwest of Riyadh City across Wadi Hanifa, using 2D electrical resistivity tomography (ERT) and ground-penetrating radar (GPR) techniques. For this purpose, four ERT profiles were measured on the downstream side of the El-Elb Dam using the Syscal Pro Switch-72 resistivity meter. In addition, a GPR survey using a 400-MHz antenna and a SIR-3000 instrument was conducted along five profiles above the stilling basins on the downstream side of the dam and one radar profile was measured outside the stilling basins area across the course of the wadi. The resultant geophysical data were interpreted with the aid of information from a field-based structural and stratigraphic evaluation of the outcropped bedrock on the banks of the wadi course. The analysis of the inverted ERT and filtered radar sections revealed several resistivity and electromagnetic reflection anomalies that are identified laterally and vertically across the measured sections. These anomalies indicate the presence of fractures and karst features affected the limestone bedrock in the dam site. These near-surface karstified and fractured strata represent a critical hazard to the structural safety of the dam.     

基于自然电位方法的土壤水分入渗过程监测

自然电位方法可以应用于水文地球物理领域获取土壤水分入渗特征.通过连续监测野外试验场的自然电位场,获取连续的自然电位平面图,直接获知了土壤水分在平面上的主要入渗区域,同时借助自然电位成像方法,揭示了土壤水分的两个主要通道,并探测出两个地下孔洞的位置,最终在地质雷达剖面图的辅助下证实了上述结论,为水文地球物理领域提供了新的技术手段和方法.     

The application of GPR and electrical resistivity tomography as useful tools in detection of sinkholes in the Cheria Basin (northeast of Algeria)

Sinkhole collapse is one of the main limitations in the development of karst areas, especially where bedrock is covered by unconsolidated material. Studies of sinkhole formation have shown that sinkholes are likely to develop in cutter (enlarged joint) zones as a result of subterranean erosion by flowing groundwater. Ground-penetrating radar (GPR) and electrical resistivity imaging or tomography (RESTOM) are well suited to mapping sinkholes because of the ability of these two techniques for detecting voids and discriminating subtle resistivity variations. Nine GPR profiles and two-dimensional electrical resistivity tomography have been applied, with relative success, to locate paleo-collapses and cavities, and to detect and characterize karst at two sinkhole sites near Cheria City where limestone is covered by about 10 m of clayey soils. The survey results suggest that GPR and RESTOM are ideal geophysical tools to aid in the detection and monitoring of sinkholes and other subsurface cavities.     

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground-penetrating radar

Across the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large-scale storms, while providing important habitats and recreational opportunities. Understanding their long-term evolution is vital in managing their condition in a changing climate. Recently, ground-penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre-scale resolution images of near-surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo-environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground-truth control. Nonetheless, GPR data interpretation can be non-intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi-automated GPR feature-extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo-landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty-organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR-interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.     

A GPR study of subsidence-creep-fault processes in Morelia, Michoacán, Mexico

In 1983, inhabitants of the City of Morelia, Michoacán, Mexico, began to observe a series of differential settlements causing damages to constructions along linear trends parallel to a system of regional faults. The same phenomenon occurs in others cities of the Mexican Volcanic Belt (MVB), such as Celaya, Aguascalientes, and Querétaro, and is linked to a structurally controlled subsidence, caused by groundwater withdrawal, and the presence of geological faults. We define this subsidence type as Subsidence-Creep-Fault Processes (SCFP), based on the necessary elements for their generation, and we studied them through geophysical and geotechnical techniques. In Morelia, the geophysical investigations have been carried out using ground-penetrating radar (GPR). GPR profiles, perpendicular to the axis of the surface fault generated by the SCFP were carried out. The common-offset single-fold profiling was used, with a central frequency of 50 MHz. In all cases it has been possible to visualize a fault plane dividing two blocks, the presence of synthetic and antithetic faults, influence zones from 20 m to 40 m, and a maximum “net throw” of 4 m. Exploration trenches followed the same direction of the profiles obtained with GPR (perpendicular to the axis of the surface fault). These trenches exposed a fault plane dividing two blocks with different lithology, generating a maximum “net throw” of 4.40 m; as well they help in the determination of influence zones that varied from 14 m to 40 m.     

探地雷达宽角反射图形拟合方法

探地雷达共中心点(CMP)法仅用于局部速度测试,为使用宽角反射法进行剖面探测,可固定一个雷达天线而将另一个天线沿测线移动,从而获取整个地质界面的反射信号.本文提出宽角反射剖面探测方法,基于平直岩层界面导出其反射信号的关系式,并开发出解释软件,可精确计算岩层几何参数和层速度,从而进行地质分层.     

Identification of cross-valley faults in the Maynardville Limestone,Oak Ridge Reservation,Tennessee, using seismic refraction tomography

First arrival times from P-wave refraction and reflection seismic surveys along Bear Creek Valley on the Oak Ridge Reservation, Tennessee, were inverted to produce refraction tomographic velocity images showing seismic velocity variations within thinly mantled karstic bedrock to a depth of approximately 20 m. Inverted velocities are consistent with two distinct bedrock groups: the Nolichucky Shale (2,730–5,150 m/s) and Maynardville Limestone (3,940–7,575 m/s). Low-velocity zones (2,700–4,000 m/s) in the tomographic images correspond to previously inferred cross-valley strike-slip faults; in places, these faults create permeability barriers that offset or block groundwater flowing along Bear Creek Valley. These faults may also force groundwater contaminants, such as dense non-aqueous phase liquids, to migrate laterally or downward, spreading contamination throughout the groundwater system. Other, previously unmapped cross-valley faults may also be visible in the tomographic images. Borehole logs suggest the low-velocity values are caused by low rigidity fractured and vuggy rock, water zones, cavities and collapse features. Surface streams, including Bear Creek, tend to lie directly above these low-velocity zones, suggesting fault and fracture control of surface drainage, in addition to the subsurface flow system. In some cases, fault zones are also associated with bedrock depressions and thicker accumulations of unconsolidated sediment.     

Geophysical investigation using resistivity and GPR methods: a case study of a lubricant oil waste disposal area in the city of Ribeirão Preto,São Paulo,Brazil

Geophysics has been shown to be effective in identifying areas contaminated by waste disposal, contributing to the greater efficiency of soundings programs and the installation of monitoring wells. In the study area, four trenches were constructed with a total volume of about 25,000 m³. They were almost totally filled with re-refined lubricating oil waste for approximately 25 years. No protection liners were used in the bottoms and laterals of the disposal trenches. The purpose of this work is to evaluate the potential of the resistivity and ground penetrating radar (GPR) methods in characterizing the contamination of this lubricant oil waste disposal area in Ribeirão Preto, SP, situated on the geological domain of the basalt spills of the Serra Geral Formation and the sandstones of the Botucatu Formation. Geophysical results were shown in 2D profiles. The geophysical methods used enabled the identification of geophysical anomalies, which characterized the contamination produced by the trenches filled with lubricant oil waste. Conductive anomalies (smaller than 185 Ωm) immediately below the trenches suggest the action of bacteria in the hydrocarbons, as has been observed in several sites contaminated by hydrocarbons in previously reported cases in the literature. It was also possible to define the geometry of the trenches, as evidenced by the GPR method. Direct sampling (chemical analysis of the soil and the water in the monitoring well) confirmed the contamination. In the soil analysis, low concentrations of several polycyclic aromatic hydrocarbons (PAHs) were found, mainly naphthalene and phenanthrene. In the water samples, an analysis verified contamination of the groundwater by lead (Pb). The geophysical methods used in the investigation provided an excellent tool for environmental characterization in this study of a lubricant oil waste disposal area, and could be applied in the study of similar areas.     

井间地震层析成像分辨率研究

通过大量数值模型试验,并根据射线密度、射线正交性和点扩散函数,对井间地震层析成像的分辨率问题进行了研究。结果表明,井间地震层析成像的分辨率在空间上是可能性变的,垂直方向的分辨率好于水平方向,而且高速异常体和低速异常体的图像形态和分辨率是非一致的。高速异常体波形CT的图像水平方向最小分辨率距离为1/2波长,垂直分辨率最小分辨率距离为3/10波长;而低速异常体的波形CT图像的水平方向的最小分辨率约为2/5波长,垂向最小分辨率约1/5个波长;走时CT的图像水平方向最小分辨率为3个像元,垂直方向的最小分辨率距离为2个像元;当井深与井间距之比为2时,图像分辨率最高。     

探地雷达测量近地表含水量的研究

为充分验证探地雷达方法用于测量近地表含水量的有效性和可靠程度,基于探地雷达的反射波法,根据介电常数、电磁波传播速度和含水量三者之间的关系,着重分析了探地雷达反射波法用于测量含水量的相关理论和特点,设计了适用于高频雷达信号的速度分析算法,建立了用于分析该问题的局部含水层模型。利用FDTD方法模拟分析了不同测量方式下,局部含水层的探地雷达响应,借助速度分析理论及Topp公式,反演得到了研究区域的含水量。     

Ground penetrating radar imaging and time-domain modelling of the infiltration of diesel fuel in a sandbox experiment

Ground penetrating radar (GPR) is a non-destructive method which, over the past 10 years, has been successfully used not only to estimate the water content of soil, but also to detect and monitor the infiltration of pollutants on sites contaminated by light non-aqueous phase liquids (LNAPL). We represented a model water table aquifer (72 cm depth) by injecting water into a sandbox that also contains several buried objects. The GPR measurements were carried out with shielded antennae of 900 and 1200 MHz, respectively, for common mid point (CMP) and constant offset (CO) profiles. We extended the work reported by Loeffler and Bano by injecting 100 L of diesel fuel (LNAPL) from the top of the sandbox. We used the same acquisition procedure and the same profile configuration as before fuel injection. The GPR data acquired on the polluted sand did not show any clear reflections from the plume pollution; nevertheless, travel times are very strongly affected by the presence of the fuel and the main changes are on the velocity anomalies. We can notice that the reflection from the bottom of the sandbox, which is recorded at a constant time when no fuel is present, is deformed by the pollution. The area close to the fuel injection point is characterized by a higher velocity than the area situated further away. The area farther away from the injection point shows a low velocity anomaly which indicates an increase in travel time. It seems that pore water has been replaced by fuel as a result of a lateral flow. We also use finite-difference time-domain (FDTD) numerical GPR modelling in combination with dielectric property mixing models to estimate the volume and the physical characteristics of the contaminated sand.     

Combined ground-penetrating radar (GPR) and electrical resistivity applications exploring groundwater potential zones in granitic terrain

Frequent failures of monsoons have forced to opt the groundwater as the only source of irrigation in non-command areas. Groundwater exploration in granitic terrain of dry land agriculture has been a major concern for farmers and water resource authorities. The hydrogeological complexities and lack of understanding of the aquifer systems have resulted in the failure of a majority of the borehole drillings in India. Hence, a combination of geophysical tools comprising ground-penetrating radar (GPR), multielectrode resistivity imaging (MERI), and vertical electrical sounding (VES) has been employed for pinpointing the groundwater potential zones in dry land agricultural of granitic terrain in India. Results obtained and verified with each other led to the detection of a saturated fracture within the environs. In GPR scanning, a 40-MHz antenna is used with specifications of 5 dielectric constant, 600 scans/nS, and 40 m depth. The anomalies acquired on GPR scans at various depths are confirmed with low-resistivity ranges of 27–50?Ω m at 23 and 27 m depths obtained from the MERI. Further, drilling with a down-the-hole hammer was carried out at two recommended sites down to 50–70 m depth, which were complimentary of VES results. The integrated geophysical anomalies have good agreement with the drilling lithologs validating the MERI and GPR data. The yields of these bore wells varied from 83 to 130 l/min. This approach is possible and can be replicated by water resource authorities in thrust areas of dry land environs of hard rock terrain around the world.     

电阻加载蝶形天线的性能研究

由于探地雷达系统大多采用超宽带窄脉冲信号,因此要求天线具有极宽的工作带宽。利用矩量法对线性离散电阻加载蝶形天线进行分析,结果表明,线性离散电阻加载在使天线获得宽带特性的同时,并使辐射的脉冲信号拖尾变得平稳,因此能够满足探地雷达系统的需要。