The presence of sedimentary layers in the earth�s subsurface results in seismic anisotropy, which makes wave velocity dependent on the propagation angle. This phenomenon gives rise to errors in seismic imaging. Among these errors are the mispositioning of migrated events and failure to retain energy during dip-moveout. Most of hydrocarbon reservoirs are defined as anisotropic media. Anisotropy is necessary not only to avoid distortions in imaging but also provides valuable information about lithology and fracture networks. To consider the influence of seismic anisotropy, an anisotropic wave equation needs to be employed. The main topics, which are discussed in this chapter, are to incorporate anisotropic effects in seismic modeling and Kirchhoff depth imaging for imaging complex structures. Two algorithms are developed, namely, (1) vertical transverse isotropy (VTI) and tilted transverse isotropy (TTI) wave modeling and (2) VTI Kirchhoff depth imaging. In the first part, a new TTI pseudo-acoustic wave equation is suggested for anisotropic forward modeling. In the second part, a VTI fast-marching eikonal solver is constructed for calculating travel times. An anelliptic VTI wave equation, which uses a nonlinear approximation, is utilized to provide the P-wave velocity information. In this study, synthetic data and a real dataset are applied to test the effectiveness of the algorithm. The spectrum comparison confirmed that the VTI algorithm produces images with higher amplitude around 30 more than isotropic condition and thus better resolution. © 2022 Elsevier Inc. All rights reserved.