A rock physics and seismic tomography study to characterize the structure of the Campi Flegrei caldera

The Campi Flegrei (CF) caldera experiences dramatic ground deformations unsurpassed anywhere in the world. The source responsible for this phenomenon is still debated. With the aim of exploring the structure of the caldera as well as the role of hydrothermal fluids on velocity changes, a multidisciplinary approach dealing with 3-D delay-time tomography and rock physics characterization has been followed. Selected seismic data were modeled by using a tomographic method based on an accurate finite-difference travel-time computation which simultaneously inverts P-wave and S-wave first-arrival times for both velocity model parameters and hypocenter locations. The retrieved P-wave and S-wave velocity images as well as the deduced Vp/Vs images were interpreted by using experimental measurements of rock physical properties on CF samples, to take into account steam/water phase transition mechanisms affecting P-wave and S-wave velocities. Also, modelling of petrophysical properties for site-relevant rocks constrains the role of overpressured fluids on velocity. A flat and low Vp/Vs anomaly lies at 4 km depth under the city of Pozzuoli. Earthquakes are located at the top of this anomaly. This anomaly implies the presence of fractured over-pressured gas-bearing formations and excludes the presence of melted rocks. At shallow depth, a high Vp/Vs anomaly located at 1 km suggests the presence of rocks containing fluids in the liquid phase. Finally, maps of the Vp*Vs product show a high Vp*Vs horse-shoe shaped anomaly located at 2 km depth. It is consistent with gravity data and well data and might constitute the on-land remainder of the caldera rim, detected below sea level by tomography using active source seismic data. For a more exhaustive description of the utilized methodologies, of synthetic tests for spatial resolution and uncertainty assessment and, the interpretation of results, the reader may refer to the paper Vanorio et al. (2005)

Three-dimensional tomography and rock properties of the Larderello-Travale

In a geothermal area, a detailed knowledge of the three-dimensional velocity structures aids the managementof the field and the further development of the geothermal source. Here,we present a high-resolution study of the three-dimensional S-wave velocity structures from microearthquake travel times for the Larderello-Travale geothermal field, Italy.We have also deduced the Vp/Vs and Vp ×Vs parameters for this area toemphasize the deep variations in the physical rock properties due to fluid content and porosity. Furthermore, effective porousmedium modelling has been performed for site-relevant lithologies, to improve our interpretation of the results in terms of rock physics signatures. This has allowed us to estimate the variation range of the seismological parameters investigated, as well as their sensitivity for suitable rock under specific physical conditions. LowVp/Vs anomalies, arising froma lower Vp compared to Vs, dominate the geothermal field of Larderello-Travale. These have been interpreted as due to steam-bearing formations. On the contrary, analysis of Vp ×Vs images provides information on the relative changes in rock porosity at depth. Comparison of tomographic section images with previously interpreted seismic lines suggests that the reflective ‘K-horizon’ delineates a transition between zones that have different porosities or crack gatherings. The ‘K-horizon’ also lies on low Vp/Vs anomalies, which suggests a steam saturation zone, despite the reduced porosity at this depth

A 3D velocity model for earthquake location in Campi Flegrei area: application to the 1982-84 uplift event

The uplift crisis of the 1982-1984 in the Campi Flegrei area underlined the importance of seismic surveillance for this volcanic caldera. One of the key elements for an effective seismic network is to make use of a reliable velocity model for earthquake location. In the present work we will discuss criteria for the construction and validation of a new 3D P-wave velocity model for earthquake location in the Campi Flegrei area built from the integration of two high-resolution 3D tomographic images of the region. The model is used for locating a group of earthquakes from the uplift event of the 1982-1984

Time-lapse sonic logs reveal patchy CO2 saturation in-situ

Based on time-lapse sonic and neutron porosity logs from the Nagaoka CO2 sequestration experiment, a P-wave velocity-saturation relation at reservoir depth is retrieved. It does not coincide with either of the end-member models of uniform and patchy saturation but falls in between even if realistic error estimates for the host rock properties are considered. Assuming a random distribution of CO2 patches it is shown that the mechanism of wave-induced flow can be evoked to explain this velocity-saturation relation. Characteristic CO2 patch size estimates range from 1 to 5 mm. Such mesoscopic heterogeneity can be responsible for attenuation and dispersion in the well logging frequency band

