The invention relates to a fracture time-frequency representation based matching pursuit frequency division method, which comprises the steps of S1, calculating a complex seismic trace by Hilbert transform on the basis of a fracture network medium model, S2, calculating the instantaneous envelope, the instantaneous phase and the instantaneous frequency of the complex seismic trace, S3, finding an envelope maximum value and a corresponding time position, and acquiring an estimated value of Morlet wavelet amplitude, S4, calculating the instantaneous phase and the instantaneous frequency at the corresponding time, and acquiring estimated values of a Morlet wavelet phase angle and the main frequency of fracture medium optimization matching, S5, subtracting a real part of the determined Morlet wavelet from the original seismic trace, and S6, repeating the steps S1-S5 until the left seismic trace amplitude is less than a given threshold. The beneficial effect lies in that the fracture time-frequency representation based matching pursuit frequency division method is good in adaptivity, high in noise resistance, and good in time and frequency resolution controllability.