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.     

ADS项目场地稳定性评价中高密度电阻率法的综合应用

以高密度电阻率法为主，综合其他方法对"十二五"国家大科学装置项目——ADS项目的场地稳定性进行评价。ADS项目场址位于广东省惠州市惠东县，通过在场区布设多条高密度电法剖面，进行二维反演后绘制电阻率剖面图，结合已有的钻孔和地质资料对其下断层的分布特征进行推断解释，查明场区内的断层产状、风化层厚度、水文地质条件等；并对厂址附近断层在未来地震作用下可能放大震害的问题进行研究，为拟建场区装置的设计和修建提供相关的技术支持。这项综合手段的评估研究可为同类国家重大项目选址时提供参考。     

高密度电法在水库坝址区勘察中的应用

近年来高密度电法被广泛应用于重大水利工程场地的工程地质调查、坝基及桥墩选址、采空区及地裂缝探测等诸多工程勘察领域。本文利用高密度电法对某水库大坝坝址区的覆盖层厚度、隐伏断层等地质结构进行探测,并在结合钻探资料和区域地质资料综合分析的基础上,取得了理想的勘察结果。查明水库区存在隐伏的断层和覆盖层厚度,为水库坝址选择提供重要的地球物理依据。      

2D and 3D resistivity inversion of Schlumberger vertical electrical soundings in Wadi El Natrun,Egypt: A case study

The Wadi El Natrun area is characterized by a very complicated geological and hydrogeological system. 45 vertical electrical soundings (Schlumberger array) were measured in the study area to elucidate the peculiarity of this unique regime, specifically the nature of waterless area. 2D and 3D resistivity inversion based on the finite element technique and regularization method were applied on the data set. 2D and 3D model resolution was investigated through the use of the Depth and Volume of Investigation Indexes. A very good matching was found between the zones of high resistivity, the waterless area, and the non-productive wells. The low resistivity zones (corresponding to Lower Pliocene clay) were also identified. The middle resistivity fresh water aquifer zones were recognized. Available results can assist in the aquifer management by selecting the most productive zone of groundwater.     

Imaging a 3D geological structure from HEM, airborne magnetic and ground ERT data in Kalat-e-Reshm area, Iran

A set of geophysical data collected in an area in Iran are analyzed to check the validity of a geological map that was prepared in connection to a mineral prospecting project and also to image the spatial electrical resistivity distribution. The data set includes helicopter electromagnetic (HEM), airborne magnetic and ground electrical resistivity measurement. Occam approach was used to invert the HEM data to model the resistivity using a layered earth model with fixed thicknesses. The algorithm is based on a nonlinear inverse problem in a least-squares sense.The algorithm was tested on a part of an HEM dataset acquired with a DIGHEM helicopter EM system at Kalat-e-Reshm, Semnan in Iran. The area contains a resistive porphyry andesite that is covered by Eocene sedimentary units. The results are shown as resistivity sections and maps confirming the existence of an arc like resistive structure in the survey area. The resistive andesite seems to be thicker than it is indicated in the geological maps. The results are compared with the reduced to the pole (RTP) airborne magnetic anomaly field data as well as with two ground resistivity profiles. We found reasonable correlations between the HEM 1D resistivity models and 2D models from electrical resistivity tomography (ERT) inversions. A 3D visualization of the 1D models along all flight lines provided a useful tool for the study of spatial variations of the resistivity structure in the investigation area.     

Gas hydrate saturation at Site C0002, IODP Expeditions 314 and 315, in the Kumano Basin,Nankai trough

The degree of gas hydrate saturation at Integrated Ocean Drilling Program (IODP) Site C0002 in the Kumano Basin, Nankai Trough, was estimated from logging‐while‐drilling logs and core samples obtained during IODP Expeditions 314 and 315. Sediment porosity data necessary for the calculation of saturation were obtained from both core samples and density logs. Two forms of the Archie equation (‘quick‐look’ and ‘standard’) were used to calculate gas hydrate saturation from two types of electrical resistivity log data (ring resistivity and bit resistivity), and a three‐phase Biot‐type equation was used to calculate gas hydrate saturation from P‐wave velocity log data. The gas hydrate saturation baseline calculated from both resistivity logs ranges from 0% to 35%, and that calculated from the P‐wave velocity log ranges from 0% to 30%. High levels of gas hydrate saturation (>60%) are present as spikes in the ring resistivity log and correspond to the presence of gas hydrate concentrations within sandy layers. At several depths, saturation values obtained from P‐wave velocity data are lower than those obtained from bit resistivity data; this discrepancy is related to the presence of free gas at these depths. Previous research has suggested that gas from deep levels in the Kumano Basin has migrated up‐dip towards the southern and seaward edge of the basin near Site C0002. The high saturation values and presence of free gas at site C0002 suggest that a large gas flux is flowing to the southern and seaward edge of the basin from a deeper and/or more landward part of the Kumano Basin, with the southern edge of the Kumano Basin (the location of site C0002) being the main area of fluid accumulation.     

