The present disclosure provides a method for optimizing bunch distance of fractured horizontal wells of shale gas, which relates to the technical field of oil exploration. The method comprises first establishing a stress field distribution model for a single fracture; then establishing an induced stress distribution model of segmented single-bunch fracturing for a horizontal well; later establishing an induced stress distribution model of segmented multi-bunch fracturing for a horizontal well; last optimizing fracturing parameters and fracture distance according to the distribution pattern of the induced stress difference.The method considers the stress barrier, stress interference effects, and the variation of the effective net pressure during the synchronous expansion of fractures, so the calculation model is more in line with the actual working conditions, has higher precision, and can provide more accurate theoretical guidance for the optimization design of segmented multi-bunch fracturing of a horizontal well.