Resistivity imaging survey to delineate subsurface seepage of hydrocarbon contaminated water at Karbala Governorate,Iraq

This research is an attempt to accomplish a 3-D resistivity imaging survey, which was carried out near a water well contaminated with hydrocarbon materials in Karbala governorate. Two-dimensional resistivity imaging measurements were collected along four parallel profiles, using a Wenner array with electrode spacing of 1 m. The RES3DINV program was used to invert the apparent resistivity data. The results displayed a resistivity distribution of the subsurface in a three-dimensional volume. Thus, both the horizontal and vertical extents of the contaminated zone were displayed. This technique revealed a low resistivity zone at depth ranges from 3 to 6 m in the investigation area, but the seepage starts at depth ranges between 2 and 3 m and continues down depth (may be to the groundwater level). This low resistivity zone is the most likely location for a subsurface seepage of contaminated water. It is clear that the sufficient measurement points along 2-D lines in a small area can increase the 3-D imaging resolution, and nearly real 3-D imaging can be achieved, when the size of subsurface anomaly compared with the electrode spacing (a) of the Wenner array is taken into consideration.     

http://www.sciencedirect.com/science/article/pii/S1674987112001545

The presence of natural voids and cavities in subsurface karstic limestones causes severe problem for civil engineering and environmental management.The presence of such features hinders the extension of urbanizations particularly in the new metropolitan.The eastern part of Saudi Arabia contains various types of karstic limestone,sinkholes,solution cavities and voids.In this context,geophysical methods particularly electrical resistivity technique is used as a cost-effective solution for investigating subsurface caves,voids,and shallow weathered zones.2-D electrical resistivity data sets have been acquired along seven profiles in the new urbanization at Al Hassa area.Data processing has been carried out taking into consideration the response of synthetic models,which simulates physical models of the most common karstic features in the area.The results are very useful to determine the extension of shallow weathered zones and to locate different cavities underneath them.The hard limestone bedrock can also be detected and traced along the surveyed profiles.     

高密度电阻率法二维成像方法及其应用

高密度电阻率法在工程勘察和环境监测中是一种很重要的勘探手段，高密度电阻率法二维成像理论与方法的研究在理论和实际应用上都有着很重要的意义。随着地球物理反演和成像技术研究的兴起和发展，层析成像技术在高密度电阻率法中的应用，提高了高密度电阻率法的勘探质量和精度。本文献电阻率法层析成像技术的前提条件及高密度电阻率法二维成像理论与方法进行了粗浅探讨，最后给出了一些应用实例。     

地下溶洞调查物探技术探讨

本文依据岩溶区岩土体的物性差异,利用波阻抗和电阻率差异,对比分析了高密度电法、瞬变电磁法、地质雷达法等几种物探方法,最终确定了高密度电法、地震映像法及声波测井和钻孔CT法等地表和地下物探手段相结合的方法来探测岩溶分布与发育情况。     

台湾中西部沉积盆地三维密度结构及其重力异常

台湾中西部沉积盆地是台湾省海陆油气勘探的重点区域。本文利用石油钻井资料和重力观测数据,依据Kriging插值及台湾中国石油公司深度-密度经验公式,计算了该地区沉积地层三维密度及剩余密度结构,并依据三维有限元解泊松方程方法,通过求解重力位计算出三维剩余密度体产生的剩余重力异常。密度结构显示,在0~6km深度范围,自更新统到先中新统的沉积地层大致分为2.0~2.30g/cm3、2.30~2.45g/cm3、2.45~2.60g/cm3、2.60~2.70g/cm3等4个密度层。其中,台中盆地陆缘区中下部2~4km深度存在较大沉积凹陷,凹陷区密度比周围低0.05~0.10g/cm3.重力计算结果表明,本区沉积地层产生的剩余重力异常在-20~+15mGal之间,深部凹陷对应负异常低值圈闭,其形态和量值与实测重力异常一致。      

三维地电断面激发极化法有限元数值模拟

用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.      

色拉哈断裂及邻区音频大地电磁三维阵列探测

为了能在川藏铁路色拉哈段隧道选址过程中最大限度地规避地质灾害体,本文针对色拉哈断裂及邻区开展了音频大地电磁三维阵列探测研究,获得了研究区可靠的三维精细电性结构.结合地表活动断裂调查结果,对区内三维电性分布特征进行了综合解释.研究表明,色拉哈断裂及邻区地表至500 m,整体表现为高阻特征,局部受断裂控制的区域为低阻特征;500 m以下电阻率下降显著,推测是断裂带在500 m以下的深部构造富含水体所致;此外,位于色拉哈断裂北侧的木格措南阶区以及南侧的断裂交汇区东侧发育两处高导体C1和C2,推测为区内含水的构造软岩,在设计隧道路线的过程中,应尽量避让.      