Relation of in situ resistivity to water content in salt rocks1

The investigation of water in salt-rock formations is of particular relevance to underground nuclear waste repositories. In the Asse salt-mine (Germany) a study into the relationship of in situ resistivity to water content has been made. Measurements were carried out in older rock-salt using an electrode array in boreholes, an electrode profile in a drift and small resistivity sensors in and around a drift seal. Further measurements were made on moist zones in a contact area of younger rock-salt and carnallitite and also in older rock-salt with anhydrite bands using electrode profiles in the drifts. The resistivities range from 10²Ωm to 10⁶Ωm. Corresponding probes have water contents from 0.01% to 1.3%. A definite relationship between resistivity and water content is revealed which can be described by Archie's law using a cementation factor m of 1.9. Porosities are between 0.08% and 1.4% and the saturations vary considerably. An explicit influence of saturation on resistivity cannot be discovered using the present data. The results enable us to estimate the in situ water content and the order of the in situ porosity using resistivity surveys at different scales. This increases significantly the safety of a nuclear repository site.     

Site response obtained from array techniques applied to the seismic noise: Two examples in Italy

The vertical component of the seismic noise has been recorded in two different sites near the towns of Mercato S. Severino and Benevento in Southern Italy by a small aperture array, in order to investigate the characteristics of the noise propagation and to study the site response. Three different array techniques have been applied in the two investigated sites: Beam Forming, High Resolution and Spatial Correlation methods. We used two simple array geometry for localising possible noise sources and estimating local shallow structure using ambient noise. The cross shaped array results effective for determining the phase velocity of waves in the case when the noise is from a single localised source; the circular array, on the other hand, is successfully used when the noise sources are distributed. The main results are: the analysis of a coherent component of the noise recorded in the two sites, interpreted as Rayleigh waves, results in reasonable velocity models; the noise recorded in the M.S.S. Plain is a space stationary signal, while the noise at the Benevento site is possibly produced by a stable noise source located close to the array. Due to this evidence, the correlation method does not yield satisfactory results when applied to the Benevento site. The 2–6 Hz spectral peaks of the noise recorded in the M.S.S. Plain can be interpreted as due to a site effect, considering the satisfactory agreement of the noise spectrum at those frequencies with the theoretical transfer function computed on the basis of the velocity model deduced from the Rayleigh waves dispersion analysis.     

The main structure and tectonic features of the Chernorud-Mukhor site on the western shore of Lake Baikal from TEM and SP measurements

Transient electromagnetic (TEM), self-potential (SP) and geoelectrical mapping measurements were carried out at the Chernorud-Mukhor site in the Priolkhonje area on the western shore of the Lake Baikal. All measurements were made along several profiles across the main strike of the regional Primorsky fault. TEM measurements were carried out in a time range from a few tens of microseconds to several tens of milliseconds. The most important result of the 1D modelling of TEM soundings is the discovery of nearly horizontal boundaries that divide high resistive overlying and well conducting underlying rocks. The resistivity of the former is in the range from 100 Ωm to 1000 Ωm, while the resistivity of the latter varies from less than 1 Ωm to several tens of Ωm. This good conductive zone could also be verified by geoelectrical mapping using Schlumberger array (AB/2=100 m). Due to high conductivity of the underlying rocks only the upper boundary of the conductive layer could be determined by TEM soundings. A regional SP anomaly with amplitude of about −450 mV has also been observed above the low resistivity zone indicating the electron nature of its conductance. Geologically, the conductive zone is represented by a graphite-bearing layer within the region of archean rocks. Since that layer extends over a large area, it may be used as a key in studying structures and tectonics of the Priolkhonje area. A 1D TEM geoelectric section shows a wide, gently sloping syncline as a probable base structure of the Chernorud-Mukhor site. Neotectonic faults divide the syncline into vertically displaced blocks that form a wide complicated graben with a total amplitude of about 250 m.     

