A simulated reservoir model, based on the permeability fractal model and three-dimensional (3D) Gaussian filter, was established to account for in-layer interlayer heterogeneity so that result conforms law of geological statistics. Combined with an embedded discrete fracture method (EDFM), a multiscale system established, forming numerical simulation multiphase flow in horizontal wells heteroge...