Using 2D inversion of magnetic resonance soundings to locate a water-filled karst conduit

A new methodology for magnetic resonance sounding (MRS) data acquisition and interpretation was developed for locating water-filled karst cavities. This methodology was used to investigate the Ouysse karst system in the Poumeyssens shaft in the Causse de Gramat (France). A new 2D numerical MRS response model was designed for improved accuracy over the previous 1D MRS approach. A special survey performed by cave divers confirmed the accuracy of the MRS results. Field results demonstrated that in favourable conditions (a low EM noise environment and a relatively shallow, large target) the MRS method, used with a coincident transmitter/receiver loop, can be an effective tool for locating a water-filled karst conduit. It was shown numerically that because an a priori orientation of the MRS profile with the karst conduit is used in the inversion scheme (perpendicular for instance), any error in this assumption introduces an additional error in locating the karst. However, the resulting error is within acceptable limits when the deviation is less than 30°. The MRS results were compared with an electrical resistivity tomography (ERT) survey. It was found that in Poumeyssens, ERT is not able to locate the water-filled karst. On the other hand, ERT provides additional information about heterogeneities in the limestone.     

Mapping of quick clay by electrical resistivity tomography under structural constraint

Geotechnical projects usually rely on traditional sounding and drilling investigations. Drilling only provides point information and the geology needs to be interpolated between these points. Near surface geophysical methods can provide information to fill those gaps. Norwegian case studies are presented to illustrate how two-dimensional electrical resistivity tomography (ERT) can be used to accurately map the extent of quick clay deposits. Quick clay may be described as highly sensitive marine clay that changes from a relatively stiff condition to a liquid mass when disturbed. Quick clay slides present a geo-hazard and therefore layers of sensitive clay need to be mapped in detail. They are usually characterized by higher resistivity than non-sensitive clay and ERT is therefore a suitable approach to identify their occurrence. However, our experience shows that ERT cannot resolve this small resistivity contrast near large anomalies such as a bedrock interface. For this reason, a constrained inversion of ERT data was applied to delineate quick clay extent both vertically and laterally. As compared to the conventional unconstrained inversions, the constrained inversion models exhibit sharper resistivity contrasts and their resistivity values agree better with in situ measurements.     

Two-dimensional, regularised inversion of VLF data

Very low frequency electromagnetic (EM) methods using VLF transmitters have found many applications in subsurface geophysical investigations. Surface measurements involving both the vertical component of the magnetic field (VLF-EM or VLF-Z) and of the apparent resistivity (VLF-R) are increasingly common. Although extensive VLF data sets have been successfully used for mapping purposes, modelling and interpretation techniques which asess the third (i.e. depth) dimension appear limited.Given a profile of VLF-R measurements the main purpose of the present study is to demonstrate an automatic method for the construction of a resistivity cross-section. The technique used is one of a new generation of regularised inversion methods. These techniques attempt to overcome the problem of equivalence/non-uniqueness in EM sounding data by constructing the resistivity distribution with the minimum amount of structure that fits the data.VLF data represent a special case of plane-wave EM sounding in that they conform, in practice, to a single-frequency technique. This fact imposes a limitation in the amount of vertical resolution that we can expect using such data. In the case of two-dimensional modelling and inversion, resolution through the cross-section is a resultant attribute from both vertical and lateral resistivity gradients within the subsurface. In order to provide insight into the practical application of regularised inversion techniques to VLF data, both synthetic and field examples are considered. Both sets of examples are primarily concerned with VLF data applied to near-surface fault mapping where the main aim is to assess the location, dip and depth extent of conductive subsurface features.     

Marine Electromagnetic Studies of Seafloor Resources and Tectonics

The past decade has been a period of rapid growth for marine electromagnetic (EM) methods, predominantly due to the industrial adoption and promotion of EM as a valuable tool for characterizing offshore hydrocarbon reservoirs. This growth is illustrated by a database of marine EM publications spanning from the early developments in the 1960’s to the present day; while over 300 peer-reviewed papers on marine EM have been published to date, more than half of these papers have been published within the last decade. This review provides an overview of these recent developments, covering industrial and academic use of marine EM for resource exploration and tectonic investigations, ranging from acquisition technology and modeling approaches to new physical and geological insights learned from recent data sets.     

Spatially distributed observations in constraining inundation modelling uncertainties

