Imaging fracture distributions of the Al-Khuff Formation outcrops using GPR and ERT geophysical techniques,Al-Qassim area,Saudi Arabia

Khuff Formation is of utmost importance in Saudi Arabia for oil and gas reservoir although it is composed mainly of limestone. This reason refers to the existence of intensive fractures which play a vital role in the increase in porosity and permeability of this formation. The fracture pattern in the study area was verified through 2D and 3D ground penetrating radar (GPR)-defined and electrical resistivity tomography (ERT)-defined surveys. In this respect, ten of 2D GPR surveys were collected along an intersected grid of profiles covering the study area while ERT data were collected along three profiles of the GPR grid. The results were interpreted in light of the field-based structural and stratigraphic assessment of the outcropping rocks. The analysis of the inverted ERT and filtered GPR sections revealed the presence of fractures. Several resistivity and electromagnetic reflection anomalies were laterally and vertically identified across the measured sections clarifying fractures that extend to a depth of 24 m in the limestone. Most fractures are oriented vertical to sub-vertical dipping both east-west and north-south.     

3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal

Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape and dimension of contamination associated with landfills.     

Electrical resistivity tomography (ERT) of a coastal carbonate aquifer (Port-Miou, SE France)

As part of a larger regional research program “KarstEAU”, the authors have applied electrical resistivity tomography (ERT) techniques to characterize heterogeneities in the Port-Miou coastal karst aquifer (Cassis, SE France). Field surveys were carried out on intensely fractured and karstified Urgonian carbonates. Extensive research has characterized macro- and micro-scale geology of the Port-Miou area and particularly underground water-filled conduits and fault/fracture and karst systems within a former quarry. The authors applied 2D ERT along two surface profiles of length 420 and 595 m to test capability for delineating subsurface conduits and possibly relationship between conduit and fault/fracture/karst orientation; and 3D ERT with a dense 120 electrode array at 1 m spacing (11 × 10 m) was applied over an area of the quarry that had been profiled using 3D georadar and which has had intensive nearby structural geological interpretation. The 2D profiling imaged several underground conduits at depths to >50 m below ground surface and below sea level, including possibly the main Port Miou submarine spring and smaller springs. The 2D profiling within the quarry provided a better understanding of the connectivity between major fractures and faults on the quarry walls and secondary springs along the coast supporting flow of the secondary springs along interpreted fracture orientations. In addition, 2D inversion-derived conductivity models indicate that high resistivity zones above sea-level are associated with non-saturated zones and low resistivity anomalies in the non-saturated zone are associated with residual clays in paleokarsts. A partitioned lower resistivity zone below sea-level can be associated with a higher porosity/permeability zone with fractures and karstic features. Inversion models of the dense 3D ERT data indicate a higher resistivity volume within the larger surveyed block. The survey characterized the non-saturated zone and the ERT resistivities are correlated with karst features interpreted by 3D georadar and visible in the inferior wall of the quarry.     

Evolution of the Moxizhuang Oil Field, Central Junggar Basin, Northwest China

Current oil saturation in the Moxizhuang (莫西庄) Oil Field in central Janggar (准噶尔) basin was evaluated by logging interpretation and measured on core samples, and the paleo-oil saturation in both the pay zones and water zones was investigated by graln-containing-oil inclusion (GOI) analysis.The pay zones in this field have low oil saturation and display low resistivity and small contrast between pay zones and water zones, and are classified as low-porosity, low oil saturation, and low resistivity reservoirs. Both the current low oil-saturation pay zones and the water zones above 4 365 m have high GOI values (up to 38%), suggesting high paleo-oil saturation. The significant difference between current oil saturation from both logging interpretation and core sample measurement and paleo-oil saturation indicated by GOI analysis suggests that this low oil-saturation field evolved from a high oil-saturation pool. Lateral re-migration and spill of formally trapped oil owing to changes in structural configuration since Neogene was the most plausible mechanism for oil loss in the Moxizhuang Oil Field.The combined effects of differential accumulation in the charge phase and the differential re-migration and spill of accumulated oil in Neogene are responsible for the complicated correlation between residual oil saturation and porosity/permeability of the reservoir sandstones and the distribution of low oil-saturation pay zones and paleo-oil zones (current water zones).     

