Horizontal transverse isotropy anisotropic parameter inversion via azimuthal seismic velocity anisotropy and its application to anisotropic 3D in-situ stress estimation

Anisotropic parameter inversion in horizontal transverse isotropy (HTI) media plays an important role in predicting the fracture density and the anisotropic in-situ stress for unconventional reservoirs. The current industry practice is to use the azimuthal PP-wave reflection coefficient to estimate the HTI anisotropic parameters. Based on the linear slip theory, we develop an innovative approach that uses azimuthal P-wave phase velocity to calculate HTI anisotropic parameters, which presents superiority over the conventional azimuthal PP-wave reflection coefficient inversion. Specifically, we first verify that the azimuthal P-wave phase velocity is for the HTI elliptical fitting than the azimuthal PP-wave reflection coefficient using the analytical formulations. Second, we sort the prestack wide-azimuth data into offset vector tile sectors and perform amplitude-variation-with-offset inversion at each azimuth. Third, we apply an elliptical fitting to the obtained azimuthal P-wave phase velocities to estimate the HTI anisotropic parameters, fracture density, and fracture direction. Fourth, based on the HTI mechanical earth model, we formulate a cost-effective 3D in-situ stress estimation method using the obtained elastic parameters and fracture compliance. Finally, field examples from the Zhaotong area, China, demonstrate that the estimated fracture density and anisotropic in-situ stress are more accurate and have higher resolution compared with conventional methods. The dominant stress regime in the study area is a strike-slip faulting regime with a governing orientation of northeast–southwest and presents good alignments with well logs, which demonstrates the reliability and accuracy of our method for predicting fracture density and anisotropic in-situ stress.