首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
The integrated use of electrical resistivity tomography (ERT) and ground penetrating radar (GPR) measurements, and in particular the joint analysis of 2D and 3D data, can represent a valid solution for target identification at complex archaeological sites. A good example, in this respect, is given by the case study of a Phoenician–Punic necropolis in the archaeological site of Nora, in southern Sardinia (Italy), where GPR and ERT measurements were collected before site excavation. In this specific case, the mix of soil and air in the buried chambers, as well as the orientation and the complex spatial distribution of these structures into the sandstone bedrock, generated a number of anomalies difficult to interpret only using 2D results. Only the integration of all GPR and ERT data in a 3D view, and the comparison with archaeological evidence after the excavation, allowed a solid interpretation of geophysical anomalies visible in the 2D sections. Overall, this case study demonstrates the efficiency of the combined use of GPR and ERT acquisitions and shows how, in general, only the joint analysis of 2D data and in a 3D view can help the interpretation of the real distribution of the buried archaeological remains at similar archaeological complex sites.  相似文献   

2.
An integrated geophysical survey has been conducted at the Tarragona’s Cathedral (Catalonia, NE Spain) with the aim to confirm the potential occurrence of archaeological remains of the Roman Temple dedicated to the Emperor Augustus. Many hypotheses have been proposed about its possible location, the last ones regarding the inner part of the Cathedral, which is one of the most renowned temples of Spain (twelfth century) evolving from Romanesque to Gothic styles. A geophysical project including electrical resistivity tomography (ERT) and ground probing radar (GPR) was planned over 1 year considering the administrative and logistic difficulties of such a project inside a cathedral of religious veneration. Finally, both ERT and GPR have been conducted during a week of intensive overnight surveys that provided detailed information on subsurface existing structures. The ERT method has been applied using different techniques and arrays, ranging from standard Wenner–Schlumberger 2D sections to full 3D electrical imaging with the advanced Maximum Yield Grid array. Electrical resistivity data were recorded extensively, making available many thousands of apparent resistivity data to obtain a complete 3D image after a full inversion. In conclusion, some significant buried structures have been revealed providing conclusive information for archaeologists. GPR results provided additional information about shallowest structures. The geophysical results were clear enough to persuade religious authorities and archaeologists to conduct selected excavations in the most promising areas that confirmed the interpretation of geophysical data. In conclusion, the significant buried structures revealed by geophysical methods under the cathedral were confirmed by archaeological digging as the basement of the impressive Roman Temple that headed the Provincial Forum of Tarraco, seat of the Concilium of Hispania Citerior Province.  相似文献   

3.
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.  相似文献   

4.
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.  相似文献   

5.
This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).  相似文献   

6.
Electrical resistivity mapping and electrical resistivity profiling are powerful instruments for investigating archaeological structures. Interpretation of geoelectrical data is complicated by near-surface anomalies and the characteristics of the applied electrode arrays. Averaging Wenner α and Wenner β data as an alternative method of focused imaging is presented to overcome these problems. The mechanism of focused imaging is explained using the sensitivity distribution of the combined arrays. Various methods of imaging geoelectrical data are examined with synthetic and field data. In electrical resistivity mapping, inversion of the data is unnecessary when using focused imaging. In electrical resistivity profiling, focused imaging gives a first idea about the subsurface resistivity distribution without achieving the quality obtainable by inversion.  相似文献   

7.
A 3D electrical resistivity imaging survey is presented in this paper. The objective was to investigate an underground wastewater system at the University of Malaya, Malaysia. Apparent resistivity data were collected along ten parallel lines using a Wenner-Schlumberger configuration; electrode cables were oriented in the x-direction with 3 m spacing. Roll-along measurements using a line spacing of 3 m were carried out covering a grid of 20 × 10 electrodes. All data sets were merged into a single data file in order to perform a 3D inversion. Two different 3D least squares algorithms, based on the robust inversion method and the smoothness-constrained technique, were used for the inversion of the apparent resistivity data. Both the horizontal and vertical extents of the anomalous zones found by inversion are displayed. The results indicate the superiority of the robust inversion method over the smoothness-constrained technique at this site. The results are in sufficient accordance with previously known information about the investigation area. The results show that 3D electrical resistivity imaging surveys, in combination with an appropriate 3D inversion method, can be highly useful for engineering and archaeological investigations as well as for environmental applications.  相似文献   

