A methodology is provided to optimize the dynamic connectivity of a discrete fracture network (DFN) model of a subsurface reservoir against observed reservoir production measures using streamlines and machine learning. Adjustment of discrete fracture network properties of the reservoir is made locally and minimizes computer processing time spent in history matching. An iterative workflow identifies history match issues between measured and predicted or simulated water cut of reservoir produced fluids. Streamline analysis quantifies injector-producer communication and identifies reservoir grid block bundles that dominate dynamic response. A genetic algorithm updates discrete fracture network properties of the reservoir model to improve dynamic history match response.