The uniform splitting of droplets on heterogeneous surfaces without external stimuli has gained increasing attention due to its wide and flexible applications. To precisely control and predict the droplet splitting, the 3D multiphase lattice Boltzmann model (LBM) is used to study the mechanism of the impact splitting on a hydrophilic surface with a superhydrophobic strip. A parallel implementation of the 3D multiphase LBM is presented to improve the performance by using an effective kernel fusion to reduce the off-chip memory data traffic. Two splitting modes with single or double liquid bridge are contrasted and analyzed. The results show that the size of the air bubble trapped during the impacting process leads to various splitting modes. By calculating the unbalanced Young's force in the three-phase contact lines, unbalanced surface tension at the superhydrophobic/hydrophilic surface is observed. In addition, we find that the increase in the initial velocity of the drop can not only shorten the splitting time but also lead to a change in splitting mode. The increase in the wettability contrast of the surface can stably shorten the splitting time. However the splitting time did not linearly shorten with the increase in strip width. Our findings would advance the understanding of drop impact splitting and possibly inspire new guidance for the design of functional interfaces. The present study also provides an efficient parallel implementation to simulate the 3D droplet motion on heterogeneous surfaces.