Seismic noise attenuation plays an important role in seismic interpretation. The empirical mode decomposition, synchrosqueezing wavelet transform, variational mode decomposition, etc., are often applied trace by trace. Multivariate empirical mode decomposition, multivariate synchrosqueezing wavelet transform, and multivariate variational mode decomposition were proposed for lateral continuity consideration. Due to large input data, mini-batch multivariate variational mode decomposition is proposed in this paper. The proposed method takes advantages both of variational mode decomposition and multivariate variational mode decomposition. This proposed method firstly segments the input data into a series of smaller ones with no overlapping and then applies multivariate variational mode decomposition to these smaller ones. High frequency-domain noise is filtered through sifting. Finally, the denoised smaller ones are concatenated to form components (or intrinsic mode functions) of the input signal. Synthetic and field data experiments validate the proposed method with different batch sizes and achieve higher signal-to-noise ratio than the variational mode decomposition method.