高密度电法在工程勘察中的应用实例

介绍了高密度电法在覆盖层勘探、滑坡体勘探、土料场分层、采空区勘探、断层调查的几个应用实例。经工程实例表明,高密度电法可在工程勘察中取得良好的应用效果。     

Geophysical assessment of subsurface coastal sediments and their engineering implications

Surface geophysics and a priori information were employed to delineate the subsurface geology at Idi-oro in Abijo, Ibeju Lekki area of Lagos, Nigeria for foundation investigation purpose. Resistivity measurement was conducted using 1-D and 2-D resistivity probing techniques. The resistivity measurements were made with ABEM tetrameter model SAS 1000 system. The 1-D vertical electrical resistivity sounding data were obtained using the Schlumberger electrode array while the 2-D resistivity data were obtained using the dipole–dipole array. The interpreted results revealed three to five subsurface geological layers. This is made up of the top soil with resistivity values that vary from 132.4 to over 2,313.5 Ω?m and thickness values that range from 0.3 to 4.8 m, the fine sand with resistivity values that vary from 221.0 to 3,032.7 Ω?m and thickness values that range from 0.4 to 5.5 m, the medium sand with resistivity values that vary from 202.8 to 1,247.7 Ω?m and thickness values that range from 4.9 to 58.4 m. On the other hand, the clayey sand has the resistivity values that vary from 146.1 to 1,744.0 Ω?m and thickness values that vary from 2.2 to 26.3 m, while the coarse sand has resistivity values that vary from 238.3 to 14,313.9 Ω?m but with no thickness value because the current terminated in this layer. The resistivity data correlate well with borehole logs. On the whole, it is concluded that the investigated area has competent sand layer that can support medium to giant engineering structures with resistivity values that vary from 202 to 14,314 Ω?m and thickness values that vary from 0.8 to 58.4 m.     

The application of GPR and electrical resistivity tomography as useful tools in detection of sinkholes in the Cheria Basin (northeast of Algeria)

Sinkhole collapse is one of the main limitations in the development of karst areas, especially where bedrock is covered by unconsolidated material. Studies of sinkhole formation have shown that sinkholes are likely to develop in cutter (enlarged joint) zones as a result of subterranean erosion by flowing groundwater. Ground-penetrating radar (GPR) and electrical resistivity imaging or tomography (RESTOM) are well suited to mapping sinkholes because of the ability of these two techniques for detecting voids and discriminating subtle resistivity variations. Nine GPR profiles and two-dimensional electrical resistivity tomography have been applied, with relative success, to locate paleo-collapses and cavities, and to detect and characterize karst at two sinkhole sites near Cheria City where limestone is covered by about 10 m of clayey soils. The survey results suggest that GPR and RESTOM are ideal geophysical tools to aid in the detection and monitoring of sinkholes and other subsurface cavities.     

Estimation of 2-D and 3-D Fracture Densities from 1-D Data - Experimental and Field Results

2-D and 3-D densities of fractures are commonly used in mining safety design, natural gas and oil production in fractured reservoirs, and the characterization of subsurface flow and transportation systems in fractured rocks. However, many field data sets are collected in 1-D frequency (f) (e.g., scanlines and borehole data). We have developed an ARC/ INFO-based technology to calculate fracture frequency and densities for a given fracture network. A series of numerical simulations are performed in order to determine the optimal orientation of a scanline, along which the maximum fracture frequency of a fracture network can be obtained. We calculated the frequency (f) and densities (both D1 and D2) of 36 natural fracture trace maps, and investigated the statistical relationship between fracture frequency and fracture density D1, i.e. D1=1.340f+ 0.034. We derived analytical solutions for converting dimensional density (D1) to non-dimensional densities (D2 and D3) assuming that fracture length distribution f     

Hydrogeophysical investigation and estimation of groundwater potentials of the Lower Palaeozoic to Precambrian crystalline basement rocks in Keffi area, north-central Nigeria, using resistivity methods

