A level-set method is used for analyzing the behaviors of gas bubbles in two fluids incompressible viscous flow domain. The governing equations are solved by using a finite volume method. The numerical results are verified by comparing with the experimental and other computational results. Computations for the deformations and motions of one or multi-bubbles in the flow domain with the initial undisturbed free interface are conducted. It can be seen that numerical results for different surface tension and density ratio arise very different behaviors of bubbles. When bubbles rise near the free interface, the free interface gives some great influence on the behaviors of bubbles. The present results computed by a level-set method give useful information about the properties of bubble motions and deformations.