Abstract: Gas microchannel flows was simulated by implicit lattice Boltzmann equation with wall equilibrium boundary condition and second-order Knudsen-number-relaxation-time relation.The results of nonlinearity of pressure and slip velocity at wall agree well with Arkilic’s analytical solution,which could verify the lattice Boltzmann method in the slip flow region.