Abstract: Owing to the limitations of both the field environment for seismic data acquisition and the performance of instruments,the seismic signals collected in seismic exploration are inevitably mixed with strong noise,thus greatly affecting the subsequent processing and interpretation.In recent years,multi-scale geometric analysis methods have become an important topic in noise suppression owing to their unique advantages.This study proposed suppressing the seismic noise using a non-local mean (NLM) algorithm in the Shearlet domain.First,the non-subsampled Shearlet transform (NSST) was performed for seismic signals.Then,the decomposed coefficient subset was further processed using the NLM method,and the weight function was improved by using eight Sobel operators to approximate the omnidirectional structure.Finally,the inverse Shearlet transform was performed for the coefficients to obtain the denoised seismic signals.Experimental results show that this combined algorithm can effectively suppress the random noise and preserve the weak events,thus showing high practicability in the seismic data processing.