The objective of this project were to: (1) Simulate parts of the reservoir using the conventional dual-porosity, dual-permeability approach, generate progressively complex reservoir models and use the production history match results to quantify formation damage; (2) Examine the effect of the numerical aspects such as grid sizes and fracture properties on the simulations results; (3) Evaluate the effect of fluid thermodynamic properties on production; and (4) Study effect of inclusion of fractures on variability in production from stochastically generated reservoir models.