Joint PP, PS1, and PS2 AVA Inversion of HTI Media

Abstract

Sedimentary strata formed by orientational vertical fractures are approximately equivalent to transversely isotropic media with a horizontal symmetry axis (HTI). Parameter inversion for the isotropic background and fracture parameters can provide valuable information for predicting fractured reservoirs. Following Bakulin’s modeling theory, we parameterized the HTI medium with the first- and second-order Lamé parameters of the isotropic background, host rock density, fracture density, and fracture azimuth. We calculated the fracture azimuth and separated split shear waves using the Alford rotation. Based on the exact reflection coefficients, we develop a joint P-to-P reflected wave (PP), PS1, and PS2 amplitude variation with angle (AVA) inversion method in the depth domain by using the second-order Levenberg–Marquardt (LM) algorithm to invert the remaining four parameters. Numerical model tests revealed that our method was more accurate compared to the traditional LM method and works for both weak/strong contrasts and anisotropy magnitudes. We used real data from the walkaway vertical seismic profiling (VSP) as an example to propose a set of strategies, thereby improving the practicality of our inversion method. To protect the fracture information contained in the wavefields, we adopted vector-feature-preserving data processing procedures to preprocess the field VSP data and separate PP, PS1, and PS2 reflections. We utilized P-wave first breaks to successfully discriminate and match these three wave modes from the same fractured formation in the depth domain; these were subsequently transformed to AVA gathers with ray tracing. We synthetically incorporated the time delays and layer velocities to build the initial models and to constrain the inversion process.