A coupled level-set/volume-of-fluid method, under the consideration of the rheological characteristics of a fluid, is employed to investigate numerical coalescence deformation of bubble pairs generated at two adjacent nozzles in carboxymethyl cellulose aqueous solutions. The satisfactory agreement between numerical results and experimental measurements proves the validity of this approach in predicting the surface evolution of bubbles. Simulated results show that the bubble coalescence process involves four stages of independent growth, rapid mergence, radial expansion, and vertical stretching. The various effects of surfactant concentration, gas flow rate, nozzle spacing, and nozzle diameter on the aspect ratio depend greatly on each coalescence period.