8.
Integrated electrical resistance tomography (ERT) and short-offset transient electromagnetic (TEM) measurements were carried out to investigate a geothermal area in the Main Central Thrust (MCT) zone of Garhwal Himalayan region, India. The study area is located around Helang on either side of Alaknanda River and it is dotted with hot water springs with water temperature of 45°–55 °C emerging at the surface.To assess the geothermal potential and its lateral and vertical extension in and around the hot water springs in the study area, 7 ERT profiles and 21 TEM stations on 7 profiles were established around the hot water spring and at far distant locations. The 2D inversion of ERT data indicates a low resistivity (< 50 Ωm) zone in the vicinity of hot springs, which appears to be associated with an underground water channel through the fractured rock. The bedrock resistivity is very high (> 1000 Ωm) whereas the resistivity of the weathered near surface soil at a far distant location from the hot spring is low (< 100 Ωm) again. A common feature of all TEM data is the sign reversal observed at roughly 10 μs. The consistent sign reversal in all TEM data indicates the existence of the multi-dimensionality of the geoelectrical structure. Therefore, the TEM data were treated by using the SLDM (Spectral Lanczos Decomposition Method) 2D/3D forward modeling code based on the finite difference algorithm. The resistivity structure obtained from ERT data was used as an input for the modeling of TEM data. Based on the joint analysis of the ERT and TEM data it can be inferred that geothermal anomalies associated with the hot spring in the MCT zone are a local feature appearing as a low resistivity zone (< 50 Ωm) at shallow depth (< 100 m) in the vicinity of the hot spring region.  相似文献   

9.
A 3D ERT study of solute transport in a large experimental tank   总被引：2，自引：0，他引：2  
A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m3. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.  相似文献   

10.
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.  相似文献   

11.
基于非结构网格的电阻率三维带地形反演   总被引：6，自引：3，他引：3       下载免费PDF全文
吴小平  刘洋  王威 《地球物理学报》2015,58(8):2706-2717
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.  相似文献   

12.
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.  相似文献   

13.
Seismic noise attenuation is very important for seismic data analysis and interpretation, especially for 3D seismic data. In this paper, we propose a novel method for 3D seismic random noise attenuation by applying noncausal regularized nonstationary autoregression (NRNA) in f–x–y domain. The proposed method, 3D NRNA (f–x–y domain) is the extended version of 2D NRNA (f–x domain). f–x–y NRNA can adaptively estimate seismic events of which slopes vary in 3D space. The key idea of this paper is to consider that the central trace can be predicted by all around this trace from all directions in 3D seismic cube, while the 2D f–x NRNA just considers that the middle trace can be predicted by adjacent traces along one space direction. 3D f–x–y NRNA uses more information from circumjacent traces than 2D f–x NRNA to estimate signals. Shaping regularization technology guarantees that the nonstationary autoregression problem can be realizable in mathematics with high computational efficiency. Synthetic and field data examples demonstrate that, compared with f–x NRNA method, f–x–y NRNA can be more effective in suppressing random noise and improve trace-by-trace consistency, which are useful in conjunction with interactive interpretation and auto-picking tools such as automatic event tracking.  相似文献   

14.
In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.  相似文献   

15.
The purpose of geophysical electrical surveys is to determine the subsurface resistivity distribution by making measurements on the ground surface. From these measurements, the true resistivity of the subsurface can be estimated. The ground resistivity is related to various geological parameters, such as the mineral and fluid content, porosity and degree of water saturation in the rock. Electrical resistivity surveys have been used for many decades in hydrogeological, mining and geotechnical investigations. More recently, they have been used for environmental surveys. To obtain a more accurate subsurface model than is possible with a simple 1-D model, a more complex model must be used. In a 2-D model, the resistivity values are allowed to vary in one horizontal direction (usually referred to as the x direction) but are assumed to be constant in the other horizontal (the y) direction. A more realistic model would be a fully 3-D model where the resistivity values are allowed to change in all three directions. In this research, a simulation of the cone penetration test and 2D imaging resistivity are used as tools to simulate the distribution of hydrocarbons in soil.  相似文献   