The performance of two modelling approaches for predicting floodplain inundation is tested using observed flood extent and 26 distributed floodplain level observations for the 1997 flood event in the town of Usti nad Orlici in the Czech Republic. Although the one‐dimensional hydrodynamic model and the integrated one‐ and two‐dimensional model are shown to perform comparably against the flood extent data, the latter shows better performance against the distributed level observations. Comparable performance in predicting the extent of inundation is found to be primarily as a result of the urban reach considered, with flood extent constrained by road and railway embankments. Uncertainty in the elevation model used in both approaches is shown to have little effect on the reliability in predicting flood extent, with a greater impact on the ability in predicting the distributed level observations. These results show that reliability of flood inundation modelling in urban reaches, where flood risk assessment is of more interest than in more rural reaches, can be improved greatly if distributed observations of levels in the floodplain are used in constraining model uncertainties. Copyright © 2005 John Wiley & Sons, Ltd.     

Contribution of the geophysical methods in characterizing the water leakage in Afamia B dam, Syria

The objective of this work was to characterize and understand the origin of the water losses problem in the Afamia B dam located at Al-Ghab basin in Syria. The survey involved various geophysical methods including a superficial electromagnetic (EM), an electrical sounding (ES) and electrical resistivity tomography (ERT). Special considerations were focused on the lake of the dam, which is currently suffering serious water infiltration and leakage through its bedrocks. The application of the ERT technique was more effective and convenient. The analysis of the integrated results of the above mentioned methods revealed a vertical serious leakage which could take place in certain locations via fractures and faults that hit the main valley and pass through the dam lake. Additionally, the presence of an alternating lithological heterogeneity between permeable and impermeable layers may lead to infiltration through the geological formations of the basin. These processes are most likely causing hydraulic connections between the Neogene superficial deposits and the underlying Cretaceous fractured and karstified carbonates rocks. Consequently, such hydraulic connections may lead to water losses and leakage throughout the dam basin. Eventually, the presence of tens of drilled wells which penetrate the Neogene and Cretaceous formations, and the existence of numerous archeological graves and subsurface tunnels complicate the situation in the study area and cause more water leakage.     

Geophysical study of the Peinan Archaeological Site,Taiwan

The Peinan archaeological site is the most intact Neolithic village with slate coffin burial complexes in Taiwan. However, the area that potentially contains significant ancient remains is covered by dense vegetation. No reliable data show the distribution of the ancient village, and no geophysical investigation has been performed at this site. To evaluate various geophysical methods under the geological setting and surface condition of the site, the physical properties of the remains were measured and four geophysical methods involving magnetic, electromagnetic (EM), electrical resistivity tomography (ERT), and ground-penetrating radar (GPR) were tested along three parallel profiles. The results imply that the EM and magnetic methods are much cost-effective and suitable for investigating the entire area. GPR and ERT methods can provide high resolution subsurface image, which are much suitable for subsequently detail investigation.The EM and magnetic surveys were thus conducted over the entire Peinan Cultural Park to understand the distribution of the ancient building remains at the Peinan site. The results of this study were verified by subsequent excavations, which indicate that the EM survey was successful in delineating the majority of the ancient village because the basements of building are highly resistive in comparison to the background sediment. The results of this investigation suggest that the ancient village was broadly distributed over the eastern part of the Peinan Culture Park and extended to the southeast.     

基于Matlab平台GUI的地震走时层析成像快速实现

地震走时层析成像是地球物理反演中的成熟方法之一,已在许多领域得到广泛应用,并取得了良好的效果。本文围绕地震走时层析成像GUI的快速实现,首先简要介绍了地震走时成像方法的原理和方法,并基于Matlab软件平台,利用GUI开发环境编写了二维横向非均匀地震走时层析成像图形用户界面（GUI）,实现了方便快捷地反演地震测深资料的速度结构并以人机交互方式快捷灵活地将图像显示出来;结合实际观测资料处理结果表明,本文编制的GUI软件能够有效地实现地震走时层析成像结果的窗口化、图像化和快捷化。     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

River–floodplain hydrology of an embanked lowland Chalk river and initial response to embankment removal

Abstract

Rivers have been channelized, deepened and constrained by embankments for centuries to increase agricultural productivity and improve flood defences. This has decreased the hydrological connectivity between rivers and their floodplains. We quantified the hydrological regime of a wet grassland meadow prior to and after the removal of river embankments. River and groundwater chemistry were also monitored to examine hydrological controls on floodplain nutrient status. Prior to restoration, the highest river flows (～2 m³ s^?1) were retained by the embankments. Under these flow conditions the usual hydraulic gradient from the floodplain to the river was reversed so that subsurface flows were directed towards the floodplain. Groundwater was depleted in dissolved oxygen (mean: 0.6 mg O₂ L^?1) and nitrate (mean: 0.5 mg NO₃ ^?-N L^?1) relative to river water (mean: 10.8 mg O₂ L^?1 and 6.2 mg NO₃ ^?-N L^?1, respectively). Removal of the embankments has reduced the channel capacity by an average of 60%. This has facilitated over-bank flow which is likely to favour conditions for improved flood storage and removal of river nutrients by floodplain sediments.

