Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.     

地球物理信号处理技术的研究及进展

从数字信号处理的原理、应用、发展过程及前景等方面对数字信号处理技术作了概述，重点总结了应用地球物理方法中和的信号处理方法和技术。对探地雷达技术中常用的数据处理技术作了详细阐述，并简要介绍了作者研制的信号处理软件。     

地质雷达在公路滑坡加固治理中的应用

地质雷达能快速、准确地探测地下潜在的地质缺陷体的空间位置、形状大小及电性质等其他一些地质参数,可以为滑坡的加固设计、施工提供可靠的资料,简述了地质雷达的勘探原理和方法,并以实例说明了地质雷达在滑坡治理中的应用方法及效果。     

Nonlinear data processing method for the signal enhancement of GPR data

An alternative data processing procedure is proposed in this paper for the purpose of enhancing the signal/noise (S/N) ratio of ground penetrating radar (GPR) data. The processing methodology is achieved by performing the logarithmic transform in conjunction with the ensemble empirical mode decomposition (EEMD), a new nonlinear data analysis method in signal processing. The synthetic model study and field example indicate that the logarithmic transform is effective in alleviating the attenuation problem. Additionally, the spectrogram obtained from Hilbert-Huang transform (HHT) shows that the decomposition sensitivity of the EEMD method is greatly improved with the aid of the logarithmic transform. This new method allows us to extract the signal components from noisy GPR data efficiently. The success of this study suggests a possible nonlinear analysis application in future GPR investigation, particularly in the filter design and gain correction.     

The use of GPR and VES in delineating a contamination plume in a landfill site: a case study in SE Brazil

This paper presents the results of the application of the Ground Penetrating Radar (GPR) method, or Georadar, in outlining a zone of contamination due to solid residues at the waste burial site of Rio Claro in the state of São Paulo, SE Brazil. A total of eight GPR profiles with 50- and 100-MHz antennae were surveyed. Six profiles were located within the landfill site and the remaining two were outside. The main objective of the GPR survey was to evaluate the side extension of contamination. A Vertical Electric Sounding (VES) survey was performed at four points within the site in order to investigate the vertical extent of the contamination plume and to define the bottom of the landfill. Two additional VESs were done outside the landfill with the purpose of determining the top of the ground water table and the geoelectric stratigraphy of the background. From the interpretation of the GPR profiles, it was possible to locate the top of the contamination plume and to infer that it was migrating laterally beyond the limits of the waste disposal site. This was observed along the profile situated close to the highway SP-127, which was about 20 m from the limit of the site. The signature of the contaminant appears as a discontinuous reflector that is believed to be a shallow ground water table. The discontinuity is marked by a shadow zone, which is characteristic of conductive contaminant residues. The contamination did not move far enough to reach a sugar cane plantation located at approximately 100 m from the border of the site. In the regions free from contamination, the ground water table was mapped at approximately 10 m of depth, and it was characterized by a strong and continuous reflector. The radar signal penetrated deep enough and enabled the identification of a second reflector at approximately 14 m deep, interpreted as the contact between the Rio Claro and the Corumbataí formations. The contact is marked by the presence of gravel characterized by ferruginous concretes, which cause the strong amplitude reflection in the GPR profile. Within the landfill site, the quantitative interpretation of the VES results showed the contamination zone. The base of the landfill varies between 11 and 15 m deep. Outside the landfill site, the VES results showed no indication of contamination and allowed the determination of the top of the ground water table and the contact between the Rio Claro and the Corumbataí formations. The results of GPR and VES showed a good agreement and the integrated interpretations were supported by local geology and information from several boreholes, about 17 m depth, on average. The bottom of the landfill reaches a maximum of 14.5 m depth.     

Changes of GPR spectra due to the presence of hydrocarbon contamination in the ground

