Sample averaging is a widely used method in machine learning to rid the random influence. Therefore, in the system with the same physical process, an averaging operator on the unit hypersphere is constructed with the idea of Riemannian gradient. Combined with the improved banded independent component analysis in the problem of blind deconvolution for seismic data, a blind deconvolution method based on averaging over unit hypersphere is proposed. Simulation result shows the effectiveness and superiority of the method.