Editor Z.W. Kundzewicz; Associate editor K. Heal

Citation Clilverd, H.M., Thompson, J.R., Heppell, C.M., Sayer, C.D., and Axmacher, J.C., 2013. River–floodplain hydrology of an embanked lowland Chalk river and initial response to embankment removal. Hydrological Sciences Journal, 58 (3), 627–650.     

Mining,Environmental, Petroleum,and Engineering Industry Applications of Electromagnetic Techniques in Geophysics

Electromagnetic (EM) techniques are extremely important as a direct detection geophysical tool utilized in the base metal industry. They were developed in countries such as Canada, whose thin conductive weathering overburden did not hamper the penetration of EM signals and enabled exploration to depths on the order of 300 m. As a result, EM techniques were used widely in North America and Scandinavia for many years before they became common in countries with a thick conductive overburden, such as Australia. The 1980s and 1990s have seen the use of EM methods move from anomaly finding to mapping, as well as the development of better, faster and more accurate computer modelling algorithms. A review of EM papers, for the years 1998 to 2002, showed that most dealt with EM techniques as mapping tools. Airborne, ground and marine EM techniques are still being developed, as are data processing and interpretation software. The advent of robust 2-D and 3-D computer modelling and inversion algorithms has led to the acceptance of EM methods as a mapping tool for many environmental and petroleum industry applications, a trend which is expected to increase.     

Comparing efficacy of different taxonomic resolutions and surrogates in detecting changes in soft bottom assemblages due to coastal defence structures

Sandy shores on the West coast of the North Adriatic Sea are extensively protected by different types of defence structures to prevent coastal erosion. Coastal defence schemes modify the hydrodynamic regime, the sediment structure and composition thus affecting the benthic assemblages. This study examines the effectiveness in detecting changes in soft bottom assemblages caused by coastal defence structures by using different levels of taxonomic resolution, polychaetes and/or bivalves as surrogates and different data transformations. A synoptic analyses of three datasets of subtidal benthic macrofauna used in studies aimed at assessing the impact of breakwaters along the North Adriatic coast has been done. Analyses of similarities and correlations between distance matrices were done using matrices with different levels of taxonomic resolution, and with polychaetes or bivalves data alone. Lentidium mediterraneum was the most abundant species in all datasets. Its abundance was not consistently related to the presence of defence structures. Moreover, distribution patterns of L. mediterraneum were masking the structure of the whole macrofaunal assemblages. Removal of L. mediterraneum from the datasets allowed the detection of changes in benthic assemblages due to coastal defences. Analyses on different levels of taxonomic resolution showed that the level of family maintained sufficient information to detect the impacts of coastal defence structures on benthic assemblages. Moreover, the outcomes depended on the transformation used. Patterns of distribution of bivalves, used as surrogates, showed low correlations with the patterns of the total macrofaunal species assemblages. Patterns of polychaetes, if identified to the species or genus level showed higher correlations with the whole dataset. However, the identification of polychaetes to species and genus level is as costly as the identification of all macrobenthic taxa at family level.This study provided additional evidences that taxonomic sufficiency is a useful tool in environmental monitoring, also in investigations on the impacts of coastal defence structures on subtidal macrofauna. The use of coarser taxonomic level, being time-efficient, would allow improving sampling designs of monitoring programs by increasing replication in space and time and by allowing long term monitoring studies.     

ELECTROMAGNETIC PROSPECTING FOR GROUNDWATER IN PRECAMBRIAN TERRAINS IN THE REPUBLIC OF UPPER VOLTA*

