The present disclosure provides a TOC prediction method for shale gas reservoirs, including: determining by well-seismic calibration a top interface Ttop and a bottom interface Tbottom of the shale gas reservoirs, and performing layer tracking in the entire area; converting the pre-stack CRP gather into angle gather seismic data; performing spectral shaping processing on the pre-stack migration pure wave seismic data; establishing an initial model, and then performing pre-stack simultaneous inversion to obtain P-wave impedance, S-wave impedance, primary-to-shear wave velocity ratio and density data volume; obtaining a TOC inversion volume A through a post-stack inversion; obtaining a TOC inversion volume B by calculation; then adding the well data for correction, and finally determining a planar distribution law of TOC content. The method can eliminate the multiple solutions of pre-stack inversion and improve the accuracy of TOC content prediction of shale.