Geoelectrical and hydrogeochemical studies for delineating seawater intrusion in the outlet of Wadi Ham, UAE

The Quaternary aquifer of Wadi Ham, UAE has been overexploited during the last two decades to meet the increasing water demands. As a result, the dynamic balance between freshwater and seawater has been disturbed and the quality of the groundwater has deteriorated. In this paper, a 2D earth resistivity survey was conducted in Wadi Ham in the area between Fujairah and Kalba to delineate the seawater intrusion. Existing monitoring wells were used to measure the horizontal and vertical variations in water salinity and thus to improve the interpretation of earth resistivity imaging data. Results of vertical electrical soundings and chemical analyses of collected water samples were used to obtain an empirical relationship between the inferred earth resistivity and the amount of total dissolved solids. This relationship was used along with the true resistivity sections resulting from the inversion of 2D resistivity data to identify three zones of water-bearing formation (fresh, brackish, and salt-water zones). Along the four 2D resistivity profiles, the depth to the fresh-brackish interface exceeded 50 m at the western part of the area and was in the order of 10 m or less in the eastern side near the shoreline. Depth to the brackish-saline water interface reached about 70 m in the western side and was in the order of 20 m in the eastern side. The thickness of the fresh water zone decreases considerably in the farming areas toward Kalba and thus the degree of seawater intrusion increases.     

Coastal-zone shallow-aquifer characterization study using geoelectrical and geochemical methods in Zirconium Complex,Atomic Energy,Pazhayakayal, Thoothukudi,India

Groundwater is a very important component of water resources in coastal aquifers in Thoothukudi. It has been established that the groundwaters in the coastal zone of the Zirconium Complex, Pazhayakayal, Thoothukudi district, Tamilnadu, India, are subjected to wave and tidal impact. The groundwater quality was studied by hydrogeological methods, 2D electrical resistivity imaging (ERI) techniques (six profiles), 11 Wenner vertical electrical sounding (VES), and well log analysis. Also, nine geochemical water samples were taken from the study area. The 2D ERI and VES surveys were carried out using WGMD-4 Ltd., Chennai, resistivity meter, multicore cable, and multielectrodes with Wenner array. The collected resistivity data were interpreted using the Res2DINV software. The research shows that the groundwaters are the result of the paleoriver flow along the Tamirabarani Channel in the western area and of the seawater intrusion in the eastern area. The fresh water is characterized by resistivity of about 10–100 Ohm?m in the study area. The resistivity of 10–50 Ohm?m indicates that the subsurface section is made up of sand, clay, and caliche. Resistivity values of more than 200 Ohm?m are specific to sand dunes. The very low resistivity (<5 Ohm?m) layer might be due to the seawater intrusion in the study area. Six water samples from the well drilled in the coastal area were analyzed, which made it possible to determine the concentrations of major and trace elements in the groundwaters. These data were used to establish the seawater intrusion and coastal environment characteristics in the study area.     

Use of time-domain electromagnetic (TDEM) method to investigate seawater intrusion in the Lebna coastal aquifer of eastern Cap Bon,Tunisia

A multi-method geophysical survey, using time-domain electromagnetic (TDEM) fields, is conducted on the north-eastern coastal aquifer of Cap Bon, Tunisia. The aim of the study is to investigate the inland extent saltwater intrusion, to estimate the fresh-saltwater transition zone, and to characterize the lithologic spatial variability of the underground layers. Moreover, the TDEM is used to determine both fluid resistivity and bulk resistivity which are required to calculate the porosity by applying the Archie’s equation. The survey of TDEM profiles is carried out along the left and right overbanks of the Lebna river from the coast until the Lebna dam. The interpretation and the analysis of these profiles allow us to better identify the flux exchanges between the river and the underlying aquifer. On the basis of the 1D and 2D longitudinal and transversal resistivity profiles as well as the existing boreholes data and groundwater quality, it seems that the freshwater-saltwater interface is moving forward inland. Furthermore, electromagnetic results combined to the Archie’s law show that the porosity values are at the ranges of 12–36 and 5–31% for the Quaternary outcropping and Pliocene formation, respectively.     

