3D probability tomography: theoretical developments and applications to high-resolution geophysical prospecting

The theory of the 3D multipole probability tomography method (3D GPT) to image source poles, dipoles, quadrupoles and octopoles, of a geophysical vector or scalar field dataset is developed. A geophysical dataset is assumed to be the response of an aggregation of poles, dipoles, quadrupoles and octopoles. These physical sources are used to reconstruct without a priori assumptions the most probable position and shape of the true geophysical buried sources, by determining the location of their centres and critical points of their boundaries, as corners, wedges and vertices. This theory, then, is adapted to the geoelectrical, gravity and self potential methods. A few synthetic examples using simple geometries and three field examples are discussed in order to demonstrate the notably enhanced resolution power of the new approach. At first, the application to a field example related to a dipole–dipole geoelectrical survey carried out in the archaeological park of Pompei is presented. The survey was finalised to recognize remains of the ancient Roman urban network including roads, squares and buildings, which were buried under the thick pyroclastic cover fallen during the 79 AD Vesuvius eruption. The revealed anomaly structures are ascribed to wellpreserved remnants of some aligned walls of Roman edifices, buried and partially destroyed by the 79 AD Vesuvius pyroclastic fall. Then, a field example related to a gravity survey carried out in the volcanic area of Mount Etna (Sicily, Italy) is presented, aimed at imaging as accurately as possible the differential mass density structure within the first few km of depth inside the volcanic apparatus. An assemblage of vertical prismatic blocks appears to be the most probable gravity model of the Etna apparatus within the first 5 km of depth below sea level. Finally, an experimental SP dataset collected in the Mt. Somma-Vesuvius volcanic district (Naples, Italy) is elaborated in order to define location and shape of the sources of two SP anomalies of opposite sign detected in the northwestern sector of the surveyed area. The modelled sources are interpreted as the polarization state induced by an intense hydrothermal convective flow mechanism within the volcanic apparatus, from the free surface down to about 3 km of depth b.s.l..

Geophysical prospecting for the layman,

cover-title.

Physical, Geochemical, Geomechanical, Geophysical and Hydrodynamic Characterization and Anisotropy Assessment of Corvio Sandstone

Corvio sandstone is a ~20 m thick unit (Corvio Formation) that appears in the top section of the Frontada Formation (Campoó Group; Lower Cretaceous) located in Northern Spain in the southern margin of the Basque-Cantabrian Basin. Up to 228 plugs were cored from four 0.3 x 0.2 x 0.5 m blocks of Corvio sandstone, to perform a comprehensive characterization of the physical, mineralogical, geomechanical, geophysical and hydrodynamic properties of this geological formation, and the anisotropic assessment of the most relevant parameters. Here we present the first data set obtained on 53 plugs which covers (i) basic physical and chemical properties including density, porosity, specific surface area and elementary analysis (XRF - CHNS); (ii) the curves obtained during unconfined and confined strengths tests, the tensile strengths, the calculated static elastic moduli and the characteristic stress levels describing the brittle behaviour of the rock; (iii) P- and S-wave velocities (and dynamic elastic moduli) and their respective attenuation factors Qp and Qs, electrical resistivity for a wide range of confining stress; and (iv) permeability and transport tracer tests. Furthermore, the geophysical, permeability and transport tests were additionally performed along the three main orthogonal directions of the original blocks, in order to complete a preliminary anisotropic assessment of the Corvio sandstone.

Cosa: arqueologia de un enclave comercial mediterráneo. Campaña de 2006

En el marco del proyecto “La ciutat romana de Cosa: arqueologia d’un enclau comercial mediterrani” , autorizado y apoyado por la Soprintendenza Archeologica por la Toscana, entre los dias 4 y 22 de septiembre de 2006 se ha realizado la segunda campaña de intervenciones arqueológicas en la ciudad romana de Cosa (Ansedonia, prov. Grosseto, Itàlia), colonia latina fundada en el 273 aC a unos 120 km. al norte de Roma. De acuerdo a los resultados obtenidos en la campanya del 2005 (localización/verificación de los límites precisos de la ínsula O-P/4-5 mediante la aplicación de técnicas de prospección geofísica completadas con la limpieza, registro y documentación arqueológica de las estructuras localizadas ) los trabajos del 2006 se han orientado hacia la identificación de la organización interna de la dicha ínsula tomando como referente el criptopórtico situado en el extremo N.E., el cual parece constituir el límite de una estructura singular (privada o pública ) estratégicamente ubicada en relación al fórum i a la Via Sacra. El trabajo de campo ha consistido en un intenso decapage con la finalidad de delimitar unidades de habitación complejas funcionalmente definidas y así la articulación existente entre ellas; en este sentido se ha podido documentar evidencias del espacio porticado superpuesto al criptopórtico así como, paralelamente a la calle 5 y en dirección a la Via Sacra, parte de habitaciones algunas de las cuales conservaban restos del pavimiento original, en un caso de mosaico. Paralelamente, se ha realizado el análisis en laboratorio de los materiales recuperados los cuales, aún procediendo de nivel superficial, empiezan a proporcionar datos sobre los diferentes momentos de ocupación de la zona, básicamente tardorepublicanos y augustales.

