This paper presents a new three-dimensional fully coupled poroelastic numerical model to simulate pressure transient response of naturally fractured reservoirs. One of the main applications of the new approach is to improve the reservoir characterization by decreasing the uncertainties associated with subsurface fracture map and to understand the interaction between fracture and matrix. The model is based on a hybrid methodology to simulate fluid flow by combining a single continuum and discrete fracture network approaches. To decrease the uncertainty associated with subsurface fracture map, a different realizations of the discrete fracture systems are generated based on image logs, seismic, geological features and the tectonic history of the studied reservoir. An iterative loops have been used between the reservoir static model and dynamic well-test data to optimize the generation of subsurface fracture maps. At each loop, a new discrete fracture system is used and the simulated pressure transient data are compared with the available field data. The process is repeated till the matching between the simulated and the measured pressure transient data is achieved using the most appropriate fracture realization. The numerical model is validated against an analytical pressure transient solution for dual-porosity system, and then applied to a real case taken from fractured basement reservoir in offshore Southern Vietnam. The developed hybrid methodology used to simulate fluid flow and rock deformation to improve the characterization of the fractured basement by simulating the pressure transient testing. The results predicted by the presented model are in a good agreement with field data; also the model was able to predict the fractures network distribution and fractures inherent properties around the tested well.