Abstract-In this study we present a step-by-step theoretical modelling approach, using established seismic wave propagation theories in anisotropic media, to generate unique anisotropic reflection patterns observed from three-dimensional pure-mode pressure (3D-PP), full-azimuth and full-offset seismic reflection data acquired over a naturally fractured tight carbonate field, onshore Texas, USA. Our aim is to gain an insight into the internal structures of the carbonate reservoir responsible for the observed anisotropic reflection patterns. From the generated model we were able to establish that the observed field seismic reflection patterns indicate azimuthal anisotropy in the form of crack induced shearwave splitting and variation in P-wave velocity with offset and azimuth. Amplitude variation with azimuth (AVAZ) analysis also confirmed multi-crack sets induced anisotropy which is characteristic of orthorhombic symmetry, evident as multiple bright and dim-amplitude azimuth directions as well as complete reversal of bright-amplitude to dim-amplitude azimuth direction as the angle of incidence increases from near (B15°) to mid (C30°) offsets. Finally, we fitted the generated P-wave velocity into an ellipse to determine the intensity and orientation (N26E) of the open crack set as well as the direction of the minimum in situ stress axis (N116E) within the reservoir. The derived information served as an aid for the design of horizontal well paths that would intercept open fractures and ensure production optimization of the carbonate reservoir, which was on production decline despite reservoir studies that indicate un-depleted reserves.