RESISTIVITY TYPE CURVES OVER OUTCROPPING VERTICAL DYKE—II *

The theoretical horizontal resistivity profiles Over an outcropping vertical dyke with two-electrode and Schlumberger electrode systems are discussed. The two-electrode array seems very useful in locating the highly conducting thick or thin vein, while the Schlumberger (gradient) array is suited to detect the conducting vein of low and moderately high resistivity contrast and as well as the resistive vein of all widths and resistivity contrasts. Besides, the inline gradient array has a distinction of establishing a clue to evaluate the resistivity contrast of the vein.     

Mineshaft imaging using surface and crosshole 3D electrical resistivity tomography: A case history from the East Pennine Coalfield, UK

Hidden mineshafts located in urban areas are a significant problem across much of the industrialized world. Electrical resistivity tomography (ERT) is a technique that can detect and characterize hidden mine entries by exploiting resistivity contrasts between the shaft and surrounding materials, resulting from either compositional or structural differences. A case study is presented in which both surface and crosshole 3D ERT surveys are used to image a hidden backfilled mineshaft at a built environment site, situated on Carboniferous Lower Coal Measures strata in the UK.Backfilled shafts generally present the greatest challenge for detection using geophysical methods, as contrasts between the fill and bedrock are typically low compared to air or water-filled conditions. Nevertheless, the shaft in this case was identified by both the surface and crosshole 3D surveys. The shaft appeared as a strongly resistive anomaly relative to background materials, which we interpreted as resulting from the disturbed fabric of the fill materials rather than any significant compositional differences. The study highlighted the respective strengths and weaknesses of the surface and crosshole ERT methodologies for this type of problem. The surface survey, which covered a non-rectangular area to accommodate irregular boundaries and other physical obstructions, provided a relatively rapid means of investigating the study site. However, this method had a limited depth of investigation and was constrained in its coverage by the locations of buildings. By contrast, the 3D crosshole method was able to image the shaft to the level of the deepest borehole electrodes. Although crosshole ERT is too expensive to be used for large-scale mineshaft surveys, this study clearly demonstrates its suitability for targeted investigations where surface methods cannot be deployed, such as scanning beneath surface structures or in situations where it is essential for resolution to be maintained with depth.     

Deep resistivity structure of the Dikili-Bergama region, west Anatolia, revealed by two- dimensional inversion of vertical electrical sounding data

Within the framework of the National Marine Geological and Geophysical Program, we re‐examined deep vertical electrical sounding (VES) data. The data, measured in 1968 by the General Directorate of Mineral Research and Exploration (MTA) of Turkey with the aim of exploring the deep resistivity structure of the Dikili–Bergama region, focus on the geothermal potential. The geoelectrical resistivity survey was conducted using a Schlumberger array with a maximum electrode half‐spacing of 4.5 km. The two‐dimensional (2D) inversion was utilized to interpret the VES data that were collected along 15‐ to 30‐km profiles. The 2D resistivity–depth cross‐sections obtained show very low resistivity values near the Dikili and Kaynarca hot springs. The 2D inversion results also indicate the presence of fault zones striking nearly N–S and E–W, and fault‐bounded graben‐horst structures that show promising potential for geothermal field resources. The 2D gravity model, which is in good agreement with the density variation of the region, supports the resistivity structure revealed by 2D inversion. The lithology information obtained from the borehole near Kaynarca also confirms the results of the resistivity interpretation and the density model.     

Challenges of using electrical resistivity method to locate karst conduits—A field case in the Inner Bluegrass Region, Kentucky

Conduits serve as major pathways for groundwater flow in karst aquifers. Locating them from the surface, however, is one of the most challenging tasks in karst research. Geophysical methods are often deployed to help locate voids by mapping variations of physical properties of the subsurface. Conduits can cause significant contrasts of some physical properties that can be detected; other subsurface features such as water-bearing fractures often yield similar contrasts, which are difficult to distinguish from the effects of the conduits. This study used electrical resistivity method to search for an unmapped karst conduit that recharges Royal Spring in the Inner Bluegrass karst region, Kentucky, USA. Three types of resistivity techniques (surface 2D survey, quasi-3D survey, and time-lapse survey) were used to map and characterize resistivity anomalies. Some of the major anomalies were selected as drilling targets to verify the existence of the conduits. Drilling near an anomaly identified by an electrical resistivity profile resulted in successful penetration of a major water-filled conduit. The drilling results also suggest that, in this study area, low resistivity anomalies in general are associated with water-bearing features. However, differences in the anomaly signals between the water-filled conduit and other water-bearing features such as water-filled fracture zones were undistinguishable. The electrical resistivity method is useful in conduit detection by providing potential drilling targets. Knowledge of geology and hydrogeology about the site and professional judgment also played important roles in locating the major conduit.     

