In the framework of a nonrelativistic potential quark model (NRPQM) for heavy quark system, we investigate the mass spectrum of the $P$-wave tetraquark states of $ccbar{c}bar{c}$ and $bbbar{b}bar{b}$. The Hamiltonian contains a linear confinement potential and parameterized one-gluon-exchange potential which includes a Coulomb type potential and spin-dependent potentials. The full-heavy tetraquark system is solved by a harmonic oscillator expansion method. With the same parameters fixed by the charmonium and bottomonium spectra, we obtained the full spectra for the $S$ and $P$-wave heavy tetraquark states. We find that the narrow structure around 6.9 GeV recently observed at LHCb in the di-$J/psi$ invariant mass spectrum can be naturally explained by the $P$-wave $ccbar{c}bar{c}$ states. Meanwhile, the observed broad structure around $6.2sim 6.8$ GeV can be consistently explained by the $S$-wave states around 6.5 GeV predicted in our previous work. Some contributions from those suppressed low-lying $P$-wave states around 6.7 GeV are also possible. Other decay channels are implied in such a scenario and they can be investigated by future experimental analysis. Considering the large discovery potential at LHCb, we give our predictions of the $P$-wave $bbbar{b}bar{b}$ states which can be searched for in the future.