A shale-gas reservoir with a multiple-fractured horizontal well (MFHW) is divided into two regions: The inner region is defined as stimulated reservoir volume (SRV), which is interconnected by the fracture network after fracturing, while the outer region is called unstimulated reservoir volume (USRV), which has not been stimulated by fracturing. Considering an arbitrary interface boundary between SRV and USRV, a composite model is presented for MFHWs in shale-gas reservoirs, which is based on multiple flow mechanisms, including adsorption/desorption, viscous flow, diffusive flow, and stress sensitivity of natural fractures. The boundary-element method (BEM) is applied to solve the production of MFHWs in shale-gas reservoirs. The accuracy of this model is validated by comparing its production solution with the result derived from an analytical method and the reservoir simulator. Furthermore, the practicability of this model is validated by matching the production history of the MFHW in a shale-gas reservoir. The result shows that the model in this work is reliable and practicable. The effects of relevant parameters on production curves are analyzed, including Langmuir volume, Langmuir pressure, hydraulic-fracture width, hydraulic-fracture permeability, natural-fracture permeability, matrix permeability, diffusion coefficient, stress-sensitivity coefficient, and the shape of the SRV. The model presented here can be used for production analysis for shale-gas-reservoir development. Recently, three kinds of numerical methods including the finite-element method (FEM) (Fan et al. 2015), the finite-difference method (FDM) (Moinfar et al. 2013), and the BEM (Kikani and Horne 1992; Zhao et al. 2016; Idorenyin and Shirif 2017) have been used to solve the production/pressure solution of the MFHW in composite reservoirs. The BEM, which is based on the fundamental solution that derives from Green's function and satisfies the governing equation, is effective in solving the production/pressure of the well in complex-boundary reservoirs (Cheng 2003). In comparison with the FEM and the FDM, the BEM can solve more efficiently the production/pressure performance of the MFHW in composite reservoirs with irregular boundaries. The reason is that the BEM needs only to discretize the boundaries into a small number of boundary elements, and solve the production/pressure solution of the boundary elements, but FEM and FDM need to divide the reservoir into numerous grids and perform an extensive calculation to solve the pressure of all grids at each timestep. The BEM has been successfully used to analyze the production/pressure performance of vertical wells, horizontal wells, verticalfractured wells, and MFHWs in reservoirs with arbitrary boundaries (Kikani and Horne 1992; Sato and Horne 1993a,b; Wang and Zhang 2009; Zhao et al. 2016). Zhao et al. (2016) used the BEM to analyze the pressure transient of MFHWs in tight gas reservoirs with an arbitrary outer boundary, and the solutions fit well with the solutions derived from the semianalytica...