Naturally Fractured Carbonate Reservoir Characterization: A Case Study of a Mature High-Pour Point Oil Field in Hungary

An integrated technical study was conducted for a field development project in West-Hungary. This study offers a better solution for estimating petrophysical properties and fracture facies vertically along the well and laterally for 3D static and dynamic models of naturally fractured reservoirs in carbonate rocks. More than 30 wells with 40 years of production history were used in order to build reliable static and dynamic models. The fracture class/facies plays essential role in spatial distribution of petrophysical properties during 3D reservoir modeling. It was defined by integrating the conventional logs, image logs, drilling parameters, and production or well test data. Three fracture facies are defined as macro-fracture (including permeable sub-seismic fault), micro-fracture, and hostrock. Subsequently the fracture-class’s spatial distribution is guided by seismic attributes of faultlikelihood combined with geological concept of fault and damage zone. As a result, the established fracture classes along the wells are validated by static and dynamic subsurface data. Spherical self-organizing map (SOM) was also implemented for predicting the fracture location in wells having limited subsurface data. Moreover, fracture lateral distribution follows the distribution of the fault zone of fault core, high-damage zone, low-damage zone, and host-rock. The higher the fault displacement the wider the damage zone and fault core formed. Macrofractures and micro-fractures frequently appear around fault core and high damage zone. While only microfractures are dominantly present in the low damaged zones. In contrast, the unfractured class is dominantly distributed in host rock area. Also, the lithologis considered in distributing the fracture class because the rock mechanic properties and number of fractures are strongly controlled by rock compositions. Once the fracture class is distributed, porosity, permeability, and water saturation are modelled in the 3D geocellular model. Finally, this fracture class also plays a role as a rock typing for reservoir simulation. The saturation height model is built using the fracture class distribution resulting the initialization, history matching process, and production forecast from 20 wells are showing excellent quality. As a novelty, this study offers a better understanding of fracture distribution and accelerates the history matching process with a more confident result of production forecast. In the absence of advanced technologies like image logs and production logging (PLT) measurements, this study still effectively assists us to recognize the fracture presence and its quality in both well-depth interval and 3D spatial space, and successfully guided us in proposing a new infill drilling with strong confidence and delivering on the high-end of expected results.