Locating oil spill response equipment and facilities: A cluster analysis approach

The response to marine accidental spillage of oil requires careful advance planning to ensure that the impact of an oil spill is minimized. In this paper a binary integer mathematical model based on the cluster analysis approach is developed to aid in determining the optimal sites for locating oil spill response equipment and facilities within a semienclosed waterway. The application of the model has been illustrated through a numerical example using hypothetical but representative data to determine optimal sites for locating oil spill response equipment and facilities in the Arabian Gulf region.     

Application of 2D and 3D geoelectrical resistivity imaging for engineering site investigation in a crystalline basement terrain,southwestern Nigeria

Orthogonal set of 2D geoelectrical resistivity field data, consisting of six parallel and five perpendicular profiles, were collected in an investigation site using the conventional Wenner array. Seven Schlumberger soundings were also conducted on the site to provide ID layering information and supplement the orthogonal 2D profiles. The observed 2D apparent resistivity data were first processed individually and then collated into 3D data set which was processed using a 3D inversion code. The 3D model resistivity images obtained from the inversion are presented as horizontal depth slices. Some distortions observed in the 2D images from the inversion of the 2D profiles are not observed in the 2D images extracted from the 3D inversion. The survey was conducted with the aim of investigating the degree of weathering and fracturing in the weathered profile, and thereby ascertaining the suitability of the site for engineering constructions as well as determining its groundwater potential.     

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.     

Electrical resistivity mapping of oil spills in a coastal environment of Lagos,Nigeria

Three years after the oil spillage and pipeline explosion that claimed about 100 human lives at Ijegun Community of Lagos–Nigeria, a combination of carefully designed 2D Electrical Resistivity Profilling and Vertical Electrical Sounding methods was deployed to map and characterise the subsurface around the contaminated site. Data acquired were processed, forward modelled and tomographically inverted to obtain the multi-dimensional resistivity distribution of subsurface. The results of the study revealed high resistivity structures that indocate the presence of contaminant (oil plumes) of different sizes and shapes around the oil leakage site. These high resistivity structures are absent in the tomograms and resistivity-depth slices computed for Iyana—a linear settlement not affected by oil spillage. The five geo-electric layers and the resistivities delineated in the area are the top soil layer, 220–670 Ωm; clayey sand layer, 300–1072 Ωm; top sand layer, 120–328 Ωm; mudstone/shale layer, 25–116 Ωm and the bottom sand layer, 15–69 Ωm. The base of the first four geo-electric layers corresponds to 3.9, 8.4, 27.2 and 34.6 m respectively. The two groundwater aquifers delineated correspond to the third and fifth geo-electric layers. The top aquifer has been infiltrated by oil plumes. The depth penetrated by the oil plume decreases from 32 m to about 24 m across the survey profiles from the two ends. It was concluded that the contaminant plumes from the oil spillage are yet to be completely degraded as at the time of the study. It is recommended that the contaminated site be remediated to remove or reduce the contaminant oil in the subsurface.     

Investigation of a shallow alluvial aquifer using geoelectrical methods: a case from Turkey

Near-surface geophysical methods are commonly used to solve a wide class of geological, engineering and environmental problems. In this study, a geoelectrical survey was performed to investigate an alluvial aquifer. The study area is located in the southwest of the Çubukluda? graben, situated in the south of ?zmir, Turkey. The geophysical studies included the electrical resistivity imaging and self-potential (SP) methods. The resistivity data were acquired along eight profiles in the northern part of the study area by a Wenner-Schlumberger electrode configuration and the data processing was achieved by a tomographic inversion technique. The SP data were collected by gradient technique along 16 profiles. Total field values were calculated for each profile by addition of the successive gradient values, then a total field SP map was obtained. The water-saturated zone in the northern part of the study area was clearly revealed by the electrical resistivity imaging and the SP survey yielded useful information on the subsurface fluid movement.     

