In our previous paper (Kanagawa et al., J. Fluid Sci. Tech., 5, 2010), we have proposed a systematic method for derivation of various types of nonlinear wave equations for plane waves in bubbly liquids. The method makes use of an asymptotic expansion with multiple scales in terms of a small wave amplitude as an expansion parameter and a set of scaling relations of physical parameters, based on basic equations of two-fluid model of bubbly flows. In this paper, we extend the method so as to handle a weakly diffracted ultrasound beam in a quiescent liquid containing a number of spherical gas bubbles distributed with a weak nonuniformity. Because of the high expandability of the original method, the extension can be accomplished by adding a scaling relation of the diameter of the beam to the original set of scaling relations. As a result, we derive a generalized Khokhlov—Zabolotskaya—Kuznetsov (KZK) equation [or a generalized Kadomtsev—Petviashvili (KP) equation] for a long wave and low frequency case.