Delineation of a quick clay zone at Smørgrav,Norway, with electromagnetic methods under geotechnical constraints

In many coastal areas of North America and Scandinavia, post-glacial clay sediments have emerged above sea level due to iso-static uplift. These clays are often destabilised by fresh water leaching and transformed to so-called quick clays as at the investigated area at Smørgrav, Norway. Slight mechanical disturbances of these materials may trigger landslides. Since the leaching increases the electrical resistivity of quick clay as compared to normal marine clay, the application of electromagnetic (EM) methods is of particular interest in the study of quick clay structures.For the first time, single and joint inversions of direct-current resistivity (DCR), radiomagnetotelluric (RMT) and controlled-source audiomagnetotelluric (CSAMT) data were applied to delineate a zone of quick clay. The resulting 2-D models of electrical resistivity correlate excellently with previously published data from a ground conductivity metre and resistivity logs from two resistivity cone penetration tests (RCPT) into marine clay and quick clay. The RCPT log into the central part of the quick clay identifies the electrical resistivity of the quick clay structure to lie between 10 and 80 Ω m. In combination with the 2-D inversion models, it becomes possible to delineate the vertical and horizontal extent of the quick clay zone. As compared to the inversions of single data sets, the joint inversion model exhibits sharper resistivity contrasts and its resistivity values are more characteristic of the expected geology. In our preferred joint inversion model, there is a clear demarcation between dry soil, marine clay, quick clay and bedrock, which consists of alum shale and limestone.     

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.     

Electrical resistivity and seismic refraction tomography to detect buried cavities

Near‐surface cavities can pose serious hazards to human safety, especially in highly urbanized town centres. The location of subsurface voids, the estimation of their size and the evaluation of the overburden thickness are necessary to assess the risk of collapse. In this study, electrical resistivity tomography (ERT) and seismic refraction tomography data are integrated in a joint interpretation process for cavity location in the city of Rome. ERT is a well established and widely employed method for cavity detection. However, additional information provided by seismic refraction tomography is capable of eliminating some potential pitfalls in resistivity data interpretation. We propose that the structure of the cavities defined by ERT can be used as a base to optimize seismic refraction tomography investigations within the framework of a joint interpretation process. Data integration and the insertion of a priori information are key issues for reducing the uncertainties associated with the inversion process and for optimizing both acquisition procedures and computation time. Herein, the two geophysical methods are tested on both synthetic and real data and the integration of the results is found to be successful in detecting isolated cavities and in assessing their geometrical characteristics. The cavity location inferred by geophysical non‐invasive methods has been subsequently confirmed by direct inspection.     

采空区地震与电法探测数据的联合反演研究及应用

综合地球物理技术在采空区的探测中发挥了重要作用.目前通常采用单方法反演、仅对不同方法反演结果进行对比解释的综合勘探方式,单方法反演的多解性严重降低了其探测精度.如何提高采空区的探测精度,对采空区进行有效探测一直被认为是地球物理技术面临的首要难题.为了提高地震与电法技术的探测精度,基于交叉梯度联合反演理论,设计了地震初至折射走时数据和高密度电法数据的联合反演算法流程,对采空区理论模型和野外实际数据进行了联合反演处理.结果发现通过两者的联合反演,不仅可以提高采空区电阻率反演模型的成像效果,而且能够获得地震单方法反演难以成像的采空区低速异常体,从而提高了地震与电法技术对采空区的探测精度.表明地震与电法探测数据联合反演是一种提高采空区探测精度的有效方法.     

Two-dimensional joint inversion of radiomagnetotelluric and direct current resistivity data

An algorithm for the two-dimensional (2D) joint inversion of radiomagnetotelluric and direct current resistivity data was developed. This algorithm can be used for the 2D inversion of apparent resistivity data sets collected by multi-electrode direct current resistivity systems for various classical electrode arrays (Wenner, Schlumberger, dipole-diplole, pole-dipole) and radiomagnetotelluric measurements jointly. We use a finite difference technique to solve the Helmoltz and Poisson equations for radiomagnetotelluric and direct current resistivity methods respectively. A regularized inversion with a smoothness constrained stabilizer was employed to invert both data sets. The radiomagnetotelluric method is not particularly sensitive when attempting to resolve near-surface resistivity blocks because it uses a limited range of frequencies. On the other hand, the direct current resistivity method can resolve these near-surface blocks with relatively greater accuracy. Initially, individual and joint inversions of synthetic radiomagnetotelluric and direct current resistivity data were compared and we demonstrated that the joint inversion result based on this synthetic data simulates the real model more accurately than the inversion results of each individual method. The developed 2D joint inversion algorithm was also applied on a field data set observed across an active fault located close to the city of Kerpen in Germany. The location and depth of this fault were successfully determined by the 2D joint inversion of the radiomagnetotelluric and direct current resistivity data. This inversion result from the field data further validated the synthetic data inversion results.     