Assessment of aquifer vulnerability to industrial waste water using resistivity measurements. A case study, along El-Gharbyia main drain, Nile Delta, Egypt

1D resistivity sounding and 2D resistivity imaging surveys were integrated with geological and hydrochemical data to assess the aquifer vulnerability and saltwater intrusion in the north of Nile Delta, Egypt. In the present study, the El-Gharbyia main drain was considered as a case study to map the sand bodies within the upper silt and clay aquitard. Twenty Schlumberger soundings and six 2D dipole-dipole profiles were executed along one profile close to the western side of the main drain. In addition, 14 groundwater samples and 4 surface water samples from the main drain were chemically analyzed to obtain the major and trace elements concentrations.The results from the resistivity and hydrochemical data were used to assess the protection of the groundwater aquifer and the potential risk of groundwater pollution. The inverted resistivities and thicknesses of the layers above the aquifer layer were used to estimate the integrated electrical conductivity (IEC) that can be used for quantification of aquifer vulnerability. According to the aquifer vulnerability assessment of an underlying sand aquifer, the southern part of the area is characterized by high vulnerability zone with slightly fresh to brackish groundwater and resistivity values of 11-23 Ω.m below the clay cap. The resistivity sections exhibit some sand bodies within the clay cap that lead to increase the recharging of surface waste water (650 mg/l salinity) and flushing the upper part of underlying saltwater aquifer. The region in the north has saltwater with resistivity less than 6 Ω.m and local vulnerable zones within the clay cap. The inverted 2D dipole-dipole profiles in the vulnerable zones, in combination with drilling information have allowed the identification of subsoil structure around the main drain that is highly affected by waste water.     

3D inversion of DC data using artificial neural networks

In this paper, we investigate the applicability of artificial neural networks in inverting three-dimensional DC resistivity imaging data. The model used to produce synthetic data for training the artificial neural network (ANN) system was a homogeneous medium of resistivity 100 Ωm with an embedded anomalous body of resistivity 1000 Ωm. The different sizes for anomalous body were selected and their location was changed to different positions within the homogeneous model mesh elements. The 3D data set was generated using a finite element forward modeling code through standard 3D modeling software. We investigated different learning paradigms in the training process of the neural network. Resilient propagation was more efficient than any other paradigm. We studied the effect of the data type used on neural network inversion and found that the use of location and the apparent resistivity of data points as the input and corresponding true resistivity as the output of networks produces satisfactory results. We also investigated the effect of the training data pool volume on the inversion properties. We created several synthetic data sets to study the interpolation and extrapolation properties of the ANN. The range of 100–1000 Ωm was divided into six resistivity values as the background resistivity and different resistivity values were also used for the anomalous body. Results from numerous neural network tests indicate that the neural network possesses sufficient interpolation and extrapolation abilities with the selected volume of training data. The trained network was also applied on a real field dataset, collected by a pole-pole array using a square grid (8 ×8) with a 2-m electrode spacing. The inversion results demonstrate that the trained network was able to invert three-dimensional electrical resistivity imaging data. The interpreted results of neural network also agree with the known information about the investigation area.     

Delineation of Near-Surface Structure in the Southern Part of 15th of May City, Cairo, Egypt Using Geological, Geophysical and Geotechnical Techniques

The integration of geological, geophysical, and geotechnical interpretation at the southern part of 15th of May City, have been used to evaluate the subsurface stratigraphy, especially the clay layer which may cause serious danger to construction. Those techniques have been used to delineate the subsurface structures as normal faults, which play a critical role on the stability of buildings. Geological setting of the area has been evaluated through the construction of a geological map from different geological sections and samples obtained from more than 30 observation points. Geophysical tools such as vertical electrical soundings (12 VESs), 2-D dipole–dipole array (7 sections), P-wave shallow seismic refraction (31 profiles) and multiple channel analysis of surface waves (31 MASW profiles) have been carried out to image the subsurface situation. Geotechnical evaluation using 26 boreholes, samples, laboratory tests and geotechnical parameters has been done at the area of interest. The geological setting demonstrates that the city had been constructed on the second and third members of Qurn Formation (Upper Eocene) composed of argillaceous limestone, marl and shale. Two normal faults are passing through the area were observed. The resistivity (VES and dipole–dipole) and seismic (P-waves and MASW) results reflect the presence of the two normal faults cross the study area, affecting the obtained section of marl, clayey marl and limestone layers. The geotechnical information indicate the presence of the normal faults and the existence of clay layer with swelling ability reaching 140%, which may cause cracks in the upper layers and/or subsidence.     