16.
An important problem of marble-quarry management is assessing the quality and the homogeneity of quarry blocks before excavation. In this study, we decided to image the limestone, which we studied in a marble quarry, in terms of layer thickness, discontinuities and cavities using the ground-penetrating-radar (GPR) method. The method was successfully applied to detect and map the fractures with the cavities in a marble layer according to depth in the Ankara City Polatli Town (Turkey) region, which represents upper Miocene–Pliocene lacustrine carbonate rocks.This paper is based on interactive transparent 3D visualisation of the 2D GPR profiles to determine changes in layer thickness and discontinuities. In addition, this paper indicates the importance of the appropriate opacity-function construction to obtain transparent 3D visualisation. Firstly we acquired and processed parallel 2D GPR profile data, then we assigned two different amplitude–colour ranges using a limited number of colours to determine the layer thickness and its discontinuities separately. We obtained a 3D volume using parallel 2D GPR data and displayed a limited amplitude range by arranging an opacity function. Therefore, we obtained transparent 3D blocks for thickness and discontinuities, and we formulated an interactive 3D display to image the horizontal, vertical and inclined discontinuities and their directions in the x–y plane versus depth. The GPR results were compared with the petrographical investigation on the basis of textural and mineralogical compositions. The vesicular textures within carbonate platform were supported by the GPR results.  相似文献   

17.
Global optimization with very fast simulated annealing (VFSA) in association with joint inversion is performed for 1D earth structures. The inherent problems of equivalence and suppression in electromagnetic (EM) and direct current (DC) resistivity methods are studied. Synthetic phase data from multifrequency sounding using a horizontal coplanar coil system and synthetic apparent resistivity data from Schlumberger DC resistivity measurements are inverted individually and jointly over different types of layered earth structures. Noisy data are also inverted. The study reveals that global optimization of individual data sets cannot solve inherent equivalence or suppression problems. Joint inversion of EM and DC measurements can overcome the problem of equivalence very well. However, a suppression problem cannot be solved even after combination of data sets. This study reveals that the K-type earth structure is easiest to resolve while the A-type is the most difficult. We also conclude that the equivalence associated with a thin resistive layer can be resolved better than that for a thin conducting layer.  相似文献   

18.
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.  相似文献   

19.
Geodynamic activity in the area of West Bohemia is typified by the occurrence of earthquake swarms, Quaternary volcanism and high flux of mantle-derived CO2. The highest swarm activity occurs beneath the eastern edge of the Cheb basin, which is delineated by the NW-SE trending morphologically pronounced Mariánské Lázn?? Fault (MLF) controlling the formation of the basin. The previous trenching survey across the MLF zone has identified several fault strands with possible Quaternary activity. In this paper we present the results of the geophysical survey focused to trace the faults signatures in geophysical sections and to build an image of near surface tectonics. The method of electric resistivity tomography (ERT) along two profiles parallel to the trench identified a strong resistivity contrast between the bodies of sandy gravels in the middle and conductive clayey sands to the west and weathered crystalline basement to the east. The 2-D ground penetration radar (GPR) sections show direct correlation of reflections with lithological boundaries identified in the trench. As expected, the GPR signal amplitudes increase with the resistivities found in the ERT sections. Two of the four faults identified in the trench are indicated in the resistivity and GPR sections. A 3-D GPR measurement has identified a spot of high amplitudes elongated parallel to the MLF trend, which coincides with the high resistivity body. To improve the signal-to-noise ratio of the time slices we stacked the GPR time slices within vertically homogeneous blocks. This provided a contrast image of the sand-gravel body including its boundaries in three dimensions. The detailed analysis of the 3-D GPR cube revealed additional fault that limits the highly reflective sands and appears to be offset by another younger fault. Our results suggest a complex fault pattern in the studied area, which deserves a further study.  相似文献   

20.
  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号