Geoelectrical survey over perched aquifers in the northern part of Upper Sakarya River Basin,Türkiye

In this study, a groundwater exploration survey was conducted using the DC Resistivity (DCR) method in a hydrogeological setting containing a perched aquifer. DCR data were gathered and an electrical tomography section was recovered using conventional four-electrode instruments with a Schlumberger array and a two-dimensional (2D) inversion scheme. The proposed scheme was tested over a synthetic three-dimensional (3D) subsurface model before deploying it in a field situation. The proposed method indicated that gathering data with simple four-electrode instruments at stations along a line and 2D inversion of datasets at multiple stations can recover depth intervals of the studied aquifer in the hydrogeological setting even if it has a 3D structure. In this study, 2D inversion of parallel profiles formed a pseudo-3D volume of the subsurface resistivity structures and mapped out multiple resistive (>25 ohm·m) bodies at shallow (between 50–100 m) and deep sections (>150 m). In general, the proposed method is convenient to encounter geological units that have limited vertical and spatial extensions in any direction and presents resistivity contrast from groundwater-bearing geologic materials.     

Imaging of wastewater percolation in heterogeneous soil using electrical resistivity tomography (ERT): a case study at east of Tenth of Ramadan City,Egypt

The objective of this study is to evaluate the effectiveness of DC resistivity surveys for imaging the wastewater percolation around the stabilization ponds in the Tenth of Ramadan City, the desert fringes of East Nile Delta, Egypt. Detailed resistivity surveys, including DC soundings and electrical resistivity tomography (ERT), were carried out along several profiles. Furthermore, synthetic modeling of ERT was designed to optimize the survey configurations and interpretation of the results. A 2D modeling of smoothness-constrained least-squares inversion scheme was applied to delineate the possible wastewater infiltration zones from oxidation ponds. Because the geoelectrical interpretation has a degree of non-uniqueness, the resistivity inversion was constrained using borehole lithological information and soil sample laboratory measurements. The DC inversion results indicate decreasing resistivity down to a depth of 15 m around waste disposal sites. The inferred soil zone close to the oxidation ponds was a mixture of sand, silt and clay. Moreover, the clay minerals were characterized by moderate swelling that could have reduced the vertical infiltration speed, causing wastewater seepage, especially around unlined disposal sites and open surface drains. Accordingly, the medium-to-low resistivity values can be attributed to wastewater leakage in clayey sand soil. Because the area slopes generally toward the northeast, the surface seepage was dominant in the shallow impermeable sandy clay subsoil. Therefore, measuring soil parameters is a complementary method to optimize resistivity interpretation, with potential for mitigating environmental hazards from wastewater leakage around disposal ponds.     

Delineation of shallow resistivity structure in the city of Burdur,SW Turkey by vertical electrical sounding measurements

