This paper examines the \(L^p(\mathbb {R}^n)\) boundedness of the k-th order commutator of parabolic singular integrals. When the kernel function \(\Omega \) satisfies the condition \(\begin{aligned} \sup \limits _{\xi \in S^{n-1}}\int _{S^{n-1}}|\Omega (y) |\left( \log \frac{1}{|\langle \xi ,y\rangle |}\right) ^{1+\gamma }d\sigma (y)<\infty \end{aligned}\) with \(\gamma >k\) , we demonstrate the boundedness within the range \(\begin{aligned} \frac{2(\gamma +1)}{2(\gamma +1)-(k+1)}<p<\frac{2(\gamma +1)}{k+1}. \end{aligned}\) This work improves several previously established results.