Modeling method for rock heterogeneities and multiple hydraulic fractures propagation based on homogenization approach and PHF-LSM