Feasibility of Seismic Characterization of Multiple Fracture Sets

Estimation of parameters of multiple fracture sets is often required for successful exploration and development of naturally fractured reservoirs. The goal of this paper is to determine the maximum number of fracture sets of a certain rheological type which, in principle, can be resolved from seismic data. The main underlying assumption is that an estimate of the complete effective stiffness tensor has been obtained, for example, from multiazimuth, multicomponent surface seismic and vertical seismic profiling (VSP) data. Although typically only a subset of the stiffness elements (or some of their combinations) may be available, this study helps to establish the limits of seismic fracture-detection algorithms. The number of uniquely resolvable fracture systems depends on the anisotropy of the host rock and the rheology and orientation of the fractures. Somewhat surprisingly, it is possible to characterize fewer vertical fracture sets than dipping ones, even though in the latter case the fracture dip has to be found from the data. For the simplest, rotationally invariant fractures embedded in either isotropic or transversely isotropic with a vertical symmetry axis (VTI) host rock, the stiffness tensor can be inverted for up to two vertical or four dipping fracture sets. In contrast, only one fracture set of the most general (microcorrugated) type, regardless of its orientation, is constrained by the effective stiffnesses. These results can be used to guide the development of seismic fracture-characterization algorithms that should address important practical issues of data acquisition, processing, and inversion for particular fracture models.

[1]  Ilya Tsvankin,et al.  Anisotropic parameters and P-wave velocity for orthorhombic media , 1997 .

[2]  V. Grechka,et al.  Multicomponent stacking-velocity tomography for transversely isotropic media , 2002 .

[3]  Michael Schoenberg,et al.  Elastic Properties of Rocks With Multiple Sets of Fractures , 1989 .

[4]  Michael E. Jones,et al.  P-wave and S-wave azimuthal anisotropy at a naturally fractured gas reservoir, Bluebell‐Altamont Field, Utah , 1999 .

[5]  A. Bakulin,et al.  Estimation of fracture parameters from reflection seismic data—Part II: Fractured models with orthorhombic symmetry , 2000 .

[6]  A. Bakulin,et al.  Estimation of Azimuthal Anisotropy And Fracture Parameters From Multiazimuthal Walkawa YVSP In the Presence of Lateral Heterogeneity , 2000 .

[7]  V. Grechka,et al.  Inversion of multicomponent, multiazimuth, walkaway VSP data for the stiffness tensor , 2003 .

[8]  G. Fryer,et al.  Reduction of the general fracture compliance matrix Z to only five independent elements , 1991 .

[9]  V. Grechka,et al.  Joint Inversion of P- And PS-waves In Orthorhombic Media: Theory And a Physical-modeling Study , 1997 .

[10]  Francis Muir,et al.  A calculus for finely layered anisotropic media , 1989 .

[11]  A. Bakulin,et al.  Estimation of fracture parameters from reflection seismic data—Part I: HTI model due to a single fracture set , 2000 .

[12]  S. Leaney,et al.  Short note: Polarization and slowness component inversion for TI anisotropy , 2000 .

[13]  R. Michelena,et al.  Fracture detection in a carbonate reservoir using a variety of seismic methods , 1999 .

[14]  M. Schoenberg.  Elastic wave behavior across linear slip interfaces , 1980 .

[15]  A. Bakulin,et al.  Estimation of fracture parameters from reflection seismic data—Part III: Fractured models with monoclinic symmetry , 2000 .

[16]  M. Schoenberg,et al.  Elastic wave propagation in media with parallel fractures and aligned cracks , 1988 .

[17]  M. Schoenberg,et al.  REFLECTION OF ELASTIC WAVES FROM PERIODICALLY STRATIFIED MEDIA WITH INTERFACIAL SLIP , 1983 .

[18]  A. Bakulin,et al.  Seismic inversion for the parameters of two orthogonal fracture sets in a VTI background medium , 2002 .

[19]  M. Schoenberg,et al.  Azimuth‐dependent tuning of seismic waves reflected from fractured reservoirs , 1999 .