The boundary integral method (BIM) is used to simulate the 3-D gas bubble, generated within the two bubble pulsation periods in proximity to a free surface in an inviscid, incompressible and irrotational flow. The present method is well validated by comparing the calculated shapes of the bubble and the free surface with both the experimental results and the numerical ones obtained by the Axisymmetric BIM code. The expansion, the collapse of the gas bubble and the further evolution of the rebounding non-spherical bubble are simulated. The various variation patts of the free surface spike and the bubble centroid for different standoff distances, the buoyancy parameters and the strength parameters are obtained to reveal the nonlinear interaction between the bubble and the free surface. The amplitude of the second maximum bubble volume and the four typical patts of the bubble jet and the free surface spike are examined in the context of the standoff distance. The large buoyancy is used to elevate the spray dome rather than the free surface spike.