Terahertz (THz) pulse imaging can be used for biomedicine, material, and food security. However, image quality is lower attributed to the THz time-domain spectroscopy system. Such as, noises and lower space resolution are presented in some THz images. To improve the THz image quality, we proposed a novel method which combining the simulation evolutionary with the symmetric fourth order partial differential equation. The image edge is first detected by using the simulation evolutionary. Then the symmetric fourth order partial differential equation (PDE) is applied on the non-edge image for noise reduction. Finally, the de-noise image is combined with the edge image obtained from the edge detection step. Experiments on four different THz images prove that the proposed method can preserve the edge information during noise reduction.