The city of Burdur, which is built on an alluvium aquifer, is located in one of the most seismically active zones in southwestern Turkey. The soil properties in the study site are characterized by unconsolidated and water-saturated sediments including silty, clayey and sandy units, and shallow groundwater level is the other characteristic of the site. Thus, the city is under soil liquefaction risk during a large earthquake. A resistivity survey including 189 vertical electrical sounding (VES) measurements was carried out in 2000 as part of a multi-disciplinary project aiming to investigate settlement properties in Burdur city and its vicinity. In the present study, the VES data acquired by using a Schlumberger array were re-processed with 1D and 2D inversion techniques to determine liquefaction potential in the study site. The results of some 1D interpretations were compared to the data from several wells drilled during the project. Also, the groundwater level map that was previously obtained by hydrological studies was extended toward north by using the resistivity data. 2D least-squares inversions were performed along nine VES profiles. This provided very useful information on vertical and horizontal extends of geologic units and water content in the subsurface. The study area is characterized by low resistivity distribution (<150 Ωm) originating from high fluid content in the subsurface. Lower resistivity (3–30 Ωm) is associated with the Quaternary and the Tertiary lacustrine sediments while relatively high resistivity (40–150 Ωm) is related to the Quaternary alluvial cone deposits. This study has also shown that the resistivity measurements are useful in the estimation of liquefaction risk in a site by providing information on the groundwater level and the fluid content in the subsurface. Based on this, we obtained a liquefaction hazard map for the study area. The liquefaction potential was classified by considering the resistivity distributions from 2D inversion of the VES profiles, the types of the sediments and the extended groundwater level map. According to this map, the study area was characterized by high liquefaction hazard risk.     

Geoelectrical investigation of oil contaminated soils in former underground fuel base: Borne Sulinowo,NW Poland

The survey has been carried out in the area of 0.23 km² of the former military underground fuel base. The oil derivative products were observed in excavations and the laboratory tests confirmed the occurrence of hydrocarbons (>C12) in soils. The purpose of the survey was to determine the spatial extent of the contamination. The studied area is covered by postglacial sediments: sands, gravels and till. The first water table was observed at a depth of 10–12 m. The detailed electromagnetic measurements with Geonics EM31-MK2 conductivity meter were performed in the whole area of the former fuel base. Obtained results were elaborated statistically and the map of apparent electrical conductivity to a depth of 6 m was created. Many local low conductivity anomalies were observed. The measurements with Geonics EM34-3XL were performed along one A–A′ profile and 1D electromagnetic modelling along with this profile was calculated to obtain the electrical conductivity cross-section to a depth of 30 m. Two-dimensional electrical resistivity imaging measurements were carried out along the same profile and the resistivity cross-section to a depth of 20 m was performed. Both conducivity and resistivity cross-sections show anomalous zones. The zones correlate with oil contaminated zones very well.     

兴蒙、吉黑地区岩石圈电性结构特征

采用具有国际先进水平的二维连续自动反演技术对收集的MT资料进行了二维反演和综合解释。得到了不同地区MT剖面的二维电性结构断面。通过二维反演断面发现了松辽盆地的东西边界较陡,而南部边界较缓的电性结构;发现在火山地区12～30km深度普遍存在低阻物体,为解释火山地区深部存在岩浆囊提供了依据。另外,于桦南—饶河重新实施的MT剖面,其二维反演断面也较清楚地揭示出:剖面西段为高阻特征且具有稳定的岩石圈厚度(80～90km),该区段对应佳木斯地块;剖面中部具有明显的电性梯度带,该梯度带为确定佳木斯地块的东界位置及其深部结构形态提供了依据;剖面东段则揭示了佳木斯地块以东地区浅部为逆冲推覆体,深部为多个高阻块体与低阻条带相间的电性结构。     

Hydrogeophysical investigations in a typical khondalitic terrain to delineate the kaolinised layer using resistivity imaging

Investigation for high yielding water wells in the khondalitic terrain (graneti ferrous silliminite gnesiss) is mostly faced with the problem of identification of the extent of the depth of kaolinisation of the aquifer. The traditional Vertical Electrical Sounding survey, Seismic Refraction survey and Very Low Frequency Electromagnetic survey could not identify the kaolinisation of the aquifer in the present investigations. The Two Dimensional (2D) Resistivity and Induced Polarization (IP) Imaging surveys are attempted for the identification of kaolinised layer and depth of kaolinisation. Number of 2D Resistivity and IP Imaging profiles were conducted near Chipurupally in Vizianagaram district of Andhra Pradesh, India along successful and failed wells located within short distances. Resistivity and IP measurements were carried out using an ABEM SAS 4000 Terrameter. The resistivity and I.P. images have provided a clear view of the thickness of the highly weathered zone (kaolinised zone) at both successful and failed wells. The highly weathered zone is identified with the resistivity values below 25 ohm.m. The depth of highly weathered material at failed well is extended about 8–10 m more deeper than the successful wells at some places to as much as 20 m more deep at some other places. This extended deeper kaolinisation of the aquifer is responsible for failure of wells. Layers having resistivities between 25–65 Ohm.m are identified as aquifer layers which are composed of moderately weathered and fractured khondalitic suit of rocks (Garnti ferrous sillimanite/biotite gneiss). Layers with resistivities greater than 65 Ohm.m are interpreted to have basement characteristics belonging to the granite gneiss. Interestingly IP imaging has not provided any greater insights in delineating the kaolinistion of the aquifer when compared to resistivity Imaging, in fact resitivity imaging has shown greater depths of kaolinisation than IP Imaging.     

