This paper presents a random noise reduction method based on ensemble empirical mode decomposition (EEMD) and wavelet threshold filtering. Firstly, we have conducted spectrum analysis and analyzed the frequency band range of effective signals and noise. Secondly, we make use of EEMD method on seismic signals to obtain intrinsic mode functions (IMFs) of each trace. Then, wavelet threshold noise reduction method is used on the high frequency IMFs of each trace to obtain new high frequency IMFs. Finally, reconstruct the desired signal by adding the new high frequency IMFs on the low frequency IMFs and the trend item together. When applying our method on synthetic seismic record and field data we can get good results.