Polynomial amplitude versus azimuth inversion in horizontally transverse isotropic media, as tested on fractured coal seams in the Surat Basin

Abstract

A new technique of amplitude versus azimuth (AVAZ) seismic inversion in horizontally transverse isotropic (HTI) media is presented. AVAZ is an effective method of characterising anisotropic variation within individual reflectors as well as characterising fractures. The compressional wave reflectivity equation in HTI media has been reformulated into a parabolic form that allows for fast and efficient inversion. The isotropic component of the azimuthal reflectivity has been separated precisely from the anisotropic component and the anisotropic component has been decoupled exactly into its constituent elliptic and anelliptic components. The exact isotropic, elliptic and anelliptic amplitude versus offset (AVO) gradient equations in HTI media are presented herein and the amount of error associated with previous approximations is also defined under the assumption of weak anisotropy. A method of calculating Thomsen’s weak anisotropy parameters using these AVO gradient terms is then outlined. Compared with the elliptic method, there is reduced error incorporated in the new AVAZ method and the error relationships of this method are compared with the Fourier method. Data from an open file 3D wide azimuth seismic survey in the Surat Basin were inverted to demonstrate the effectiveness of the techniques which are presented herein. Seismic amplitudes from six azimuthal stacks were extracted over two horizons and inverted around a well where full-wave sonic and density logs were acquired. For both horizons, the error between the inverted anisotropy parameters from seismic and the inverted anisotropy parameters from wire line logs were found to be less than 5% for both horizons.