A new method for geoelectrical investigations underwater

A new electrical method is proposed for determining the apparent resistivity of multi-earth layers located underwater. The method is based on direct current geoelectric sounding principles. A layered earth model is used to simulate the stratigraphic target. The measurement array is of pole-pole type; it is located underwater and is orientated vertically. This particular electrode configuration is very useful when conventional electrical methods cannot be used, especially if the water depth becomes very important. The calculated apparent resistivity shows a substantial quality increase in the measured signal caused by the underwater targets, from which little or no response is measured using conventional surface electrode methods. In practice, however, different factors such as water stratification, underwater streams or meteorological conditions complicate the interpretation of the field results. A case study is presented, where field surveys carried out on Lake Geneva were interpreted using the calculated apparent resistivity master-curves.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2007

A prior long-term and complex evaluation of the already available data on the geophysical prospecting during the first season work carried out at 2006, at the archaeological site of Tchinguiz Tepe of Termez, took place to decide the strategy to follow during the campaign of 2007. This previous evaluation of the information, on one hand, leaded to the decision to increase the geophysical prospecting at Tchinguiz Tepe, on the other hand, to decide the exact location of areas where the archaeological interventions.would carry out. The main objective at the beginning of this new season was to crosscheck the reliabilityof the measurements and, at the same time, to establish the unknown up to the present archaeologicaland chronological sequence of Tchinguiz Tepe. Meanwhile, the geophysical prospecting also wasextended to the outskirts of the city were the localisation of an unknown up to now Buddhist Monasterywas possible.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2008

This year the development of our project can be divided into two main clearly different parts, on one hand the laboratory work, where the sampled ceramic individuals has been prepared and analyzed and the elaboration of the data obtained during the excavation of 2007 has been finished and, on the other hand, the field work developed at the archaeological site during this specific year (2008).In the mark of the analytical work a significant number of ceramic individuals (144) from the different stratigraphical units from various areas of the excavation of Termez and Tchinguiz Tepe sampled duringthe new and previous field works has been archaeometrically characterized. This specific material includedindividuals dated into the Hellenistic and Sassanian period, which has been confirmed by C14dating upon organic samples.At the same time, in the mark of the field work of 2008 the archaeological record, already started tobe under study during the excavation of 2007, has been completed and two new archaeological recordshave been registered on of which is located in the area of Tchinguiz Tepe. For the archeological studythe information of the previous geophysical prospecting has been indisputably taken into considerationand the same methodology has been applied to crosscheck the latter archaeological results.

Seismic attenuation due to wave-induced fluid flow at microscopic and macroscopic scales

Wave-induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro-cracks or broken grain contacts, causes the stiffness of the so-called modified dry frame to be complex-valued and frequency-dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave-induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2007

A prior long-term and complex evaluation of the already available data on the geophysical prospecting during the first season work carried out at 2006, at the archaeological site of Tchinguiz Tepe of Termez, took place to decide the strategy to follow during the campaign of 2007. This previous evaluation of the information, on one hand, leaded to the decision to increase the geophysical prospecting at Tchinguiz Tepe, on the other hand, to decide the exact location of areas where the archaeological interventions.would carry out. The main objective at the beginning of this new season was to crosscheck the reliabilityof the measurements and, at the same time, to establish the unknown up to the present archaeologicaland chronological sequence of Tchinguiz Tepe. Meanwhile, the geophysical prospecting also wasextended to the outskirts of the city were the localisation of an unknown up to now Buddhist Monasterywas possible.

Preliminary report of the first season work of the International Pluridisciplinary Archaeological Expedition to Bactria 2008

