The principle focus of this work is on proposing a new model to evaluate the production performance of complexly-fractured well in shale reservoirs for conducting rate transient analysis (RTA) and pressure transient analysis (PTA). The model is established on the trilinear-flow idealization presented by Brown et al. (2009) for fractured horizontal well and associated with continuum geometry, fractal characterization and anomalous diffusion. The coupled model could take into account the non-uniform distribution of fracture network within stimulated reservoir volume (SRV) and property heterogeneity among hydraulic fractures in the configuration of multiple fractured horizontal well (MFHW).The second focus is put on developing a novel semi-analytical solution to fill the gap between trilinear-flow model for single-fracture hypothesis and actual MFHW configuration. The new solution could account for the effect of fracture-production interaction and property heterogeneity among fractures. Calculative results show that the type, sequence and duration of flow regimes are determined by fractal characteristics, associated anomalous diffusion, fracture number/spacing, fracture conductivity, etc. Hence, it is possible to observe the feature of long-period linear flow trends in the production performance of fractured well and explain the presentation of new flow regimes which are generally not identified by conventional trilinear-flow model. In addition, approximate solutions for corresponding flow regimes have also been proposed to interpret the production performance using straight forward method. Furthermore, a field example from Sichuan Shale in China is presented to demonstrate the practical application of the new model for interpreting production data analysis of MFHW. The model provides us with new knowledge and insight into understanding production characteristics and allows us to correctly predict well performance in fractured shale reservoirs