By developing the multiple scales method, we analytically study the dynamics properties of gap soliton of BoseEinstein condensate in optical lattices. It is shown that the gap soliton will appear at Brillouin zone edge of linear band spectrum of the condensates when the interatomic interaction strength is larger than the lattice depth. Moreover, the density of gap soliton starts to be relatively small, while it increases with time and becomes stable.