Horizontal well drilling and hydraulic fracturing technologies play an essential role in improving gas recovery from shale reservoirs. However, they are expensive technologies and resource intensive production strategies. Optimization design of horizontal well spacing and fracture stage placement helps in striking a balance between gas production and economic benefits. However, previous researches relied mainly on numerical simulation technique which is computationally expensive and time-consuming. To reduce the computational burden, a novel multi-fidelity support vector regression (MFSVR) surrogate model assisted horizontal well spacing and fracture stage placement integrated optimization method, namely WSF-MFSVR, is proposed in this study. In the WSF-MFSVR method, both low-fidelity (LF) and high-fidelity (HF) numerical simulation models were applied to establish the multi-fidelity (MF) surrogate model so as to lessen the computational burden and guarantee its quality. In order to enhance the evaluation accuracy, the particle swarm optimization (PSO) algorithm was adopted to find the optimal hyper-parameters of the MFSVR model. Two cases with different wells and fracture types based on the shale gas reservoir with Barnett shale properties were employed to verify the WSF-MFSVR method. The results indicated that a combination of 300 H F and 3500 LF samples was the most suitable for establishing the MFSVR model to approximate the numerical simulation model. In terms of computational efficiency, the WSF-MFSVR method was about 50 times faster than the HF numerical simulation model-based method. Furthermore, the relative hyper-area difference (RHD) and overall spread (OS) of the WSF-MFSVR method were similar to that of the HF numerical simulation model-based method. However, the RHD and OS of the WSF-MFSVR method were superior to that of the LF numerical simulation model-based and single-fidelity support vector regression model-assisted methods. The data of these indexes quantitatively showed the superior convergence and diversity of the final optimal solutions obtained by WSF-MFSVR method.