Fractured reservoirs contain about 85 and 90 percent of oil and gas resources respectively in Iran. A comprehensive study and investigation of fractures as the main factor affecting fluid flow or perhaps barrier seems necessary for reservoir development studies. High degrees of heterogeneity and sparseness of data have incapacitated conventional deterministic methods in fracture network modeling. Recently, simulated annealing (SA) has been applied to generate stochastic realizations of spatially correlated fracture networks by assuming that the elastic energy of fractures follows Boltzmann distribution. Although SA honors local variability, the objective function of geometrical fracture modeling is defined for homogeneous conditions. In this study, after the introduction of SA and the derivation of the energy function, a novel technique is presented to adjust the model with highly heterogeneous data for a fractured field from the southwest of Iran. To this end, the regular object-based model is combined with a grid-based technique to cover the heterogeneity of reservoir properties. The original SA algorithm is also modified by being constrained in different directions and weighting the energy function to make it appropriate for heterogeneous conditions. The simulation results of the presented approach are in good agreement with the observed field data.