A practical approach to develop a proper anisotropic rock physics model for media comprising of multiple fracture sets in the absence of sophisticated laboratory/wireline data

For economical production from a fractured reservoir, a characteristic analysis of the fracture parameters like its density and orientation within the reservoir is essential to improve the fluid flow during extraction. This study deals with the development of a proper anisotropic rock physics model for a media with multiple fracture sets to study the spatial distribution of important fracture parameters i.e., fracture density and orientation in the absence of sophisticated laboratory/wireline and pre-stack seismic data. The crest of hydrocarbon producing fault-bounded Balkassar Anticline in Northern Potwar, Upper Indus Basin, Pakistan is selected as a case study representing a potential zone for development of fractures at reservoir level (Sakesar Limestone). The methodology consists of the interpretation of 3D post-stack seismic and conventional wireline log data to demarcate the reservoir containing fractures. The Ant-tracking discrete fracture network (DFN) attribute is applied on 3D post-stack seismic data to obtain an initial estimate about the presence of fracture corridors and their orientations. Based on this initial estimate, a proper rock physics model has been developed utilizing inverse Gassmann relations, T-matrix approximation, and Brown and Korringa relations. The output from the developed rock physics model has been displayed in the form of 13 effective independent elastic stiffness constants (monoclinic symmetry–representing media comprising of multiple fracture sets) as a function of fracture densities and azimuthal fracture orientations. A clear decreasing trend in effective elastic stiffness constants with increasing fracture densities can be observed. Similarly, a periodic trend of effective elastic stiffness constants with fracture orientations can be observed. These trends are more or less expected, but they would have been difficult to quantify without a proper rock physics model. The use of independent effective elastic constants for the generation of synthetic seismic amplitude versus angle and azimuth (AVAZ) data and its correlation with observed seismic AVAZ data in a geostatistical sense has been discussed.