River terrace sand and gravel deposit reserve estimation using three-dimensional electrical resistivity tomography for bedrock surface detection

We describe the application of 3D electrical resistivity tomography (ERT) to the characterisation and reserve estimation of an economic fluvial sand and gravel deposit. Due to the smoothness constraints used to regularise the inversion, it can be difficult to accurately determine the geometry of sharp interfaces. We have therefore considered two approaches to interface detection that we have applied to the 3D ERT results in an attempt to provide an accurate and objective assessment of the bedrock surface elevation. The first is a gradient-based approach, in which the steepest gradient of the vertical resistivity profile is assumed to correspond to the elevation of the mineral/bedrock interface. The second method uses an intrusive sample point to identify the interface resistivity at a location within the model, from which an iso-resistivity surface is identified that is assumed to define the interface. Validation of these methods has been achieved through direct comparison with observed bedrock surface elevations that were measured using real-time-kinematic GPS subsequent to the 3D ERT survey when quarrying exposed the bedrock surface. The gradient-based edge detector severely underestimated the depth to bedrock in this case, whereas the interface resistivity method produced bedrock surface elevations that were in close agreement with the GPS-derived surface. The failure of the gradient-based method is attributed to insufficient model sensitivity in the region of the bedrock surface, whereas the success of the interface resistivity method is a consequence of the homogeneity of the mineral and bedrock, resulting in a consistent interface resistivity. These results highlight the need for some intrusive data for model validation and for edge detection approaches to be chosen on the basis of local geological conditions.     

High resolution 3D ERT to help GPR data interpretation for researching archaeological items in a geologically complex subsurface

Muro Leccese (Lecce) contains one the most important Messapian archaeological sites in southern Italy.The archaeological interest of the site arises from the discovery of the remains of Messapian walls, tombs, roads, etc. (4th–2nd centuries BC) in the neighbourhood. The archaeological remains were found at about 0.3 m depth.At present the site belongs to the municipality, which intends to build a new sewer network through it. The risk of destroying potentially interesting ancient archaeological structures during the works prompted an archaeological survey of the area. The relatively large dimensions of the area (almost 10,000 m²), together with time and cost constraints, made it necessary to use geophysical investigations as a faster means to ascertain the presence of archaeological items. Since the most important targets were expected to be located at a soil depth of about 0.3 m, a ground-penetrating radar (GPR) survey was carried out in an area located near the archaeological excavations. Unfortunately the geological complexity did not allow an easy interpretation of the GPR data.Therefore a 3D electrical resistivity tomography (ERT) scan was conducted in order to resolve these interpretation problems.A three-way comparison of the results of the dense ERT measurements parallel to the x axis, the results of the measurements parallel to the y axis and the combined results was performed.Subsequently the synthetic model approach was used to provide a better characterization of the resistivity anomalies visible on the ERT field data.The 3D inversion results clearly illustrate the capability to resolve in view of quality 3D structures of archaeological interest. According to the presented data the inversion models along one direction (x or y) seems to be adequate in reconstructing the subsurface structures.Naturally field data produce good quality reconstructions of the archaeological features only if the x-line and y-line measurements are considered together. Despite the increased computational time required by the 3D acquisition and 3D inversion schemes, good quality results can be produced.     

Assessment of the reliability of 2D inversion of apparent resistivity data

The reliability of inversion of apparent resistivity pseudosection data to determine accurately the true resistivity distribution over 2D structures has been investigated, using a common inversion scheme based on a smoothness‐constrained non‐linear least‐squares optimization, for the Wenner array. This involved calculation of synthetic apparent resistivity pseudosection data, which were then inverted and the model estimated from the inversion was compared with the original 2D model. The models examined include (i) horizontal layering, (ii) a vertical fault, (iii) a low‐resistivity fill within a high‐resistivity basement, and (iv) an upfaulted basement block beneath a conductive overburden. Over vertical structures, the resistivity models obtained from inversion are usually much sharper than the measured data. However, the inverted resistivities can be smaller than the lowest, or greater than the highest, true model resistivity. The substantial reduction generally recorded in the data misfit during the least‐squares inversion of 2D apparent resistivity data is not always accompanied by any noticeable reduction in the model misfit. Conversely, the model misfit may, for all practical purposes, remain invariant for successive iterations. It can also increase with the iteration number, especially where the resistivity contrast at the bedrock interface exceeds a factor of about 10; in such instances, the optimum model estimated from inversion is attained at a very low iteration number. The largest model misfit is encountered in the zone adjacent to a contact where there is a large change in the resistivity contrast. It is concluded that smooth inversion can provide only an approximate guide to the true geometry and true formation resistivity.     

