A new 2D numerical model to predict the underwater acoustic propagation is obtained by exploring the potential of the Partition of Unity Method (PUM) enriched with plane waves. The aim of the work is to obtain sound pressure level distributions when multiple operational noise sources are present, in order to assess the acoustic impact over the marine fauna. The model takes advantage of the suitability of the PUM for solving the Helmholtz equation, especially for the practical case of large domains and medium frequencies. The seawater acoustic absorption and the acoustic reflectance of the sea surface and sea bottom are explicitly considered, and perfectly matched layers (PML) are placed at the lateral artificial boundaries to avoid spurious reflexions. The model includes semi-analytical integration rules which are adapted to highly oscillatory integrands with the aim of reducing the computational cost of the integration step. In addition, we develop a novel strategy to mitigate the ill-conditioning of the elemental and global system matrices. Specifically, we compute a low-rank approximation of the local space of solutions, which in turn reduces the number of degrees of freedom, the CPU time and the memory footprint. Numerical examples are presented to illustrate the capabilities of the model and to assess its accuracy.