Based on the perturbation theory and Beoulli equation,equations of aspherical oscillation of two interacting bubbles are derived.This system is then used for the numerical investigation of the deformation of the two bubbles surfaces in a spherical ultrasound field in liquid.We find that the details of the aspherical oscillation of two bubbles are shown by the analysis of a2(t) and b2(t) that describe the surface deformation of bubbles 1 and 2,respectively.