Inversion of differences in frequency components of azimuthal seismic data for indicators of oil-bearing fractured reservoirs based on an attenuative cracked model

Abstract

Based on a model of attenuative cracked rock, we have derived a simplified and frequency-dependent stiffness matrix associated with (1) a rock volume containing aligned and partially saturated cracks and (2) a new indicator of oil-bearing fractured reservoirs, which is related to pressure relaxation in cracked rocks and influenced by fluid viscosity and saturation. Starting from the mathematical form of a perturbation in this stiffness matrix across a reflecting interface separating two attenuative cracked media, we set up a linearized P-wave to P-wave reflection coefficient as an azimuthally and frequency-dependent function of dry rock elastic properties, dry fracture weaknesses, and the new indicator. By varying this reflection coefficient with azimuthal angle, we derive a further expression referred to as the quasidifference in elastic impedance, or [Formula: see text], which is primarily affected by the dry fracture weaknesses and the new indicator. An inversion approach is established to use differences in frequency components of seismic amplitudes to estimate these weaknesses and the indicator based on the derived [Formula: see text]. In synthetic inversion tests, we determine that the approach produces interpretable parameter estimates in the presence of data with a moderate signal-to-noise ratio (S/N). Testing on a real data set suggests that reliable fracture weakness and indicator are generated by the approach; fractured and oil-bearing reservoirs are identified through a combination of the dry fracture weakness and the new indicator.