The location of hydrocarbon contamination in the ground using the GPR method is based mainly on information taken from reflected signals. In the cases investigated in Polish contaminated sites, such signals were very seldom recorded. A long time after spillage, contamination takes the form of plumes with different size and distribution, which depends on geological and hydraulic properties of the ground. In this paper, it is shown that the set of hydrocarbon plumes should be described with a stochastic model, and such plumes may generate the scattered waves which cause changes in the power spectra. It has been observed that the power spectra of GPR signals over contaminated areas are quite different from such spectra over clear ones. These differences were discussed in this paper on the basis of theoretical analysis, numerical modelling and the results of GPR terrain surveys.     

基于速度估计的改进型线性变换有限差分偏移在探地雷达中的应用

为了解决传统线性变换有限差分偏移算法对变速介质的高敏感性问题,对加入了速度变量的标量波动方程进行线性变换,在横向非均匀介质条件下,详细推导了包含混合因子项的线性变换差分偏移方程,实现了横向非均匀介质中陡倾角反射界面的正确归位,同时改进型算法可以灵活地选择混合因子项,这样就增强了原算法对复杂地质问题的适应能力,在此基础上结合速度估计算法,减小多次尝试输入法所产生的误差,提出一种基于速度估计的改进型线性变换有限差分偏移算法,并用该算法对实测探地雷达剖面进行了偏移处理,通过对比偏移前后的雷达资料得知,当横向非均匀介质中同时存在多次反射质点和主体反射层时,改进的偏移算法能有效地提高雷达剖面的偏移精度.     

The effect of roughness configuration on velocity profiles in an artificial channel

In order to determine the effect of bed roughness on velocity distribution, we used seven different configurations of bed roughness, with 16 test runs of varying discharge and slope for each configuration. For each run, one-dimensional velocity profiles were measured at 1 cm vertical increments over the crest of the roughness element, and at intervals of 4·25 cm downstream. Results indicate that velocity profile shape remains fairly constant for a given slope and roughness configuration as discharge increases. As slope increases, the profiles become less linear, with a much larger near-bed velocity gradient and a more pronounced velocity peak close to 0·6 flow depth at the measurement point immediately downstream from the roughness element. The zone of large near-bed velocity gradients increases in both length and depth as roughness concentration decreases, up to a length/height ratio of about 9, at which point maximum flow resistance occurs. Longitudinal roughness elements do not create nearly as much flow resistance as do transverse elements. Rates of velocity increase suggest that roughness elements spaced at a length/height ratio of about 9 are most effective at creating flow resistance over a range of discharges in channels with steeper slopes. © 1998 John Wiley & Sons, Ltd.     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

高密度电法在地下防渗墙检测中的应用

对高密度电法在防渗墙无损检测中的应用方法进行了对比研究,通过实际工程应用结果,对高密度电法四极削面及三极剖面两种方法进行对比研究;对于符合检测条件的防渗墙,采用三极剖面法,可以得到比较满意的检测效果;此方法快速、低耗、成果直观、灵敏度高、适应性强,便于工程实际应用,具有广阔的应用前景和推广价值。     

基于Pade逼近的Cole-Cole频散介质GPR有限元正演

针对Cole-Cole频散介质中的复介电常数是jω的分数次幂函数,传统的时域有限元法难以离散及计算时间域分数阶导数,本文采用Pade逼近算法将含有时间分数阶导数的Cole-Cole频散介质电磁波方程推导为一组整数阶辅助微分方程,提出了一种适用于Cole-Cole频散介质的GPR有限元正演模拟算法.在复数伸展坐标系下,通过在频率域Cole-Cole频散介质电磁波方程中引入2个中间变量,并将其变换到时间域,从而以变分形式将PML边界条件加载到Cole-Cole频散介质GPR有限元方程组中,并给出了详细的求解公式.在此基础上,编制了基于Pade逼近的Cole-Cole频散介质GPR有限元正演程序,利用该程序对均匀模型进行计算,并与解析解进行对比,验证了本文构建的GPR有限元正演算法的正确性和有效性.设计了一个复杂Cole-Cole频散介质GPR模型,利用本文构建的GPR有限元正演算法进行模拟并与非频散介质模型的模拟结果进行对比,分析了电磁波在Cole-Cole频散介质中传播衰减增强、子波延伸,分辨率降低等传播规律,有助于实测雷达资料更可靠、更准确的解释.模拟结果表明,基于Pade逼近的GPR有限元正演算法可用于复杂Cole-Cole频散介质结构模拟,且具有较高的计算精度.

     

