Horizontal wells have been used to increase reservoir recovery, especially in unconventional reservoirs, and hydraulic fracturing has been applied to further extend the contact with the reservoir to increase the efficiency of development. In the past, many models, analytical or numerical, were developed to describe the flow behavior in horizontal wells with fractures. Source solution is one of the analytical/semi-analytical approaches. To solve fractured well problems, source methods were advanced from point sources to volumetric source, and pressure change inside fractures was considered in the volumetric source method. This study aims at developing a method that can predict horizontal well performance and the model can also be applied to horizontal wells with multiple fractures in complex natural fracture networks. The method solves the problem by superposing a series of slab sources under transient or pseudosteady-state flow conditions. The principle of the method comprises the calculation of semi-analytical response of a rectilinear reservoir with closed outer boundaries. A statistically assigned fracture network is used in the study to represent natural fractures based on the spacing between fractures and fracture geometry. The multiple dominating hydraulic fractures are then added to the natural fracture system to build the physical model of the problem. Each of the hydraulic fractures is connected to the horizontal wellbore, and the natural fractures are connected to the hydraulic fractures through the network description. Each fracture, natural or hydraulically induced, is treated as a series of slab sources. The analytical solution of superposed slab sources provides the base of the approach, and the overall flow from each fracture and the effect between the fractures are modeled by applying superposition principle to all of the fractures. It is assumed that hydraulic fractures are the main fractures that connect with the wellbore and the natural fractures are branching fractures which only connect with the main fractures. The fluid inside of the branch fractures flows into the main fractures, and the fluid of the main fracture from both the reservoir and the branch fractures flows to the wellbore. Predicting well performance in a complex fracture network system is extremely challenged. The statistical nature of natural fracture networks changes the flow characteristic from that of a single linear fracture. Simply using the single fracture model for individual fracture, and then adding the flow from each fracture for the network could introduce significant error. This study provides a semi-analytical approach to estimate well performance in a complex fracture network system.