This year the development of our project can be divided into two main clearly different parts, on one hand the laboratory work, where the sampled ceramic individuals has been prepared and analyzed and the elaboration of the data obtained during the excavation of 2007 has been finished and, on the other hand, the field work developed at the archaeological site during this specific year (2008).In the mark of the analytical work a significant number of ceramic individuals (144) from the different stratigraphical units from various areas of the excavation of Termez and Tchinguiz Tepe sampled duringthe new and previous field works has been archaeometrically characterized. This specific material includedindividuals dated into the Hellenistic and Sassanian period, which has been confirmed by C14dating upon organic samples.At the same time, in the mark of the field work of 2008 the archaeological record, already started tobe under study during the excavation of 2007, has been completed and two new archaeological recordshave been registered on of which is located in the area of Tchinguiz Tepe. For the archeological studythe information of the previous geophysical prospecting has been indisputably taken into considerationand the same methodology has been applied to crosscheck the latter archaeological results.

CR1Dmod: A Matlab program to model 1D complex resistivity effects in electrical and electromagnetic surveys

We have constructed a forward modelling code in Matlab, capable of handling several commonly used electrical and electromagnetic methods in a 1D environment. We review the implemented electromagnetic field equations for grounded wires, frequency and transient soundings and present new solutions in the case of a non-magnetic first layer. The CR1Dmod code evaluates the Hankel transforms occurring in the field equations using either the Fast Hankel Transform based on digital filter theory, or a numerical integration scheme applied between the zeros of the Bessel function. A graphical user interface allows easy construction of 1D models and control of the parameters. Modelling results are in agreement with other authors, but the time of computation is less efficient than other available codes. Nevertheless, the CR1Dmod routine handles complex resistivities and offers solutions based on the full EM-equations as well as the quasi-static approximation. Thus, modelling of effects based on changes in the magnetic permeability and the permittivity is also possible.

An Integrated Hydro geological Study of the Muvattupuzha River Basin, Kerala, India

The present investigation on the Muvattupuzha river basin is an integrated approach based on hydrogeological, geophysical, hydrogeochemical parameters and the results are interpreted using satellite data. GIS also been used to combine the various spatial and non-spatial data. The salient finding of the present study are accounted below to provide a holistic picture on the groundwaters of the Muvattupuzha river basin. In the Muvattupuzha river basin the groundwaters are drawn from the weathered and fractured zones. The groundwater level fluctuations of the basin from 1992 to 2001 reveal that the water level varies between a minimum of 0.003 m and a maximum of 3.45 m. The groundwater fluctuation is affected by rainfall. Various aquifer parameters like transmissivity, storage coefficient, optimum yield, time for full recovery and specific capacity indices are analyzed. The depth to the bedrock of the basin varies widely from 1.5 to 17 mbgl. A ground water prospective map of phreatic aquifer has been prepared based on thickness of the weathered zone and low resistivity values (<500 ohm-m) and accordingly the basin is classified in three phreatic potential zones as good, moderate and poor. The groundwater of the Muvattupuzha river basin, the pH value ranges from 5.5 to 8.1, in acidic nature. Hydrochemical facies diagram reveals that most of the samples in both the seasons fall in mixing and dissolution facies and a few in static and dynamic natures. Further study is needed on impact of dykes on the occurrence and movement of groundwater, impact of seapages from irrigation canals on the groundwater quality and resources of this basin, and influence of inter-basin transfer of surface water on groundwater.

Investigación geofísica de los paleovalles pleistocenos del NE de Cataluña: el ejemplo de la Riera de Tossa (Costa Brava meridional)

The geometry of Riera de Tossa (Costa Brava, Catalonia) Pleistocene paleovalley has been studied by two geophysical prospecting methods: vertical electric sounding (VES-DC) and seismic refraction. The results show the existence of a sudden slope change of the thalweg that permits to differentiate two sections where erosion processes operated in unequal intensity. Based on geological and geomorphological available data, it can be inferred that this slope rupture is caused by the litological control imposed by the alineation of resistant materials, which limited the headward erosion extend of the creek during the Upper Pleistocene

Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional- and shear-wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse-echo method. The measurements were made both in vacuum-dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin-impregnated polished thin sections, X-ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional- and shear-wave velocities (V-p and V-s, respectively) decrease with increasing porosity and that V-p decreases approximately twice as fast as V-s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore-structure-dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot-Gassmann velocity values are greater than the measured velocity values due to the rock-fluid interaction. This is not accounted for in the Biot-Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time-average relationships overestimated the measured velocities even more than the Biot model.