Two-axis underwater channel often exists in deep ocean. Because of the coupling between surface channel and SOFAR channel, sound propagation in the two-axis underwater channel is complex and so its calculations of acoustic fields are difficult. The beam-displacement ray-mode (BDRM) theory is a normal mode method for propagation modeling in the common horizontally stratified shallow water. We improve the theory, proposing a new method for computing the upper boundary reflection coefficient, and apply it to calculate the acoustic fields of two-axis underwater channel. Transmission losses in the two-axis underwater channel are calculated by the BDRM theory. The results are in good agreement with the KRAKEN code and the computational speed excels those of the other methods.