Magnetotelluric resistivity modeling for 3D characterization of geothermal reservoirs in the Western side of Mt. Aso, SW Japan

Geothermal reservoirs are usually located at a depth range of 2 to 5 km, so to efficiently utilize such resources an advanced prospecting method is needed to detect these deep geologic structures. This study aimed to three-dimensionally characterize geothermal reservoirs by a combination of Magnetotelluric (MT) survey, inversion analysis of apparent resistivity, and interpolation of the resistivity data obtained. The western side of Mt. Aso crater, southwest Japan, was chosen as the case study area. Three hot springs exist there and a fault is assumed to go in the direction connecting them. A MT survey was carried out at 26 sites and the data processed by a remote reference method to reduce artificial noises. Based on skewness and Mohr circle analyses of the impedance tensor, the local geologic structure at each site could be approximated as horizontally layered and therefore, a one-dimensional inversion analysis was applied to the MT raw data. The resultant resistivity column data were then interpolated by the three-dimensional optimization principle method. The resistivity distributions obtained clarified continuous conductors with especially low resistivity (less than 10 Ω·m) at the hot springs along the fault. Because the resistivity decreases largely with an abundance of clay minerals, the conductors could be considered to correspond with the cap rocks. Thus, two geothermal reservoirs, whose shapes were estimated to be pillar, were detected under the cap rocks in an elevation range of − 1000 to − 3000 m. By comparing the resistivity distributions with the temperature distributions based on fluid-flow calculations at a steady state using FEM, the validity of the location and dimension of the estimated reservoirs were confirmed.     

Analysis of Influence Characteristics of Site Ground Motion in Ya'an Area,Sichuan

为分析雅安地区场地结构引起的震动特性,选取该地区重大安评项目中的实际测试钻孔进行场地地震反应分析,采用等效线性化方法,并以NGA强震记录为输入,开展场地反应计算,在全面分析该地区场地地震特征的基础上,以峰值速度为参数指标,研究场地条件引起的峰值速度放大效应。研究结果表明,雅安地区场地卓越周期约为0.2 s,自振周期为0.2 s的建筑结构破坏相对严重;随着基岩峰值加速度的增大,峰值速度呈增大趋势,峰值速度与峰值加速度变化趋势并非完全一致,分析场地效应时应考虑多种地震动参数的变化规律;不同地震输入下峰值速度放大系数为1～1.8;峰值速度放大系数k_v在不同强度基岩输入下,大致呈随等效剪切波速的增大而减小的趋势,软土场地更易引起场地峰值速度的放大,这种现象在芦山地震中较为明显。     

Electromagnetic mapping of electrical conductivity beneath the Columbia basalts1

Sedimentary rocks beneath the Columbia River Basalt Group are recognized as having potential for oil and gas production, but the overlying layered basalts effectively mask seismic reflections from the underlying sediments. Four electromagnetic (EM) methods have been applied on profiles crossing Boylston Ridge, a typical east–west trending anticline of the Yakima Fold Belt, in an attempt to map the resistivity interface between the basalts and the sediments and to map variations in structure and resistivity within the sediments. The EM surveys detected strong variations in resistivity within the basalts, and in particular the continuous magnetotelluric array profiling (EMAP) revealed resistivity lows beneath the surface anticlines. These low resistivity zones probably coincide with fracturing in the core of the anticlines and they appear to correlate well with similar zones of low seismic velocity observed on a nearby seismic profile. The controlled-source EM surveys (in-loop transient, long-offset transient, and variable-offset frequency-domain) were designed in anticipation of relatively uniform high resistivity basalts, and were found to have been seriously distorted by the intrabasalt conductors discovered in the field. In particular, the resistivity sections derived from 1D inversions were found to be inconsistent and misleading. The EMAP survey provided the most information about the subsurface resistivity distribution, and was certainly the most cost-effective. However, both controlled-source and EMAP surveys call for accurate 2D or 3D inversion to accommodate the geological objectives of this project.