Reliable characterization of subsurface fractures information within the carbonate's resources plays a key role in reservoir description in terms of the storage and movement of fluid. Commonly, anisotropy information classifies heterogeneity in the rock fabric due to fractures, pores, and local stress variations. Amplitude Variation with Azimuth investigation examines azimuthal variations in the seismic amplitude information to provide useful information about fractures and anisotropy of a horizontally transverse isotropic (HTI).This paper presents modeling, analysis, and interpretation of Azimuthal amplitude variation with angles (AVAz) from 3D seismic and calibrate results with the wells data to extract fractures information in terms of fractures magnitude and orientations. We employed two most used data-driven alternate approaches 1) near-offset Rüger's equation and 2) azimuthal Fourier Coefficient (FC) for analyzing the AVAz responses in seismic data to estimate the fractures density and their orientations, and results are calibrated with the bore hole image data (FMI). These approaches proved practically convenient to quantify the azimuthal amplitude variations with two sets of AVAZ attributes. The results of AVAz modeling and analysis required good quality of input data. Before AVAZ modeling and attributes analysis, seismic data is pre-conditioned, and the data quality of the 3D seismic dataset is improved. Several qualitative and quantitative QC measures (e.g., comparison of angle stacks, AVAz analysis, cross-plots of correlation coefficients, and time shifts) are employed to ascertain the improvement of the health of the seismic data. Azimuthal AVA modeling helped to understand the azimuthal amplitude response due to anisotropy caused by fractures/stress. The predicted seismic response from the models is used to assess the sensitivity of AVAZ with the parameters of the anisotropic model. After the AVAZ modeling, the outcomes of the AVAZ attributes are generated using to 3D seismic and calibrated with the well data set to identify the anomaly of anisotropic gradient (Bani) and second Fourier coefficients (R2) by relating to fractures density or stress magnitude. These attributes helped to identify the amplitude variations with offset and azimuth in the deep inhomogeneous tight carbonate reservoir interval (Patala) of DS hydrocarbon field discovered in the north Potwar Basin, onshore Pakistan. High score of Bani or R2 observed at well location at the tight Paleocene carbonate (Patala formation) reservoirs exhibited the presence of fractures which is confirmed from the image logs (e.g., FMI).The results demonstrated that the studied formations exhibit significant heterogeneity due to the faults/fractures and stresses associated with the complex geological structure of the basin as well as variations in the rock fabric. Hence, accurate analysis of the fractures information from seismic data is quite useful not only in understanding the reservoir heterogeneity to optimize the production, but also in predicting the completion and development strategies that may be most effective.