The interaction of bubbles is the key to understand gas-liquid bubbling flow.Two-dimensional axis-symmetry computational fluid dynamics simulations on the interactive bubbles were performed with VOF method,which was validated by experimental work.It is testified that several different bubble interactive behaviors could be acquired under different conditions.Firstly,for large bubbles(d: 4,6,8,10 mm),the trailing bubble ris-ing velocity and aspect ratio have negative correlations with liquid viscosity and surface tension.The influences of viscosity and surface tension on leading bubble are negligible.Secondly,for smaller bubbles(d: 1,2 mm),the re-sults are complicated.The two bubbles tend to move together due to the attractive force by the wake and the po-tential repulsive force.Especially for high viscous or high surface tension liquid,the bubble pairs undergo several times acceleration and deceleration.In addition,bubble deformation plays an important role during bubble inter-action which cannot be neglected.