Recent research has established that the dominant P-wave attenuation mechanism in reservoir rocks at seismic frequencies is because of wave-induced fluid flow ͑mesoscopic loss͒. The P-wave induces a fluid-pressure difference at mesoscopic-scale inhomogeneities ͑larger than the pore size but smaller than the wavelength, typically tens of cen-timeters͒ and generates fluid flow and slow ͑diffusion͒ Biot waves ͑continuity of pore pressure is achieved by energy conversion to slow P-waves, which diffuse away from the interfaces͒. In this context, we consider a periodically stratified medium and investigate the amount of attenuation ͑and velocity dispersion͒ caused by different types of heterogeneities in the rock properties, namely, porosity, grain and frame moduli, permeability, and fluid properties. The most effective loss mechanisms result from porosity variations and partial saturation, where one of the fluids is very stiff and the other is very compliant, such as, a highly permeable sandstone at shallow depths, saturated with small amounts of gas ͑around 10% saturation͒ and water. Grain-and frame-moduli variations are the next cause of attenuation. The relaxation peak moves towards low frequencies as the ͑background͒ permeability decreases and the viscosity and thickness of the layers increase. The analysis indicates in which cases the seismic band is in the relaxed regime, and therefore, when the Gassmann equation can yield a good approximation to the wave velocity.

We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO 2 sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO 2 , methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.

No abstract

Crystal orientation fabric (COF) describes the intrinsic anisotropic nature of ice and is an important parameter for modeling glacier flow. We present the results of three‐component active‐source seismic observations from the Whillans Ice Stream (WIS), a fast‐flowing ice stream in West Antarctica. Surface‐wave dispersion analysis, ray tracing, and traveltime inversion of compressional (P) and shear (S) waves reveal the presence of transversely isotropic ice with a vertical axis of symmetry (VTI) beneath approximately 65 m of isotropic firn. The ice stream is characterized by weak anisotropy, involving an average ice thickness of approximately 780 m. The analysis indicates that about 95% of the ice mass is anisotropic, and the crystalline c axes span within an average broad cone angle of 73 ± 10° with respect to the vertical axis. Moreover, the mean temperature T (below the firn) estimated from seismic data is −15 ± 5°C. These data do not show evidence of englacial seismic reflectivity, which indicates the lack of abrupt changes in the COF. The presence of azimuthal anisotropy due to transversely compressive flow or fractures aligned along a preferential direction is also excluded. We suggest that the observed VTI ice structure is typical of large ice streams in regions where basal sliding and bed deformation dominate over internal glacial deformation.

ABSTRACT. Microtremor measurements and the horizontal-to-vertical spectral ratio (HVSR) technique, generally used for site effect studies as well as to determine the thickness of soft sedimentary layers, can effectively be applied to map the thickness of glaciers. In this work the radio-echo sounding, geoelectric and active seismic methods, widely employed to image the earth interior, are applied to verify the reliability of the HVSR technique in Alpine and Antarctic glacial environments. The technique has been used to analyze passive seismic data from glaciers of the Adamello and Ortles-Cevedale massifs (Italy), the Bernese Oberland Alps (Switzerland) and from the Whillans Ice Stream (West Antarctica). Comparing with the results obtained from the different geophysical imaging methods, we show that the resonance frequency in the HVSR spectra correlates well with the ice thickness at the site, in a wide range from a few tens of meters to more than 800 m. The reliability of the method mainly depends on the coupling of sensors at the glacier surface and on the basal impedance contrast. This passive seismic technique offers a logistically efficient and cost effective method to map glacier and ice-sheet thicknesses. Moreover, under certain conditions, it allows reliable estimations of the basal seismic properties.