The use of resistivity sounding and two-dimensional (2-D) resistivity imaging was investigated with the aim of delineating and estimating the groundwater potential in Keffi area. Rock types identified are mainly gneisses and granites. Twenty-five resistivity soundings employing the Schlumberger electrode array were conducted across the area. Resistivity sounding data obtained were interpreted using partial curve matching approach and 1-D inversion algorithm, RESIST version 1.0. The 2-D resistivity imaging was also carried out along two traverses using dipole–dipole array, and the data obtained were subjected to finite element method modeling using DIPRO inversion algorithm to produce a two-dimensional subsurface geological model. Interpretation of results showed three to four geoelectrical layers. Layer thickness values were generally less than 2 m for collapsed zone, and ranged from 5 to 30 m for weathered bedrock (saprolite). Two major aquifer units, namely weathered bedrock (saprolite) aquifer and fractured bedrock (saprock) aquifer, have been delineated with the latter usually occurring beneath the former in most areas. Aquifer potentials in the area were estimated using simple schemes that involved the use of three geoelectrical parameters, namely: depth to fresh bedrock, weathered bedrock (saprolite) resistivity and fractured bedrock (saprock) resistivity. The assessment delineated the area into prospective high, medium and low groundwater potential zones.     

Electrical resistivity investigation of potential cavities underlying a proposed ash disposal area

Measurements of electrical resistivity anomalies using a Wenner electrode array in a profiling mode were conducted over a 105-acre area near Crystal River, Florida Analyzed results of collected data were used to determine the likelihood of subsurface solution cavities beneath a proposed ash disposal area “Sounding” measurements with electrode spacing values extending up to 100 feet and Lee electrode measurements were used to identify and locate specific anomaly maxima identified anomalies on the order of 200 chms-feet contrasted with expected values of 300 to 400 ohms-feet and were attributed to lower resistivities inherent in water-filled solution cavities below a shallow water table Prominent cavity complexes at 5 to 15 feet depth were predicted in areas of low resistivity anomalies in the otherwise geoelectrically uniform, Eocene Ocala limestone Shallow test borings confirmed the presence of a significant concentration of shallow cavities in the predicted areas This study demonstrated the feasibility of electrical resistivity measurements as a reconnaissance tool to economically determine subsurface cavities and to aid in developmental and environmental planning     

高密度电法勘探数据三维处理技术及城市勘察中的应用

在以往的高密度数据反演解释过程中,将地质结构简化成二维,通过地电断面图、切片图进行分析.实际上,地质体具有三维地电结构,二维反演难免会受到旁侧效应的影响,二维图件也不能清晰直观地反映地质体的空间分布情况.本文以某陶瓷厂空洞探测为例,利用实测资料构建三维数据场,在反演参数优化的基础上,选择基于最小二乘法的三维反演技术对数据场进行处理,将反演结果利用三维可视化技术显示出来,对采空区位置及空间形态给出结论.研究结果显示,三维处理技术能够从空间上,多手段、多视角分析物探数据,有效地消除了旁侧效应的影响; 物探结果直观清晰地呈现出空洞的分布位置、范围以及空间特征,与钻探结果吻合较好.     

The 1.80–1.74-Ga gabbro–anorthosite–rapakivi Korosten Pluton in the Ukrainian Shield: a 3-D geophysical reconstruction of deep structure

The deep structure of the gabbro–anorthosite–rapakivi granite (“AMCG-type”) Korosten Pluton (KP) in the northwestern Ukrainian Shield was studied by 3-D modelling of the gravity and magnetic fields together with previous seismic data. The KP occupies an area of ca. 12,500 km² and comprises several layered gabbro-anorthositic intrusions enveloped by large volumes of rapakivi-type granitoids. Between 1.80 and 1.74 Ga, the emplacement of mafic and associated granitoid melts took place in several pulses. The 3-D geophysical reconstruction included: (a) modelling of the density distribution in the crust using the observed Bouguer anomaly field constrained by seismic data on Moho depth, and (b) modelling of the magnetic anomaly field in order to outline rock domains of various magnetisation, size and shape in the upper and lower crust. The density modelling was referred to three depth levels of 0 to 5, 5 to 18, and 18 km to Moho, respectively. The 3-D reconstruction demonstrates close links between the subsurface geology of the KP and the structure of the lower crust. The existence of a non-magnetic body with anomalously high seismic velocity and density is documented. Most plausibly, it represents a gabbroic stock (a parent magma chamber) with a vertical extent of ca. 20 km, penetrating the entire lower crust. This stock has a half-cylindrical shape and a diameter of ca. 90 km. It appears to be connected with a crust–mantle transitional lens previously discovered by EUROBRIDGE seismic profiling. The position of the stock relative to the subsurface outlines of the KP is somewhat asymmetric. This may be due to a connection between the magmatism and sets of opposite-dipping faults initially developed during late Palaeoproterozoic collisional deformation in the Sarmatian crustal segment. Continuing movements and disturbances of the upper mantle and the lower crust during post-collisional tectonic events between 1.80 and 1.74 Ga may account for the long-lived, recurrent AMCG magmatism.     