Porosities of building limestones: using the solid density to assess data quality

A good knowledge of the volume-fraction porosity is essential in any technical work on porous materials. In construction materials the porosity is commonly measured by the Archimedes buoyancy method, from which the bulk density of the test specimen is also obtained. The porosity and the bulk density together fix the solid density of the specimen, as only two of the three quantities are independent. The solid density, although rarely discussed, is determined by the mineralogy of the specimen, and therefore can provide a valuable check on the accuracy of porosity and bulk density measurements. Our analysis of published data on calcitic limestones shows that the solid density is generally close to the ideal crystallographic density of calcite. Small deviations can often be traced to variations in mineral composition. However some published porosity–density data are inconsistent with the known mineralogy. Deviations which cannot be ascribed to composition may be assumed to arise from measurement errors. We show the value of using the solid density as a quality check on the measured porosity. We recommend that the solid density should always be calculated for this purpose when the Archimedes method is used. This check can be useful also when porosities are measured by helium pycnometry or by mercury intrusion porosimetry

Time-lapse characterization of hydrothermal seawater and microbial interactions with basaltic tephra at Surtsey Volcano

A new International Continental Drilling Program (ICDP) project will drill through the 50-yearoldedifice of Surtsey Volcano, the youngest of the Vestmannaeyjar Islands along the south coast of Iceland, to perform interdisciplinary time-lapse investigations of hydrothermal and microbial interactions with basaltic tephra. The volcano, created in 1963–1967 by submarine and subaerial basaltic eruptions, was first drilled in 1979. In October 2014, a workshop funded by the ICDP convened 24 scientists from 10 countries for 3 and a half days on Heimaey Island to develop scientific objectives, site the drill holes, and organize logistical support. Representatives of the Surtsey Research Society and Environment Agency of Iceland also participated. Scientific themes focus on further determinations of the structure and eruptive processes of the type locality of Surtseyan volcanism, descriptions of changes in fluid geochemistry and microbial colonization of the subterrestrial deposits since drilling 35 years ago, and monitoring the evolution of hydrothermal and biological processes within the tephra deposits far into the future through the installation of a Surtsey subsurface observatory. The tephra deposits provide a geologic analog for developing specialty concretes with pyroclastic rock and evaluating their long-term performance under diverse hydrothermal conditions

KG2B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 1: Measurements, pressure dependence and pore-fluid effects

Measuring the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore-fluids. For that purpose a collaborative benchmarking exercise involving 24 laboratories was organized for measuring the permeability of a single low permeability material, the Grimsel granodiorite, at a common effective confining pressure (5 MPa). The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyze the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low permeability measurements. In total 39 measurements were collected that allowed us to discuss the influence of (i) pore-fluid, (ii) measurement method, (iii) sample size and (iv) pressure sensitivity. Discarding some outliers from the bulk data set (4 out of 39) an average permeability of 1.11 × 10−18 m² with a standard deviation of 0.57 × 10−18 m² was obtained. The most striking result was the large difference in permeability for gas measurements compared to liquid measurements. Regardless of the method used, gas permeability was higher than liquid permeability by a factor approximately 2 (kgas = 1.28 × 10−18 m² compared to kliquid = 0.65 × 10−18 m²). Possible explanations are that (i) liquid permeability was underestimated due to fluid-rock interactions (ii) gas permeability was overestimated due to insufficient correction for gas slippage and/or (iii) gases and liquids do not probe exactly the same porous networks. The analysis of Knudsen numbers shows that the gas permeability measurements were performed in conditions for which the Klinkenberg correction is sufficient. Smaller samples had a larger scatter of permeability values, suggesting that their volume were below the Representative Elementary Volume. The pressure dependence of permeability was studied by some of the participating teams in the range 1–30 MPa and could be fitted to an exponential law k = ko.exp(–γPeff) with γ = 0.093 MPa−1. Good practice rules for measuring permeability in tight materials are also provided

Elastic Properties of Clay Dry Mineral Aggregates, Suspension and Sandstones

Geophysical Journal International, v. 155, n. 1, p. 319-326, 2003. http://dx.doi.org/10.1046/j.1365-246X.2003.02046.xInternational audienc