Numerical Modeling for Q-anisotropy in Fractured Media

A method for numerical modelling of three-component (3C) three-dimensional (3D) wavefields in horizontally layered attenuative medium is derived. The method is based on delta-operators, six-dimensional Stroh formalism and Haskell-Thomson matrix method. The earth model consists of several multi-fractured layers, each of which represented by anisotropic medium model with attenuation. The attenuative fractured medium model is described by Schoenberg's linear slip formalism developed further more by introduction of the complex-valued weaknesses in the stiffness matrix. The numerical modelling method is validated for simplified earth model, which is attenuative HTI medium due to single fracture set. The Q-anisotropy and velocity anisotropy revealed from simulated wavefiedls fairly match theory-predicted quantities.