Seismic reflection coefficient inversion method based on total variation minimized constraint

The invention provides a seismic reflection coefficient inversion method based on total variation minimized constraint, which inverses post-stack seismic data in a channel-to-channel and time window-to-time window mode. The method comprises: for a single time window, first performing Fourier transform on the post-stack seismic data in the time window and seismic wavelet extracted in advance to obtain a reflection coefficient frequency domain expression; performing Fourier transform on a time domain reflection coefficient; extracting a real part and a virtual part based on reflection coefficient even-odd separation, and constructing a corresponding solution equation; employing minimized total variation constraint to solve the equation based on a traditional conjugate gradient algorithm to obtain even and odd reflection coefficients and reconstruct an original reflection coefficient in the time window; and successively performing inversion on all time windows to obtain the reflection coefficient of a single-channel seismic recording. The method successively performs inversion on seismic refection coefficients, and is in favor of increasing a seismic data resolution and improving reservoir prediction precision.