Abstract: By using the improved Keller-Miksis equation, the mutually interacting gas bubbles in a non single frequency acoustic field was investigated numerically. We assumed spherical symmetry and a coupling of the bubble oscillations, which radius were 2 μm and 5 μm respectively. The numerical results showed that the maximum radius and the secondary Bjerknes force of the coupled double bubbles were related to the amplitude, the amplitude ratio and driving frequency of the driving sound pressure. When the amplitude ratio and driving frequency were constant, the larger the driving sound pressure, the larger the maximum radius and the greater the secondary bjerknes force; When the driving sound pressure and driving frequency were constant, the larger the amplitude ratio, the larger the maximum radius and the greater the secondary bjerknes force; When the driving sound pressure and the amplitude ratio were constant, The larger the driving frequency, the smaller the maximum radius and the smaller the secondary bjerknes force.