利用实测重力垂直梯度反演长白山地区一剖面的深部构造

针对重力梯度高分辨率的特点,利用在长白山地区实测的重力垂直梯度数据,采用梯度空间参量图反演其深部构造。台阶模型试验表明: 重力梯度空间参量图能给出构造倾角和倾面的信息,结合重力梯度剖面和梯度空间参量图可以构建出地下构造的几何模型,进而对一些复杂构造进行解释。通过对比实测布格重力异常和实测重力梯度异常,重力梯度比重力异常的分辨率更高; 将梯度法应用到实测重力梯度数据的处理中,结果表明: 该方法对确定密度变化界面的水平位置和深度具有非常好的效果。     

天津市重力数据反演解释

重力数据反演是获取地下地质信息的有效手段。本文通过对天津市全区密度和电性参数的统计分析,了解天津市地层的物性分布规律。在对研究区布格重力异常数据采用归一化总水平导数垂向导数（NVDR-THDR）技术处理的基础上,对断裂构造进行了识别,对构造单元进行了划分。然后结合研究区地质资料,应用Geosoft软件中的密度界面反演模块对各构造单元内的密度界面起伏进行反演。并以GOCAD软件为平台,构建了天津市三维可视化地质模型,借助大地电磁（MT）测深点虚拟钻孔,实现了对主要地质界面和构造的标定。最后利用布格重力异常数据拟合同位置地质剖面,验证了模型的可靠性。应用上述反演成果,结合区内地质认识,推断宝坻断裂（F₁）、蓟运河断裂（F₂）、F₃、杨柳青断裂（F₄）、F₇在中生代后活动性较强,沧东断裂（F₆）和天津断裂（F₃₀）在古近纪后活动性较强。     

Integrated geophysical interpretation for the area located at the eastern part of Ismailia Canal, Greater Cairo, Egypt

The integrated geophysical interpretation for the different geophysical tools such as resistivity and gravity is usually used to define the structural elements, stratigraphic units, groundwater potentiality, and depth to the basement rocks. In the present work, gravity and resistivity data were utilized for detecting the groundwater aquifer and structural elements, as well as the upper and lower surfaces of the subsurface basaltic sheet in an area located at the eastern side of Ismailia Canal, northeastern Greater Cairo, Egypt. Two hundred and ten gravity stations were measured using an Autograv instrument through a grid pattern of 50?×?50 m. The different required corrections were carried out, such as drift, elevation, tide, and latitude corrections. The final corrected data represented by the Bouguer anomaly map were filtered using high- and low-pass filters into regional and residual gravity anomaly maps. The resulting residual gravity anomaly map was used for gravity modeling to calculate the depths to the upper and lower surfaces of the basaltic sheet. The resulting gravity models indicated that the depths to the upper surface of the basaltic sheet are ranged between 26 and 314 m, where the shallower depths were found around the southern and eastern parts. The depths to the lower surface of the basaltic sheet are varied from 86 to 338 m, and the thickness of the basaltic sheet is ranged from 24 to 127 m, where the biggest thicknesses were found around the southern and northern parts of the study area. Forty-two vertical electrical soundings (VES) were carried out using Schlumberger configuration with AB/2 spacings ranged from 1.5 to 500 m. 1D quantitative interpretation was carried out through manual and analytical interpretations. The VES data were also inverted assuming a 3D resistivity distribution. The results from the 3D resistivity inversion indicated that the subsurface section consists of sand, sandstone, and sandy–clays of Miocene deposits overlying the basalts. Such basaltic features (of Oligocene age) are underlain by Gabal Ahmar Formation of Oligocene deposits, which are composed of sand and sandstone. Therefore, two aquifers were deduced in the area. The first is the Miocene aquifer (shallower) and the other is the Oligocene aquifer (deeper).     

三维地电断面电阻率测深有限元数值模拟

用有限单元法进行了电导率分块均匀的三维点源电场电阻率测深的数值模拟.首先给出了三维构造中点源电场的边值问题、变分问题; 然后, 用有限单元法求解变分问题, 将区域剖分成六面体单元, 在单元中进行三线性函数插值, 将变分方程化为线性代数方程组; 最后解方程组, 得各节点的电位值, 进而计算出地表的视电阻率.对几例较典型的地电模型进行试算, 结果表明本方法是行之有效的      

激发极化数据二维有限元反演法在金矿勘探中的应用

给出电阻率、激发极化法数据的二维有限元反演方法和原理 ,并将此方法用于高龙金矿电法勘探资料解释中 ,解释结果与已知地质资料基本吻合。由此证明 ,利用双线性模型进行二维反演的有效性及其独特的优点。