3-D Magnetotelluric Exploration of Tenerife Geothermal System (Canary Islands,Spain)

The resistivity structure of the Tenerife geothermal system has been determined by the 3-D inversion of data from different magnetotelluric surveys. In this paper, the ocean and topography effects on the magnetotelluric data were investigated by constructing a 3-D conceptual geoelectrical model of the island. The study showed that these effects should be taken into account in order to obtain a reliable subsurface model of the island. Data from 148 sites were used during three-dimensional inversion. The most interesting feature in the final geoelectrical model of the geothermal system is a low resistivity structure (<10 Ωm) above the resistive core of the system. The low resistivity structure has been interpreted as a hydrothermal clay alteration cap typically generated in the conventional geothermal systems. The resistivity model has been correlated with a recent seismic velocity model, showing that a low resistivity structure surrounds an area with high P wave velocity and medium–high resistivity. This medium–high resistivity area can be associated with a slowly solidified magma and, therefore, with a hotter part of the geothermal system.     

Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration

In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.     

Characterizing solute transport in undisturbed soil cores using electrical and X‐ray tomographic methods

Solute transport in undisturbed soil is a complex process and detailed information on the transport characteristics is needed to provide fundamental understanding of the processes involved. X‐ray computer tomography (CT) and electrical resistivity tomography (ERT) have been used to gain information on the transport characteristics. Both methods are non‐intrusive and do not disturb the soil, in contrast to other methods. CT provides high resolution information on bulk density and macropores, while ERT provides a three‐dimensional image of the internal resistivity structure. By adding a suitable solute under steady‐state flow, the internal resistivity changes can be interpreted as a change in resident concentrations. In our experiment two cores from different field sites were investigated. The ERT measurements revealed two transport modes (one fast and one slow) in one of the cores and only one mode in the other. This was consistent with the results of transfer function modelling on the independently measured breakthrough curves (BTCs). The fast transport mode is perhaps a result of many connected macropores, detected by CT, but this could not be verified with the ERT measurements because of the coarser resolution. However, with ERT in both cases we were able to explain the observed BTC qualitatively. Copyright © 1999 John Wiley & Sons, Ltd.     

Joint elastic‐electrical effective medium models of reservoir sandstones

Improvements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic‐electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three‐phase effective medium model for sandstones with pore‐filling clay minerals based on the combined self‐consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P‐wave velocity and electrical resistivity and of the divergent trends of clean and clay‐rich sandstones at higher porosities. Using only a few well‐constrained input parameters, the new model offers a practical way to predict in situ porosity and clay content in brine saturated sandstones from co‐located P‐wave velocity and electrical resistivity data sets.     

GPR to constrain ERT data inversion in cavity searching: Theoretical and practical applications in archeology

I used theoretical forward models to show that a cavity embedded in a stratified sedimentary sequence can induce an equivalence problem in the ERT data inversion. Conductive top soil increases the misfit between the ground feature and the ERT model. The misfit depends on array and stratigraphy sequences. The latter induce an equivalence problem that manifests itself as wrong cavity depth positioning. The misfit is greater in the data acquired with Schlumberger array than with dipole–dipole.The ambiguity of ERT data inversion problems was tested in the detection of cavities linked to an 8th–6th century B.C. Sabine tomb, 3 m wide × 3 m long × 2 m high, excavated from a shaly gray volcanic ash (cinerite) layer covered by semi-lithoid tuff and top soil layers. In the real study I reduced the ambiguity in the inverse problem of ERT data using a priori information on geometry and resistivity of the cavity. The constrains were carried out from georadar data acquired with 80 and 200 MHz antenna. I demonstrate that this procedure has a practical application in cavity detection, and is a key to the reduction of the uncertainty inherent in the inversion process of ERT data.     

Geophysical investigation of tumuli by means of surface 3D Electrical Resistivity Tomography

