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 rising 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 results are complicated. The two bubbles tend to move together due to the attractive force by the wake and the potential 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 interaction which cannot be neglected. (C) 2019 The Chemical Industry and Engineering Society of China, and Chemical Industry Press Co., Ltd. All rights reserved.