Abstract: The traditional motor fault diagnosis method based on the audio signal is insufficient to obtain the feature information of the motor audio signal and the fault diagnosis precision is not high. In order to solve the above problems, a mine motor fault diagnosis method based on optimized variational mode decomposition (VMD) and convolutional neural network CNN bidirectional long short-term memory (BiLSTM) is proposed. The whale algorithm (WOA) optimized VMD is used to decompose the motor audio signal to address the issues of modal aliasing and endpoint effects. The motor audio signal is decomposed into K intrinsic mode functions (IMF). After Pearson correlation coefficient screening, the 13-dimensional static MFCC feature parameters of the main IMF component are extracted. In order to obtain the dynamic features of the signal, the first and second-order difference coefficients of the 13-dimensional static MFCC are extracted to form a 39-dimensional feature vector. By combining dynamic and static features, the performance of fault diagnosis can be improved. In order to improve the precision of fault diagnosis, a BiLSTM layer is introduced into the CNN. The CNN extracts local features of the audio signal in the spatial dimension. The BiLSTM preserves bidirectional time series information of the audio signal in the temporal dimension. It captures long-distance dependencies of the audio signal, thereby maximizing the preservation of global and local features. The experimental results show the following points. ① Each IMF component of VMD decomposition has an independent center frequency and uniform distribution, and exhibits sparsity in the frequency domain. It can effectively avoid modal aliasing problems. In IMF solving, VMD decomposition avoids endpoint effects in empirical mode decomposition (EMD) and ensemble empirical mode decomposition (EEMD) through mirror extension. ② The fault diagnosis accuracy based on 13-dimensional static MFCC features is 97.5%. The fault diagnosis accuracy based on 39-dimensional dynamic and static MFCC features is 1.11% higher than that based on 13-dimensional static MFCC features. ③ The accuracy of the diagnostic model based on CNN-BiLSTM reaches 98.61%, which is 5.83%, 4.17%, and 3.89% higher than the current universal diagnostic models CNN, BiLSTM, and CNN-LSTM, respectively.