古龙南地区低阻油层胶结指数和饱和度指数影响因素实验

采用变胶结指数m和饱和度指数n的阿尔奇公式仍是目前定量评价低阻油层饱和度的主要方法之一。然而,由于不同地区低阻油层的成因不同,从而使得低阻油层m和n的影响因素随地区不同而发生变化。因此,针对不同地区低阻油层,有必要研究其阿尔奇公式m和n的影响因素及变化规律。依据古龙南地区葡萄花低电阻率油层划分标准以及考虑岩性、物性、电性变化,选取了低阻油层岩样,设计了岩样岩石物理实验。基于岩石物理实验数据,绘制m和n与物性参数、孔隙结构参数、黏土附加导电能力参数的交会图。经分析可知：低阻油层的m和n值小于常规油层的m和n值,并随孔渗综合指数、中值半径、平均半径、半径均值的增大而增大。低阻油层的m值随泥质体积分数、阳离子交换容量、微孔隙体积与可动流体体积之比的增大而减小,随核磁共振测量的T₂几何平均值的增加而增大;n值随孔隙度、孔隙度和泥质体积分数比值的增加而增大。低阻油层的m值与T₂几何平均值、微孔隙体积与可动流体体积之比关系较好,n值与孔渗综合指数、平均半径、最大半径、半径均值的相关性非常好。研究表明,低阻油层的m值受孔隙结构和泥的影响较大,而n受孔隙结构影响较大。可用T₂几何平均值或微孔隙体积与可动流体体积之比确定低阻油层的m值,可用孔渗综合指数、平均半径、最大半径、半径均值确定低阻油层的n值。     

The buried southern continuation of the Oaxaca-Juarez terrane boundary and Oaxaca Fault,southern Mexico: Magnetotelluric constraints

Thirty magnetotelluric soundings were made along two NW–SE profiles to the north and south of Oaxaca City in southern Mexico. The profiles crossed the N–S Oaxaca Fault and the Oaxaca-Juarez terrane boundary defined by the Juarez mylonitic complex. Dimensionality analysis of the MT data showed that the subsurface resistivity structure is 2D or 3D. The Oaxaca and correlative Guichicovi terranes consist of ca. 1–1.4 Ga granulitic continental crust overlain by Phanerozoic sedimentary rocks, characterized by high and low resistivities, respectively. The Juarez terrane consists of oceanic Mesozoic metavolcanic and metasedimentary rocks, characterized by a low to medium resistivity layer, that is approximately 10 km thick. The Oaxaca Fault is a Cenozoic aged, normal fault that reactivated the dextral and thrust Juarez mylonitic complex north of Oaxaca City: its location south of Oaxaca City is uncertain. In the southern profile, the MT data show a ca. 20–50 km wide, west-dipping, relatively low resistivity zone material that extends through the entire crust. This is inferred to be the Juarez terrane bounded on either side by the ca. 1–1.4 Ga granulites. The Oaxaca Fault is imaged only by a major electrical resistivity discontinuity (low to the west, high to the east) along both the western border of the Juarez mylonitic complex (northern profile) and the San Miguel de la Cal mountains (southern profile) suggesting continuity.