A sparse lattice representation lattice Boltzmann method algorithm is implemented on Graphics Processing Units (GPU) to accelerate pore scale flow simuation.Prefomance testing shows that sparse lattice representation approach grately reduces memory requirement and maintains performance under low porosity compared with basic algorithm.Overall speedup reaches two orders of magnitude compared with serial code.Various factors including collision model,float number precision,and GPU that affect computing speed of the algorithm are invesgated independently.It indicates that MRT model runs as fast as LBGK model on new generation of GPU cards.While on old GPU cards,MRT model's computing speed matchs LBGK only when using single precision float.