Tumuli are artificially erected small hills that cover monumental tombs or graves. In this work, the surface three-dimensional (3D) Electrical Resistivity Tomography (ERT) method, composed of dense parallel two-dimensional (2D) tomographies, was used to investigate the properties of the tumuli filling material and to resolve buried archaeological structures inside the tumuli.The effectiveness of the method was investigated by numerical modeling and through 3D inversion of synthetic apparent resistivity data. A resistivity model that simulates the inhomogeneous tumulus material and the tombs that are buried inside the tumulus was assumed. The Dipole–Dipole (DD), Pole–Dipole (PD), Pole–Pole (PP), Gradient (GRAD), Midpoint-Potential-Referred (MPR) and Schlumberger Reciprocal (SCR) arrays, which are suitable for multichannel resistivity instruments, were tested. The tumulus topography (pyramid or capsized cup) was incorporated into the inversion procedure through a distorted finite element mesh. The inversion procedure was based on a smoothness constrained Gauss–Newton algorithm in which the Active Constraint Balancing (ACB) method was also applied in order to enhance the least-squares resolving power and stability.Synthetic modeling showed that the different tumulus layers and the horizontal contact of the artificial tumulus material with the natural background soil were reconstructed by all of the tested electrode arrays. Generally, PD and the GRAD arrays comprise the optimum choices to investigate the subsurface properties of a tumulus and locate buried tombs. The MPR model was inferior to the GRAD model, while the DD, PP and SCR models had the poorest resolution. It was also shown that the inversion models are practically independent from the survey direction and the topography shape of the tumulus.The real field data collected employing the PD array along a small tumulus from the archaeological site of Vergina in northern Greece enhanced the synthetic modeling results. The inversion model outlined a number of archaeological structures that exhibit a high possibility to correlate with graves. Overall, this work signifies that the surface 3D ERT method can provide a valuable tool in the non-destructive archaeological exploration of tumuli.     

大地电磁的人工鱼群最优化约束反演

大地电磁的反演问题是非线性,如果采用线性反演方法容易陷入局部极小,使得反演结果非唯一性严重.本文将人工鱼群算法引入到地球物理反演之中,提出了非线性的大地电磁人工鱼群最优化反演.该方法不需要进行偏导数的求取,可以对反演的范围进行约束,以减小反演结果的非唯一性.同时我们对搜索步长进行了改进,给出适用于大地电磁反演的人工鱼群参数.大量的理论数据试算表明,人工鱼群反演算法能够较好地寻找到全局最优解.实测数据的处理结果表明,该方法可以用来处理实际资料,并且能够取得很好的应用效果.     

航空电磁拟三维模型空间约束反演

为了克服时间域航空电磁数据单点反演结果中常见的电阻率或层厚度横向突变造成数据难以解释的问题,通过引入双向约束实现航空电磁拟三维空间约束反演.除考虑沿测线方向相邻测点之间的横向约束外,同时还考虑了垂直测线方向测点在空间上的相互约束.为此,首先设计拟三维模型中固定层厚和可变层厚两种空间约束反演方案,然后通过在目标函数中引入沿测线和垂直测线方向上的模型参数约束矩阵,并使用L-BFGS算法使目标函数最小化,获得最优拟三维模型空间反演解.基于理论模型和实测数据反演,对单点反演与两种空间约束反演方案的有效性进行比较,证明本文空间约束反演算法对于噪声的压制效果好,反演的界面连续光滑,同时内存需求和反演时间少,是一种快速有效的反演策略.     

Laterally constrained inversion of helicopter-borne frequency-domain electromagnetic data

Helicopter-borne frequency-domain electromagnetic (HEM) surveys are used for fast high-resolution, three-dimensional resistivity mapping. Standard interpretation tools are often based on layered earth inversion procedures which, in general, explain the HEM data sufficiently. As a HEM system is moved while measuring, noise on the data is a common problem. Generally, noisy data will be smoothed prior to inversion using appropriate low-pass filters and consequently information may be lost.For the first time the laterally constrained inversion (LCI) technique has been applied to HEM data combined with the automatic generation of dynamic starting models. The latter is important because it takes the penetration depth of the electromagnetic fields, which can heavily vary in survey areas with different geological settings, into account. The LCI technique, which has been applied to diverse airborne and ground geophysical data sets, has proven to be able to improve the HEM inversion results of layered earth structures. Although single-site 1-D inversion is generally faster and — in case of strong lateral resistivity variations — more flexible, LCI produces resistivity — depth sections which are nearly identical to those derived from noise-free data.The LCI results are compared with standard single-site Marquardt–Levenberg inversion procedures on the basis of synthetic data as well as field data. The model chosen for the generation of synthetic data represents a layered earth structure having an inhomogeneous top layer in order to study the influence of shallow resistivity variations on the resolution of deep horizontal conductors in one-dimensional inversion results. The field data example comprises a wide resistivity range in a sedimentary as well as hard-rock environment.If a sufficient resistivity contrast between air and subsurface exists, the LCI technique is also very useful in correcting for incorrect system altitude measurements by using the altitude as a constrained inversion parameter.     

电阻率和速度随机分布的MT与地震联合反演

在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度.     

飞行高度同时反演的固定翼航空瞬变电磁-维反演

航空电磁测量记录中,不仅感生电动势测量数据有观测误差,而且高度计测量数据也有误差,直接进行常规反演往往导致反演结果不可靠,研究飞行高度数据有误差下的反演算法具有实际意义.本文以层状模型的固定翼时间域航空电磁多分量理论响应数据为例,提出了两种针对飞行高度计记录数据有误差时的正则化反演算法,一个是自适应正则化反演方法,另一...