Azimuthal Seismic Amplitude Difference Inversion for Fracture Weakness

Abstract

Fracture weakness prediction is an important task in fractured reservoir analysis. We propose a new method to use seismic amplitude differences between azimuths to estimate the normal and tangential fracture weaknesses under the assumption that the anisotropic perturbation of the reflection coefficient is mainly induced by fractures. We first derive an expression of the reflection coefficient in terms of the normal and tangential fracture weaknesses for the case of an interface separating two fractured media. Then we use the linear fitting method to get the relationship between the two fracture weaknesses, and change the variables to precondition the inversion problem. The Bayesian framework, under the hypothesis of a Cauchy distribution prior information and a Gaussian distribution likelihood function, is employed to construct the objective function, and an initial low-frequency constraint is introduced to the objective function to make the inversion more stable. The conjugate gradient algorithm is adopted to solve the inverse problem. Tests on both synthetic and real data demonstrate that the normal and tangential fracture weaknesses can be estimated reasonably in the case of seismic data containing a moderate noise, and our inversion approach appears to be a stable method for predicting the fracture weaknesses.