The “Autorité des amenegements des valées des Voltas (AVV)” is establishing new rural settlements in the Volta valleys. First, a survey of available water supplies is performed. Economic aquifers in Precambrian terrains are deep (15–50 m) and usually occur in fractured zones accompanying faults. Such zones can be identified on aerial photographs, but their precise location on the ground is virtually impossible by visual means. Because of the small size of the aquifers, a location error of 5 m can make the difference between a productive well and a dry hole. Traditionally, resistivity profiling has been used as the means of locating the fractured zones in the field. Our studies suggest that the task can be performed faster, cheaper and more accurately by VLF and EM methods. Because of the limited choice of transmitting stations reccivable in Upper Volta, the VLF method is not sufficiently sensitive to detect conductors with a strike between 45° and 105°. The results obtained with a multifrequency, horizontal-loop EM (HLEM) system were satisfactory in all investigated areas. During the 1980 field season, 35 target areas were surveyed. Of the 24 holes drilled so far, 23 are productive. The weathered layer is a source of distinctive HLEM anomalies, which are characteristic of the underlying rocks. Therefore, different interpretational procedures had to be developed for granitic and volcano-sedimentary areas. Despite the high background level of in-phase and quadrature components, which varied with thickness and conductivity of the weathered layer, aquifers could be detected at a depth greater than 30 m. Attempts were made to interpret the HLEM results quantitatively using two models: a three-layer medium and a valley discontinuity. The latter model is more realistic, but more scale modelling will have to be performed to permit development of viable interpretational procedures. Meanwhile, phasor diagrams based on drilling and resistivity sounding data can be used to estimate the aquifer depth.     

Electromagnetic sounding of geothermal zones

Electromagnetic (EM) data provide a substantial contribution to the geophysical mapping and monitoring of geothermal reservoirs. This paper presents an up-to-date picture of the achievements of EM methods for geothermal exploration as they have emerged over the last few years. It has been proved that EM sounding of geothermal zones and distant monitoring macro-parameters of the reservoirs, fluid-filled faults and other elements of the geothermal system are possible provided that modern 3-D inversion is used along with techniques that reduce the effects of industrial and geological noise. In addition, geological and petrophysical data also need to be included in the analysis.     

A Bayesian trans-dimensional approach for the fusion of multiple geophysical datasets

We propose a Bayesian fusion approach to integrate multiple geophysical datasets with different coverage and sensitivity. The fusion strategy is based on the capability of various geophysical methods to provide enough resolution to identify either subsurface material parameters or subsurface structure, or both. We focus on electrical resistivity as the target material parameter and electrical resistivity tomography (ERT), electromagnetic induction (EMI), and ground penetrating radar (GPR) as the set of geophysical methods. However, extending the approach to different sets of geophysical parameters and methods is straightforward. Different geophysical datasets are entered into a trans-dimensional Markov chain Monte Carlo (McMC) search-based joint inversion algorithm. The trans-dimensional property of the McMC algorithm allows dynamic parameterisation of the model space, which in turn helps to avoid bias of the post-inversion results towards a particular model. Given that we are attempting to develop an approach that has practical potential, we discretize the subsurface into an array of one-dimensional earth-models. Accordingly, the ERT data that are collected by using two-dimensional acquisition geometry are re-casted to a set of equivalent vertical electric soundings. Different data are inverted either individually or jointly to estimate one-dimensional subsurface models at discrete locations. We use Shannon's information measure to quantify the information obtained from the inversion of different combinations of geophysical datasets. Information from multiple methods is brought together via introducing joint likelihood function and/or constraining the prior information. A Bayesian maximum entropy approach is used for spatial fusion of spatially dispersed estimated one-dimensional models and mapping of the target parameter. We illustrate the approach with a synthetic dataset and then apply it to a field dataset. We show that the proposed fusion strategy is successful not only in enhancing the subsurface information but also as a survey design tool to identify the appropriate combination of the geophysical tools and show whether application of an individual method for further investigation of a specific site is beneficial.     

Electrical resistivity tomography survey for delineating uncharted mine galleries in West Bengal,India*

The history of subsidence, fires, flooding and other kinds of environmental hazards related to shallow coal workings in India goes back to colonial times some 300 years ago. As coal production accelerated in modern times, so did the environmental and socio-economic drawbacks related to exploitation. In the mid-1980s, a hydropneumatic sand-stowing method was developed to fill in abandoned galleries but their exact location had to be known. Unfortunately, most of these old workings are uncharted and consequently large tracts of land cannot be stabilized. A research programme making use of integrated surface, borehole and cross-hole geophysical methods was undertaken over a five-year span to try to solve this problem. Surface geophysical methods, being cheaper and faster than their cross- and downhole counterparts, were used to cover larger areas on an exploratory basis, while cross-hole methods were employed to locate more accurately one or a network of galleries to be perforated by drillhole(s) and used as a conduit for sand stowing. We report the results of one of the cross-hole geophysical methods: electrical resistivity tomography (ERT). A pole–dipole configuration is used and both cross-hole and surface-borehole methodologies are tested. Forward modelling and inversion of synthetic data making use of downhole and surface physical and geometrical parameters are presented first. This phase is followed by the inversion of real data. It is concluded that ERT is not applicable for the detection of dry voids, but is effective in a waterlogged environment which is estimated to represent 85–90% of the cases. In waterlogged galleries, ERT is applicable in both cross-hole and surface-downhole modes, the latter allowing a larger surface coverage at low cost. ERT is thus a reliable geophysical tool to image water-filled voids and an adequate technique to address environmental and geotechnical problems.     