The role of GPR techniques in determining ice cave properties: Peña Castil ice cave,Picos de Europa

The structure and ice content of ice caves are poorly understood. Ground penetrating radar (GPR) can provide useful insights but has only rarely been applied to ice caves. This paper interprets GPR images (radargrams) in terms of internal structure, stratification, compaction, thickness and volume of the ice block in the Peña Castil ice cave (Central Massif of Picos de Europa, northern Spain), providing the endokarst geometry of the ice cave in GPR data reflections. Eight radargrams were obtained by applying a shielded ground‐coupled antenna with a nominal frequency of 400 MHz. Although the radargrams do not depict the ice–basal bedrock interface, they suggest that the ice block is at least 54 m deep and similarly thick. Some curved reflection signatures suggest a potential vertical displacement in the block of ice, and thus certain dynamics in the ice body. Other images show numerous interbedded clasts and thin sediment layers imaged as banded reflections. In this particular cave a direct visual inspection of the ice stratigraphy is a difficult task but GPR provides clear reflectivity patterns of some of its internal features, making GPR a suitable instrument for this and future studies to achieve a better and broader understanding of the internal behavior of ice caves.     

曲靖非相干散射雷达功率剖面的初步观测与分析

介绍了非相干散射雷达功率剖面与电离层电子密度剖面之间的关系,给出了曲靖非相干散射雷达功率剖面的一些观测结果,结合电离层垂直探测数据,初步介绍了曲靖非相干散射雷达在夜间电离层波状扰动、低纬电离层异常区北驼峰位置、电离层电子密度日落增强和暴时变化等观测研究中的应用.本文结果丰富了对我国曲靖地区电离层空间天气特性的认识,显示了该雷达在我国低纬电离层空间天气观测研究中的良好应用前景.     

The use of regional and alluvial groundwater by riparian trees in the wet-dry tropics of northern Australia

Riparian trees play a critical role in the ecological function of rivers, yet are threatened by anthropogenic change to the hydrological cycle. Identifying the sources of water used by riparian trees can inform sustainable water policy. We used isotopic analysis complemented by measurements of plant water relations to assess water sources for riparian trees at two sites with contrasting hydrogeological processes; one with an alluvial aquifer overlaying an aquitard, and one where fault-induced preferential pathways in the aquitard allowed the flow of deeper, older groundwater from a regional aquifer to the alluvium. At both sites, plant water potential, stomatal conductance, and plant water isotope composition in the xylem sap of riparian trees were collected from two landscape positions, the riverbank and floodplain. We used a Bayesian mixing model (MixSIAR) to assess differences in the proportion of water sources for sites and landscape positions. We found that xylem water isotope values differed between the two sites in line with their hydrogeological characteristics, with trees at the regional aquifer site using water sourced from the regional groundwater and trees at the site with only an alluvial aquifer present using a mixture of water sources, with no dominant source identified. Higher plant predawn water potential values at the regional site indicated greater water availability and support the inference that plants were using more groundwater at the regional site compared to the alluvial site. Trees closer to the river had higher isotope values, indicative of surficial water sources i.e. shallow soil water and river water. Our findings show that the water sources used by riparian trees reflect local hydrogeology and resource availability. Water managers should identify and protect plant water sources to ensure maintenance of riparian trees.     

The use of mixing models and chemical geothermometers for estimating underground temperatures in geothermal systems

