The accurate calculation of the sound field is one of the most concerning issues in hydroacoustics. The one-dimensional spectral method has been used to correctly solve simplified underwater acoustic propagation models, but it is difficult to solve actual ocean acoustic fields using this model due to its application conditions and approximation error. Therefore, it is necessary to develop a direct solution method for the two-dimensional Helmholtz equation of ocean acoustic propagation without using simplified models. Here, two commonly used spectral methods, Chebyshev–Galerkin and Chebyshev–collocation, are used to correctly solve the two-dimensional Helmholtz model equation. Since Chebyshev–collocation does not require harsh boundary conditions for the equation, it is then used to solve ocean acoustic propagation. The numerical calculation results are compared with analytical solutions to verify the correctness of the method. Compared with the mature Kraken program, the Chebyshev–collocation method exhibits higher numerical calculation accuracy. Therefore, the Chebyshev–collocation method can be used to directly solve the representative two-dimensional ocean acoustic propagation equation. Because there are no model constraints, the Chebyshev–collocation method has a wide range of applications and provides results with high accuracy, which is of great significance in the calculation of realistic ocean sound fields.