Numerical simulation of hydraulic fracture propagation in heterogeneous unconventional reservoir