The influence of anisotropic wind velocity and temperature inhomogeneities on the attenuation of infrasound field intensity with increasing distance from a point source and on its altitude distribution is studied. The field is calculated as a function of receiver height and horizontal distance from the source using method of the pseudo-differential parabolic equation for the atmosphere with model realizations of anisotropic effective sound speed fluctuations. These realizations are obtained from the nonlinear shaping model for the gravity wave perturbations which produces the fluctuations with both the vertical and horizontal spectra consistent with the observed spectra. When propagating in the stratospheric and thermospheric wave guides the multiple scattering of infrasound field from the anisotropic fluctuations results in certain vertical wave number spectra and probability density functions of infrasound intensity fluctuations in the stratospheric (altitudes 30–40 km) and mesospheric layers (50–70 km). The statistical characteristics of the intensity fluctuations as a function of distance from the source (up to 2200 km) were studied. The same characteristics were obtained for the infrasound field scattered from the inhomogeneities whose vertical profile was retrieved from the infrasound signals from surface explosions detected in the shadow zone.