Numerical simulation of catastrophic flood: the case study of hypothetical failure of the Bielkowo hydro-power plant reservoir

The numerical modeling of flood wave propagation following the hypothetical breaks of the embankments of the Bielkowo hydro-power plant storage reservoir (Kolbudy II Reservoir) on the Radunia River in Poland has been presented. The results of computations were used to estimate the parameters of the flood waves, which are indispensable for the flood zone determination and mapping and then for the flood risk analysis. When estimating the reach and area of the inundation, related to the embankments failures, digital terrain model, and mathematical model of flood wave propagation are necessary. For the numerical simulations of flood, the mathematical model of free surface, two-dimensional unsteady water flow was applied. Four locations of potential breaks of the reservoir embankments were considered. The computed flood zones were presented on the flood hazard maps. The maps have been used by the local authorities and the dam owner to manage the flood risk related to hydro-power plants operations on the Radunia River. This type of research has been done for the first time for the water plant managed by the ENERGA Elektrownie Straszyn.     

Three-dimensional interpretation of TEM soundings

We describe the approach to the interpretation of electromagnetic (EM) sounding data which iteratively adjusts the three-dimensional (3D) model of the environment by local one-dimensional (1D) transformations and inversions and reconstructs the geometrical skeleton of the model. The final 3D inversion is carried out with the minimal number of the sought parameters. At each step of the interpretation, the model of the medium is corrected according to the geological information. The practical examples of the suggested method are presented.     

基于同步大地电磁时间序列依赖关系的噪声处理

本文从信号与系统的角度讨论了同步大地电磁时间序列信号之间的依赖关系,选取高信噪比的时间序列信号作为先验数据,用最小二乘法估算依赖关系;结合参考道的数据,合成本地道含噪声时段的数据,最后用合成数据替代噪声段数据,组成新数据,从而在时域中去除大地电磁噪声.西藏地区高信噪比实测数据的试算结果表明,无论电场还是磁场,信号之间的依赖关系是相对稳定的,只与先验数据的长度有关,与时间无关;虽然不同参考点之间的依赖关系不同,但都可以精确合成本地点数据,与参考点地下电性结构和参考距离无关.仿真实验显示,去噪后的信号与原始信号基本一致.实测数据处理结果表明,该方法可以有效去除强噪声干扰,抑制中高频段的近场源效应,同时保留了微弱的有效信号,保证了处理结果的正确性.最后针对方差比方法无法识别的方波噪声,提出了一种简单的平移方法,成功去除了持续时间大于窗口长度的方波噪声;将该方法与远参考技术结合,可以有效抑制近场源噪声干扰,获得光滑连续并且可信的测深资料.     

Electrical imaging of lateritic weathering mantles over granitic and metamorphic basement of eastern Senegal, West Africa

The electrical properties of several tens of metres of lateritic weathering mantle were investigated by using electrical resistivity tomography (ERT) in two basement areas of eastern Senegal. The field survey was conducted along two profiles providing continuous coverage. Colour-modulated pseudosections of apparent resistivity vs. pseudo-depth were plotted for all survey lines, giving an approximate image of the subsurface structure. In the area underlain by granitic basement, the pseudosection suggests a very inhomogeneous weathered layer in which the apparent resistivity changes more rapidly than thickness. In the second area, underlain by schists, the lateral changes in electrical properties are less pronounced than those of the granitic area. Interpretation of 2D Wenner resistivity data yielded considerable detail about the regolith, even without pit information. In both areas, the near-surface topsoil comprising undersaturated lateritic material is highly resistive. The intermediate layer with low resistivities (e.g., 20–100 Ωm) contains clays including small quantities of water. The third, highly resistive layer reflects the granitic basement. Comparison of ERT survey results with pit information shows general agreement and suggests that ERT can be used as a fast and efficient exploration tool to map the thick lateritic weathering mantle in tropical basement areas with hard rock geology.