Although fetal phonocardiogram (fPCG) signals have become a good indicator for discovered heart disease, they may be contaminated by various noises that reduce the signals quality and the final diagnosis decision. Moreover, the noise may cause the risk of the data to misunderstand the heart signal and to misinterpret it. The main objective of this paper is to effectively remove noise from the fPCG signal to make it clinically feasible. So, we proposed a novel noise reduction method based on Improved Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (ICEEMDAN), wavelet threshold and Crow Search Algorithm (CSA). This noise reduction method, named ICEEMDAN-DWT-CSA, has three major advantages. They were, (i) A better suppress of mode mixing and a minimized number of IMFs, (ii) A choice of wavelet corresponding to the study signal proven by the literature and (iii) Selection of the optimal threshold value. Firstly, the noisy fPCG signal is decomposed into Intrinsic Mode Functions (IMFs) by the (ICEEMDAN). Each noisy IMFs were decomposed by the Discrete Wavelet Transform (DWT). Then, the optimal threshold value using the (CSA) technique is selected and the thresholding function is carried out in the detail’s coefficients. Secondly, each denoised (IMFs) is reconstructed by applying the Inverse Discrete Wavelet Transform (IDWT). Finally, all these denoised (IMFs) are combined to get the denoised fPCG signal. The performance of the proposed method has been evaluated by Signal to Noise Ratio (SNR), Mean Square Error (MSE) and the Correlation Coefficient (COR). The experiment gave a better result than some standard methods.