Thomsen's parameters from p-wave measurements in a spherical sample

The aim of this paper is to understand the seismic anisotropy of the overburden shale in an oilfield in the North West Shelf of Western Australia. To this end, we first find the orientation of the symmetry axis of a spherical shale sample from measurements of ultrasonic P-wave velocities in 132 directions at the reservoir pressure. After transforming the data to the symmetry axis coordinates, we find Thomsen's anisotropy parameters d and ? using these measurements and measurements of the shear-wave velocity along the symmetry axis from a well log. To find these anisotropy parameters, we use a very fast simulated re-annealing algorithm with an objective function that contains only the measured ray velocities, their numerical derivatives and the unknown elasticity parameters. The results show strong elliptical anisotropy in the overburden shale. This approach produces smaller uncertainty of Thomsen parameter d than more direct approaches.