Abstract: Two-axis underwater channel often exists in deep ocean. Sound propagation in two-axis channel is a benchmark problem for computational methods of underwater acoustics, and some computational methods are ill-posed. In this paper, the generalized phase-integral (WKBZ) normal mode approach is extended and applied to dealing with the sound propagation in two-axis underwater channel.Calculations show that the extended WKBZ approach is effective.