Numerical modeling of hydraulic fracturing process in a heterogeneous medium through a coupling scheme between the cohesive fracture and the lattice-Boltzmann models