Estimation of Seismic Wave Attenuation from 3D Seismic Data: A Case Study of OBC Data Acquired in an Offshore Oilfield

Previous studies performed in Abu Dhabi oilfields, United Arab Emirates, revealed the direct link of seismic wave attenuation to petrophysical properties of rocks. However, all those studies were based on zero offset VSP data, which limits the attenuation estimation at one location only. This is due to the difficulty of estimating attenuation from 3D seismic data, especially in carbonate rocks. To overcome this difficulty, we developed a workflow based on the centroid frequency shift method and Gabor transform which is optimized by using VSP data. The workflow was applied on 3D Ocean Bottom Cable seismic data. Distinct attenuation anomalies were observed in highly heterogeneous and saturated zones, such as the reservoirs and aquifers. Scattering shows significant contribution in attenuation anomalies, which is unusual in sandstones. This is due to the complex texture and heterogeneous nature of carbonate rocks. Furthermore, attenuation mechanisms such as frictional relative movement between fluids and solid grains, are most likely other important causes of attenuation anomalies. The slight lateral variation of attenuation reflects the lateral homogeneous stratigraphy of the oilfield. The results demonstrate the potential of seismic wave attenuation for delineating heterogeneous zones with high fluid content, which can substantially help for enhancing oil recovery.

NUMERICAL ALGORITHM OF SEISMIC ATTENUATION ESTIMATION IN ANISOTROPIC FRACTURED POROUS FLUID-SATURATED MEDIA

In our work we investigate the effect of transport and elastic properties anisotropy on seismic attenuation due to fracture-to-fracture wave-induced fluid flow using numerical algorithm of estimation of seismic wave attenuation in anisotropic fractured porous fluid-saturated media. Algorithm is based on numerical solution of anisotropic Biot equations using finite-difference scheme on staggered grid. We perform a set of numerical experiments to model wave propagation in fractured media with anisotropic fractured-filling material providing wave-induced fluid flow within interconnected fractures. Recorded signals are used for numerical estimation of inverse quality factor. Results demonstrate the effect of fracture-filling material anisotropy on seismic wave attenuation.

Seismic wave attenuation due to fluid pressure diffusion at the mesoscopic scale: an experimental and numerical study

<p>Monitoring of the subsurface with seismic methods can be improved by better understanding the attenuation of seismic waves due to fluid pressure diffusion (FPD). In porous rocks saturated with multiple fluid phases the attenuation of seismic waves by FPD is sensitive to the mesoscopic scale distribution of the respective fluids. The relationship between fluid distribution and seismic wave attenuation could be used, for example, to assess the effectiveness of residual trapping of carbon dioxide (CO2) in the subsurface. Determining such relationships requires validating models of FPD with accurate laboratory measurements of seismic wave attenuation and modulus dispersion over a broad frequency range, and, in addition, characterising the fluid distribution during experiments. To address this challenge, experiments were performed on a Berea sandstone sample in which the exsolution of CO2 from water in the pore space of the sample was induced by a reduction in pore pressure. The fluid distribution was determined with X-ray computed tomography (CT) in a first set of experiments. The CO2 exosolved predominantly near the outlet, resulting in a heterogeneous fluid distribution along the sample length. In a second set of experiments, at similar pressure and temperature conditions, the forced oscillation method was used to measure the attenuation and modulus dispersion in the partially saturated sample over a broad frequency range (0.1 - 1000 Hz). Significant P-wave attenuation and dispersion was observed, while S-wave attenuation and dispersion were negligible. These observations suggest that the dominant mechanism of attenuation and dispersion was FPD. The attenuation and dispersion by FPD was subsequently modelled by solving Biot’s quasi-static equations of poroelasticity with the finite element method. The fluid saturation distribution determined from the X-ray CT was used in combination with a Reuss average to define a single phase effective fluid bulk modulus. The numerical solutions agree well with the attenuation and modulus dispersion measured in the laboratory, supporting the interpretation that attenuation and dispersion was due to FPD occurring in the heterogenous distribution of the coexisting fluids. The numerical simulations have the advantage that the models can easily be improved by including sub-core scale porosity and permeability distributions, which can also be determined using X-ray CT. In the future this could allow for conducting experiments on heterogenous samples.</p>

Modelling viscoelastic wave phenomenon by homogenisation of the poroelasticity equations

Poroelastic effects have been of great interest in the seismic literature as they have been identified as a major cause of wave attenuation in heterogeneous porous media. The observed attenuation in the seismic wave can be explained in part by energy loss to fluid motion in the pores. On the other hand, it is known that the attenuation is particularly pronounced in stratified structures where the scale of spatial heterogeneity is much smaller than the seismic wavelength. Understanding of poroelastic effects on seismic wave attenuation in heterogeneous porous media has largely relied on numerical experiments. In this work, we present a homogenisation technique to obtain an upscaled viscoelastic model that captures seismic wave attenuation when the sub-seismic scale heterogeneity is periodic. The upscaled viscoelastic model directly relates seismic wave attenuation to the material properties of the heterogeneous structure. We verify our upscaled viscoelastic model against a full poroelastic model in numerical experiments. Our homogenisation technique suggests a new approach for solving coupled equations of motion.