Application of various chemical geothermometers and mixing models indicate underground temperatures of 260°C, 280°C and 265°C in the Geysir, Hveravellir and Landmannalaugar geothermal fields in Iceland, respectively. Mixing of the hot water with cold water occurs in the upflow zones of all these geothermal systems. Linear relations between chloride, boron and δ¹⁸O constitute the main evidence for mixing, which is further substantiated by chloride, silica and sulphate relations in the Geysir and Hveravellir fields.A new carbonate-silica mixing model is proposed which is useful in distinguishing boiled and non-boiled geothermal waters. This model can also be used to estimate underground temperatures using data from warm springs. This model, as well as the chloride-enthalpy model and the Na-Li, and CO₂-gas geothermometers, invariably yield similar results as the quartz geothermometer sometimes also does. By contrast, the Na-K and the Na-K-Ca geothermometers yield low values in the case of boiling hot springs, largely due to loss of potassium from solution in the upflow. The results of these geothermometers are unreliable for mixed waters due to leaching subsequent to mixing.     

我国市县地球物理监测资源评估

我国地球物理监测站网主要由GNSS、定点形变、重力、地电、地磁、地下流体六大观测站网组成,由国家站、省级站和市县站构成基本监测单元。市县站观测手段多、站点数量庞大,是地球物理监测站网的重要组成部分。为摸清市县站监测资源底数、统筹利用市县监测数据,中国地震局历时2年完成了我国市县地震监测资源评估工作。基于此项工作,本文从观测环境、观测系统、数据质量和预报应用四个方面对我国市县地球物理监测资源进行简要分析,评估结果显示: 37.8%的测项评估为合格,可直接纳入国家地球物理数据库管理,47.9%的测项评估为基本合格,可纳入国家市县数据库管理,14.3%的测项评估为不合格,需要停测或者优化改造。针对评估结果,简要分析了市县监测运维工作存在的主要问题并提出具体可行的解决建议,为提升市县站监测质量提供参考。     

The use of seismic reflection profiles to probe the crystalline basement in Britain

In an attempt to assess the feasibility of a future programme of seismic reflection profiling to investigate major structural features deep within the crust of Britain, existing commercial offshore data across shallow basement have been investigated together with lines on land, fired conventionally with explosives, in which record lengths have been extended to at least 12 s two-way time. Profiles crossing Lewisian basement, a granite batholith, and Caledonian and Variscan structures have been examined. The seismic information has been combined with gravity and magnetic surveys, bore-hole logs and surface geology to give self-consistent interpretational models. Features such as intrusive bodies and thrust faults can be discerned within the basement. These are clearly visible on the seismic record section obtained along the MOIST line, fired along the north coast of Scotland especially for deep reflections, and, for the first time, normal-incidence reflections from the Moho are clearly observed.     

On the use of resistance envelopes to identify the controls on slope stability in the tropics

In tropical and sub-tropical slopes, soil suction may in certain circumstances play a role in maintaining slope stability. Resistance envelope methods are outlined that provide a means of assessing the threshold soil water conditions for stability. In addition, this technique enables the likely failure depth to be identified. Application to slopes in St. Lucia, West Indies, show the accordance of resistance envelope predictions with stability analysis results. A methodology for the geomorphological investigation of stability processes in tropical slopes is proposed.     

北京边界层大气污染物的垂直廓线监测与分析

基于扫描差分光学吸收光谱（DOAS）系统,于2007年8月27日~9月4日期间对北京市朝阳区大气污染物SO2,HCHO,O3和NO2的垂直分布进行了连续监测,并对污染物垂直廓线进行了分析,详细探讨了SO2垂直分布特征以及夜间NO对O3的滴定作用.结果表明,SO2浓度通常没有明显垂直分层分布特征,但在清晨风速较低时呈现负梯度变化.研究发现夜间O3、NO和NO2之间存在稳态作用,表明城市区域O3的滴定主要来自于地面NO的直接排放,显示出NO的滴定作用对夜间O3、NO和NO2的垂直变化起到重要作用.