Abstract-Generation of naturally fractured reservoir subsurface fracture maps and prediction its production potential are considered complex process due to insufficient data available such as bore hole images, core data and proper reservoir simulator model.To overcome such shortcomings, the industry has relied on geo-statistical analyses of hard and soft data, which are often referred to as static data. This paper presents an integrated workflow that models and predicting fractured reservoirs performance, through the use of gradient-based inversion techniques and discrete fracture network modelling (DFN), which-through the inversion of well test data (i.e., dynamic data)-aims to optimise fracture properties and then predicting of the reservoir production potential. The first step in the workflow is to identify flow contributing fracture sets by analysing available core descriptions, borehole images, conventional log data and production data. Once the fracture sets are identified, the fracture intensity is statistically populated in the inter-well space. In the second step, 3D block-based permeability tensors are calculated based on flow through discrete fractures, and the fracture intensity is then propagated away from the wellbore, i.e., by relating to permeability tensors with fracture intensity. In the final step (fracture optimisation), the fracture properties are computed by DFN modelling, which includes distribution, orientation and geometry in different realisations. Fluid flow is simulated in these discrete fractures to estimate pressure change and pressure derivatives. The production rate associated with drill stem test that performed within this reservoir area has been successfully simulated using the optimised subsurface fracture map that has been generated from the first step.